|
| |
|
|
A101786
|
|
G.f. satisfies: A(x) = 1 + x*A(x)/(1 - 2*x^2*A(x)^2).
|
|
0
|
|
|
|
1, 1, 1, 3, 9, 25, 77, 247, 801, 2657, 8969, 30635, 105785, 368745, 1295493, 4582767, 16309953, 58357313, 209798289, 757461011, 2745281705, 9984464761, 36428252541, 133293594343, 489028250465, 1798543861537, 6629635284505
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
0,4
|
|
|
COMMENTS
|
Formula may be derived using the Lagrange Inversion theorem (cf. A049124).
|
|
|
LINKS
|
Table of n, a(n) for n=0..26.
|
|
|
FORMULA
|
a(n) = Sum_{k=0..[(n-1)/2]} C(n-k-1, k)*C(n, 2*k)*2^k/(2*k+1) for n>0, with a(0)=1. G.f.: A(x) = [ series reversion of x*(1-2*x^2)/(1+x-2*x^2) ]/x.
|
|
|
EXAMPLE
|
Generated from Fibonacci polynomials (A011973) and
coefficients of odd powers of 1/(1-x):
a(1) = 1*1/1
a(2) = 1*1/1 + 0*1*2/3
a(3) = 1*1/1 + 1*3*2/3
a(4) = 1*1/1 + 2*6*2/3 + 0*1*2^2/5
a(5) = 1*1/1 + 3*10*2/3 + 1*5*2^2/5
a(6) = 1*1/1 + 4*15*2/3 + 3*15*2^2/5 + 0*1*2^3/7
a(7) = 1*1/1 + 5*21*2/3 + 6*35*2^2/5 + 1*7*2^3/7
a(8) = 1*1/1 + 6*28*2/3 + 10*70*2^2/5 + 4*28*2^3/7 + 0*1*2^4/9
This process is equivalent to the formula:
a(n) = Sum_{k=0..[(n-1)/2]} C(n-k-1,k)*C(n,2*k)*2^k/(2*k+1).
|
|
|
PROG
|
(PARI) {a(n)=if(n==0, 1, sum(k=0, (n-1)\2, binomial(n-k-1, k)*binomial(n, 2*k)*2^k/(2*k+1)))}
|
|
|
CROSSREFS
|
Cf. A011973, A001045, A049124, A101785.
Sequence in context: A001189 A212352 A198180 * A217995 A192371 A192465
Adjacent sequences: A101783 A101784 A101785 * A101787 A101788 A101789
|
|
|
KEYWORD
|
nonn
|
|
|
AUTHOR
|
Paul D. Hanna, Dec 16 2004
|
|
|
STATUS
|
approved
|
| |
|
|