login
a(n) = Sum_{p|n, p prime} p^sopf(n/p).
2

%I #4 Feb 05 2024 11:14:28

%S 0,1,1,4,1,17,1,4,27,57,1,41,1,177,368,4,1,251,1,153,2530,2169,1,41,

%T 3125,8361,27,561,1,5568,1,4,178478,131361,94932,275,1,524649,1596520,

%U 153,1,37514,1,8313,6686,8389137,1,41,823543,78157,129145076,32937,1,251

%N a(n) = Sum_{p|n, p prime} p^sopf(n/p).

%t a[n_] := Sum[p, {p, Select[Divisors[n], PrimeQ]}]; Table[DivisorSum[n, #^a[n/#] &, PrimeQ[#] &], {n, 80}]

%Y Cf. A008472 (sopf), A369744, A369911.

%K nonn,easy

%O 1,4

%A _Wesley Ivan Hurt_, Feb 05 2024