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 18 results. Next

A095661 Fifth column (m=4) of (1,3)-Pascal triangle A095660.

Original entry on oeis.org

3, 13, 35, 75, 140, 238, 378, 570, 825, 1155, 1573, 2093, 2730, 3500, 4420, 5508, 6783, 8265, 9975, 11935, 14168, 16698, 19550, 22750, 26325, 30303, 34713, 39585, 44950, 50840, 57288, 64328, 71995, 80325, 89355, 99123, 109668, 121030, 133250, 146370
Offset: 0

Views

Author

Wolfdieter Lang, Jun 11 2004

Keywords

Comments

If Y is a 3-subset of an n-set X then, for n>=6, a(n-6) is the number of 4-subsets of X having at most one element in common with Y. - Milan Janjic, Nov 23 2007
Row 3 of the convolution array A213550. - Clark Kimberling, Jun 20 2012

Crossrefs

Partial sums of A006503.
Cf. A213550.

Programs

Formula

G.f.: (3-2*x)/(1-x)^5.
a(n) = (n+12)*binomial(n+3, 3)/4 = 3*b(n)-2*b(n-1), with b(n) := binomial(n+4, 4); cf. A000332.
a(n) = Sum_{k=1..n} Sum_{i=1..k} i*(n-k+3), with offset 1. - Wesley Ivan Hurt, Sep 25 2013

A095663 Eighth column (m=7) of (1,3)-Pascal triangle A095660.

Original entry on oeis.org

3, 22, 92, 288, 750, 1716, 3564, 6864, 12441, 21450, 35464, 56576, 87516, 131784, 193800, 279072, 394383, 547998, 749892, 1012000, 1348490, 1776060, 2314260, 2985840, 3817125, 4838418, 6084432, 7594752, 9414328, 11594000, 14191056
Offset: 0

Views

Author

Wolfdieter Lang, Jun 11 2004

Keywords

Comments

If Y is a 3-subset of an n-set X then, for n>=9, a(n-9) is the number of 7-subsets of X having at most one element in common with Y. - Milan Janjic, Nov 23 2007

Crossrefs

Seventh column: A095662. Ninth column: A095664.

Formula

G.f.: (3-2*x)/(1-x)^8.
a(n)= binomial(n+6, 6)*(n+21)/7 = 3*b(n)-2*b(n-1), with b(n):=binomial(n+7, 7); cf. A000580.

A095662 Seventh column (m=6) of (1,3)-Pascal triangle A095660.

Original entry on oeis.org

3, 19, 70, 196, 462, 966, 1848, 3300, 5577, 9009, 14014, 21112, 30940, 44268, 62016, 85272, 115311, 153615, 201894, 262108, 336490, 427570, 538200, 671580, 831285, 1021293, 1246014, 1510320, 1819576, 2179672, 2597056, 3078768, 3632475
Offset: 0

Views

Author

Wolfdieter Lang, Jun 11 2004

Keywords

Comments

If Y is a 3-subset of an n-set X then, for n>=8, a(n-8) is the number of 6-subsets of X having at most one element in common with Y. - Milan Janjic, Nov 23 2007

Crossrefs

Sixth column: A000574. Eighth column: A095663.

Programs

  • Mathematica
    CoefficientList[Series[(3-2x)/(1-x)^7,{x,0,40}],x] (* or *) LinearRecurrence[{7,-21,35,-35,21,-7,1},{3,19,70,196,462,966,1848},40] (* Harvey P. Dale, Mar 30 2014 *)

Formula

G.f.: (3-2*x)/(1-x)^7.
a(n)= binomial(n+5, 5)*(n+18)/6 = 3*b(n)-2*b(n-1), with b(n):=binomial(n+6, 6); cf. A000579.
a(0)=3, a(1)=19, a(2)=70, a(3)=196, a(4)=462, a(5)=966, a(6)=1848, a(n)=7*a(n-1)-21*a(n-2)+35*a(n-3)-35*a(n-4)+21*a(n-5)-7*a(n-6)+a(n-7). - Harvey P. Dale, Mar 30 2014

A095664 Ninth column (m=8) of (1,3)-Pascal triangle A095660.

Original entry on oeis.org

3, 25, 117, 405, 1155, 2871, 6435, 13299, 25740, 47190, 82654, 139230, 226746, 358530, 552330, 831402, 1225785, 1773783, 2523675, 3535675, 4884165, 6660225, 8974485, 11960325, 15777450, 20615868, 26700300, 34295052, 43709380, 55303380, 69494436, 86764260, 107666559, 132835365, 162994065
Offset: 0

