A351253 a(n) = A276085(A351251(n)).

0, 1, 2, 1, 4, 5, 6, 7, 8, 7, 4, 11, 12, 13, 14, 7, 16, 17, 18, 19, 14, 19, 22, 23, 24, 13, 26, 25, 28, 29, 30, 31, 32, 31, 34, 5, 36, 37, 38, 37, 34, 41, 12, 43, 44, 37, 46, 47, 48, 19, 38, 49, 52, 53, 54, 49, 26, 55, 58, 59, 60, 61, 62, 31, 64, 65, 66, 67, 68, 67, 34, 71, 72, 73, 74, 61, 76, 47, 78, 79, 74, 79

Offset

0,3

Links

Antti Karttunen, Table of n, a(n) for n = 0..30030

Formula

a(n) = A276085(A351251(n)) = A276085(A276086(n) / A324198(n)).

a(n) = n - A351254(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)); };
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A351251(n) = { my(u=A276086(n)); (u/gcd(n, u)); };
A351253(n) = A276085(A351251(n));

Crossrefs

Cf. A276085, A276086, A324198, A351251, A351254.

Cf. also A351233.

Sequence in context: A058359 A261608 A351250 * A110332 A052947 A159287

Adjacent sequences: A351250 A351251 A351252 * A351254 A351255 A351256

Keyword

nonn,easy

Author

Antti Karttunen, Feb 05 2022

Status

approved