login
A323903
a(n) = A002487(A122111(n)).
6
1, 1, 1, 2, 1, 2, 1, 3, 4, 2, 1, 3, 1, 2, 4, 3, 1, 4, 1, 3, 4, 2, 1, 3, 8, 2, 7, 3, 1, 4, 1, 5, 4, 2, 8, 8, 1, 2, 4, 3, 1, 4, 1, 3, 7, 2, 1, 5, 14, 12, 4, 3, 1, 9, 8, 3, 4, 2, 1, 8, 1, 2, 7, 5, 8, 4, 1, 3, 4, 12, 1, 6, 1, 2, 18, 3, 14, 4, 1, 5, 9, 2, 1, 8, 8, 2, 4, 3, 1, 9, 14, 3, 4, 2, 8, 5, 1, 16, 7, 6, 1, 4, 1, 3, 18
OFFSET
1,4
FORMULA
a(n) = A002487(A122111(n)) = A002487(A322865(n)).
a(p) = 1 for all primes p.
PROG
(PARI)
A002487(n) = { my(a=1, b=0); while(n>0, if(bitand(n, 1), b+=a, a+=b); n>>=1); (b); }; \\ From A002487
A064989(n) = {my(f); f = factor(n); if((n>1 && f[1, 1]==2), f[1, 2] = 0); for (i=1, #f~, f[i, 1] = precprime(f[i, 1]-1)); factorback(f)};
A122111(n) = if(1==n, n, prime(bigomega(n))*A122111(A064989(n)));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 09 2019
STATUS
approved