%I#29 2021年4月17日06:28:29
%S 1,10,55340253020008166870142858012501280111111340100000044,
%电话:90909091008333341852076923076924071428578557190666666680272,
%电话:625000006251280588235294117648055555556112370526315789473684220
%N 10种颜色的N珠项链数量。
%H Eric Weisstein的数学世界,<a href=“http://mathworld.wolfram.com/Necklace.html“>项链</a>
%H<a href=“/index/Ne#项链”>项链相关序列的索引条目</a>
%F a(n)=(1/n)*和{d|n}φ(d)*10^(n/d)=A054617(n)/n,n>0。
%F G.F:1-总和{n>=1}φ(n)*log(1-10*x^n)/n.-Herbert Kociemba,2016年11月2日
%F a(0)=1;a(n)=(1/n)*和{k=1..n}10^gcd(n,k).-_伊利亚·古特科夫斯基,2021年4月17日
%总资产=1+10*x+55*x^2+340*x^3+2530*x^4+20008*x^5+166870*x^6+。。。
%p with(combstruct):A:=[N,{N=循环(并集(Z$10))},未标记]:seq(计数(A,大小=N),N=0..19);#_Zerinvary Lajos,2007年12月5日
%t mx=40;系数列表[系列[1-Sum[EulerPhi[i]日志[1-10*x^i]/i,{i,1,mx}],{x,0,mx}],x](*_Herbert Kociemba_,2016年11月2日*)
%t k=10;前缀[表[DivisorSum[n,EulerPhi[#]k^(n/#)&]/n,{n,1,30}](*_Robert A.Russell_,2018年9月21日*)
%o(PARI)a(n)=如果(n==0,1,1/n*sumdiv(n,d,eulerphi(d)*10^(n/d));\\_阿尔图格·阿尔坎,2018年9月21日
%A075195的Y列10。
%Y参考A054617。
%K nonn公司
%0、2
%A _N.J.A.Sloane,2000年4月16日
%E编辑:克里斯蒂安·G·鲍尔,2002年9月7日
%E a(0)由_Herbert Kociemba修正,2016年11月2日
|