A117671 -id:A117671 - cp's OEIS frontend

A025174 a(n) = binomial(3n-1, n-1).

0, 1, 5, 28, 165, 1001, 6188, 38760, 245157, 1562275, 10015005, 64512240, 417225900, 2707475148, 17620076360, 114955808528, 751616304549, 4923689695575, 32308782859535, 212327989773900, 1397281501935165, 9206478467454345, 60727722660586800, 400978991944396320

Offset: 0

Wouter Meeussen, Emeric Deutsch

Number of standard tableaux of shape (2n-1,n). Example: a(2)=5 because in the top row we can have 123, 124, 125, 134, or 135. - Emeric Deutsch, May 23 2004
Number of peaks in all generalized {(1,2),(1,-1)}-Dyck paths of length 3n.
Positive terms in this sequence are the numbers k such that k and 2k are consecutive terms in a row of Pascal's triangle. 1001 is the only k such that k, 2k, and 3k are consecutive terms in a row of Pascal's triangle. - J. Lowell, Mar 11 2023

			L.g.f.: L(x) = x + 5*x^2/2 + 28*x^3/3 + 165*x^4/4 + 1001*x^5/5 + 6188*x^6/6 + ...
where G(x) = exp(L(x)) satisfies G(x) = 1 + x*G(x)^3, and begins:
exp(L(x)) = 1 + x + 3*x^2 + 12*x^3 + 55*x^4 + 273*x^5 + ... + A001764(n)*x^n + ...

B. C. Berndt, Ramanujan's Notebooks Part I, Springer-Verlag, see Entry 14, Corollary 1, p. 71.

Robert Israel, Table of n, a(n) for n = 0..1190
Paul Barry, On the Central Antecedents of Integer (and Other) Sequences, J. Int. Seq., Vol. 23 (2020), Article 20.8.3.
D. Kruchinin and V. Kruchinin, A Generating Function for the Diagonal T2n,n in Triangles, Journal of Integer Sequence, Vol. 18 (2015), article 15.4.6.
W. Mlotkowski and K. A. Penson, Probability distributions with binomial moments, arXiv preprint arXiv:1309.0595 [math.PR], 2013.
Emanuele Munarini, Shifting Property for Riordan, Sheffer and Connection Constants Matrices, Journal of Integer Sequences, Vol. 20 (2017), Article 17.8.2.
Stepan Orevkov, Asymptotics of the number of lattice triangulations of rectangles of width 4 and 5, arXiv:2412.17065 [math.CO], 2024. See p. 16.
N. J. Wildberger and Dean Rubine, A Hyper-Catalan Series Solution to Polynomial Equations, and the Geode, Amer. Math. Monthly (2025). See section 12.

Cf. A001764 (binomial(3n,n)/(2n+1)), A117671 (C(3n+1,n+1)), A004319, A005809, A006013, A013698, A045721, A117671, A165817, A224274, A236194.

Cf. A000108, A001045, A005156, A224274.

[Binomial(3*n-1,n-1): n in [0..30]]; // Vincenzo Librandi, Nov 12 2014

with(combinat):seq(numbcomp(3*i,i), i=0..20); # Zerinvary Lajos, Jun 16 2007

Table[ GegenbauerC[ n, n, 1 ]/2, {n, 0, 24} ]
Join[{0},Table[Binomial[3n-1,n-1],{n,20}]] (* Harvey P. Dale, Oct 19 2022 *)
nmax=23; CoefficientList[Series[(2+HypergeometricPFQ[{1/3,2/3},{1/2,1},27x/4])/3-1,{x,0,nmax}],x]Range[0,nmax]! (* Stefano Spezia, Dec 31 2024 *)

vector(30, n, n--; binomial(3*n-1, n-1)) \\ Altug Alkan, Nov 04 2015

