#此b文件由Antti Karttunen在GP/PARI计算器版本2.9.4(发布)amd64运行linux(x86-64/GMP-6.1.2内核)64位版本中计算,编译时间:2017年12月19日,gcc版本7.3.0(Ubuntu 7.3.0-1ubuntu1),三、先得汉斯哈维曼为A156552首批万项A156552首万项A156552首批万项A156552首万项条款准备的因子分解文件与:wget https://oeis.org/A156552/A156552.txt 的V15655252sigs=readvec(“A15655552.txt”);;\\然后把它作为一个平价向量读入。 #a23243(n)=如果(n<=2,n-1,我的(prsig=2,n-1,我的(prsig=prsig=2,n)我的(prsig=1,我的(prsig=prV156552Sigs[n],ps=prsig[1],es=prsig[2]);prod(i=1,1,#ps,((ps[i]^(1+es[i])))-1)/(ps[i]-1 1))); #A064989(n)={我的(f);f=因子(n);如果(n>1&&f[1,1]=2,f[1,2]=0),f[1,2]=0;为(i=1,#f ~,f[i,1]=预充(f[i,1]—1)预充(f[i,1]—1));因素背包(f)}(f)金(f)};f=因子(n)n)n);如果(n>1)f(n)f[#A156552(n)=如果(1==n,0,if(!(n%2),1+(2*A156552(n/2)),2*A156552(A064989(n)); #A324731(n)=(2*A156552(n))-->