|
程序
|
(PARI)
分配(234567890);
A003961号(n) =my(f=因子(n));对于(i=1,#f~,f[i,1]=下一素数(f[i、1]+1));factorback(f);\\使用的代码米歇尔·马库斯
A252746print(up_to_n)={my(s,i=0,n=0);对于(n=0,up_to_n,if(0==n,s=1;lev=vector(1);lev[1]=2,oldlev=lev;lev=向量(2*长度(oldlev));s=0;对于(i=0.,(2^n)-1,lev[i+1]=if(!(i%2),A003961号(旧列夫[(i\2)+1]),2*旧列夫((i\ 2)+1);s+=如果((i%2),(lev[i+1]>lev[i]),0));写入(“b252746.txt”,n,“”,s);};
A252746打印(23);\\术语a(0)。。用该程序计算了a(23)。
|