login
A351455
a(n) = A064989(A001615(A003961(n))).
2
1, 1, 2, 2, 1, 2, 2, 4, 6, 1, 5, 4, 4, 2, 2, 8, 3, 6, 2, 2, 4, 5, 6, 8, 5, 4, 18, 4, 1, 2, 17, 16, 10, 3, 2, 12, 10, 2, 8, 4, 7, 4, 2, 10, 6, 6, 8, 16, 14, 5, 6, 8, 6, 18, 5, 8, 4, 1, 29, 4, 13, 17, 12, 32, 4, 10, 4, 6, 12, 2, 31, 24, 3, 10, 10, 4, 10, 8, 10, 8, 54, 7, 12, 8, 3, 2, 2, 20, 25, 6, 8
OFFSET
1,3
FORMULA
Multiplicative with a(p^e) = A064989((q+1)*q^(e-1)), where q = nextPrime(p) = A151800(p).
a(n) = A003557(n) * A351450(n).
PROG
(PARI)
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
A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A064989(n) = { my(f = factor(n>>valuation(n, 2))); for(i=1, #f~, f[i, 1] = precprime(f[i, 1]-1)); factorback(f); };
CROSSREFS
Coincides with A326042 on squarefree numbers (A005117, and apparently on no other numbers).
Cf. also A351441.
Sequence in context: A253193 A027869 A156748 * A075445 A216612 A194447
KEYWORD
nonn,mult
AUTHOR
Antti Karttunen, Feb 12 2022
STATUS
approved