login
A354873
a(n) = A258851(n) / gcd(A056239(n), A258851(n)).
3
1, 1, 2, 1, 7, 1, 4, 3, 11, 1, 5, 1, 3, 19, 8, 1, 33, 1, 32, 13, 7, 1, 52, 5, 25, 9, 22, 1, 53, 1, 16, 37, 31, 41, 14, 1, 35, 11, 14, 1, 73, 1, 64, 87, 41, 1, 64, 7, 85, 55, 19, 1, 135, 29, 116, 31, 49, 1, 136, 1, 53, 15, 32, 23, 107, 1, 32, 73, 117, 1, 204, 1, 61, 35, 54, 79, 127, 1, 208, 27, 67, 1, 47, 43, 71, 22
OFFSET
2,3
COMMENTS
Denominator of fraction A056239(n) / A258851(n).
FORMULA
a(n) = A258851(n) / A354871(n) = A258851(n) / gcd(A056239(n), A258851(n)).
PROG
(PARI)
A056239(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i, 2] * primepi(f[i, 1]))); }
A258851(n) = (n*sum(i=1, #n=factor(n)~, n[2, i]*primepi(n[1, i])/n[1, i])); \\ From A258851
A354873(n) = { my(u=A258851(n)); (u/gcd(u, A056239(n))); };
\\ Or:
A354873(n) = denominator(A056239(n)/A258851(n));
CROSSREFS
Cf. A056239, A258851, A278510, A354871, A354872 (numerators).
Sequence in context: A155992 A178211 A305564 * A345301 A160417 A117044
KEYWORD
nonn,frac
AUTHOR
Antti Karttunen, Jun 11 2022
STATUS
approved