Views

Author

Wolfdieter Lang, Jun 11 2004

Keywords

Comments

If Y is a 3-subset of an n-set X then, for n>=10, a(n-10) is the number of 8-subsets of X having at most one element in common with Y. - Milan Janjic, Nov 23 2007

Crossrefs

Eighth column: A095663. Tenth column: A095665.

Programs

  • PARI
    x='x+O('x^66); Vec((3-2*x)/(1-x)^9) \\ Joerg Arndt, May 11 2013

Formula

a(n)= binomial(n+7, 7)*(n+24)/8 = 3*b(n)-2*b(n-1), with b(n):=binomial(n+8, 8); cf. A000581.
G.f.: (3-2*x)/(1-x)^9.

A095665 Tenth column (m=9) of (1,3)-Pascal triangle A095660.

Original entry on oeis.org

3, 28, 145, 550, 1705, 4576, 11011, 24310, 50050, 97240, 179894, 319124, 545870, 904400, 1456730, 2288132, 3513917, 5287700, 7811375, 11347050, 16231215, 22891440, 31865925, 43826250, 59603700, 80219568, 106919868, 141214920
Offset: 0

Views

Author

Wolfdieter Lang, Jun 11 2004

Keywords

Comments

If Y is a 3-subset of an n-set X then, for n>=11, a(n-11) is the number of 9-subsets of X having at most one element in common with Y. - Milan Janjic, Nov 23 2007

Crossrefs

Ninth column: A095663.

Programs

  • Mathematica
    Table[Binomial[n+8,8] (n+27)/9,{n,0,30}] (* or *) LinearRecurrence[{10,-45,120,-210,252,-210,120,-45,10,-1},{3,28,145,550,1705,4576,11011,24310,50050,97240},30] (* Harvey P. Dale, Oct 13 2017 *)

Formula

a(n)= binomial(n+8, 8)*(n+27)/9 = 3*b(n)-2*b(n-1), with b(n):=binomial(n+9, 9); cf. A000582.
G.f.: (3-2*x)/(1-x)^10.

A007559 Triple factorial numbers (3*n-2)!!! with leading 1 added.

Original entry on oeis.org

1, 1, 4, 28, 280, 3640, 58240, 1106560, 24344320, 608608000, 17041024000, 528271744000, 17961239296000, 664565853952000, 26582634158080000, 1143053268797440000, 52580450364682240000, 2576442067869429760000, 133974987529210347520000, 7368624314106569113600000
Offset: 0

Views

Author

Keywords

Comments

a(n) is the number of increasing quaternary trees on n vertices. (See A001147 for ternary and A000142 for binary trees.) - David Callan, Mar 30 2007
a(n) is the product of the positive integers k <= 3*n that have k modulo 3 = 1. - Peter Luschny, Jun 23 2011
See A094638 for connections to differential operators. - Tom Copeland, Sep 20 2011
Partial products of A016777. - Reinhard Zumkeller, Sep 20 2013
For n > 2, a(n) is a Zumkeller number. - Ivan N. Ianakiev, Jan 28 2020
a(n) is the number of generalized permutations of length n related to the degenerate Eulerian numbers (see arXiv:2007.13205), cf. A336633. - Orli Herscovici, Jul 28 2020

Examples

			G.f. = 1 + x + 4*x^2 + 28*x^3 + 280*x^4 + 3640*x^5 + 58240*x^6 + ...
a(3) = 28 and a(4) = 280; with top row of M^3 = (28, 117, 108, 27), sum = 280.
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n)= A035469(n, 1), n >= 1, (first column of triangle A035469(n, m)).
Cf. A107716. - Gary W. Adamson, Oct 22 2009
Cf. A095660. - Gary W. Adamson, Jul 19 2011
Subsequence of A007661. A007696, A008548.
a(n) = A286718(n,0), n >= 0.
Row sums of A336633.

