%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