OFFSET
0,2
PROG
(PARI)
allocatemem(234567890);
A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ Using code of Michel Marcus
A252741print(up_to_n) = { my(s, i=0, n=0); for(n=0, up_to_n, if(0 == n, s = 1, if(1 == n, s = 2; lev = vector(1); lev[1] = 2, oldlev = lev; lev = vector(2*length(oldlev)); s = 1; for(i = 0, (2^(n-1))-1, lev[i+1] = if((i%2), A003961(oldlev[(i\2)+1]), 2*oldlev[(i\2)+1]); s *= lev[i+1]))); write("b252741.txt", n, " ", s/n!)); };
A252741print(6); \\ Produces the same terms.
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Dec 21 2014
STATUS
approved