n = 2的
的60个数字的循环
[(0,),(1,),(2,),(3,),(4,),(5,),(6,),(7,)
,(8,),(9,),(10,),(11,),(12,),(13,),(14,),(15,
),(16,),(17,),(18,),(19,),(20,),(21,),(22,[1,
7]),(23,),(24,),(25,),(26,),(27,),(28,),(29,)
,(30,),(31,),(32,),(33,),(34,),(35,),(36,),(3
7,),(38,),(39,),(40,),(41,),(42,),(43,),(44,[
3,0]),(45,),(46,),(47,),(48,),(49,),(50,),(51,[4,9
]),(52,),(53,),(54,),(55,),(56,),(57,),(58,),
(59,)]
的20个数字的循环
[(0,),(1,),(2,),(3,),(4,),(5,),(6,),(7,)
,(8,),(9,),(10,),(11,),(12,),(13,),(14,),(15,
),(16,),(17,),(18,),(19,)]
的3个
, ,
的12个
, , , , , , , , , , ,
的4个
, , ,
n = 2的齐了正好100个
以后只罗列最小值及其循环个数
n = 3
1
1
124
31
124
31
4
31
124
31
124
62
62
62
62
2
31
31
31
31
共20个序列
使用了三个haskell函数
let ff l = iterate (\(i, l) -> (i + 1, (tail l) ++ [(sum l) \`mod\` 10]) (0, l)
计算指定起始序列的后续序列， 无限长度
let cc l = take 2 \$ filter (\(_, x) -> x == l) ff l
得到由起始序列生成的序列的周期， 第二个(i, l)的i就是
let check l n = sort \$ map(\(_, x) -> foldl(\a, b -> 10 * a + b)0, x) \$ take n \$ ff l
得到起始序列l生成的序列的每项的转成10进制数字，按大小排列， 共n个， n由cc l得到，以便检查某个起始序列是否包含在已知循环
:lol
还可以让程序自己产生所有的不重合起始序列及其循环周期
但程序一行长度就可能太长了
在交互里有点不方便
这个总结貌似不错
http://mathmu.cn/doc/PolyFacZp.html
$ x^6-x^5-x^4-x^3-x^2-x-1 mod 5= (x^4+3*x^3+3*x^2+2*x+3)*(x+2)*(x+4) $
根据shshsh_0510上面的结论，而且如果$x^4+3x^3+3x^2+2x+3$在$F_5$中不可约，那么设这个方程在$F_{5^4}$中的解为$x_1,x_2,x_3,x_4$,在加上解$3,1$,我们有n=6时候的通项公式为
$a_1x_1^n+a_2x_2^n+a_3x_3^n+a_4x_4^n+a_5*3^n+a_6$
其中$5^4-1=624$是$x_1,x_2,x_3,x_4$的周期，也就是n=6时，数列的周期必然是624的因子。
上面只考虑了模5的情况，还要考虑模2的情况，也就是对${x^7+1}/{x+1}$在$F_2$上进行因子分解。
呵呵，刚想问你怎么不管F2了？
上面的分解是用maple，x^4+3x^3+3x^2+2x+3 应该是不可约的了。
对于x^n+x^n-1+...+1 的分解 在f2上 当n为素数幂-1时，可以分解为分园多项式，然后有个复杂的定理可以运用默比乌斯函数容易的对分园多项式进行分解。
但这个x^n-x^n-1-...-1，而且是F5上，而且是任意n，应该是没简单推导方法了。
现在好像还没有多项式算法吧。
Factor(x^7+1) mod 2=(x+1)*(x^3+x+1)*(x^3+x^2+1)
于是，方程在扩域 F(2^3)上完全分解，其根r的周期为7的约数，周期为lcd(624,7)=4368的约数