%我#29 2021年5月8日23:04:58

%编号13,787287098742568040768964072101962770117827746413785812040，

%电话：16292367218419415066889402329808512240811265925044，

%电话：341239286758115241588538022965570

%N全13移位下长度为N的轨道数（其周期点由A001022计数）。

%GF（13）上n次一元不可约多项式的个数_罗伯特·伊斯雷尔（Robert Israel），2015年1月7日

%C带有n个13种颜色珠子的林登单词数（非周期项链）_安德鲁·霍罗伊，2017年12月10日

%H Indranil Ghosh，n的表，n=1..100的a（n）</a>

%H Yash Puri和Thomas Ward，<a href=“http://www.fq.math.ca/Scanned/39-5/puri.pdf“>卢卡斯序列特有的动力学性质，《斐波纳契季刊》，第39卷，第5期（2001年11月），第398-402页。

%H Y.Puri和T.Ward，<a href=“http://www.cs.uwaterloo.ca/journals/JIS/VOL4/WARD/short.html“>周期轨道的算法和增长</a>，J.Integer Seqs.，第4卷（2001年），#01.2.1。

%H T.病房，<a href=“http://www.mth.uea.ac.uk/~h720/research/files/integersequences.html“>完全可实现的序列</a>

%F a（n）=（1/n）*和{d|n}mu（d）13^（n/d）。

%F G.F:Sum_{k>=1}mu（k）*log（1/（1-13*x^k））/k.-_Ilya Gutkovskiy_，2019年5月19日

%e a（2）=78，因为在整个13移位中有169个周期2点和13个固定点，所以必须有长度为2的（169-13）/2=78个轨道。

%pf:=n->add（numtheory:-mobius（d）*13^（n/d），d=numtheori:-除数（n））/n；

%p序列（f（n），n=1..100）；#_罗伯特·伊斯雷尔（Robert Israel），2015年1月7日

%t a[n_]：=（1/n）*和[MoebiusMu[d]*13^（n/d），{d，除数[n]}]；表[a[n]，{n，20}]（*_Indranil Ghosh，2017年3月26日*）

%o（PARI）a（n）=sumdiv（n，d，moebius（d）*13^（n/d））/n；\\_Michel Marcus，2015年1月7日

%o（Python）

%o来自症状输入除数，mobius

%o打印（[sum（mobius（d）*13**（n//d）for d in divisors（n））//n for n in range（1,21）]）#_Indranil Ghosh_，2017年3月26日

%A074650的Y列13。

%Y参考A001022。

%K nonn公司

%O 1,1号机组

%A Thomas Ward，2001年3月21日