Programs

  • GAP
    List([0..20], n-> Product([0..n-1], k-> 3*k+1 )); # G. C. Greubel, Aug 20 2019
  • Haskell
    a007559 n = a007559_list !! n
    a007559_list = scanl (*) 1 a016777_list
    -- Reinhard Zumkeller, Sep 20 2013
    
  • Magma
    b:= func< n | (n lt 2) select n else (3*n-2)*Self(n-1) >;
    [1] cat [b(n): n in [1..20]]; // G. C. Greubel, Aug 20 2019
    
  • Maple
    A007559 := n -> mul(k, k = select(k-> k mod 3 = 1, [$1 .. 3*n])): seq(A007559(n), n = 0 .. 17); # Peter Luschny, Jun 23 2011
    # second Maple program:
    b:= proc(n) option remember; `if`(n<1, 1, n*b(n-3)) end:
    a:= n-> b(3*n-2):
    seq(a(n), n=0..20);  # Alois P. Heinz, Dec 18 2024
  • Mathematica
    a[ n_] := If[ n < 0, 1 / Product[ k, {k, - 2, 3 n - 1, -3}],
      Product[ k, {k, 1, 3 n - 2, 3}]]; (* Michael Somos, Oct 14 2011 *)
    FoldList[Times,1,Range[1,100,3]] (* Harvey P. Dale, Jul 05 2013 *)
    Range[0, 19]! CoefficientList[Series[((1 - 3 x)^(-1/3)), {x, 0, 19}], x] (* Vincenzo Librandi, Oct 08 2015 *)
  • Maxima
    a(n):=if n=1 then 1 else (n)!*(sum(m/n*sum(binomial(k,n-m-k)*(-1/3)^(n-m-k)* binomial (k+n-1,n-1),k,1,n-m),m,1,n)+1); /* Vladimir Kruchinin, Aug 09 2010 */
    
  • PARI
    {a(n) = if( n<0, (-1)^n / prod(k=0,-1-n, 3*k + 2), prod(k=0, n-1, 3*k + 1))}; /* Michael Somos, Oct 14 2011 */
    
  • PARI
    my(x='x+O('x^33)); Vec(serlaplace((1-3*x)^(-1/3))) /* Joerg Arndt, Apr 24 2011 */
    
  • Sage
    def A007559(n) : return mul(j for j in range(1,3*n,3))
    [A007559(n) for n in (0..17)]  # Peter Luschny, May 20 2013
    

Formula