G.f.: z*g^2/(1-3*z*g^2), where g=g(z) is given by g=1+z*g^3, g(0)=1, that is, (in Maple command) g := 2*sin(arcsin(3*sqrt(3*z)/2)/3)/sqrt(3*z). - Emeric Deutsch, May 22 2003
a(n) = Sum_{k=0..n} ((3k+1)/(2n+k+1))C(3n, 2n+k)*A001045(k). - Paul Barry, Oct 07 2005
Hankel transform of a(n+1) is A005156(n+1). - Paul Barry, Apr 14 2008
G.f.: x*B'(x)/B(x) where B(x) is the g.f. of A001764. - Vladimir Kruchinin Feb 03 2013
D-finite with recurrence: 2*n*(2*n-1)*a(n) -3*(3*n-1)*(3*n-2)*a(n-1)=0. - R. J. Mathar, Feb 05 2013
Logarithmic derivative of A001764; g.f. of A001764 satisfies G(x) = 1 + x*G(x)^3. - Paul D. Hanna, Jul 14 2013
G.f.: (2*cos((1/3)*arcsin((3/2)*sqrt(3*x)))-sqrt(4-27*x))/(3*sqrt(4-27*x)). - Emanuele Munarini, Oct 14 2014
a(n) = Sum_{k=1..n} binomial(n-1,n-k)*binomial(2*n,n-k). - Vladimir Kruchinin, Nov 12 2014
a(n) = [x^n] C(x)^n for n >= 1, where C(x) = (1 - sqrt(1 - 4*x))/(2*x) is the generating function for A000108 (Ramanujan). - Peter Bala, Jun 24 2015
From Peter Bala, Nov 04 2015: (Start)
Without the initial term 0, the o.g.f. equals f(x)*g(x)^2, where f(x) is the o.g.f. for A005809 and g(x) is the o.g.f. for A001764. g(x)^2 is the o.g..f for A006013. More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(3*n + k,n). Cf. A045721 (k = 1), A004319 (k = 3), A236194 (k = 4), A013698 (k = 5), A165817 (k = -1), A117671 (k = -2). (End)
G.f.: ( 2F1(1/3,2/3;1/2;27*x/4)-1)/3. - R. J. Mathar, Jan 27 2020
O.g.f. without the initial term 0, in the form g(x)=(2*cos(arcsin((3*sqrt(3)*sqrt(x))/2)/3)/sqrt(4-27*x)-1)/(3*x), satisfies the following algebraic equation: 1+(9*x-1)*g(x)+x*(27*x-4)*g(x)^2+x^2*(27*x-4)*g(x)^3=0. - Karol A. Penson, Oct 11 2021
O.g.f. equals f(x)/(1 - 2*f(x)), where f(x) = series reversion (x/(1 + x)^3) = x + 3*x^2 + 12*x^3 + 55*x^4 + 273*x^5 + ... is the o.g.f. of A001764 with the initial term omitted. Cf. A224274. - Peter Bala, Feb 03 2022
Right-hand side of the identities (1/2)*Sum_{k = 0..n} (-1)^(n+k)*C(x*n,n-k)*C((x+2)*n+k-1,k) = C(3*n-1,n-1) and (1/3)*Sum_{k = 0..n} (-1)^k* C(x*n,n-k)*C((x-3)*n+k-1,k) = C(3*n-1,n-1), both valid for n >= 1 and x arbitrary. - Peter Bala, Feb 28 2022
a(n) ~ 2^(-2*n)*3^(3*n)/(2*sqrt(3*n*Pi)). - Stefano Spezia, Apr 25 2024
a(n) = Sum_{k = 0..n-1} binomial(2*n+k-1, k) = Sum_{k = 0..n-1} (-1)^(n+k+1)* binomial(3*n, k). - Peter Bala, Jul 21 2024
E.g.f.: (2 + hypergeom([1/3, 2/3], [1/2, 1], 27*x/4))/3 - 1. - Stefano Spezia, Dec 31 2024

A045721 a(n) = binomial(3*n+1,n).

Original entry on oeis.org

1, 4, 21, 120, 715, 4368, 27132, 170544, 1081575, 6906900, 44352165, 286097760, 1852482996, 12033222880, 78378960360, 511738760544, 3348108992991, 21945588357420, 144079707346575, 947309492837400, 6236646703759395

Offset: 0

Emeric Deutsch

Number of leaves in all noncrossing rooted trees on n nodes on a circle.
Number of standard tableaux of shape (n-1,1^(2n-3)). - Emeric Deutsch, May 25 2004
a(n) = number of Dyck (2n-3)-paths with exactly one descent of odd length. For example, a(3) counts all 5 Dyck 3-paths except UDUDUD. - David Callan, Jul 25 2005
a(n+2) gives row sums of A119301. - Paul Barry, May 13 2006
a(n) is the number of paths avoiding UU from (0,0) to (3n,n) and taking steps from {U,H}. - Shanzhen Gao, Apr 15 2010
Central coefficients of triangle A078812. - Vladimir Kruchinin, May 10 2012
Row sums of A252501. - L. Edson Jeffery, Dec 18 2014

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Milan Janjic, Two Enumerative Functions
D. Kruchinin and V. Kruchinin, A Method for Obtaining Generating Function for Central Coefficients of Triangles, Journal of Integer Sequence, Vol. 15 (2012), article 12.9.3.
W. Mlotkowski and K. A. Penson, Probability distributions with binomial moments, arXiv preprint arXiv:1309.0595 [math.PR], 2013.
Index entries for sequences related to rooted trees

Cf. A252501, A263134 (partial sums), A001764, A004319, A005809, A006013, A013698, A025174, A117671, A165817, A236194.

[Binomial(3*n+1, n): n in [0..20]]; // Vincenzo Librandi, Aug 07 2014

[seq( binomial(3*n+1,n),n=0..40)]; # N. J. A. Sloane, Jun 09 2007

Table[Binomial[3 n + 1, n], {n, 0, 20}] (* Vincenzo Librandi, Aug 07 2014 *)

a(n)=binomial(3*n+1,n) \\ Charles R Greathouse IV, Mar 18 2014

