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
Links
- Antti Karttunen, Table of n, a(n) for n = 0..30030
Programs
-
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));