a(n) = Product_{k=0..n-1} (3*k + 1).
a(n) = (3*n - 2)!!!.
a(n) = A007661(3*n-2).
E.g.f.: (1-3*x)^(-1/3).
a(n) ~ sqrt(2*Pi)/Gamma(1/3)*n^(-1/6)*(3*n/e)^n*(1 - (1/36)/n - ...). - Joe Keane (jgk(AT)jgk.org), Nov 22 2001
a(n) = 3^n*Pochhammer(1/3, n).
a(n) = Sum_{k=0..n} (-3)^(n-k)*A048994(n, k). - Philippe Deléham, Oct 29 2005
a(n) = n!*(1+Sum_{m=1..n} (m/n)*Sum_{k=1..n-m} binomial(k, n-m-k)*(-1/3)^(n-m-k)*binomial(k+n-1, n-1)), n>1. - Vladimir Kruchinin, Aug 09 2010
From Gary W. Adamson, Jul 19 2011: (Start)
a(n) = upper left term in M^n, M = a variant of Pascal (1,3) triangle (Cf. A095660); as an infinite square production matrix:
1, 3, 0, 0, 0,...
1, 4, 3, 0, 0,...
1, 5, 7, 3, 0,...
...
a(n+1) = sum of top row terms of M^n. (End)
a(n) = (-2)^n*Sum_{k=0..n} (3/2)^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
G.f.: 1/Q(0) where Q(k) = 1 - x*(3*k+1)/( 1 - x*(3*k+3)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 21 2013
G.f.: G(0)/2, where G(k)= 1 + 1/(1 - x*(3*k+1)/(x*(3*k+1) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 26 2013
E.g.f.: E(0)/2, where E(k)= 1 + 1/(1 - x*(3*k+1)/(x*(3*k+1) + (k+1)/E(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 26 2013
Let D(x) = 1/sqrt(1 - 2*x) be the e.g.f. for the sequence of double factorial numbers A001147. Then the e.g.f. A(x) for the triple factorial numbers satisfies D( Integral_{t=0..x} A(t) dt ) = A(x). Cf. A007696 and A008548. - Peter Bala, Jan 02 2015
O.g.f.: hypergeom([1, 1/3], [], 3*x). - Peter Luschny, Oct 08 2015
a(n) = 3^n * Gamma(n + 1/3)/Gamma(1/3). - Artur Jasinski, Aug 23 2016
a(n) = (-1)^n / A008544(n), 0 = a(n)*(+3*a(n+1) -a(n+2)) +a(n+1)*a(n+1) for all n in Z. - Michael Somos, Sep 30 2018
D-finite with recurrence: a(n) +(-3*n+2)*a(n-1)=0, n>=1. - R. J. Mathar, Feb 14 2020
Sum_{n>=1} 1/a(n) = (e/9)^(1/3) * (Gamma(1/3) - Gamma(1/3, 1/3)). - Amiram Eldar, Jun 29 2020

Extensions

Better description from Wolfdieter Lang

A055998 a(n) = n*(n+5)/2.

Original entry on oeis.org

0, 3, 7, 12, 18, 25, 33, 42, 52, 63, 75, 88, 102, 117, 133, 150, 168, 187, 207, 228, 250, 273, 297, 322, 348, 375, 403, 432, 462, 493, 525, 558, 592, 627, 663, 700, 738, 777, 817, 858, 900, 943, 987, 1032, 1078, 1125, 1173, 1222, 1272
Offset: 0

Views

Author

Barry E. Williams, Jun 14 2000

Keywords

Comments

If X is an n-set and Y a fixed (n-3)-subset of X then a(n-3) is equal to the number of 2-subsets of X intersecting Y. - Milan Janjic, Aug 15 2007
Bisection of A165157. - Jaroslav Krizek, Sep 05 2009
a(n) is the number of (w,x,y) having all terms in {0,...,n} and w=x+y-1. - Clark Kimberling, Jun 02 2012
Numbers m >= 0 such that 8m+25 is a square. - Bruce J. Nicholson, Jul 26 2017
a(n-1) = 3*(n-1) + (n-1)*(n-2)/2 is the number of connected, loopless, non-oriented, multi-edge vertex-labeled graphs with n edges and 3 vertices. Labeled multigraph analog of A253186. There are 3*(n-1) graphs with the 3 vertices on a chain (3 ways to label the middle graph, n-1 ways to pack edges on one of connections) and binomial(n-1,2) triangular graphs (one way to label the graphs, pack 1 or 2 or ...n-2 on the 1-2 edge, ...). - R. J. Mathar, Aug 10 2017
a(n) is also the number of vertices of the quiver for PGL_{n+1} (see Shen). - Stefano Spezia, Mar 24 2020
Starting from a(2) = 7, this is the 4th column of the array: natural numbers written by antidiagonals downwards. See the illustration by Kival Ngaokrajang and the cross-references. - Andrey Zabolotskiy, Dec 21 2021

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, p. 193.

Crossrefs

a(n) = A095660(n+1, 2): third column of (1, 3)-Pascal triangle.
Row n=2 of A255961.

Programs

Formula

G.f.: x*(3-2*x)/(1-x)^3.
a(n) = A027379(n), n > 0.
a(n) = A126890(n,2) for n > 1. - Reinhard Zumkeller, Dec 30 2006
a(n) = A000217(n) + A005843(n). - Reinhard Zumkeller, Sep 24 2008
If we define f(n,i,m) = Sum_{k=0..n-i} binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-m-j), then a(n) = -f(n,n-1,3), for n >= 1. - Milan Janjic, Dec 20 2008
a(n) = A167544(n+8). - Philippe Deléham, Nov 25 2009
a(n) = a(n-1) + n + 2 with a(0)=0. - Vincenzo Librandi, Aug 07 2010
a(n) = Sum_{k=1..n} (k+2). - Gary Detlefs, Aug 10 2010
a(n) = A034856(n+1) - 1 = A000217(n+2) - 3. - Jaroslav Krizek, Sep 05 2009
Sum_{n>=1} 1/a(n) = 137/150. - R. J. Mathar, Jul 14 2012
a(n) = 3*n + A000217(n-1) = 3*n - floor(n/2) + floor(n^2/2). - Wesley Ivan Hurt, Jun 15 2013
a(n) = Sum_{i=3..n+2} i. - Wesley Ivan Hurt, Jun 28 2013
a(n) = 3*A000217(n) - 2*A000217(n-1). - Bruno Berselli, Dec 17 2014
a(n) = A046691(n) + 1. Also, a(n) = A052905(n-1) + 2 = A055999(n-1) + 3 for n>0. - Andrey Zabolotskiy, May 18 2016
E.g.f.: x*(6+x)*exp(x)/2. - G. C. Greubel, Apr 05 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(2)/5 - 47/150. - Amiram Eldar, Jan 10 2021
From Amiram Eldar, Feb 12 2024: (Start)
Product_{n>=1} (1 - 1/a(n)) = -5*cos(sqrt(33)*Pi/2)/(4*Pi).
Product_{n>=1} (1 + 1/a(n)) = 15*cos(sqrt(17)*Pi/2)/(2*Pi). (End)

A000285 a(0) = 1, a(1) = 4, and a(n) = a(n-1) + a(n-2) for n >= 2.

Original entry on oeis.org

1, 4, 5, 9, 14, 23, 37, 60, 97, 157, 254, 411, 665, 1076, 1741, 2817, 4558, 7375, 11933, 19308, 31241, 50549, 81790, 132339, 214129, 346468, 560597, 907065, 1467662, 2374727, 3842389, 6217116, 10059505, 16276621, 26336126, 42612747, 68948873, 111561620, 180510493, 292072113, 472582606
Offset: 0

Views

Author

Keywords

Comments

a(n-1) = Sum_{k=0..ceiling((n-1)/2)} P(4;n-1-k,k), n >= 1, with a(-1)=3. These are the sums over the SW-NE diagonals in P(4;n,k), the (4,1) Pascal triangle A093561. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs. Also SW-NE diagonal sums in the Pascal (1,3) triangle A095660.
In general, for a Fibonacci sequence beginning with 1,b we have a(n) = (2^(-1-n)*((1-sqrt(5))^n*(1+sqrt(5)-2b) + (1+sqrt(5))^n*(-1+sqrt(5)+2b)))/sqrt(5). In this case we have b=4. - Herbert Kociemba, Dec 18 2011
Pisano period lengths: 1, 3, 8, 6, 20, 24, 16, 12, 24, 60, 5, 24, 28, 48, 40, 24, 36, 24, 18, 60, ... - R. J. Mathar, Aug 10 2012
a(n) = number of independent vertex subsets (i.e., the Merrifield-Simmons index) of the tree obtained from the path tree P_{n-1} by attaching two pendant edges to one of its endpoints (n >= 2). Example: if n=3, then we have the star tree with edges ab, ac, ad; it has 9 independent vertex subsets: empty, a, b, c, d, bc, cd, bd, bcd.
For n >= 2, the number a(n-1) is the dimension of a commutative Hecke algebra of type D_n with independent parameters. See Theorem 1.4 and Corollary 1.5 in the link "Hecke algebras with independent parameters". - Jia Huang, Jan 20 2019
For n >= 1, a(n) is the number of edge covers of the tadpole graph T_{3,n-1} with T_{3,0} interpreted as just the cycle C_3. Example: If n=2, we have C_3 and P_1 joined by a bridge, which is just the triangle with a pendant, and this graph has 5 edge covers. In general, because of the path portion of the graph, the number of edge covers of T{3,n-1} satisfies the same recurrence as Fibonacci sequence and it starts with 4,5. - Feryal Alayont, Aug 27 2023
Eswarathasan (1978) called these numbers "pseudo-Fibonacci numbers", and proved that 1, 4, and 9 are the only squares in this sequence. If the recurrence is extended to negative indices, then there is only one more square, a(-9) = 81. Eswarathasan (1979) proved that none of the terms (even with negative indices) are twice a square. - Amiram Eldar, Mar 09 2024
For n>2, a(n) + (-1)^ceiling(n/2) is the number of ways to tile this strip of length n-1, with a central staircase, using unit squares and dominoes:
_
||
|||_|||_|||_|||. - Greg Dresden and Runhe Zhang, Sep 07 2024

Examples

			G.f. = 1 + 4*x + 5*x^2 + 9*x^3 + 14*x^4 + 23*x^5 + 37*x^6 + 60*x^7 + ...
		

References

  • Richard E. Merrifield and Howard E. Simmons, Topological Methods in Chemistry, Wiley, New York, 1989. pp. 131.
  • Joe Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 224.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Essentially the same as A104449, which only has A104449(0)=3 prefixed.
Cf. A090888, A101220, A109754, A091157 (subsequence of primes).

Programs

  • GAP
    F:=Fibonacci;; List([0..40], n-> F(n+2) +2*F(n) ); # G. C. Greubel, Nov 08 2019
  • Haskell
    a000285 n = a000285_list !! n
    a000285_list = 1 : 4 : zipWith (+) a000285_list (tail a000285_list)
    -- Reinhard Zumkeller, Apr 28 2011
    
  • Magma
    a0:=1; a1:=4; [GeneralizedFibonacciNumber(a0,a1,n): n in [0..30]]; // Bruno Berselli, Feb 12 2013
    
  • Maple
    with(combinat):a:=n->2*fibonacci(n)+fibonacci(n+2): seq(a(n), n=0..34);
  • Mathematica
    LinearRecurrence[{1,1},{1,4},40] (* or *) Table[(3*LucasL[n]- Fibonacci[n])/2,{n,40}] (* Harvey P. Dale, Jul 18 2011 *)
    a[ n_]:= Fibonacci[n] + LucasL[n+1]; (* Michael Somos, May 28 2014 *)
  • Maxima
    a[0]:1$ a[1]:4$ a[n]:=a[n-1]+a[n-2]$ makelist(a[n],n,0,30); /* Martin Ettl, Oct 25 2012 */
    
  • PARI
    Vec((1+3*x)/(1-x-x^2)+O(x^40)) \\ Charles R Greathouse IV, Nov 20 2012
    
  • Sage
    f=fibonacci; [f(n+2) +2*f(n) for n in (0..40)] # G. C. Greubel, Nov 08 2019
    

Formula

G.f.: (1+3*x)/(1-x-x^2). - Simon Plouffe in his 1992 dissertation
Row sums of A131775 starting (1, 4, 5, 9, 14, 23, ...). - Gary W. Adamson, Jul 14 2007
a(n) = 2*Fibonacci(n) + Fibonacci(n+2). - Zerinvary Lajos, Oct 05 2007
a(n) = ((1+sqrt(5))^n - (1-sqrt(5))^n)/(2^n*sqrt(5)) + (3/2)* ((1+sqrt(5))^(n-1) - (1-sqrt(5))^(n-1))/(2^(n-2)*sqrt(5)). Offset 1. a(3)=5. - Al Hakanson (hawkuu(AT)gmail.com), Jan 14 2009
a(n) = 3*Fibonacci(n+2) - 2*Fibonacci(n+1). - Gary Detlefs, Dec 21 2010
a(n) = A104449(n+1). - Michael Somos, Apr 07 2012
From Michael Somos, May 28 2014: (Start)
a(n) = A101220(3, 0, n+1).
a(n) = A109754(3, n+1).
a(k) = A090888(2, k-1), for k > 0.
a(-1 - n) = (-1)^n * A013655(n).
a(n) = Fibonacci(n) + Lucas(n+1), see Mathematica field. (End)
11*Fibonacci(n+1) = a(n+3) - a(n-2) = 3*a(n-1) + 2*a(n). - Manfred Arens and Michel Marcus, Jul 14 2014
a(n) = (9*F(n) + F(n-3))/2. - J. M. Bergot, Jul 15 2017
a(n-1) = 3 * A000045(n) + A000045(n+1). - R. J. Mathar, Feb 14 2024
E.g.f.: (cosh(x/2) + sinh(x/2))*(5*cosh(sqrt(5)*x/2) + 7*sqrt(5)*sinh(sqrt(5)*x/2))/5. - Stefano Spezia, Dec 31 2024

A228196 A triangle formed like Pascal's triangle, but with n^2 on the left border and 2^n on the right border instead of 1.

Original entry on oeis.org

0, 1, 2, 4, 3, 4, 9, 7, 7, 8, 16, 16, 14, 15, 16, 25, 32, 30, 29, 31, 32, 36, 57, 62, 59, 60, 63, 64, 49, 93, 119, 121, 119, 123, 127, 128, 64, 142, 212, 240, 240, 242, 250, 255, 256, 81, 206, 354, 452, 480, 482, 492, 505, 511, 512, 100, 287, 560, 806, 932, 962, 974, 997, 1016, 1023, 1024
Offset: 1

Views

Author

Boris Putievskiy, Aug 15 2013

Keywords

Comments

The third row is (n^4 - n^2 + 24*n + 24)/12.
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 04 2013

Examples

			The start of the sequence as a triangular array read by rows:
   0;
   1,  2;
   4,  3,  4;
   9,  7,  7,  8;
  16, 16, 14, 15, 16;
  25, 32, 30, 29, 31, 32;
  36, 57, 62, 59, 60, 63, 64;
		

Crossrefs

Cf. We denote Pascal-like triangle with L(n) on the left border and R(n) on the right border by (L(n),R(n)). A007318 (1,1), A008949 (1,2^n), A029600 (2,3), A029618 (3,2), A029635 (1,2), A029653 (2,1), A037027 (Fibonacci(n),1), A051601 (n,n) n>=0, A051597 (n,n) n>0, A051666 (n^2,n^2), A071919 (1,0), A074829 (Fibonacci(n), Fibonacci(n)), A074909 (1,n), A093560 (3,1), A093561 (4,1), A093562 (5,1), A093563 (6,1), A093564 (7,1), A093565 (8,1), A093644 (9,1), A093645 (10,1), A095660 (1,3), A095666 (1,4), A096940 (1,5), A096956 (1,6), A106516 (3^n,1), A108561(1,(-1)^n), A132200 (4,4), A134636 (2n+1,2n+1), A137688 (2^n,2^n), A160760 (3^(n-1),1), A164844(1,10^n), A164847 (100^n,1), A164855 (101*100^n,1), A164866 (101^n,1), A172171 (1,9), A172185 (9,11), A172283 (-9,11), A177954 (int(n/2),1), A193820 (1,2^n), A214292 (n,-n), A227074 (4^n,4^n), A227075 (3^n,3^n), A227076 (5^n,5^n), A227550 (n!,n!), A228053 ((-1)^n,(-1)^n), A228074 (Fibonacci(n), n).
Cf. A000290 (row 1), A153056 (row 2), A000079 (column 1), A000225 (column 2), A132753 (column 3), A118885 (row sums of triangle array + 1), A228576 (generalized Pascal's triangle).

Programs

  • GAP
    T:= function(n,k)
        if k=0 then return n^2;
        elif k=n then return 2^n;
        else return T(n-1,k-1) + T(n-1,k);
        fi;
      end;
    Flat(List([0..12], n-> List([0..n], k-> T(n,k) ))); # G. C. Greubel, Nov 12 2019
  • Maple
    T:= proc(n, k) option remember;
          if k=0 then n^2
        elif k=n then 2^k
        else T(n-1, k-1) + T(n-1, k)
          fi
        end:
    seq(seq(T(n, k), k=0..n), n=0..10); # G. C. Greubel, Nov 12 2019
  • Mathematica
    T[n_, k_]:= T[n, k] = If[k==0, n^2, If[k==n, 2^k, T[n-1, k-1] + T[n-1, k]]]; Table[T[n, k], {n,0,10}, {k,0,n}]//Flatten (* G. C. Greubel, Nov 12 2019 *)
    Flatten[Table[Sum[i^2 Binomial[n-1-i, n-k-i], {i,1,n-k}] + Sum[2^i Binomial[n-1-i, k-i], {i,1,k}], {n,0,10}, {k,0,n}]] (* Greg Dresden, Aug 06 2022 *)
  • PARI
    T(n,k) = if(k==0, n^2, if(k==n, 2^k, T(n-1, k-1) + T(n-1, k) )); \\ G. C. Greubel, Nov 12 2019
    
  • Python
    def funcL(n):
       q = n**2
       return q
    def funcR(n):
       q = 2**n
       return q
    for n in range (1,9871):
       t=int((math.sqrt(8*n-7) - 1)/ 2)
       i=n-t*(t+1)/2-1
       j=(t*t+3*t+4)/2-n-1
       sum1=0
       sum2=0
       for m1 in range (1,i+1):
          sum1=sum1+funcR(m1)*binomial(i+j-m1-1,i-m1)
       for m2 in range (1,j+1):
          sum2=sum2+funcL(m2)*binomial(i+j-m2-1,j-m2)
       sum=sum1+sum2
    
  • Sage
    @CachedFunction
    def T(n, k):
        if (k==0): return n^2
        elif (k==n): return 2^n
        else: return T(n-1, k-1) + T(n-1, k)
    [[T(n, k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Nov 12 2019
    

Formula

T(n,0) = n^2, n>0; T(0,k) = 2^k; T(n, k) = T(n-1, k-1) + T(n-1, k) for n,k > 0. [corrected by G. C. Greubel, Nov 12 2019]
Closed-form formula for general case. Let L(m) and R(m) be the left border and the right border of Pascal like triangle, respectively. We denote binomial(n,k) by C(n,k).
As table read by antidiagonals T(n,k) = Sum_{m1=1..n} R(m1)*C(n+k-m1-1, n-m1) + Sum_{m2=1..k} L(m2)*C(n+k-m2-1, k-m2); n,k >=0.
As linear sequence a(n) = Sum_{m1=1..i} R(m1)*C(i+j-m1-1, i-m1) + Sum_{m2=1..j} L(m2)*C(i+j-m2-1, j-m2), where i=n-t*(t+1)/2-1, j=(t*t+3*t+4)/2-n-1, t=floor((-1+sqrt(8*n-7))/2); n>0.
Some special cases. If L(m)={b,b,b...} b*A000012, then the second sum takes form b*C(n+k-1,j). If L(m) is {0,b,2b,...} b*A001477, then the second sum takes form b*C(n+k,n-1). Similarly for R(m) and the first sum.
For this sequence L(m)=m^2 and R(m)=2^m.
As table read by antidiagonals T(n,k) = Sum_{m1=1..n} (2^m1)*C(n+k-m1-1, n-m1) + Sum_{m2=1..k} (m2^2)*C(n+k-m2-1, k-m2); n,k >=0.
As linear sequence a(n) = Sum_{m1=1..i} (2^m1)*C(i+j-m1-1, i-m1) + Sum_{m2=1..j} (m2^2)*C(i+j-m2-1, j-m2), where i=n-t*(t+1)/2-1, j=(t*t+3*t+4)/2-n-1, t=floor((-1+sqrt(8*n-7))/2).
As a triangular array read by rows, T(n,k) = Sum_{i=1..n-k} i^2*C(n-1-i, n-k-i) + Sum_{i=1..k} 2^i*C(n-1-i, k-i); n,k >=0. - Greg Dresden, Aug 06 2022

Extensions

Cross-references corrected and extended by Philippe Deléham, Dec 27 2013

A028329 Twice central binomial coefficients.

Original entry on oeis.org

2, 4, 12, 40, 140, 504, 1848, 6864, 25740, 97240, 369512, 1410864, 5408312, 20801200, 80233200, 310235040, 1202160780, 4667212440, 18150270600, 70690527600, 275693057640, 1076515748880, 4208197927440, 16466861455200, 64495207366200, 252821212875504, 991837065896208
Offset: 0

Views

Author

Keywords

Comments

Central elements in the even-Pascal triangle A028326.
If Y is a 3-subset of an 2n-set X then, for n>=3, a(n-1) is the number of (n+1)-subsets of X having at least two elements in common with Y. - Milan Janjic, Dec 16 2007
a(n) denotes the number of ways one can reach the (n,n) point in an n X n grid via the point (n-1, n-1) starting from (0,0) when moving right and up is allowed [From Avik Roy (avik_3.1416(AT)yahoo.co.in), Jan 29 2009]
It appears that a(n-1) is also the number of quivers in the mutation class of twisted types BD_n and CD_n for n >= 3. - Christian Stump, Nov 03 2010
This is the case m = n+1 in the Catalan's formula (2m)!*(2n)!/(m!*(m+n)!*n!) - see Umberto Scarpis in References. - Bruno Berselli, Apr 27 2012
From Ran Pan, Feb 01 2016: (Start)
a(n) is the number of North-East paths from (0,0) to (n+1,n+1) that bounce off the diagonal y = x an even number of times. Details can be found in Section 4.2 in Pan and Remmel's link.
a(n) is the number of North-East paths from (0,0) to (n+1,n+1) that cross the diagonal y = x an even number of times. Details can be found in Section 4.3 in Pan and Remmel's link. (End)

References

  • Umberto Scarpis, Sui numeri primi e sui problemi dell'analisi indeterminata in Questioni riguardanti le matematiche elementari, Nicola Zanichelli Editore (1924-1927, third Edition), page 11.

Crossrefs

Bisection of A047073, A063886.
First differences of A054113.

Programs

  • Magma
    [2*(n+1)*Catalan(n): n in [0..30]]; // G. C. Greubel, Jul 13 2024
    
  • Maple
    seq(add(binomial(2*n,n),k=1..2),n=0..23); # Zerinvary Lajos, Dec 14 2007
  • Mathematica
    Table[2Binomial[2n,n],{n,0,30}] (* Harvey P. Dale, Aug 08 2011 *)
  • PARI
    a(n)=2*binomial(2*n,n)
    
  • SageMath
    [2*binomial(2*n,n) for n in range(31)] # G. C. Greubel, Jul 13 2024

Formula

G.f.: 2/sqrt(1 - 4*x).
a(n) = 2*A000984(n).
a(n) = 2 * binomial(2*n, n).
a(n) = A100320(n) = A095660(2*n,n) for n > 0. - Reinhard Zumkeller, Apr 08 2012
G.f.: G(0), where G(k)= 1 + 1/(1 - 2*x*(2*k + 1)/(2*x*(2*k + 1) + (k + 1)/ G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 07 2013
a(n) = binomial(2*n+2, n+1) - A162551(n). - Ran Pan, Feb 01 2016
D-finite with recurrence: n*a(n) + 2*(-2*n+1)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
E.g.f.: 2*exp(2*x)*BesselI(0, 2*x). - Stefano Spezia, May 11 2024

Extensions

Edited by Michael Somos, Sep 13 2003
Showing 1-10 of 18 results. Next