login
A346252
a(n) = A342002(A328622(n)).
4
0, 1, 2, 7, 1, 5, 2, 9, 16, 47, 11, 37, 4, 13, 22, 59, 17, 49, 1, 7, 13, 41, 8, 31, 3, 11, 19, 53, 14, 43, 2, 11, 20, 61, 13, 47, 24, 83, 142, 389, 107, 319, 38, 111, 184, 473, 149, 403, 17, 69, 121, 347, 86, 277, 31, 97, 163, 431, 128, 361, 4, 15, 26, 73, 19, 59, 34, 103, 172, 449, 137, 379, 48, 131, 214, 533, 179, 463, 27, 89
OFFSET
0,3
FORMULA
a(n) = A342001(A346102(n)) = A342002(A328622(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]));
A342001(n) = (A003415(n) / A003557(n));
A346102(n) = { my(m=1, p=2); while(n, m *= (p^((((p%2)+1)*n)%p)); n = n\p; p = nextprime(1+p)); (m); };
(PARI)
A342002(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= p^(e>0); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); };
A346252(n) = A342002(A328622(n)); \\ Rest of program in A328622.
CROSSREFS
KEYWORD
nonn,base,look
AUTHOR
Antti Karttunen, Jul 11 2021
STATUS
approved