r: =proc(n)选项记住;局部q;
`如果`(n<10,n,irem(n,10,'q')*10^(长度(n)-1)+r(q))
结束时间:
b: =proc(n)选项记忆;局部k;
如果n=0,则为0
b(n-1)+1中k的else
当isprime(k)或2<>加(i[2],i=ifactors(k)[2])时
做od;k个
fi(菲涅耳)
结束时间:
a: =proc(n)选项记忆;局部m,k;
m: =r(b(n));
对于k从4到m do
如果irem(m,k)=0且不是isprime(k)且
添加(i[2],i=ifactors(k)[2])=2,然后返回kfi
od;0
结束时间:
序列(a(n),n=1..100);#Alois P.Heinz,2012年3月28日