A122399
a(n) = Sum_{k=0..n} k^n * k! * Stirling2(n,k).
Original entry on oeis.org
1, 1, 9, 211, 9285, 658171, 68504709, 9837380491, 1863598406805, 450247033371451, 135111441590583909, 49300373690091496171, 21495577955682021043125, 11037123350952586270549531, 6591700149366720366704735109
Offset: 0
E.g.f.: A(x) = 1 + x + 9*x^2/2! + 211*x^3/3! + 9285*x^4/4! + 658171*x^5/5! + ...
such that
A(x) = 1 + (exp(x)-1) + (exp(2*x)-1)^2 + (exp(3*x)-1)^3 + (exp(4*x)-1)^4 + ...
The e.g.f. is also given by the series:
A(x) = 1/2 + exp(x)/(1+exp(x))^2 + exp(4*x)/(1+exp(2*x))^3 + exp(9*x)/(1+exp(3*x))^4 + exp(16*x)/(1+exp(4*x))^5 + exp(25*x)/(1+exp(5*x))^6 + ...
or, equivalently,
A(x) = 1/2 + exp(-x)/(1+exp(-x))^2 + exp(-2*x)/(1+exp(-2*x))^3 + exp(-3*x)/(1+exp(-3*x))^4 + exp(-4*x)/(1+exp(-4*x))^5 + exp(-5*x)/(1+exp(-5*x))^6 + ...
-
a := n -> add(k^n*k!*combinat[stirling2](n,k),k=0..n); # Max Alekseyev, Feb 01 2007
-
Flatten[{1,Table[Sum[k^n*k!*StirlingS2[n,k],{k,0,n}],{n,1,20}]}] (* Vaclav Kotesovec, Jun 21 2013 *)
-
{a(n)=polcoeff(sum(m=0, n, m^m*m!*x^m/prod(k=1, m, 1-m*k*x+x*O(x^n))), n)}
for(n=0, 20, print1(a(n), ", ")) \\ Paul D. Hanna, Jan 05 2013
-
{a(n)=n!*polcoeff(sum(k=0, n, (exp(k*x +x*O(x^n)) - 1)^k), n)}
for(n=0,25,print1(a(n),", ")) \\ Paul D. Hanna, Oct 26 2014
-
/* From e.g.f. infinite series: */
\p100 \\ set precision
{A=Vec(serlaplace(sum(n=0, 500, 1.*exp(n^2*x +O(x^26))/(1 + exp(n*x +O(x^26)))^(n+1)) ))}
for(n=0, #A-1, print1(round(A[n+1]), ", ")) \\ Paul D. Hanna, Oct 30 2014
A249459
a(n) = Sum_{k=0..n} k^(2*n).
Original entry on oeis.org
1, 1, 17, 794, 72354, 10874275, 2438235715, 762963987380, 317685943157892, 169842891165484965, 113394131858832552133, 92465351109879998121806, 90431265068257318469676710, 104479466717230437574945525959, 140782828210237288756752539959687
Offset: 0
-
Table[Sum[k^(2*n),{k,1,n}],{n,1,20}]
Table[n!*SeriesCoefficient[Sum[Exp[k^2*x], {k, 1, n}],{x,0,n}], {n,1,20}]
-
a(n)=n!*polcoeff(sum(k=0, n, exp(k*x+x*O(x^n))^k), n);
for(n=1, 20, print1(a(n), ", "))
-
my(N=20, x='x+O('x^N)); Vec(sum(k=0, N, (k^2*x)^k/(1-k^2*x))) \\ Seiichi Manyama, Dec 03 2021
A338040
E.g.f.: Sum_{j>=0} 4^j * (exp(j*x) - 1)^j.
Original entry on oeis.org
1, 4, 132, 11140, 1763076, 449262724, 168055179012, 86720706877060, 59029852191779076, 51241585497612147844, 55245853646893977682692, 72423868722672448652558980, 113447698393867318106045295876, 209271794145089904620369489016964
Offset: 0
-
Flatten[{1, Table[Sum[4^j * j^n * j! * StirlingS2[n, j], {j, 0, n}], {n, 1, 20}]}]
nmax = 20; CoefficientList[Series[1 + Sum[4^j*(Exp[j*x] - 1)^j, {j, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]!
-
{a(n)=local(X=x+x*O(x^n)); n!*polcoeff(sum(m=0, n, 4^m*(exp(m*X)-1)^m), n)}
A221077
E.g.f.: Sum_{n>=0} tanh(n*x)^n.
Original entry on oeis.org
1, 1, 8, 160, 5888, 345856, 29677568, 3502489600, 544181977088, 107675615297536, 26435436140822528, 7885689342279024640, 2809177794704769548288, 1177952320402008693538816, 574318105367992485583781888, 322156963576521588458420961280, 206009256195720974104252003647488
Offset: 0
E.g.f.: A(x) = 1 + x + 8*x^2/2! + 160*x^3/3! + 5888*x^4/4! + 345856*x^5/5! +...
where
A(x) = 1 + tanh(x) + tanh(2*x)^2 + tanh(3*x)^3 + tanh(4*x)^4 + tanh(5*x)^5 +...
-
nmax = 20; CoefficientList[Series[1 + Sum[Tanh[k*x]^k, {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, May 31 2022 *)
Join[{1}, Table[Sum[2^n * k^n * Sum[(-1)^j * Binomial[k, j] * Sum[(-1)^m * Binomial[j + m - 1, m] * StirlingS2[n, m] * m! / 2^m, {m, 1, n}], {j, 0, k}], {k, 0, n}], {n, 1, 20}]] (* Vaclav Kotesovec, Jun 01 2022 *)
-
{a(n)=local(X=x+x*O(x^n), Egf); Egf=sum(m=0, n, tanh(m*X)^m); n!*polcoeff(Egf, n)}
for(n=0,20,print1(a(n),", ") )
-
{a(n)=local(X=x+x*O(x^n), Egf); Egf=sum(m=0, n, (exp(2*m*X)-1)^m/(exp(2*m*X)+1)^m); n!*polcoeff(Egf, n)}
for(n=0,20,print1(a(n),", ") )
A221078
E.g.f.: Sum_{n>=0} tan(n*x)^n.
Original entry on oeis.org
1, 1, 8, 164, 6400, 404176, 37541888, 4814990144, 815074508800, 176018678814976, 47223034903789568, 15407438848482919424, 6007522256082907955200, 2758698201106509138251776, 1473586749521302260021198848, 905915791153129699969076117504
Offset: 0
E.g.f.: A(x) = 1 + x + 8*x^2/2! + 164*x^3/3! + 6400*x^4/4! + 404176*x^5/5! +...
where
A(x) = 1 + tan(x) + tan(2*x)^2 + tan(3*x)^3 + tan(4*x)^4 + tan(5*x)^5 +...
-
nmax = 20; CoefficientList[Series[1 + Sum[Tan[k*x]^k, {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, May 31 2022 *)
Join[{1}, Table[Sum[(-1)^((n-k)/2) * 2^n * k^n * Sum[(-1)^j * Binomial[k, j] * Sum[(-1)^m * Binomial[j + m - 1, m] * StirlingS2[n, m] * m! / 2^m, {m, 1, n}], {j, 0, k}], {k, 0, n}], {n, 1, 20}]] (* Vaclav Kotesovec, Jun 01 2022 *)
-
{a(n)=local(X=x+x*O(x^n), Egf); Egf=sum(m=0, n, tan(m*X)^m); n!*polcoeff(Egf, n)}
for(n=0,20,print1(a(n),", ") )
A245322
E.g.f.: Sum_{n>=0} sin(n*x)^n.
Original entry on oeis.org
1, 1, 8, 161, 6016, 360421, 31628288, 3823725821, 609263681536, 123729353398441, 31195066498285568, 9560281195915697081, 3500145542231863853056, 1508772905238685631514061, 756360258034794813559144448, 436312320288025061112662937941, 286966475921556619941746443288576
Offset: 0
-
nmax=20; Flatten[{1,Rest[CoefficientList[Series[Sum[Sin[k*x]^k, {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]!]}]
Flatten[{1,Table[Sum[(-1)^k * (n-2*k)^n * 2^(2*k-n) * Sum[Binomial[n-2*k,j] * (-1)^j * (n-2*k-2*j)^n,{j,0,n-2*k}],{k,0,n/2}],{n,1,20}]}]
-
{a(n)=n!*polcoeff(sum(k=0, n, sin(k*x+x*O(x^n))^k), n)}
for(n=0, 20, print1(a(n), ", "))
A195415
E.g.f.: Sum_{n>=1} tanh(n*x)^n = Sum_{n>=1} a(n)*4^(n-1)*x^n/n!.
Original entry on oeis.org
1, 2, 10, 92, 1351, 28982, 855100, 33214232, 1642999501, 100843185962, 7520379392890, 669760178257172, 70211429619908851, 8558006664633638942, 1200128210993564085880, 191861070874818576596912, 34685967730611200643509401, 7041037426518318365605795922
Offset: 1
E.g.f.: A(x) = x + 8*x^2/2! + 160*x^3/3! + 5888*x^4/4! + 345856*x^5/5! +...
or, equivalently,
A(x) = x + 2*4*x^2/2! + 10*4^2*x^3/3! + 92*4^3*x^4/4! + 1351*4^4*x^5/5! +...
where
A(x) = tanh(x) + tanh(2*x)^2 + tanh(3*x)^3 + tanh(4*x)^4 + tanh(5*x)^5 +...
-
seq(coeff(n!/4^(n-1)*series(add(tanh(n*x)^n, n = 1..100), x, 101), x, n), n = 1..100); # Peter Bala, May 29 2022
-
nmax = 20; Rest[CoefficientList[Series[Sum[Tanh[k*x]^k, {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]! / 4^(Range[0, nmax] - 1)] (* Vaclav Kotesovec, May 31 2022 *)
-
{a(n)=local(X=x+x*O(x^n), Egf); Egf=sum(m=1, n, tanh(m*X)^m); n!/4^(n-1)*polcoeff(Egf, n)}
A249489
a(n) = [x^n/n!] Sum_{k=0..n} cosh(k*x)^k.
Original entry on oeis.org
1, 0, 9, 0, 12070, 0, 126447741, 0, 5100496997940, 0, 562605048135059545, 0, 138523311740417986721274, 0, 66543520389763227261554370645, 0, 56664734898911130799849838608991176, 0, 79610326854782816434044397510470501877041
Offset: 0
-
Flatten[{1,Table[Sum[Sum[Binomial[k,j] * k^n*(k-2*j)^n / 2^k,{j,0,k}],{k,0,n}],{n,1,20}]}]
-
{a(n)=n!*polcoeff(sum(k=0, n, cosh(k*x+x*O(x^n))^k), n)}
for(n=0, 20, print1(a(n), ", "))
-
{a(n) = sum(k=0,n,sum(j=0,k, binomial(k, j) * k^n*(k-2*j)^n / 2^k ))}
for(n=0, 20, print1(a(n), ", ")) \\ Paul D. Hanna, Oct 15 2018, using Vaclav's formula.
A249698
E.g.f.: Sum_{n>=0} Product_{k=1..n} cosh(k*x).
Original entry on oeis.org
1, 0, 6, 0, 2426, 0, 7553776, 0, 90192976308, 0, 2939813898295990, 0, 213701821328573755046, 0, 30292525174041077292043440, 0, 7609302838629919155170452856136, 0, 3152886110080180503361685427596189430, 0, 2038143533263759863560759054752335955960482
Offset: 0
-
Table[n!*SeriesCoefficient[Sum[Product[Cosh[k*x],{k,1,j}],{j,0,n}],{x,0,n}],{n,0,20}]
nn=20; tab = ConstantArray[0,nn]; tab[[1]] = Series[Cosh[x],{x,0,nn}]; Do[tab[[k]] = Series[tab[[k-1]]*Cosh[k*x],{x,0,nn}],{k,2,nn}]; Flatten[{1,Table[kk!*Sum[Coefficient[tab[[k]],x^kk],{k,1,kk}],{kk,1,nn}]}] (* more efficient *)
-
{a(n)=local(X=x+x*O(x^n), Egf); Egf=sum(m=0, n, prod(k=1, m, cosh(k*X))); n!*polcoeff(Egf, n)}
for(n=0, 20, print1(a(n), ", "))
A320419
E.g.f.: Sum_{n>=0} 2^n * sinh(n*x)^n.
Original entry on oeis.org
1, 2, 32, 1298, 98816, 12116642, 2181373952, 541793612978, 177515752718336, 74174630255081282, 38495436789222735872, 24292625097918019749458, 18317925825330618728185856, 16266073932645598088605425122, 16800468023465020621665905672192, 19969924961381649826994229325322738
Offset: 0
E.g.f.: A(x) = 1 + 2*x + 32*x^2/2! + 1298*x^3/3! + 98816*x^4/4! + 12116642*x^5/5! + 2181373952*x^6/6! + 541793612978*x^7/7! + ...
such that
A(x) = 1 + 2*sinh(x) + 4*sinh(2*x)^2 + 8*sinh(3*x)^3 + 16*sinh(4*x)^4 +...
or, equivalently,
A(x) = 1 + exp(x)*(1 - exp(-2*x)) + exp(4*x)*(1 - exp(-4*x))^2 + exp(9*x)*(1 - exp(-6*x))^3 + exp(16*x)*(1 - exp(-8*x))^4 + exp(25*x)*(1 - exp(-10*x))^5 + ...
RELATED SERIES.
A(log(1+x)) = 1 + 2*x + 15*x^2 + 201*x^3 + 3807*x^4 + 93103*x^5 + 2788528*x^6 + 98816388*x^7 + 4043274742*x^8 + ... + A319466(n)*x^n + ...
A(-log(1-x)) = 1 + 2*x + 17*x^2 + 233*x^3 + 4457*x^4 + 109599*x^5 + 3294200*x^6 + 117023348*x^7 + 4796944724*x^8 + ... + A319947(n)*x^n + ...
-
{a(n) = n! * polcoeff(sum(k=0, n, 2^k * sinh(k*x + x*O(x^n))^k ), n)}
for(n=0, 30, print1(a(n), ", "))
Showing 1-10 of 10 results.
Comments