login
A351254
a(n) = A276085(gcd(n, A276086(n))).
6
0, 0, 0, 2, 0, 0, 0, 0, 0, 2, 6, 0, 0, 0, 0, 8, 0, 0, 0, 0, 6, 2, 0, 0, 0, 12, 0, 2, 0, 0, 0, 0, 0, 2, 0, 30, 0, 0, 0, 2, 6, 0, 30, 0, 0, 8, 0, 0, 0, 30, 12, 2, 0, 0, 0, 6, 30, 2, 0, 0, 0, 0, 0, 32, 0, 0, 0, 0, 0, 2, 36, 0, 0, 0, 0, 14, 0, 30, 0, 0, 6, 2, 0, 0, 30, 6, 0, 2, 0, 0, 0, 30, 0, 2, 0, 0, 0, 0, 60, 2, 6
OFFSET
0,4
LINKS
FORMULA
a(n) = A276085(A324198(n)) = A276085(gcd(n, A276086(n))).
a(n) = n - A351253(n).
PROG
(PARI)
A002110(n) = prod(i=1, n, prime(i));
A276085(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*A002110(primepi(f[k, 1])-1)); };
A324198(n) = { my(m=1, p=2, orgn=n); while(n, m *= (p^min(n%p, valuation(orgn, p))); n = n\p; p = nextprime(1+p)); (m); };
CROSSREFS
Cf. also A351234.
Cf. A324583 (positions of the rest of zeros after a(0)=0).
Sequence in context: A176192 A028633 A316898 * A356303 A238429 A332671
KEYWORD
nonn,easy
AUTHOR
Antti Karttunen, Feb 05 2022
STATUS
approved