

A345310


a(n) = Sum_{pn, p prime} p^gcd(p,n/p).


0



0, 2, 3, 4, 5, 5, 7, 4, 27, 7, 11, 7, 13, 9, 8, 4, 17, 29, 19, 9, 10, 13, 23, 7, 3125, 15, 27, 11, 29, 10, 31, 4, 14, 19, 12, 31, 37, 21, 16, 9, 41, 12, 43, 15, 32, 25, 47, 7, 823543, 3127, 20, 17, 53, 29, 16, 11, 22, 31, 59, 12, 61, 33, 34, 4, 18, 16, 67, 21, 26, 14, 71, 31
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,2


COMMENTS

a(p) = Sum_{pp} p^gcd(p,p/p) = p^1 = p, for p prime.


LINKS



EXAMPLE

a(18) = Sum_{p18} p^gcd(p,18/p) = 2^gcd(2,9) + 3^gcd(3,6) = 2^1 + 3^3 = 29.


MATHEMATICA

Table[Sum[k^GCD[k, n/k] (PrimePi[k]  PrimePi[k  1]) (1  Ceiling[n/k] + Floor[n/k]), {k, n}], {n, 100}]


CROSSREFS



KEYWORD

nonn


AUTHOR



STATUS

approved



