login
A348991
a(n) = A333791(A276086(n)).
1
0, 0, 0, 1, 0, 4, 0, 3, 2, 12, 8, 39, 0, 18, 12, 67, 48, 214, 0, 93, 62, 342, 248, 1089, 0, 468, 312, 1717, 1248, 5464, 0, 5, 4, 20, 16, 65, 2, 30, 24, 109, 90, 346, 12, 155, 124, 554, 460, 1751, 62, 780, 624, 2779, 2310, 8776, 312, 3905, 3124, 13904, 11560, 43901, 0, 40, 32, 153, 128, 492, 16, 219, 178, 788, 664, 2495
OFFSET
0,6
PROG
(PARI)
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A332993(n) = if(1==n, n, n + A332993(n/vecmin(factor(n)[, 1])));
A332994(n) = if(1==n, n, n + A332994(n/vecmax(factor(n)[, 1])));
A333791(n) = (A332993(n)-A332994(n));
CROSSREFS
Cf. A060735 (gives the positions of other zeros after the initial a(0)=0).
Sequence in context: A096303 A262695 A021252 * A195775 A363728 A362562
KEYWORD
nonn,base,look
AUTHOR
Antti Karttunen, Nov 07 2021
STATUS
approved