(平价)
A064989号(n) ={my(f);f=因子(n);如果((n>1&&f[1,1]==2),f[1,2]=0;对于(i=1,#f~,f[i,1]=预素数(f[i,1]-1));factorback(f)};
邮编:A156552(n) =if(1==n,0,if(!(n%2),1+(2*邮编:A156552(n/2)),2*邮编:A156552(A064989号(n) ));
A323243(n) =如果(1==n,0,西格玛(邮编:A156552(n) );
(平价)
\\对于计算术语a(n),n>~4000,请使用Hans Havermann的因子分解文件https://oeis.org/邮编:A156552/a156552.txt
v156552sigs=readvec(“a156552.txt”);\\首先将其作为PARI向量读入。
A323243(n) =如果(n<=2,n-1,my(prsig=v156552sigs[n],ps=prsig[1],es=prsig[2]);prod(i=1,#ps,((ps[i]^(1+es[i])-1)/(ps[i]-1));\\n然后播放sigma
\\安蒂·卡尔图宁2019年3月15日
|