login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A205815 a(n) = Sum_{k=1..n} binomial(n,k) * sigma(n,k) * 2^(n-k). 4

%I #8 Oct 08 2016 12:26:33

%S 1,17,136,1585,16986,282338,4784900,101750689,2359918963,62200943002,

%T 1792160567088,56765070059074,1946195069937314,72080471103535786,

%U 2862427829603121696,121449533922041845569,5480386857784931063958,262149577935595804303451

%N a(n) = Sum_{k=1..n} binomial(n,k) * sigma(n,k) * 2^(n-k).

%C Here sigma(n,k) equals the sum of the k-th powers of the divisors of n.

%F Logarithmic derivative of A205814.

%F a(n) = Sum_{d|n} ((d+2)^n - 2^n).

%F a(n) ~ exp(2) * n^n. - _Vaclav Kotesovec_, Oct 08 2016

%e L.g.f.: L(x) = x + 17*x^2/2 + 136*x^3/3 + 1585*x^4/4 + 16986*x^5/5 +...

%e Exponentiation yields the g.f. of A205814:

%e exp(L(x) = 1 + x + 9*x^2 + 54*x^3 + 482*x^4 + 4239*x^5 + 55561*x^6 +...

%e Illustration of terms.

%e a(2) = 2*sigma(2,1)*2 + 1*sigma(2,2)*1 = 2*3*2 + 1*5*1 = 17;

%e a(3) = 3*sigma(3,1)*4 + 3*sigma(3,2)*2 + 1*sigma(3,3)*1 = 3*4*4 + 3*10*2 + 1*28*1 = 136;

%e a(4) = 4*sigma(4,1)*8 + 6*sigma(4,2)*4 + 4*sigma(4,3)*2 + 1*sigma(4,3)*1 = 4*7*8 + 6*21*4 + 4*73*2 + 1*273*1 = 1585.

%t Table[Sum[Binomial[n, k] * DivisorSigma[k, n] * 2^(n-k), {k, 1, n}], {n, 1, 20}] (* _Vaclav Kotesovec_, Oct 08 2016 *)

%o (PARI) {a(n)=sum(k=1, n, binomial(n, k)*sigma(n, k)*2^(n-k))}

%Y Cf. A205814 (exp), A205812.

%K nonn

%O 1,2

%A _Paul D. Hanna_, Feb 01 2012

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 13:42 EDT 2024. Contains 371254 sequences. (Running on oeis4.)