|
|
A025565
|
|
a(n) = T(n,n-1), where T is array defined in A025564.
|
|
15
|
|
|
1, 2, 4, 10, 26, 70, 192, 534, 1500, 4246, 12092, 34606, 99442, 286730, 829168, 2403834, 6984234, 20331558, 59287740, 173149662, 506376222, 1482730098, 4346486256, 12754363650, 37461564504, 110125172682, 323990062452, 953883382354
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
a(n+1) is the number of UDU-free paths of n upsteps (U) and n downsteps (D), n>=0. - David Callan, Aug 19 2004
Also the number of different integer sets { k_1, k_2, ..., k_(i+1) } with Sum_{j=1..i+1} k_j = i and k_j >= 0, see the "central binomial coefficients" (A000984), without all sets in which any two successive k_j and k_(j+1) are zero. See the partition problem eq. 3.12 on p. 19 in my dissertation below. - Eva Kalinowski, Oct 18 2012
|
|
LINKS
|
|
|
FORMULA
|
G.f.: x*sqrt((1+x)/(1-3*x)).
Sum_{i=0..n} Sum_{j=0..i} (-1)^(n-i)*a(j)*a(i-j) = 3^n. - Mario Catalani (mario.catalani(AT)unito.it), Jul 02 2003
a(1) = 1, a(n) = M(n-1) + Sum_{k=1..n-1} M(k-1)*a(n-k) with M=A001006, the Motzkin Numbers. - Reinhard Zumkeller, Mar 30 2012
D-finite with recurrence: (-n+1)*a(n) +2*(n-1)*a(n-1) +3*(n-3)*a(n-2)=0. - R. J. Mathar, Dec 02 2012
G.f.: G(0), where G(k) = 1 + 4*x*(4*k+1)/( (1+x)*(4*k+2) - x*(1+x)*(4*k+2)*(4*k+3)/(x*(4*k+3) + (1+x)*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 26 2013
a(n) = n*hypergeom([2-n, 1/2-n/2, 1-n/2], [2, -n], 4). - Peter Luschny, Jul 12 2016
a(n) = (-1)^n*2*hypergeom([3/2, 2-n], [2], 4) for n > 1. - Peter Luschny, Jan 30 2017
|
|
EXAMPLE
|
G.f. = x + 2*x^2 + 4*x^3 + 10*x^4 + 26*x^5 + 70*x^6 + 192*x^7 + 534*x^8 + ...
|
|
MAPLE
|
seq( add(binomial(i-2, k)*(binomial(i-k, k+1)), k=0..floor(i/2)), i=1..30 ); # Detlef Pauly (dettodet(AT)yahoo.de), Nov 09 2001
# Alternatively:
a := n -> `if`(n=1, 1, 2*(-1)^n*hypergeom([3/2, 2-n], [2], 4)):
|
|
MATHEMATICA
|
T[_, 0] = 1; T[1, 1] = 2; T[n_, k_] /; 0 <= k <= 2n := T[n, k] = T[n-1, k-2] + T[n-1, k-1] + T[n-1, k]; T[_, _] = 0;
a[n_] := T[n-1, n-1];
|
|
PROG
|
(Haskell)
a025565 n = a025565_list !! (n-1)
a025565_list = 1 : f a001006_list [1] where
f (x:xs) ys = y : f xs (y : ys) where
y = x + sum (zipWith (*) a001006_list ys)
(Sage)
def A():
a, b, n = 1, 1, 1
yield a
while True:
yield a + b
n += 1
a, b = b, ((3*(n-1))*a+(2*n-1)*b)//n
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|