A347160 Sum of 4th powers of distinct prime divisors of n that are <= sqrt(n).

0, 0, 0, 16, 0, 16, 0, 16, 81, 16, 0, 97, 0, 16, 81, 16, 0, 97, 0, 16, 81, 16, 0, 97, 625, 16, 81, 16, 0, 722, 0, 16, 81, 16, 625, 97, 0, 16, 81, 641, 0, 97, 0, 16, 706, 16, 0, 97, 2401, 641, 81, 16, 0, 97, 625, 2417, 81, 16, 0, 722, 0, 16, 2482, 16, 625, 97, 0, 16, 81, 3042

Offset

1,4

Links

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

Formula

G.f.: Sum_{k>=1} prime(k)^4 * x^(prime(k)^2) / (1 - x^prime(k)).

Mathematica

Table[DivisorSum[n, #^4 &, # <= Sqrt[n] && PrimeQ[#] &], {n, 1, 70}]
nmax = 70; CoefficientList[Series[Sum[Prime[k]^4 x^(Prime[k]^2)/(1 - x^Prime[k]), {k, 1, nmax}], {x, 0, nmax}], x] // Rest

Crossrefs

Cf. A001159, A005065, A005068, A063962, A097974, A098002, A347143, A347158, A347159.

Sequence in context: A326852 A070570 A347158 * A005077 A059060 A331140

Adjacent sequences: A347157 A347158 A347159 * A347161 A347162 A347163

Keyword

nonn

Author

Ilya Gutkovskiy, Aug 20 2021

Status

approved