cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 15 results. Next

A026300 Motzkin triangle, T, read by rows; T(0,0) = T(1,0) = T(1,1) = 1; for n >= 2, T(n,0) = 1, T(n,k) = T(n-1,k-2) + T(n-1,k-1) + T(n-1,k) for k = 1,2,...,n-1 and T(n,n) = T(n-1,n-2) + T(n-1,n-1).

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 1, 3, 5, 4, 1, 4, 9, 12, 9, 1, 5, 14, 25, 30, 21, 1, 6, 20, 44, 69, 76, 51, 1, 7, 27, 70, 133, 189, 196, 127, 1, 8, 35, 104, 230, 392, 518, 512, 323, 1, 9, 44, 147, 369, 726, 1140, 1422, 1353, 835, 1, 10, 54, 200, 560, 1242, 2235, 3288, 3915, 3610, 2188
Offset: 0

Views

Author

Keywords

Comments

Right-hand columns have g.f. M^k, where M is g.f. of Motzkin numbers.
Consider a semi-infinite chessboard with squares labeled (n,k), ranks or rows n >= 0, files or columns k >= 0; number of king-paths of length n from (0,0) to (n,k), 0 <= k <= n, is T(n,n-k). - Harrie Grondijs, May 27 2005. Cf. A114929, A111808, A114972.

Examples

			Triangle starts:
  [0] 1;
  [1] 1, 1;
  [2] 1, 2,  2;
  [3] 1, 3,  5,   4;
  [4] 1, 4,  9,  12,   9;
  [5] 1, 5, 14,  25,  30,  21;
  [6] 1, 6, 20,  44,  69,  76,   51;
  [7] 1, 7, 27,  70, 133, 189,  196,  127;
  [8] 1, 8, 35, 104, 230, 392,  518,  512,  323;
  [9] 1, 9, 44, 147, 369, 726, 1140, 1422, 1353, 835.
		

References

  • Harrie Grondijs, Neverending Quest of Type C, Volume B - the endgame study-as-struggle.
  • A. Nkwanta, Lattice paths and RNA secondary structures, in African Americans in Mathematics, ed. N. Dean, Amer. Math. Soc., 1997, pp. 137-147.

Crossrefs

Reflected version is in A064189.
Row sums are in A005773.
T(n,n) are Motzkin numbers A001006.
Other columns of T include A002026, A005322, A005323.

