A331601 a(n) = A002487(A241909(n)).
1, 1, 1, 2, 1, 4, 1, 3, 2, 8, 1, 7, 1, 14, 4, 3, 1, 4, 1, 11, 8, 22, 1, 9, 2, 64, 3, 43, 1, 18, 1, 5, 14, 110, 4, 9, 1, 162, 22, 47, 1, 34, 1, 127, 7, 440, 1, 13, 2, 12, 64, 191, 1, 8, 8, 97, 110, 1002, 1, 23, 1, 752, 11, 5, 14, 112, 1, 1249, 162, 16, 1, 17, 1, 610, 4, 897, 4, 220, 1, 111, 3, 4882, 1, 121, 22, 5494, 440, 281, 1, 26, 8, 7623, 1002
Offset: 1
Keywords
Links
Programs
-
PARI
A002487(n) = { my(a=1, b=0); while(n>0, if(bitand(n, 1), b+=a, a+=b); n>>=1); (b); }; \\ From A002487 A241909(n) = if(1==n||isprime(n),2^primepi(n),my(f=factor(n),h=1,i,m=1,p=1,k=1); while(k<=#f~, p = nextprime(1+p); i = primepi(f[k,1]); m *= p^(i-h); h = i; if(f[k,2]>1, f[k,2]--, k++)); (p*m)); A331601(n) = A002487(A241909(n));