a(n) is asymptotic to c/sqrt(n)*(27/4)^n with c=0.73... - Benoit Cloitre, Jan 27 2003
G.f.: gz^2/(1-3zg^2), where g=g(z) is given by g=1+zg^3, g(0)=1, i.e. (in Maple command) g := 2*sin(arcsin(3*sqrt(3*z)/2)/3)/sqrt(3*z). - Emeric Deutsch, May 22 2003
a(n+2) = C(3n+1,n) = Sum_{k=0..n} C(3n-k,n-k). - Paul Barry, May 13 2006
a(n+2) = C(3n+1,2n+1) = A078812(2n,n). - Paul Barry, Nov 09 2006
G.f.: A(x)=(2*cos(asin((3^(3/2)*sqrt(x))/2)/3)* sin(asin((3^(3/2)* sqrt(x))/2)/3))/(sqrt(3)*sqrt(1-(27*x)/4)*sqrt(x)). - Vladimir Kruchinin, Jun 10 2012
From Peter Luschny, Sep 04 2012: (Start)
O.g.f.: hypergeometric2F1([2/3, 4/3], [3/2], x*27/4).
a(n) = (n+1)*hypergeometric2F1([-2*n, -n], [2], 1).  (End)
D-finite with recurrence 2*n*(2*n+1)*a(n) - 3*(3*n-1)*(3*n+1)*a(n-1) = 0. - R. J. Mathar, Feb 05 2013
a(n) = Sum_{r=0..n} C(n,r) * C(2*n+1,r). - J. M. Bergot, Mar 18 2014
From Peter Bala, Nov 04 2015: (Start)
a(n) = Sum_{k = 0..n} 1/(2*k + 1)*binomial(3*n - 3*k,n - k)*binomial(3*k, k).
O.g.f. equals f(x)*g(x), where f(x) is the o.g.f. for A005809 and g(x) is the o.g.f. for A001764. More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(3*n + k,n). Cf. A025174 (k = 2), A004319 (k = 3), A236194 (k = 4), A013698 (k = 5), A165817 (k = -1), A117671 (k = -2). (End)
a(n) = [x^n] 1/(1 - x)^(2*(n+1)). - Ilya Gutkovskiy, Oct 10 2017
O.g.f.: (i/24)*((4*sqrt(4 - 27*z) + 12*i*sqrt(3)*sqrt(z))^(2/3) - (4*sqrt(4 - 27*z) - 12*i*sqrt(3)*sqrt(z))^(2/3))*sqrt(3)*8^(1/3)*sqrt(4 - 27*z)/(sqrt(z)*(-4 + 27*z)), where i = sqrt(-1). - Karol A. Penson, Dec 13 2023
a(n-1) = (1/(4*n))*binomial(2*n, n)^2 * (1 - 3*((n - 1)/(n + 1))^3 + 5*((n - 1)*(n - 2)/((n + 1)*(n + 2)))^3 - 7*((n - 1)*(n - 2)*(n - 3)/((n + 1)*(n + 2)*(n + 3)))^3 + - ...) for n >= 1. Cf. A112029. - Peter Bala, Aug 08 2024

Simpler definition from Ira M. Gessel, May 26 2007. This change means that most of the offsets in the comments will now need to be changed too.

A165817 Number of compositions (= ordered integer partitions) of n into 2n parts.

Original entry on oeis.org

1, 2, 10, 56, 330, 2002, 12376, 77520, 490314, 3124550, 20030010, 129024480, 834451800, 5414950296, 35240152720, 229911617056, 1503232609098, 9847379391150, 64617565719070, 424655979547800, 2794563003870330, 18412956934908690, 121455445321173600

Offset: 0

Thomas Wieder, Sep 29 2009

Number of ways to put n indistinguishable balls into 2*n distinguishable boxes.

Number of rankings of n unlabeled elements for 2*n levels.

			Let [1,1,1], [1,2] and [3] be the integer partitions of n=3.
Then [0,0,0,1,1,1], [0,0,0,0,1,2] and [0,0,0,0,0,3] are the corresponding partitions occupying 2*n = 6 positions.
We have to take into account the multiplicities of the parts including the multiplicities of the zeros.
Then
  [0,0,0,1,1,1] --> 6!/(3!*3!) = 20
  [0,0,0,0,1,2] --> 6!/(4!*1!*1!) = 30
  [0,0,0,0,0,3] --> 6!/(5!*1!) = 6
and thus a(3) = 20+30+6=56.
a(2)=10, since we have 10 ordered partitions of n=2 where the parts are distributed over 2*n=4 boxes:
  [0, 0, 0, 2]
  [0, 0, 1, 1]
  [0, 0, 2, 0]
  [0, 1, 0, 1]
  [0, 1, 1, 0]
  [0, 2, 0, 0]
  [1, 0, 0, 1]
  [1, 0, 1, 0]
  [1, 1, 0, 0]
  [2, 0, 0, 0].

Alois P. Heinz, Table of n, a(n) for n = 0..400
W. Mlotkowski and K. A. Penson, Probability distributions with binomial moments, arXiv preprint arXiv:1309.0595 [math.PR], 2013.

Cf. A000079, A001700, A059481, A081204, A001764, A004319, A006013, A005809, A013698, A025174, A045721, A117671, A236194.

[Binomial(3*n-1, n): n in [0..30]]; // Vincenzo Librandi, Aug 07 2014

for n from 0 to 16 do
a[n] := 9*sqrt(3)*GAMMA(n+5/3)*GAMMA(n+4/3)*27^n/(Pi*GAMMA(2*n+3))
end do;

Table[Binomial[3 n - 1, n], {n, 0, 20}] (* Vincenzo Librandi, Aug 07 2014 *)

vector(30, n, n--; binomial(3*n-1, n)) \\ Altug Alkan, Nov 04 2015

from math import comb
def A165817(n): return comb(3*n-1,n) if n else 1 # Chai Wah Wu, Oct 11 2023

def A165817(n):
    return rising_factorial(2*n,n)/falling_factorial(n,n)
[A165817(n) for n in (0..22)]   # Peter Luschny, Nov 21 2012

