Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A347232
a(n) = A346485(A276086(n)), where A346485 is Möbius transform of A342001.
4
0, 1, 1, 3, 1, 1, 1, 5, 6, 14, 4, 4, 1, 1, 2, 2, 0, 0, 1, 1, 2, 2, 0, 0, 1, 1, 2, 2, 0, 0, 1, 7, 8, 20, 6, 6, 10, 34, 44, 92, 24, 24, 6, 6, 12, 12, 0, 0, 6, 6, 12, 12, 0, 0, 6, 6, 12, 12, 0, 0, 1, 1, 2, 2, 0, 0, 4, 4, 8, 8, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 2, 2, 0, 0, 4, 4, 8, 8
OFFSET
0,4
FORMULA
a(n) = A346485(A276086(n)).
PROG
(PARI)
A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A003557(n) = (n/factorback(factorint(n)[, 1]));
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A342001(n) = (A003415(n) / A003557(n));
A346485(n) = sumdiv(n, d, moebius(n/d)*A342001(d));
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Antti Karttunen, Aug 26 2021
STATUS
approved