A359733 a(n) = (1/2) * Sum_{d|n} (2*d)^(n/d).
1, 4, 7, 20, 21, 88, 71, 296, 373, 1084, 1035, 5084, 4109, 16496, 20787, 67728, 65553, 286516, 262163, 1070180, 1189937, 4194568, 4194327, 17760824, 16827241, 67109228, 72150655, 269503660, 268435485, 1104603808, 1073741855, 4303389216, 4476371181
Offset: 1
Programs
-
Mathematica
a[n_] := DivisorSum[n, (2*#)^(n/#) &] / 2; Array[a, 33] (* Amiram Eldar, Aug 14 2023 *)
-
PARI
a(n) = sumdiv(n, d, (2*d)^(n/d))/2;
-
PARI
my(N=40, x='x+O('x^N)); Vec(sum(k=1, N, k*x^k/(1-2*k*x^k)))
Formula
G.f.: Sum_{k>0} k * x^k / (1 - 2 * k* x^k).