|
|
A202878
|
|
Expansion of e.g.f.: exp(16*x/(1-x)) / sqrt(1-x^2).
|
|
2
|
|
|
1, 16, 289, 5776, 126025, 2972176, 75186241, 2027520784, 57988974481, 1751546371600, 55668326576641, 1855807478279056, 64713593898036889, 2354701531657512976, 89209297718289390625, 3512141211682081889296, 143435878498076017059361
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
FORMULA
|
a(n) = ( Sum_{k=0..floor(n/2)} 4^(n-2*k)/2^k * n!/((n-2*k)!*k!) )^2.
a(n) ~ n^n*exp(8*sqrt(n)-8-n)/2 * (1+22/(3*sqrt(n))). - Vaclav Kotesovec, May 23 2013
D-finite with recurrence: a(n) = (n+15)*a(n-1) + (n-1)*(n+15)*a(n-2) - (n-1)*(n-2)^2*a(n-3). - Vaclav Kotesovec, May 23 2013
|
|
EXAMPLE
|
E.g.f.: A(x) = 1 + 16*x + 289*x^2/2! + 5776*x^3/3! + 126025*x^4/4! + ...
where A(x) = 1 + 4^2*x + 17^2*x^2/2! + 76^2*x^3/3! + 355^2*x^4/4! + 1724^2*x^5/5! + ... + A202879(n)^2*x^n/n! + ...
|
|
MATHEMATICA
|
CoefficientList[Series[Exp[16*x/(1-x)]/Sqrt[1-x^2], {x, 0, 20}], x]* Range[0, 20]! (* Vaclav Kotesovec, May 23 2013 *)
|
|
PROG
|
(PARI) {a(n)=n!*polcoeff(exp(16*x/(1-x)+x*O(x^n))/sqrt(1-x^2+x*O(x^n)), n)}
(PARI) {a(n)=n!^2*polcoeff(exp(4*x+x^2/2+x*O(x^n)), n)^2}
(Magma) R<x>:=PowerSeriesRing(Rationals(), 30); Coefficients(R!(Laplace( Exp(16*x/(1-x))/Sqrt(1-x^2) ))); // G. C. Greubel, Jun 22 2022
(SageMath)
P.<x> = PowerSeriesRing(QQ, prec)
return P( exp(16*x/(1-x))/sqrt(1-x^2) ).egf_to_ogf().list()
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|