|
%I #5 Aug 20 2021 22:52:35
%S 0,0,0,8,0,8,0,8,27,8,0,35,0,8,27,8,0,35,0,8,27,8,0,35,125,8,27,8,0,
%T 160,0,8,27,8,125,35,0,8,27,133,0,35,0,8,152,8,0,35,343,133,27,8,0,35,
%U 125,351,27,8,0,160,0,8,370,8,125,35,0,8,27,476,0,35,0,8,152
%N Sum of cubes of distinct prime divisors of n that are <= sqrt(n).
%F G.f.: Sum_{k>=1} prime(k)^3 * x^(prime(k)^2) / (1 - x^prime(k)).
%t Table[DivisorSum[n, #^3 &, # <= Sqrt[n] && PrimeQ[#] &], {n, 1, 75}]
%t nmax = 75; CoefficientList[Series[Sum[Prime[k]^3 x^(Prime[k]^2)/(1 - x^Prime[k]), {k, 1, nmax}], {x, 0, nmax}], x] // Rest
%Y Cf. A001158, A005064, A005067, A063962, A097974, A098002, A280375, A347157, A347160.
%K nonn
%O 1,4
%A _Ilya Gutkovskiy_, Aug 20 2021
|