a(n) = 9*sqrt(3)*GAMMA(n+5/3)*GAMMA(n+4/3)*27^n/(Pi*GAMMA(2*n+3)).
a(n) = binomial(3*n-1, n);
Let denote P(n) = the number of integer partitions of n,
p(i) = the number of parts of the i-th partition of n,
d(i) = the number of different parts of the i-th partition of n,
m(i,j) = multiplicity of the j-th part of the i-th partition of n.
Then one has:
a(n) = Sum_{i=1..P(n)} (2*n)!/((2*n-p(i))!*(Prod_{j=1..d(i)} m(i,j)!)).
a(n) = rf(2*n,n)/ff(n,n), where rf is the rising factorial and ff the falling factorial. - Peter Luschny, Nov 21 2012
G.f.: A(x) = x*B'(x)/B(x), where B(x) satisfies B(x)^3-2*B(x)^2+B(x)=x, B(x)=A006013(x). - Vladimir Kruchinin, Feb 06 2013
G.f.: A(x) = sqrt(3*x)*cot(asin((3^(3/2)*sqrt(x))/2)/3)/(sqrt(4-27*x)). - Vladimir Kruchinin, Mar 18 2015
a(n) = Sum_{k=0..n} binomial(n-1,n-k)*binomial(2*n,k). - Vladimir Kruchinin, Oct 06 2015
From Peter Bala, Nov 04 2015: (Start)
The o.g.f. equals f(x)/g(x), where f(x) is the o.g.f. for A005809 and g(x) is the o.g.f. for A001764. More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(3*n + k,n). Cf. A045721 (k = 1), A025174 (k = 2), A004319 (k = 3), A236194 (k = 4), A013698 (k = 5) and A117671 (k = -2). (End)
a(n) = [x^n] 1/(1 - x)^(2*n). - Ilya Gutkovskiy, Oct 03 2017
a(n) = A059481(2n,n). - Alois P. Heinz, Oct 17 2022
From Peter Bala, Feb 14 2024: (Start)
a(n) = (-1)^n * binomial(-2*n, n).
a(n) = hypergeom([1 - 2*n, -n], [1], 1).
The g.f. A(x) satisfies A(x/(1 + x)^3) = 1/(1 - 2*x).
Sum_{n >= 0} a(n)/9^n = (1 + 4*cos(Pi/9))/3.
Sum_{n >= 0} a(n)/27^n = (3 + 4*sqrt(3)*cos(Pi/18))/9.
Sum_{n >= 0} a(n)*(2/27)^n = (2 + sqrt(3))/3. (End)
From Peter Bala, Sep 16 2024: (Start)
a(n) = Sum_{k = 0..n} binomial(n+k-1, k)*binomial(2*n-k-1, n-k).
More generally, a(n) = Sum_{k = 0..n} (-1)^k*binomial(x*n, k)*binomial((x+3)*n-k-1, n-k) for arbitrary x.
a(n) = (2/3) * Sum_{k = 0..n} (-1)^k*binomial(x*n+k-1, k)*binomial((x+3)*n, n-k) for n >= 1 and arbitrary x. (End)
G.f.: 1/(3-2*g) where g = 1+x*g^3 is the g.f. of A001764. - Seiichi Manyama, Aug 17 2025

a(0) prepended and more terms from Alois P. Heinz, Apr 04 2012

A004319 a(n) = binomial(3*n, n - 1).

Original entry on oeis.org

1, 6, 36, 220, 1365, 8568, 54264, 346104, 2220075, 14307150, 92561040, 600805296, 3910797436, 25518731280, 166871334960, 1093260079344, 7174519270695, 47153358767970, 310325523515700, 2044802197953900, 13488561475572645, 89067326568860640, 588671286046028640

Offset: 1

N. J. A. Sloane

Milton Abramowitz and Irene A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.

Seiichi Manyama, Table of n, a(n) for n = 1..1000
Milan Janjic, Two Enumerative Functions.
Milton Abramowitz and Irene A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Emanuele Munarini, Shifting Property for Riordan, Sheffer and Connection Constants Matrices, Journal of Integer Sequences, Vol. 20 (2017), Article 17.8.2.
Franck Ramaharo, Statistics on some classes of knot shadows, arXiv:1802.07701 [math.CO], 2018.

Cf. A045721 (k=1), A025174 (k=2), A236194 (k=4), A013698 (k=5), A165817 (k=-1), A117671 (k=-2).

Cf. A000139, A001764, A005809, A006013, A236194, A001791, A004331, A004343, A004356, A004369, A004382.

A004319 := proc(n)
binomial(3*n,n-1);
end proc: # R. J. Mathar, Aug 10 2015

Table[Binomial[3n, n - 1], {n, 20}] (* Harvey P. Dale, Sep 21 2011 *)

a(n):=sum((binomial(3*i-1,2*i-1)*binomial(3*n-3*i-3,2*n-2*i-2))/(2*n-2*i-1),i,1,n-1)/2; /* Vladimir Kruchinin, May 15 2013 */

vector(30, n, binomial(3*n, n-1)) \\ Altug Alkan, Nov 04 2015

