A354873 a(n) = A258851(n) / gcd(A056239(n), A258851(n)).

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

Author

Antti Karttunen, Jun 11 2022

Keywords

• nonn
• frac

Comments

Denominator of fraction A056239(n) / A258851(n).

Links

• Antti Karttunen, Table of n, a(n) for n = 2..16384
• Index entries for sequences computed from indices in prime factorization

Crossrefs

Cf. A056239, A258851, A278510, A354871, A354872 (numerators).

Programs

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));

Formula

a(n) = A258851(n) / A354871(n) = A258851(n) / gcd(A056239(n), A258851(n)).