b: =proc(n,i)选项记忆;局部p;
p: =`if`(i<1,1,ithprime(i));
`如果`(n=0或i<1,1,max(b(n,i-1),
seq(p^j*b(n-p^j,i-1),j=1..ilog[p](n))
结束时间:
g: =n->b(n,`if`(n<8,3,numtheory[pi](ceil(1.328*isqrt(n*ilog(n))))):
a: =proc(n)选项记忆;局部k;
对于来自1+`if`(n=1,0,a(n-1))的k
而不是irem(g(k),k)=0做od;k个
结束时间:
seq(a(n),n=1..70)#阿洛伊斯·海因茨2013年5月22日
|