G.f.: (g-1)/(1-3*z*g^2), where g = g(z) is given by g = 1 + z*g^3, g(0) = 1, i.e. (in Maple notation), g := 2*sin(arcsin(3*sqrt(3*z)/2)/3)/sqrt(3*z). - Emeric Deutsch, May 22 2003
a(n) = Sum_{i=0..n-1} binomial(i+2*n, i). - Ralf Stephan, Jun 03 2005
D-finite with recurrence -2*(2*n+1)*(n-1)*a(n) + 3*(3*n-1)*(3*n-2)*a(n-1) = 0. - R. J. Mathar, Feb 05 2013
a(n) = (1/2) * Sum_{i=1..n-1} binomial(3*i - 1, 2*i - 1)*binomial(3*n - 3*i - 3, 2*n - 2*i - 2)/(2*n - 2*i - 1). - Vladimir Kruchinin, May 15 2013
G.f.: x*hypergeom2F1(5/3, 4/3; 5/2; 27x/4). - R. J. Mathar, Aug 10 2015
a(n) = n*A001764(n). - R. J. Mathar, Aug 10 2015
From Peter Bala, Nov 04 2015: (Start)
With offset 0, the o.g.f. equals f(x)*g(x)^3, where f(x) is the o.g.f. for A005809 and g(x) is the o.g.f. for A001764. More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(3*n + k, n). See the cross-references. (End)
G.f.: cos(t)/(2*sqrt(1 - (27*x)/4)) - sin(t)/(sqrt(3)*sqrt(x)), where t = arcsin((sqrt(27*x))/2)/3. - Vladimir Kruchinin, May 13 2016
a(n) = [x^(2*n+1)] 1/(1 - x)^n. - Ilya Gutkovskiy, Oct 10 2017
a(n) = binomial(n+1, 2) * A000139(n). - F. Chapoton, Feb 23 2024

A013698 a(n) = binomial(3*n+2, n-1).

Original entry on oeis.org

1, 8, 55, 364, 2380, 15504, 100947, 657800, 4292145, 28048800, 183579396, 1203322288, 7898654920, 51915526432, 341643774795, 2250829575120, 14844575908435, 97997533741800, 647520696018735, 4282083008118300

Offset: 1

Joachim.Rosenthal(AT)nd.edu (Joachim Rosenthal), Emeric Deutsch

Degree of variety K_{2,n}^1. Also number of double-rises (or odd-level peaks) in all generalized {(1,2),(1,-1)}-Dyck paths of length 3(n+1).
Number of dissections of a convex (2n+2)-gon by n-2 noncrossing diagonals into (2j+2)-gons, 1<=j<=n-1.
a(n) is the number of lattice paths from (0,0) to (3n+1,n-1) avoiding two consecutive up-steps. - Shanzhen Gao, Apr 20 2010

Robert Israel, Table of n, a(n) for n = 1..1000
M. S. Ravi et al., Dynamic Pole Assignment and Schubert Calculus, SIAM J. Control Optimization, 34 (1996), 813-832, esp. p. 825.
Daniel W. Stasiuk, An Enumeration Problem for Sequences of n-ary Trees Arising from Algebraic Operads, Master's Thesis, University of Saskatchewan-Saskatoon (2018).

Cf. A013699 (q=2), A013700 (q=3), A013701 (q=4), A013702 (q=5).

A column of triangle A102537. Cf. A001764, A004319, A005809, A006013, A025174, A045721, A117671, A165817, A236194.

List([1..25], n-> Binomial(3*n+2, n-1)) # G. C. Greubel, Mar 21 2019

[Binomial(3*n+2, n-1): n in [1..25]]; // Vincenzo Librandi, Aug 10 2015

seq(binomial(3*n+2,n-1), n=0..30); # Robert Israel, Aug 09 2015

Table[Binomial[3*n+2, n-1], {n, 25}] (* Arkadiusz Wesolowski, Apr 02 2012 *)

first(m)=vector(m,n,binomial(3*n+2, n-1)); /* Anders Hellström, Aug 09 2015 */

[binomial(3*n+2, n-1) for n in (1..25)] # G. C. Greubel, Mar 21 2019

G.f.: g/((g-1)^3*(3*g-1)) where g*(1-g)^2 = x. - Mark van Hoeij, Nov 09 2011
a(n) = Sum_{k=0..n-1} binomial(2*n+k+2,k). - Arkadiusz Wesolowski, Apr 02 2012
D-finite with recurrence 2*(2*n+3)*(n+1)*a(n) -n*(67*n+34)*a(n-1) +30*(3*n-1)*(3*n-2)*a(n-2)=0. - R. J. Mathar, Feb 05 2013
a(n+1) = (3*n+5)*(3*n+4)*(3*n+3)*a(n)/((2*n+5)*(2*n+4)*n). - Robert Israel, Aug 09 2015
With offset 0, the o.g.f. equals f(x)*g(x)^5, where f(x) is the o.g.f. for A005809 and g(x) is the o.g.f. for A001764. More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(3*n + k,n). Cf. A045721 (k = 1), A025174 (k = 2), A004319 (k = 3), A236194 (k = 4), A165817 (k = -1), A117671 (k = -2). - Peter Bala, Nov 04 2015

A118113 Even Fibbinary numbers + 1; also 2*Fibbinary(n) + 1.

Original entry on oeis.org

1, 3, 5, 9, 11, 17, 19, 21, 33, 35, 37, 41, 43, 65, 67, 69, 73, 75, 81, 83, 85, 129, 131, 133, 137, 139, 145, 147, 149, 161, 163, 165, 169, 171, 257, 259, 261, 265, 267, 273, 275, 277, 289, 291, 293, 297, 299, 321, 323, 325, 329, 331, 337, 339, 341, 513, 515, 517

Offset: 0

