A337341 a(n) = A337339(n) - n.

0, 3, 10, 37, 20, 107, 54, 357, 304, 211, 74, 1001, 132, 95, 598, 3265, 164, 2795, 246, 1965, 1492, 739, 398, 9089, 1176, 1275, 7786, 349, 452, 5483, 654, 29493, 2080, 1591, 2930, 25277, 804, 2343, 3574, 17821, 884, 13571, 1062, 6801, 15268, 3739, 1358, 81965, 7272, 10755, 4462, 11653, 1688, 70259, 4086, 8765, 6556, 807

Offset: 1

Author

Antti Karttunen, Aug 24 2020

Keywords

• nonn

Comments

Superficially, it seems that A000027 offers an lower bound for A337339 as there seems to be no negative terms in this sequence, but of course it is not guaranteed.

Links

• Antti Karttunen, Table of n, a(n) for n = 1..8191
• Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
• Index entries for sequences computed from indices in prime factorization
• Index entries for sequences related to sigma(n)

Crossrefs

Cf. A000027, A337339, A337340.

Programs

PARI

A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A337341(n) = { my(h=A003961(n),s=(h^2),u=(s+1)/2); ((u/gcd(1+sigma(s), u)) - n); };

Formula

a(n) = A337339(n) - n.