%I #15 Mar 02 2022 10:02:22
%S 0,4,27,4,3125,31,823543,4,27,3129,285311670611,31,302875106592253,
%T 823547,3152,4,827240261886336764177,31,1978419655660313589123979,
%U 3129,823570,285311670615,20880467999847912034355032910567,31,3125,302875106592257
%N a(n) = Sum_{p|n, p prime} p^p.
%H Seiichi Manyama, <a href="/A351366/b351366.txt">Table of n, a(n) for n = 1..388</a>
%F G.f.: Sum_{k>=1} (prime(k) * x)^prime(k)/(1 - x^prime(k)). - _Seiichi Manyama_, Mar 02 2022
%e a(6) = 31; a(6) = Sum_{p|6, p prime} p^p = 2^2 + 3^3 = 31.
%o (PARI) my(N=40, x='x+O('x^N)); concat(0, Vec(sum(k=1, N, isprime(k)*(k*x)^k/(1-x^k)))) \\ _Seiichi Manyama_, Mar 02 2022
%Y Cf. A001221, A351365.
%K nonn
%O 1,2
%A _Wesley Ivan Hurt_, Feb 08 2022
|