Labos Elemer, Apr 13 2006

m for which binomial(3*m-2,m) (see A117671) is odd, since by Kummer's theorem that happens exactly when the binary expansions of m and 2*m-2 have no 1 bit at the same position in each, and so m odd and no 11 bit pairs except optionally the least significant 2 bits. - Kevin Ryde, Jun 14 2025

Alois P. Heinz, Table of n, a(n) for n = 0..1000

Cf. A000108, A118112, A022341.
Cf. A003714 (Fibbinary numbers), A022340 (even Fibbinary numbers).
Cf. A117671, A263190, A171791, A263075.

F:= combinat[fibonacci]:
b:= proc(n) local j;
      if n=0 then 0
    else for j from 2 while F(j+1)<=n do od;
         b(n-F(j))+2^(j-2)
      fi
    end:
a:= n-> 2*b(n)+1:
seq(a(n), n=0..70);  # Alois P. Heinz, Aug 03 2012

Select[Table[Mod[Binomial[3*k,k], k+1], {k,1200}], #>0&]

a(n) = A022340(n) + 1.
a(n) = 2*A003714(n) + 1.
Solutions to {x : binomial(3x,x) mod (x+1) != 0 } are given in A022341. The corresponding values of binomial(3x,x) mod (x+1) are given here.

New definition from T. D. Noe, Dec 19 2006

A236194 a(n) = binomial(3n+1, n-1).

Original entry on oeis.org

1, 7, 45, 286, 1820, 11628, 74613, 480700, 3108105, 20160075, 131128140, 854992152, 5586853480, 36576848168, 239877544005, 1575580702584, 10363194502115, 68248282427325, 449972009097765, 2969831763694950, 19619725782651120, 129728497393775280

Offset: 1

Bruno Berselli, Jan 20 2014

This sequence is related to A006013 by a(n)/n = A006013(n)/2.

Bruno Berselli, Table of n, a(n) for n = 1..100
Emanuele Munarini, Shifting Property for Riordan, Sheffer and Connection Constants Matrices, Journal of Integer Sequences, Vol. 20 (2017), Article 17.8.2.

Cf. A006013; A025174: C(3n-1, n-1); A117671: C(3n+1, n+1).
Second column of the triangle A159841.
Third column of the triangle A119301.
Cf. A001764, A004319, A005809, A013698, A045721, A165817.
Cf. A045721, A117671.

Magma
```
[Binomial(3*n+1,n-1): n in [1..30]];
```
Mathematica
```
Table[Binomial[3n+1, n-1], {n, 30}]
```

makelist(binomial(3*n+4,n),n,0,40); /* Emanuele Munarini, Oct 14 2014 */

vector(30, n, binomial(3*n+1, n-1)) \\ Altug Alkan, Nov 04 2015

[binomial(3*n+1,n-1) for n in range(1,31)] # G. C. Greubel, Nov 09 2022

G.f.: (sqrt(4-27*x)*cos((2/3)*arcsin((3/2)*sqrt(3*x))) + sqrt(3*x)*sin((2/3)*arcsin((3/2)*sqrt(3*x))) - sqrt(4-27*x))/(3*sqrt(4-27*x)*x^2). - Emanuele Munarini, Oct 14 2014
From Peter Bala, Nov 04 2015: (Start)
With offset 0, the o.g.f. equals f(x)*g(x)^4, where f(x) is the o.g.f. for A005809 and g(x) is the o.g.f. for A001764.
More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(3*n + k,n). Cf. A045721 (k = 1), A025174 (k = 2), A004319 (k = 3), A013698 (k = 5), A165817 (k = -1), A117671 (k = -2). (End)
a(n) = [x^n] x/(1 - x)^(2*n+3). - Ilya Gutkovskiy, Oct 10 2017
From Karol A. Penson, Mar 02 2024: (Start)
G.f.: ((sqrt(3)*sqrt(x)*i + sqrt(4 - 27*x))*(4*sqrt(4 - 27*x) - 12*i*sqrt(3)*sqrt(x))^(2/3) + (-sqrt(3)*sqrt(x)*i + sqrt(4 - 27*x))*(4*sqrt(4 - 27*x) + 12*i*sqrt(3)*sqrt(x))^(2/3) - 8*sqrt(4 - 27*x))/(24*sqrt(4 - 27*x)*x^2), where i is the imaginary unit, i=sqrt(-1).
G.f.: hypergeometric3F2([5/3,2,7/3],[5/2,3],27*x/4).
G.f. = G satisfies the algebraic equation: 1 + (7*z-1)*G + (27*z-4)*z^2*G^2 + (27*z-4)*z^4*G^3 = 0. (End)

A258993 Triangle read by rows: T(n,k) = binomial(n+k,n-k), k = 0..n-1.

Original entry on oeis.org

1, 1, 3, 1, 6, 5, 1, 10, 15, 7, 1, 15, 35, 28, 9, 1, 21, 70, 84, 45, 11, 1, 28, 126, 210, 165, 66, 13, 1, 36, 210, 462, 495, 286, 91, 15, 1, 45, 330, 924, 1287, 1001, 455, 120, 17, 1, 55, 495, 1716, 3003, 3003, 1820, 680, 153, 19, 1, 66, 715, 3003, 6435, 8008, 6188, 3060, 969, 190, 21

Offset: 1

Reinhard Zumkeller, Jun 22 2015

T(n,k) = A085478(n,k) = A007318(A094727(n),A004736(k)), k = 0..n-1;
rounded(T(n,k)/(2*k+1)) = A258708(n,k);
rounded(sum(T(n,k)/(2*k+1)): k = 0..n-1) = A000967(n).

			.  n\k |  0  1    2    3     4     5     6     7    8    9  10 11
. -----+-----------------------------------------------------------
.   1  |  1
.   2  |  1  3
.   3  |  1  6    5
.   4  |  1 10   15    7
.   5  |  1 15   35   28     9
.   6  |  1 21   70   84    45    11
.   7  |  1 28  126  210   165    66    13
.   8  |  1 36  210  462   495   286    91    15
.   9  |  1 45  330  924  1287  1001   455   120   17
.  10  |  1 55  495 1716  3003  3003  1820   680  153   19
.  11  |  1 66  715 3003  6435  8008  6188  3060  969  190  21
.  12  |  1 78 1001 5005 12870 19448 18564 11628 4845 1330 231 23  .

Reinhard Zumkeller, Rows n = 1..125 of triangle, flattened

If a diagonal of 1's is added on the right, this becomes A085478.
Essentially the same as A143858.
Cf. A007318, A004736, A094727.
Cf. A027941 (row sums), A117671 (central terms), A143858, A000967, A258708.
T(n,k): A000217 (k=1), A000332 (k=2), A000579 (k=3), A000581 (k=4), A001287 (k=5), A010965 (k=6), A010967 (k=7), A010969 (k=8), A010971 (k=9), A010973 (k=10), A010975 (k=11), A010977 (k=12), A010979 (k=13), A010981 (k=14), A010983 (k=15), A010985 (k=16), A010987 (k=17), A010989 (k=18), A010991 (k=19), A010993 (k=20), A010995 (k=21), A010997 (k=22), A010999 (k=23), A011001 (k=24), A017714 (k=25), A017716 (k=26), A017718 (k=27), A017720 (k=28), A017722 (k=29), A017724 (k=30), A017726 (k=31), A017728 (k=32), A017730 (k=33), A017732 (k=34), A017734 (k=35), A017736 (k=36), A017738 (k=37), A017740 (k=38), A017742 (k=39), A017744 (k=40), A017746 (k=41), A017748 (k=42), A017750 (k=43), A017752 (k=44), A017754 (k=45), A017756 (k=46), A017758 (k=47), A017760 (k=48), A017762 (k=49), A017764 (k=50).
T(n+k,n): A005408 (k=1), A000384 (k=2), A000447 (k=3), A053134 (k=4), A002299 (k=5), A053135 (k=6), A053136 (k=7), A053137 (k=8), A053138 (k=9), A196789 (k=10).
Cf. A165253.

Flat(List([1..12], n-> List([0..n-1], k-> Binomial(n+k,n-k) ))); # G. C. Greubel, Aug 01 2019

a258993 n k = a258993_tabl !! (n-1) !! k
a258993_row n = a258993_tabl !! (n-1)
a258993_tabl = zipWith (zipWith a007318) a094727_tabl a004736_tabl

[Binomial(n+k,n-k): k in [0..n-1], n in [1..12]]; // G. C. Greubel, Aug 01 2019

Table[Binomial[n+k,n-k], {n,1,12}, {k,0,n-1}]//Flatten (* G. C. Greubel, Aug 01 2019 *)

T(n,k) = binomial(n+k,n-k);
for(n=1, 12, for(k=0,n-1, print1(T(n,k), ", "))) \\ G. C. Greubel, Aug 01 2019

[[binomial(n+k,n-k) for k in (0..n-1)] for n in (1..12)] # G. C. Greubel, Aug 01 2019

T(n,k) = A085478(n,k) = A007318(A094727(n),A004736(k)), k = 0..n-1;
rounded(T(n,k)/(2*k+1)) = A258708(n,k);
rounded(sum(T(n,k)/(2*k+1)): k = 0..n-1) = A000967(n).

A331322 a(n) = (3*n + 1)!/(n!)^3.

Original entry on oeis.org

1, 24, 630, 16800, 450450, 12108096, 325909584, 8779605120, 236637794250, 6380456082000, 172080900531540, 4641917845743360, 125235075213284400, 3379123922914656000, 91184624634161304000, 2460769070127233057280, 66411927755894739034170, 1792432652235221330334000

Offset: 0

Peter Luschny, Jan 18 2020

nonn

Diagonal of the rational function 1 / (1 - x - y - z)^2. - Ilya Gutkovskiy, Apr 23 2025

G. C. Greubel, Table of n, a(n) for n = 0..500

Cf. A000108, A006480, A117671, A331431.

[(n+1)^2*Binomial(3*n+1,n+1)*Catalan(n): n in [0..25]]; // G. C. Greubel, Mar 22 2022

a := n -> (3*n+1)!/(n!)^3: seq(a(n), n=0..17); # Or:
hypergeom([2/3, 4/3], [1], 27*x): ser := series(%, x, 20):
seq(coeff(%, x, n), n=0..17); # Or:
a := proc(n) option remember; if n=0 then 1 else 3*(9 - n^(-2))*a(n-1) fi end:
# 4th Maple program:
W:=proc(x)sqrt(3)*MeijerG([[], [0, 0]], [[1/3, -1/3], []], x/27)/(18*Pi);end;
a:=proc(n) round(evalf[32](int(x^n*W(x),x=0..27)));end;
seq(a(n),n=0..17);
# Karol A. Penson, Jul 28 2023

Table[(3*n+1)*Binomial[3*n,n]*Binomial[2*n,n], {n,0,25}] (* G. C. Greubel, Mar 22 2022 *)

[(3*n+1)*binomial(2*n,n)*binomial(3*n,n) for n in (0..25)] # G. C. Greubel, Mar 22 2022

a(n) = [x^n] hypergeom([2/3, 4/3], [1], 27*x).
a(n) = 3*(9 - n^(-2))*a(n-1) for n > 0.
a(n) = (-1)^n*A331431(2*n, n).
a(n) = (n+1)^2*A117671(n)*A000108(n). - G. C. Greubel, Mar 22 2022
From Karol A. Penson, Jul 28 2023: (Start)
a(n) = Integral_{x=0..27} x^n*W(x) dx, where the weight function W(x) is defined on (0, 27) and it can be expressed with the Meijer G-function MeijerG as: W(x) = (sqrt(3)/(18*Pi))*MeijerG([[],[0,0]],[[-1/3,1/3],[]],x/27). The function W(x) is positive on its support (0, 27), is singular at x=0, and decreases monotonically to zero at x = 27.
The function W(x) is unique as it is the solution of the Hausdorff moment problem with the moments a(n). Due to the presence of two equal parameters (0,0) in MeijerG, it is not certain if W(x) can be represented by other known special functions. (End)
From  Peter Bala, Oct 10 2024: (Start)
a(n) = (3*n + 1)*A006480(n).
a(n-1) = 1/(8*n^3) * Sum_{k = 0..2*n} (-1)^(n+k) * k*(2*n-k)^3 * binomial(2*n, k)^3 for n >= 1.
a(n-1) = 1/(4*n^2) * Sum_{k = 0..2*n-1} (-1)^(n+k) * k^3 * binomial(2*n, k)^2 * binomial(2*n-1, k) for n >= 1. (End)

A374440 Triangle read by rows: T(n, k) = T(n - 1, k) + T(n - 2, k - 2), with boundary conditions: if k = 0 or k = 2 then T = 1; if k = 1 then T = n - 1.

Original entry on oeis.org

1, 1, 0, 1, 1, 1, 1, 2, 1, 0, 1, 3, 1, 1, 1, 1, 4, 1, 3, 2, 0, 1, 5, 1, 6, 3, 1, 1, 1, 6, 1, 10, 4, 4, 3, 0, 1, 7, 1, 15, 5, 10, 6, 1, 1, 1, 8, 1, 21, 6, 20, 10, 5, 4, 0, 1, 9, 1, 28, 7, 35, 15, 15, 10, 1, 1, 1, 10, 1, 36, 8, 56, 21, 35, 20, 6, 5, 0

Offset: 0

Peter Luschny, Jul 21 2024

Member of the family of Lucas-Fibonacci polynomials.

			Triangle starts:
  [ 0]  1;
  [ 1]  1,  0;
  [ 2]  1,  1,  1;
  [ 3]  1,  2,  1,  0;
  [ 4]  1,  3,  1,  1,  1;
  [ 5]  1,  4,  1,  3,  2,  0;
  [ 6]  1,  5,  1,  6,  3,  1,  1;
  [ 7]  1,  6,  1, 10,  4,  4,  3,  0;
  [ 8]  1,  7,  1, 15,  5, 10,  6,  1,  1;
  [ 9]  1,  8,  1, 21,  6, 20, 10,  5,  4,  0;
  [10]  1,  9,  1, 28,  7, 35, 15, 15, 10,  1, 1;

Cf. A374441.

Cf. A000032 (Lucas), A001611 (even sums, Fibonacci + 1), A000071 (odd sums, Fibonacci - 1), A001911 (alternating sums, Fibonacci(n+3) - 2), A025560 (row lcm), A073028 (row max), A117671 & A025174 (central terms), A057979 (subdiagonal), A000217 (column 3).

T := proc(n, k) option remember; if k = 0 or k = 2 then 1 elif k > n then 0
elif k = 1 then n - 1 else T(n - 1, k) + T(n - 2, k - 2) fi end:
seq(seq(T(n, k), k = 0..n), n = 0..9);
T := (n, k) -> ifelse(k = 0, 1, binomial(n - floor(k/2), ceil(k/2)) -
binomial(n - ceil((k + irem(k + 1, 2))/2), floor(k/2))):

T(n, k) = binomial(n - floor(k/2), ceiling(k/2)) - binomial(n - ceiling((k + even(k))/2), floor(k/2)) if k > 0, T(n, 0) = 1, where even(k) = 1 if k is even, otherwise 0.

Columns with odd index agree with the odd indexed columns of A374441.

cp's OEIS Frontend

A025174 a(n) = binomial(3n-1, n-1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A045721 a(n) = binomial(3*n+1,n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

Extensions

A165817 Number of compositions (= ordered integer partitions) of n into 2n parts.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Python

Sage

Formula

Extensions

A004319 a(n) = binomial(3*n, n - 1).

Views

Author

Keywords

References

Links

Crossrefs

Programs

Maple

Mathematica

Maxima

PARI

Formula

A013698 a(n) = binomial(3*n+2, n-1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A118113 Even Fibbinary numbers + 1; also 2*Fibbinary(n) + 1.

Views

Author

Keywords

Comments