login
A342917
a(n) = A001615(n) / gcd(1+n, A001615(n)), where A001615 is Dedekind psi, n * Product_{p|n, p prime} (1 + 1/p).
3
1, 1, 1, 6, 1, 12, 1, 4, 6, 18, 1, 24, 1, 8, 3, 24, 1, 36, 1, 12, 16, 36, 1, 48, 15, 14, 9, 48, 1, 72, 1, 16, 24, 54, 4, 72, 1, 20, 7, 72, 1, 96, 1, 8, 36, 72, 1, 96, 28, 30, 18, 84, 1, 108, 9, 32, 40, 90, 1, 144, 1, 32, 3, 96, 14, 144, 1, 36, 48, 144, 1, 144, 1, 38, 30, 120, 16, 168, 1, 16, 54, 126, 1, 192, 54, 44, 15, 144
OFFSET
1,4
COMMENTS
The scatter plot shows two distinct "fans" separated by a gap. Why?
LINKS
FORMULA
a(n) = A001615(n) / A342915(n) = A001615(n) / gcd(1+n, A001615(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
A342917(n) = { my(u=A001615(n)); (u/gcd(1+n, u)); };
CROSSREFS
Cf. also A160595.
Sequence in context: A216605 A342635 A107348 * A040041 A162933 A304252
KEYWORD
nonn,look
AUTHOR
Antti Karttunen, Mar 29 2021
STATUS
approved