A345301 a(n) = Sum_{p|n, p prime} p^pi(n/p).

0, 1, 1, 2, 1, 7, 1, 4, 9, 13, 1, 17, 1, 23, 52, 16, 1, 43, 1, 41, 130, 43, 1, 113, 125, 77, 81, 113, 1, 270, 1, 64, 364, 145, 968, 371, 1, 275, 898, 881, 1, 1328, 1, 377, 1354, 535, 1, 1241, 2401, 1137, 2476, 681, 1, 2699, 4456, 2913, 6922, 1053, 1, 10710, 1, 2079, 8962

Offset

1,4

Comments

If p is prime, a(p) = Sum_{p|p} p^pi(p/p) = p^0 = 1.

Links

Table of n, a(n) for n=1..63.

Example

a(12) = Sum_{p|12} p^pi(12/p) = 2^pi(6) + 3^pi(4) = 2^3 + 3^2 = 17.

Mathematica

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

Prog

(Python)
from sympy import primefactors, primepi
def A345301(n): return sum(p**primepi(n//p) for p in primefactors(n)) # Chai Wah Wu, Jun 13 2021

Crossrefs

Cf. A000720.

Sequence in context: A178211 A305564 A354873 * A160417 A117044 A280691

Adjacent sequences: A345298 A345299 A345300 * A345302 A345303 A345304

Keyword

nonn

Author

Wesley Ivan Hurt, Jun 13 2021

Status

approved