A344699 a(n) = A344697(A108951(n)).
1, 1, 1, 6, 1, 6, 1, 4, 72, 6, 1, 4, 1, 6, 72, 24, 1, 48, 1, 4, 72, 6, 1, 24, 2160, 6, 18, 4, 1, 48, 1, 16, 72, 6, 2160, 288, 1, 6, 72, 24, 1, 48, 1, 4, 18, 6, 1, 16, 5760, 288, 72, 4, 1, 108, 2160, 24, 72, 6, 1, 288, 1, 6, 18, 96, 2160, 48, 1, 4, 72, 288, 1, 64, 1, 6, 108, 4, 5760, 48, 1, 16, 2592, 6, 1, 288, 2160
Offset: 1
Keywords
Links
Crossrefs
Programs
-
PARI
A034386(n) = prod(i=1, primepi(n), prime(i)); A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) }; A001615(n) = if(1==n,n, my(f=factor(n)); prod(i=1, #f~, f[i, 1]^f[i, 2] + f[i, 1]^(f[i, 2]-1))); \\ After code in A001615 A344697(n) = { my(u=A001615(n)); (u/gcd(u,sigma(n))); }; A344699(n) = A344697(A108951(n));