|
|
A045445
|
|
Number of nonisomorphic systems of catafusenes for the unsymmetrical schemes (group C_s) with two appendages (see references for precise definition).
|
|
6
|
|
|
0, 1, 6, 29, 132, 590, 2628, 11732, 52608, 237129, 1074510, 4893801, 22395420, 102943815, 475139070, 2201301575, 10234016880, 47731093715, 223273611810, 1047265325255, 4924606035900, 23211459517120, 109642275853176, 518959629394294, 2460993383491632, 11691102386417575
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,3
|
|
COMMENTS
|
Number of 3-Motzkin paths of length n (i.e., lattice paths from (0,0) to (n,0) that do not go below the line y = 0 and consist of steps U = (1,1), D = (1,-1) and three types of steps H = (1,0)) that start with a U step. Example: a(4) = 29 because we have UDUD, UUDD, 9 UDHH paths, 9 UHDH paths and 9 UHHD paths. - Emeric Deutsch, Mar 26 2004
Here, n is the total number of hexagons in the system, which is usually denoted by h in most of the references below. In Cyvin, Brunvoll, and Cyvin (1992), Table 1, p. 28, it seems that the rooted hexagon is "distinguished", and the sequence is shifted by 1. - Petros Hadjicostas, May 26 2019
|
|
LINKS
|
Eric Weisstein's World of Mathematics, Fusene.
|
|
FORMULA
|
G.f.: (1/2)*(7*x^2 - 6*x + 1 + (3*x-1)*sqrt(5*x^2-6*x+1))/x^2. - Vladeta Jovovic, Jul 19 2001
a(n) = binomial(2n+2, n+1)/(n+2) + Sum_{k=1..n} binomial(2k, k)*binomial(n-1, k-1)*(3k-2n-3)/((n-k+1)*(k+1)) (n >= 2). - Emeric Deutsch, Mar 26 2004
Recurrence: (n-2)*(n+2)*a(n) = 3*(n-1)*(2*n-1)*a(n-1) - 5*(n-2)*(n-1)*a(n-2). - Vaclav Kotesovec, Oct 08 2012
a(n) = (2/(n+1))*Sum_{m=0..n-1} C(n+1,m)*C(2*n-2*m+2,n-m-1). - Vladimir Kruchinin Oct 18 2022
Let h(n) = hypergeom([-n-2, -n+2], [-n+1/2], -1/4) then a(n) = A002057(n-2)*h(n) = (2*(n-1)/(n+2))*CatalanNumber(n)*h(n). - Peter Luschny, Oct 23 2022
|
|
MAPLE
|
a := n -> binomial(2*n+2, n+1)/(n+2) + add(binomial(2*k, k)*binomial(n-1, k-1)*(3*k-2*n-3)/(n-k+1)/(k+1), k=1..n): 0, seq(a(n), n=2..23);
# Alternative:
a := n -> (2*(n - 1)/(n + 2))*(binomial(2*n, n) / (n + 1))*hypergeom([-n-2, -n+2], [-n + 1/2], -1/4): seq(simplify(a(n)), n = 1..26); # Peter Luschny, Oct 23 2022
|
|
MATHEMATICA
|
a[n_] = Binomial[2n+2, n+1]/(n+2) + Sum[Binomial[2k, k]*Binomial[n-1, k-1]*(3k-2n-3)/(n-k+1)/(k+1), {k, 1, n}];
Table[SeriesCoefficient[(1/2)*(7*x^2-6*x+1+(3*x-1)*Sqrt[5*x^2-6*x+1])/x^2, {x, 0, n}], {n, 1, 23}] (* Vaclav Kotesovec, Oct 08 2012 *)
|
|
PROG
|
(PARI) x='x+O('x^66); concat([0], Vec((1/2)*(7*x^2-6*x+1+(3*x-1)*sqrt(5*x^2-6*x+1))/x^2)) \\ Joerg Arndt, May 04 2013
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|