|
|
A025231
|
|
a(n) = a(1)*a(n-1) + a(2)*a(n-2) + ...+ a(n-1)*a(1) for n >= 3.
|
|
5
|
|
|
2, 3, 12, 57, 300, 1686, 9912, 60213, 374988, 2381322, 15361896, 100389306, 663180024, 4421490924, 29712558576, 201046204173, 1368578002188, 9366084668802, 64403308499592, 444739795023054, 3082969991029800
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (1 - sqrt(1 - 8*x + 4*x^2))/2. - Michael Somos, Jun 08 2000
n*a(n) = (8*n - 12)*a(n - 1) - (4*n - 12)*a(n - 2). [Richard Choulet, Dec 16 2009]
G.f.: 1 + x - G(0); G(k) = 1 + 2*x - 3*x/G(k + 1); (continued fraction, 1-step). - Sergei N. Gladkovskii, Jan 05 2012
a(n) ~ sqrt(4*sqrt(3) - 6)*(4 + 2*sqrt(3))^n/(4*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Oct 07 2012
0 = a(n)*(+16*a(n+1) - 80*a(n+2) + 16*a(n+3)) + a(n+1)*(+16*a(n+1) + 56*a(n+2) - 20*a(n+3)) + a(n+2)*(+4*a(n+2) + a(n+3)) if n>0. - Michael Somos, Apr 10 2014
a(n) = Sum_{k=0..n} C(k+1,n-k)*2^(2*k+1-n)*(-1)^(n-k)*C(2*k,k)/(k+1). - Vladimir Kruchinin, Apr 21 2023
|
|
EXAMPLE
|
G.f. = 2*x + 3*x^2 + 12*x^3 + 57*x^4 + 300*x^5 + 1686*x^6 + 9912*x^7 + ...
|
|
MATHEMATICA
|
Table[SeriesCoefficient[(1 - Sqrt[1 - 8*x + 4*x^2])/2, {x, 0, n}], {n, 1, 20}] (* Vaclav Kotesovec, Oct 07 2012 *)
|
|
PROG
|
(PARI) a(n)=polcoeff((1-sqrt(1-8*x+4*x^2+x*O(x^n)))/2, n)
(Maxima)
a(n):=sum(binomial(k+1, n-k)*2^(2*k+1-n)*(-1)^(n-k)*binomial(2*k, k)/(k+1), k, 0, n); /* _Vladimir Kruchinin, Apr 21 2023 */
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,eigen
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|