Programs

  • Haskell
    a026300 n k = a026300_tabl !! n !! k
    a026300_row n = a026300_tabl !! n
    a026300_tabl = iterate (\row -> zipWith (+) ([0,0] ++ row) $
                                    zipWith (+) ([0] ++ row) (row ++ [0])) [1]
    -- Reinhard Zumkeller, Oct 09 2013
    
  • Maple
    A026300 := proc(n,k)
       add(binomial(n,2*i+n-k)*(binomial(2*i+n-k,i) -binomial(2*i+n-k,i-1)), i=0..floor(k/2));
    end proc: # R. J. Mathar, Jun 30 2013
  • Mathematica
    t[n_, k_] := Sum[ Binomial[n, 2i + n - k] (Binomial[2i + n - k, i] - Binomial[2i + n - k, i - 1]), {i, 0, Floor[k/2]}]; Table[ t[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Robert G. Wilson v, Jan 03 2011 *)
    t[, 0] = 1; t[n, 1] := n; t[n_, k_] /; k>n || k<0 = 0; t[n_, n_] := t[n, n] = t[n-1, n-2]+t[n-1, n-1]; t[n_, k_] := t[n, k] = t[n-1, k-2]+t[n-1, k-1]+t[n-1, k]; Table[t[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Apr 18 2014 *)
    T[n_, k_] := Binomial[n, k] Hypergeometric2F1[1/2 - k/2, -k/2, n - k + 2, 4];
    Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Peter Luschny, Mar 21 2018 *)
  • PARI
    tabl(nn) = {for (n=0, nn, for (k=0, n, print1(sum(i=0, k\2, binomial(n, 2*i+n-k)*(binomial(2*i+n-k, i)-binomial(2*i+n-k, i-1))), ", ");); print(););} \\ Michel Marcus, Jul 25 2015

Formula

T(n,k) = Sum_{i=0..floor(k/2)} binomial(n, 2i+n-k)*(binomial(2i+n-k, i) - binomial(2i+n-k, i-1)). - Herbert Kociemba, May 27 2004
T(n,k) = A027907(n,k) - A027907(n,k-2), k<=n.
Sum_{k=0..n} (-1)^k*T(n,k) = A099323(n+1). - Philippe Deléham, Mar 19 2007
Sum_{k=0..n} (T(n,k) mod 2) = A097357(n+1). - Philippe Deléham, Apr 28 2007
Sum_{k=0..n} T(n,k)*x^(n-k) = A005043(n), A001006(n), A005773(n+1), A059738(n) for x = -1, 0, 1, 2 respectively. - Philippe Deléham, Nov 28 2009
T(n,k) = binomial(n, k)*hypergeom([1/2 - k/2, -k/2], [n - k + 2], 4). - Peter Luschny, Mar 21 2018
T(n,k) = [t^(n-k)] [x^n] 2/(1 - (2*t + 1)*x + sqrt((1 + x)*(1 - 3*x))). - Peter Luschny, Oct 24 2018
The n-th row polynomial R(n,x) equals the n-th degree Taylor polynomial of the function (1 - x^2)*(1 + x + x^2)^n expanded about the point x = 0. - Peter Bala, Feb 26 2023

Extensions

Corrected and edited by Johannes W. Meijer, Oct 05 2010

A091867 Triangle read by rows: T(n,k) = number of Dyck paths of semilength n having k peaks at odd height.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 1, 3, 0, 1, 3, 4, 6, 0, 1, 6, 15, 10, 10, 0, 1, 15, 36, 45, 20, 15, 0, 1, 36, 105, 126, 105, 35, 21, 0, 1, 91, 288, 420, 336, 210, 56, 28, 0, 1, 232, 819, 1296, 1260, 756, 378, 84, 36, 0, 1, 603, 2320, 4095, 4320, 3150, 1512, 630, 120, 45, 0, 1, 1585, 6633, 12760, 15015, 11880, 6930, 2772, 990, 165, 55, 0, 1
Offset: 0

Views

Author

Emeric Deutsch, Mar 10 2004

Keywords

Comments

Number of ordered trees with n edges having k leaves at odd height. Row sums are the Catalan numbers (A000108). T(n,0)=A005043(n). Sum_{k=0..n} k*T(n,k) = binomial(2n-2,n-1).
T(n,k)=number of Dyck paths of semilength n and having k ascents of length 1 (an ascent is a maximal string of consecutive up steps). Example: T(4,2)=6 because we have UdUduud, UduuddUd, uuddUdUd, uudUdUdd, UduudUdd and uudUddUd (the ascents of length 1 are indicated by U instead of u).
T(n,k) is the number of Łukasiewicz paths of length n having k level steps (i.e., (1,0)). A Łukasiewicz path of length n is a path in the first quadrant from (0,0) to (n,0) using rise steps (1,k) for any positive integer k, level steps (1,0) and fall steps (1,-1) (see R. P. Stanley, Enumerative Combinatorics, Vol. 2, Cambridge Univ. Press, Cambridge, 1999, p. 223, Exercise 6.19w; the integers are the slopes of the steps). Example: T(4,1)=4 because we have HU(2)DD, U(2)HDD, U(2)DHD and U(2)DDH, where H=(1,0), U(1,1), U(2)=(1,2) and D=(1,-1). - Emeric Deutsch, Jan 06 2005
T(n,k) = number of noncrossing partitions of [n] containing k singleton blocks. Also, T(n,k) = number of noncrossing partitions of [n] containing k adjacencies. An adjacency is an occurrence of 2 consecutive integers in the same block (here 1 and n are considered consecutive). In fact, the statistics # singletons and # adjacencies have a symmetric joint distribution.
Exponential Riordan array [e^x*(Bessel_I(0,2x)-Bessel_I(1,2x)),x]. - Paul Barry, Mar 03 2011
T(n,k) is the number of ordered trees having n edges and exactly k nodes with one child. - Geoffrey Critzer, Feb 25 2013
From Tom Copeland, Nov 04 2014: (Start)
Summing the coeff. of the partitions in A134264 for a Lagrange inversion formula (see also A249548) containing (h_1)^k = (1')^k gives this triangle, so this array's o.g.f. H(x,t) = x + t * x^2 + (1 + t^2) * x^3 ... is the inverse of the o.g.f. of A104597 with a sign change, i.e., H^(-1)(x,t) = (x-x^2) / [1 + (t-1)(x-x^2)] = Cinv(x)/[1 + (t-1)Cinv(x)] = P[Cinv(x),t-1] where Cinv(x)= x * (1-x) is the inverse of C(x) = [1-sqrt(1-4*x)]/2, an o.g.f. for the Catalan numbers A000108, and P(x,t) = x/(1+t*x) with inverse Pinv(x,t) = -P(-x,t) = x/(1-t*x). Therefore,
O.g.f.: H(x,t) = C[Pinv(x,t-1)] = C[P(x,1-t)] = C[x/(1-(t-1)x)] = {1-sqrt[1-4*x/(1-(t-1)x)]}/2 (for A091867). Reprising,
Inverse O.g.f.: H^(-1)(x,t) = x*(1-x) / [1 + (t-1)x(1-x)] = P[Cinv(x),t-1].
From general arguments in A134264, the row polynomials are an Appell sequence with lowering operator d/dt, having the umbral property (p(.,t)+a)^n=p(n,t+a) with e.g.f. = e^(x*t)/w(x), where 1/w(x)= e.g.f. of first column for the Motzkin numbers in A005043. (Mislabeled argument corrected on Jan 31 2016.)
Cf. A124644 (t-shifted polynomials), A026378 (t=-4), A001700 (t=-3), A005773 (t=-2), A126930 (t=-1) and A210736 (t=-1, a(0)=0, unsigned), A005043 (t=0), A000108 (t=1), A007317 (t=2), A064613 (t=3), A104455 (t=4), A030528 (for inverses).
(End)
The sequence of binomial transforms A126930, A005043, A000108, ... in the above comment appears in A126930 and the link therein to a paper by F. Fite et al. on page 42. - Tom Copeland, Jul 23 2016

Examples

			T(4,2)=6 because we have (ud)uu(ud)dd, uu(ud)dd(ud), uu(ud)(ud)dd, (ud)(ud)uudd, (ud)uudd(ud) and uudd(ud)(ud) (here u=(1,1), d=(1,-1) and the peaks at odd height are shown between parentheses).
Triangle begins:
   1;
   0,   1;
   1,   0,   1;
   1,   3,   0,   1;
   3,   4,   6,   0,  1;
   6,  15,  10,  10,  0,  1;
  15,  36,  45,  20, 15,  0, 1;
  36, 105, 126, 105, 35, 21, 0, 1;
  ...
		

References

  • R. Sedgewick and P. Flajolet, Analysis of Algorithms, Addison and Wesley, 1996, page 254 (first edition)

Crossrefs

Programs

  • Maple
    T := proc(n,k) if k>n then 0 elif k=n then 1 else (binomial(n+1,k)/(n+1))*sum(binomial(n+1-k,j)*binomial(n-k-j-1,j-1),j=1..floor((n-k)/2)) fi end: seq(seq(T(n,k),k=0..n),n=0..12);
    T := (n,k) -> (-1)^(n+k)*binomial(n,k)*hypergeom([-n+k,1/2],[2],4): seq(seq(simplify(T(n, k)), k=0..n), n=0..10); # Peter Luschny, Jul 27 2016
    # alternative Maple program:
    b:= proc(x, y, t) option remember; expand(`if`(x=0, 1,
          `if`(y>0, b(x-1, y-1, 0)*z^irem(t*y, 2), 0)+
          `if`(y (p-> seq(coeff(p, z, i), i=0..n))(b(2*n, 0$2)):
    seq(T(n), n=0..16);  # Alois P. Heinz, May 12 2017
  • Mathematica
    nn=10;cy = ( 1 + x - x y - ( -4x(1+x-x y) + (-1 -x + x y)^2)^(1/2))/(2(1+x-x y)); Drop[CoefficientList[Series[cy,{x,0,nn}],{x,y}],1]//Grid  (* Geoffrey Critzer, Feb 25 2013 *)
    Table[Which[k == n, 1, k > n, 0, True, (Binomial[n + 1, k]/(n + 1)) Sum[Binomial[n + 1 - k, j] Binomial[n - k - j - 1, j - 1], {j, Floor[(n - k)/2]}]], {n, 0, 11}, {k, 0, n}] // Flatten (* Michael De Vlieger, Jul 25 2016 *)

Formula

T(n, k) = [binomial(n+1, k)/(n+1)]*Sum_{j=1..floor((n-k)/2)} binomial(n+1-k, j)*binomial(n-k-j-1, j-1) for kn. G.f.=G=G(t, z) satisfies z(1+z-tz)G^2-(1+z-tz)G+1=0. T(n, k)=r(n-k)*binomial(n, k), where r(n)=A005043(n) are the Riordan numbers.
G.f.: 1/(1-xy-x^2/(1-x-xy-x^2/(1-x-xy-x^2/(1-x-xy-x^2/(1-... (continued fraction). - Paul Barry, Aug 03 2009
Sum_{k=0..n} T(n,k)*x^k = A126930(n), A005043(n), A000108(n), A007317(n), A064613(n), A104455(n) for x = -1,0,1,2,3,4 respectively. - Philippe Deléham, Dec 03 2009
Sum_{k=0..n} (-1)^(n-k)*T(n,k)*x^k = A168491(n), A099323(n+1), A001405(n), A005773(n+1), A001700(n), A026378(n+1), A005573(n), A122898(n) for x = -1, 0, 1, 2, 3, 4, 5, 6 respectively. - Philippe Deléham, Dec 03 2009
E.g.f.: e^(x+xy)*(Bessel_I(0,2x)-Bessel_I(1,2x)). - Paul Barry, Mar 10 2010
From Tom Copeland, Nov 06 2014: (Start)
O.g.f.: H(x,t) = {1-sqrt[1-4x/(1-(t-1)x)]}/2 (shifted index, as given in Copeland's comment, see comp. inverse there).
H(x,t)= x / [1-(C.+(t-1))x] = Sum_{n>=1} (C.+ (t-1))^(n-1)*x^n umbrally, e.g., (a.+b.)^2 = a_0*b_2 + 2 a_1*b1_+ a_0*b_2, where (C.)^n = C_n are the Catalan numbers (1,1,2,5,14,..) of A000108.
This shows directly that the lowering operator for the polynomials is D=d/dt, i.e., D p(n,t)= D(C. + (t-1))^n = n * (C. + (t-1))^(n-1) = n*p(n-1,t), so that the polynomials form an Appell sequence, and that p(n,0) gives a Motzkin sum, or Riordan, number A005043.
(End)
T(n,k) = (-1)^(n+k)*binomial(n,k)*hypergeom([k-n,1/2],[2],4). - Peter Luschny, Jul 27 2016

A300574 Coefficient of x^n in 1/((1-x)(1+x^3)(1-x^5)(1+x^7)(1-x^9)...).

Original entry on oeis.org

1, 1, 1, 0, 0, 1, 2, 1, 0, 0, 2, 2, 1, 0, 2, 3, 2, 0, 2, 4, 4, 0, 1, 4, 6, 2, 1, 4, 8, 4, 2, 4, 10, 6, 2, 3, 12, 10, 4, 2, 13, 14, 8, 2, 14, 18, 12, 2, 14, 22, 18, 3, 14, 26, 26, 6, 14, 29, 34, 10, 14, 32, 44, 16, 14, 34, 56, 26, 16, 34, 67, 38, 20, 34, 78, 52, 26
Offset: 0

Views

Author

Gus Wiseman, Mar 08 2018

Keywords

Comments

By Theorem 1 of Craig, the values a(n) in this list are known to be nonnegative. Combined with Theorem 2 of Seo and Yee, this shows that a(n) = |number of partitions of n into odd parts with an odd index minus the number of partitions of n into odd parts with an even index|. - William Craig, Dec 31 2021

References

  • Seunghyun Seo and Ae Ja Yee, Index of seaweed algebras and integer partitions, Electronic Journal of Combinatorics, 27:1 (2020), #P1.47. See Conjecture 1 and Theorem 2.

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[1/QPochhammer[x, -x^2], {x, 0, 100}], x]
    nmax = 100; CoefficientList[Series[Product[1/((1+x^(4*k-1))*(1-x^(4*k-3))), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jun 04 2019 *)

Formula

O.g.f.: Product_{n >= 0} 1/(1 - (-1)^n x^(2n+1)).
a(n) = Sum (-1)^k where the sum is over all integer partitions of n into odd parts and k is the number of parts not congruent to 1 modulo 4.
a(n) has average order Gamma(1/4) * exp(sqrt(n/3)*Pi/2) / (2^(9/4) * 3^(1/8) * Pi^(3/4) * n^(5/8)). - Vaclav Kotesovec, Jun 04 2019

A300866 Signed recurrence over binary strict trees: a(n) = 1 - Sum_{x + y = n, 0 < x < y < n} a(x) * a(y).

Original entry on oeis.org

1, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1, -1, 1, 1, -2, 3, -1, -3, 8, -8, 1, 14, -26, 22, 10, -59, 90, -52, -74, 238, -291, 80, 417, -930, 915, 124, -1991, 3483, -2533, -2148, 9011, -12596, 5754, 14350, -37975, 42735, -4046, -77924, 154374, -133903, -56529, 376844, -591197, 355941, 522978, -1706239
Offset: 0

Views

Author

Gus Wiseman, Mar 13 2018

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=1-Sum[a[k]*a[n-k],{k,1,(n-1)/2}];
    Array[a,40]
  • PARI
    seq(n)={my(v=vector(n)); for(n=1, n, v[n] = 1 - sum(k=1, (n-1)\2, v[k]*v[n-k])); concat([1], v)} \\ Andrew Howroyd, Aug 27 2018

A292137 G.f.: Im(1/(i*x; x)_inf), where (a; q)_inf is the q-Pochhammer symbol, i = sqrt(-1).

Original entry on oeis.org

0, 1, 1, 0, 0, 0, -1, -2, -2, -2, -2, -3, -3, -2, -2, -2, -1, 1, 2, 2, 4, 6, 7, 8, 10, 13, 14, 14, 15, 17, 17, 15, 15, 16, 14, 10, 8, 6, 1, -5, -10, -14, -21, -31, -38, -43, -53, -64, -71, -77, -86, -97, -104, -108, -115, -124, -127, -125, -127, -130, -125, -116
Offset: 0

Views

Author

Seiichi Manyama, Sep 09 2017

Keywords

Examples

			Product_{k>=1} 1/(1 - i*x^k) = 1 + (0+1i)*x + (-1+1i)*x^2 + (-1+0i)*x^3 + (-1+0i)*x^4 + (-1+0i)*x^5 + (-2-1i)*x^6 + (-1-2i)*x^7 + ...
		

Crossrefs

Programs

  • Maple
    N:= 100:
    S := convert(series( add( (-1)^n*x^(2*n+1)/(mul(1 - x^k,k = 1..2*n+1)), n = 0..N ), x, N+1 ), polynom):
    seq(coeff(S, x, n), n = 0..N); # Peter Bala, Jan 15 2021
  • Mathematica
    Im[CoefficientList[Series[1/QPochhammer[I*x, x], {x, 0, 100}], x]] (* Vaclav Kotesovec, Sep 17 2017 *)

Formula

1/(i*x; x)_inf is the g.f. for A292136(n) + i*a(n).
a(n) = Sum (-1)^((k - 1)/2) where the sum is over all integer partitions of n into an odd number of parts and k is the number of parts. - Gus Wiseman, Mar 08 2018
G.f.: Sum_{n >= 0} (-1)^n * x^(2*n+1)/Product_{k = 1..2*n+1} (1 - x^k). - Peter Bala, Jan 15 2021

A174295 Matrix inverse of A174294.

Original entry on oeis.org

1, -1, 1, -1, 0, 1, 0, -1, 0, 1, -1, 0, 0, 0, 1, 1, -2, -1, 1, 0, 1, -3, 2, 0, -2, 2, 0, 1, 6, -7, -3, 3, -3, 3, 0, 1, -15, 14, 3, -10, 7, -4, 4, 0, 1, 36, -37, -12, 19, -19, 12, -5, 5, 0, 1, -91, 90, 24, -54, 42, -30, 18, -6, 6, 0, 1, 232, -233, -67, 127, -115, 73, -43, 25, -7, 7
Offset: 0

Views

Author

Mats Granvik, Mar 15 2010

Keywords

Comments

First column is a signed version of A099323 with an additional leading 1.
First 5 rows as in A054525.

Examples

			Table begins:
  n\k|...0...1...2...3...4...5...6...7...8...9..10
  ---|--------------------------------------------
  0..|...1
  1..|..-1...1
  2..|..-1...0...1
  3..|...0..-1...0...1
  4..|..-1...0...0...0...1
  5..|...1..-2..-1...1...0...1
  6..|..-3...2...0..-2...2...0...1
  7..|...6..-7..-3...3..-3...3...0...1
  8..|.-15..14...3.-10...7..-4...4...0...1
  9..|..36.-37.-12..19.-19..12..-5...5...0...1
  10.|.-91..90..24.-54..42.-30..18..-6...6...0...1
		

Crossrefs

Programs

  • Mathematica
    t[n_, k_]:= t[n, k]= If[k<0 || k>n, 0, If[k==0 || k==n, 1, If[k==1, Mod[n, 2], t[n-1, k-1] +t[n-2, k-1] -t[n-1, k] -t[n-2, k] ]]]; (* t = A174294 *)
    M:= With[{m=30}, Table[t[n, k], {n,0,m}, {k,0,m}]];
    T:= Inverse[M];
    Table[T[[n+1, k+1]], {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Nov 25 2021 *)

Formula

Sum_{k=0..n} T(n, k) = A000007(n).
T(n, 0) = A174297(n).

A174297 First column of A174295.

Original entry on oeis.org

1, -1, -1, 0, -1, 1, -3, 6, -15, 36, -91, 232, -603, 1585, -4213, 11298, -30537, 83097, -227475, 625992, -1730787, 4805595, -13393689, 37458330, -105089229, 295673994, -834086421, 2358641376, -6684761125, 18985057351, -54022715451
Offset: 0

Views

Author

Mats Granvik, Mar 15 2010

Keywords

Comments

First 6 terms as in Mobius function A008683. Signed version of A099323 with an additional leading 1.

Crossrefs

Programs

  • Magma
    a:= func< n | n lt 2 select (-1)^n else (&+[(-1)^(k+1)*Binomial(n-2, k)*Catalan(k): k in [0..n-2]]) >;
    [a(n): n in [0..30]]; // G. C. Greubel, Nov 25 2021
    
  • Mathematica
    a[n_]:= a[n]= If[n<2, (-1)^n, Sum[(-1)^(j+1)*Binomial[n-2, j]*CatalanNumber[j], {j, 0, n-2}]]; Table[a[n], {n,0,40}] (* G. C. Greubel, Nov 25 2021 *)
  • Sage
    [1,-1]+[sum( (-1)^(j+1)*binomial(n-2,j)*catalan_number(j) for j in (0..n-2) ) for n in (2..40)] # G. C. Greubel, Nov 25 2021

Formula

a(n) = -(-3)^(n-3/2)*hypergeometric2F1([3/2, n-1],[2],4) for n > 2. - Mark van Hoeij, Jul 02 2010
a(n) = (-1)^n if n < 2 otherwise Sum_{j=0..n-2} (-1)^(j-1)*binomial(n-2, j)*Catalan(j). - G. C. Greubel, Nov 25 2021

A300863 Signed recurrence over enriched p-trees: a(n) = (-1)^(n - 1) + Sum_{y1 + ... + yk = n, y1 >= ... >= yk > 0, k > 1} a(y1) * ... * a(yk).

Original entry on oeis.org

1, 0, 2, 2, 6, 14, 34, 82, 214, 566, 1482, 4058, 10950, 30406, 83786, 235714, 658286, 1874254, 5293674, 15189810, 43312542, 125075238, 359185586, 1043712922, 3015569582, 8800146182, 25565402802, 74918274562, 218572345718, 642783954238, 1882606578002
Offset: 1

Views

Author

Gus Wiseman, Mar 13 2018

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=(-1)^(n-1)+Sum[Times@@a/@y,{y,Select[IntegerPartitions[n],Length[#]>1&]}];
    Array[a,40]

Formula

O.g.f.: (-1/(1+x) + Product 1/(1-a(n)x^n))/2.

A171651 Triangle T, read by rows : T(n,k) = A007318(n,k)*A005773(n+1-k).

Original entry on oeis.org

1, 2, 1, 5, 4, 1, 13, 15, 6, 1, 35, 52, 30, 8, 1, 96, 175, 130, 50, 10, 1, 267, 576, 525, 260, 75, 12, 1, 750, 1869, 2016, 1225, 455, 105, 14, 1, 2123, 6000, 7476, 5376, 2450, 728, 140, 16, 1, 6046, 19107, 27000, 22428, 12096, 4410, 1092, 180, 18, 1
Offset: 0

Views

Author

Philippe Deléham, Dec 14 2009

Keywords

Examples

			Triangle begins:
   1;
   2,   1;
   5,   4,  1;
  13,  15,  6, 1;
  35,  52, 30, 8, 1;
  ...
		

Crossrefs

Programs

  • Maple
    b:= proc(u, d, t) option remember; `if`(u=0 and d=0, 1/2,
          expand(`if`(u=0, 0, b(u-1, d, 2)*`if`(t=3, x, 1))
          +`if`(d=0, 0, b(u, d-1, `if`(t=2, 3, 1)))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n+1$2, 1)):
    seq(T(n), n=0..12);  # Alois P. Heinz, Apr 29 2015
    # second program:
    A171651:= (n, k)-> binomial(n,k)*add((-1)^(n-k-j)*binomial(n-k,j)*binomial(2*j+1,j+1),j=0..n-k): seq(print(seq(A171651(n, k), k=0..n)), n=0..9);  # Mélika Tebni, Dec 16 2023
  • Mathematica
    b[u_, d_, t_] := b[u, d, t] = If[u == 0 && d == 0, 1/2, Expand[If[u == 0, 0, b[u-1, d, 2]*If[t == 3, x, 1]] + If[d == 0, 0, b[u, d-1, If[t == 2, 3, 1]]]]];
    T[n_] := Function[p, Table[Coefficient[p, x, i], {i, 0, n}]][b[n+1, n+1, 1] ];
    Table[T[n], {n, 0, 12}] // Flatten (* Jean-François Alcover, May 21 2016, after Alois P. Heinz *)

Formula

Sum_{k, 0<=k<=n} T(n,k)*x^k = A168491(n), A099323(n), A001405(n), A005773(n+1), A001700(n), A026378(n+1), A005573(n), A122898(n) for x = -3, -2, -1, 0, 1, 2, 3, 4 respectively.
E.g.f. of column k: exp(x)*(BesselI(0,2*x)+BesselI(1,2*x))*x^k / k!. - Mélika Tebni, Dec 16 2023

Extensions

Corrected by Philippe Deléham, Dec 18 2009

A300864 Signed recurrence over strict trees: a(n) = -1 + Sum_{y1 + ... + yk = n, y1 > ... > yk > 0, k > 1} a(y1) * ... * a(yk).

Original entry on oeis.org

1, 1, 0, 1, 0, 0, 2, -1, 0, 4, -6, 6, 6, -24, 38, -17, -64, 188, -230, -6, 662, -1432, 1286, 1210, -6362, 10692, -5530, -18274, 57022, -74364, 174, 216703, -489544, 467860, 391258, -2256430, 3948206, -2234064, -6725362, 21920402, -29716570, 2095564, 84595798, -198418242, 197499846
Offset: 1

Views

Author

Gus Wiseman, Mar 13 2018

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=-1+Sum[Times@@a/@y,{y,Select[IntegerPartitions[n],Length[#]>1&&UnsameQ@@#&]}];
    -Array[a,40]
Showing 1-10 of 15 results. Next