|
%I #25 Jun 02 2019 04:37:20
%S 1,5,10,25,26,80,50,161,163,290,122,988,170,796,1580,2305,290,5561,
%T 362,10670,9404,5912,530,58436,16251,19258,66340,118640,842,381740,
%U 962,431105,547172,268214,509500,3534037,1370,1056880,4813052,8616326,1682
%N a(n) = Sum_{d divides n} d^(n/d + 1).
%H Seiichi Manyama, <a href="/A078308/b078308.txt">Table of n, a(n) for n = 1..6284</a>
%F G.f.: Sum_{n>0} n^2*x^n/(1-n*x^n).
%F L.g.f.: -log(Product_{ k>0 } (1-k*x^k)) = Sum_{ n>=0 } (a(n)/n)*x^n. - _Benedict W. J. Irwin_, Jul 04 2016
%p A078308 := proc(n)
%p add( d^(n/d+1),d=numtheory[divisors](n)) ;
%p end proc:
%p seq(A078308(n),n=1..10) ; # _R. J. Mathar_, Dec 14 2011
%t Table[CoefficientList[Series[-Log[Product[(1 - k x^k), {k, 1, 60}]], {x, 0, 60}],x][[n + 1]] (n), {n, 1, 60}] (* _Benedict W. J. Irwin_, Jul 04 2016 *)
%o (PARI) a(n) = sumdiv(n, d, d^(n/d+1)); \\ _Michel Marcus_, Jul 04 2016
%o (PARI) N=66; x='x+O('x^N); Vec(x*deriv(-log(prod(k=1, N, 1-k*x^k)))) \\ _Seiichi Manyama_, Jun 02 2019
%Y Cf. A055225, A006906, A022661.
%K nonn
%O 1,2
%A _Vladeta Jovovic_, Nov 22 2002
|
