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

A000581 a(n) = binomial coefficient C(n,8).

Original entry on oeis.org

1, 9, 45, 165, 495, 1287, 3003, 6435, 12870, 24310, 43758, 75582, 125970, 203490, 319770, 490314, 735471, 1081575, 1562275, 2220075, 3108105, 4292145, 5852925, 7888725, 10518300, 13884156, 18156204, 23535820, 30260340, 38608020, 48903492, 61523748, 76904685
Offset: 8

Views

Author

Keywords

Comments

Figurate numbers based on 8-dimensional regular simplex. - Jonathan Vos Post, Nov 28 2004
Just as A005712 and A000574 are described as the coefficients of x^4 and x^5 in the expansion of (1+x+x^2)^n, so should this sequence be described as the coefficients of x^3 therein. - R. K. Guy, Oct 19 2007
Product of 8 consecutive numbers divided by 8!. - Artur Jasinski, Dec 02 2007
In this sequence there are no primes. - Artur Jasinski, Dec 02 2007
a(n) = number of (n-8)-digit numbers with nondescending digits. E.g., a(9) = 9 = {1,2,3,..,9}, a(10) = 45 = {11-19, 22-29, 33-39, ..., 99} [0 is counted as a zero-digit number rather than a 1-digit number]. - Toby Gottfried, Feb 14 2012
a(n) =fallfac(n, 8)/8! = binomial(n, 8) is also the number of independent components of an antisymmetric tensor of rank 8 and dimension n >= 8 (for n = 1..7 this becomes 0). Here fallfac is the falling factorial. - Wolfdieter Lang, Dec 10 2015
Number of compositions (ordered partitions) of n+1 into exactly 9 parts. - Juergen Will, Jan 02 2016
Number of weak compositions (ordered weak partitions) of n-8 into exactly 9 parts. - Juergen Will, Jan 02 2016
Partial sums of A000580. - Art Baker, Mar 26 2019

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.
  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 196.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 7.
  • J. C. P. Miller, editor, Table of Binomial Coefficients. Royal Society Mathematical Tables, Vol. 3, Cambridge Univ. Press, 1954.
  • 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

Programs

Formula

G.f.: x^8/(1-x)^9.
a(n) = A110555(n+1,8). - Reinhard Zumkeller, Jul 27 2005
a(n) = n*(n-1)*(n-2)*(n-3)*(n-4)*(n-5)*(n-6)*(n-7)/8!. - Artur Jasinski, Dec 02 2007
Sum_{k>=8} 1/a(k) = 8/7. - Tom Edgar, Sep 10 2015
Sum_{n>=8} (-1)^n/a(n) = A001787(8)*log(2) - A242091(8)/7! = 1024*log(2) - 74432/105 = 0.9065224171... - Amiram Eldar, Dec 10 2020

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Mar 17 2000
Some formulas referring to other offsets rewritten by R. J. Mathar, Jul 07 2009
3 more terms from William Boyles, Aug 06 2015

A026729 Square array of binomial coefficients T(n,k) = binomial(n,k), n >= 0, k >= 0, read by downward antidiagonals.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 2, 1, 0, 0, 1, 3, 1, 0, 0, 0, 3, 4, 1, 0, 0, 0, 1, 6, 5, 1, 0, 0, 0, 0, 4, 10, 6, 1, 0, 0, 0, 0, 1, 10, 15, 7, 1, 0, 0, 0, 0, 0, 5, 20, 21, 8, 1, 0, 0, 0, 0, 0, 1, 15, 35, 28, 9, 1, 0, 0, 0, 0, 0, 0, 6, 35, 56, 36, 10, 1, 0, 0, 0, 0, 0, 0, 1, 21, 70, 84, 45, 11, 1, 0, 0, 0, 0
Offset: 0

Views

Author

N. J. A. Sloane, Jan 19 2003

Keywords

Comments

The signed triangular matrix T(n,k)*(-1)^(n-k) is the inverse matrix of the triangular Catalan convolution matrix A106566(n,k), n=k>=0, with A106566(n,k) = 0 if nPhilippe Deléham, Aug 01 2005
As a number triangle: unsigned version of A109466. - Philippe Deléham, Oct 26 2008
A063967*A130595 as infinite lower triangular matrices. - Philippe Deléham, Dec 11 2008
Modulo 2, this sequence becomes A106344. - Philippe Deléham, Dec 18 2008
Let {a_(k,i)}, k>=1, i=0,...,k, be the k-th antidiagonal of the array. Then s_k(n) = Sum_{i=0..k}a_(k,i)* binomial(n,k) is the n-th element of the k-th column of A111808. For example, s_1(n) = binomial(n,1) = n is the first column of A111808 for n>1, s_2(n) = binomial(n,1) + binomial(n,2) is the second column of A111808 for n>1, etc. Therefore, in cases k=3,4,5,6,7,8, s_k(n) is A005581(n), A005712(n), A000574(n), A005714(n), A005715(n), A005716(n), respectively. Besides, s_k(n+5) = A064054(n). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
As a triangle, T(n,k) = binomial(k,n-k). - Peter Bala, Nov 27 2015
For all n >= 0, k >= 0, the k-th homology group of the n-torus H_k(T^n) is the free abelian group of rank T(n,k) = binomial(n,k). See the Math Stack Exchange link below. - Jianing Song, Mar 13 2023

Examples

			Array begins
  1 0 0 0 0 0 ...
  1 1 0 0 0 0 ...
  1 2 1 0 0 0 ...
  1 3 3 1 0 0 ...
  1 4 6 4 1 0 ...
As a triangle, this begins
  1
  0 1
  0 1 1
  0 0 2 1
  0 0 1 3 1
  0 0 0 3 4 1
  0 0 0 1 6 5 1
  ...
Production array is
  0    1
  0    1   1
  0   -1   1   1
  0    2  -1   1  1
  0   -5   2  -1  1  1
  0   14  -5   2 -1  1  1
  0  -42  14  -5  2 -1  1  1
  0  132 -42  14 -5  2 -1  1  1
  0 -429 132 -42 14 -5  2 -1  1  1
  ... (Cf. A000108)
		

Crossrefs

The official entry for Pascal's triangle is A007318. See also A052553 (the same array read by upward antidiagonals).
Cf. A030528 (subtriangle for 1<=k<=n).

Programs

  • GAP
    nmax:=15;; T:=List([0..nmax],n->List([0..nmax],k->Binomial(n,k)));;
    b:=List([2..nmax],n->OrderedPartitions(n,2));;
    a:=Flat(List([1..Length(b)],i->List([1..Length(b[i])],j->T[b[i][j][1]][b[i][j][2]]))); # Muniru A Asiru, Jul 17 2018
  • Magma
    /* As triangle: */ [[Binomial(k, n-k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Nov 29 2015
    
  • Maple
    seq(seq(binomial(k,n-k),k=0..n),n=0..12); # Peter Luschny, May 31 2014
  • Mathematica
    Table[Binomial[k, n - k], {n, 0, 12}, {k, 0, n}] // Flatten (* Michael De Vlieger, Nov 28 2015 *)

Formula

As a number triangle, this is defined by T(n,0) = 0^n, T(0,k) = 0^k, T(n,k) = T(n-1,k-1) + Sum_{j, j>=0} (-1)^j*T(n-1,k+j)*A000108(j) for n>0 and k>0. - Philippe Deléham, Nov 07 2005
As a triangle read by rows, it is [0, 1, -1, 0, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 22 2006
As a number triangle, this is defined by T(n, k) = Sum_{i=0..n} (-1)^(n+i)*binomial(n, i)*binomial(i+k, i-k) and is the Riordan array ( 1, x*(1+x) ). The row sums of this triangle are F(n+1). - Paul Barry, Jun 21 2004
Sum_{k=0..n} x^k*T(n,k) = A000007(n), A000045(n+1), A002605(n), A030195(n+1), A057087(n), A057088(n), A057089(n), A057090(n), A057091(n), A057092(n), A057093(n) for n=0,1,2,3,4,5,6,7,8,9,10. - Philippe Deléham, Oct 16 2006
T(n,k) = A109466(n,k)*(-1)^(n-k). - Philippe Deléham, Dec 11 2008
G.f. for the triangular interpretation: -1/(-1+x*y+x^2*y). - R. J. Mathar, Aug 11 2015
For T(0,0) = 0, the triangle below has the o.g.f. G(x,t) = [t*x(1+x)]/[1-t*x(1+x)]. See A109466 for a signed version and inverse, A030528 for reverse and A102426 for a shifted version. - Tom Copeland, Jan 19 2016

A005712 Coefficient of x^4 in expansion of (1+x+x^2)^n.

Original entry on oeis.org

1, 6, 19, 45, 90, 161, 266, 414, 615, 880, 1221, 1651, 2184, 2835, 3620, 4556, 5661, 6954, 8455, 10185, 12166, 14421, 16974, 19850, 23075, 26676, 30681, 35119, 40020, 45415, 51336, 57816, 64889, 72590, 80955, 90021, 99826, 110409, 121810, 134070
Offset: 2

Views

Author

Keywords

Comments

a(n) = A111808(n,4) for n>3. - Reinhard Zumkeller, Aug 17 2005
If a 2-set Y and 2-set Z, having one element in common, are subsets of an n-set X then a(n-3) is the number of 5-subsets of X intersecting both Y and Z. - Milan Janjic, Oct 03 2007
Antidiagonal sums of the convolution array A213781. - Clark Kimberling, Jun 22 2012

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 78.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n)= A027907(n, 4), n >= 2 (fifth column of trinomial coefficients).

Programs

  • Magma
    I:=[1, 6, 19, 45, 90]; [n le 5 select I[n] else 5*Self(n-1)-10*Self(n-2)+10*Self(n-3)-5*Self(n-4)+Self(n-5): n in [1..40]]; Vincenzo Librandi, Jun 16 2012
    
  • Maple
    seq(binomial(n+2,n-2) + binomial(n+1,n-2) - binomial(n,n-2), n=2..50); # Zerinvary Lajos, May 16 2006
    A005712:=(-1-z+z**2)/(z-1)**5; # Conjectured (correctly) by Simon Plouffe in his 1992 dissertation.
    A005712 := n -> GegenbauerC(`if`(4A005712(n)), n=2..20); # Peter Luschny, May 10 2016
  • Mathematica
    CoefficientList[Series[(1+x-x^2)/(1-x)^5,{x,0,40}],x] (* Vincenzo Librandi, Jun 16 2012 *)
    LinearRecurrence[{5,-10,10,-5,1},{1,6,19,45,90},40] (* Harvey P. Dale, Apr 30 2015 *)
  • PARI
    Vec((x^2)*(1+x-x^2)/(1-x)^5+O(x^99)) \\ Charles R Greathouse IV, Sep 23 2012

Formula

G.f.: (x^2)*(1+x-x^2)/(1-x)^5.
a(n) = binomial(n+2,n-2) + binomial(n+1,n-2) - binomial(n,n-2). - Zerinvary Lajos, May 16 2006
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). Vincenzo Librandi, Jun 16 2012
a(n) = binomial(n,2) + 3*binomial(n,3) + binomial(n,4) (see our comment in A026729). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
a(n) = GegenbauerC(N, -n, -1/2) where N = 4 if 4Peter Luschny, May 10 2016
E.g.f.: exp(x)*x^2*(12 + 12*x + x^2)/24. - Stefano Spezia, Jul 09 2023

Extensions

More terms from Vladeta Jovovic, Oct 02 2000

A111808 Left half of trinomial triangle (A027907), triangle read by rows.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 1, 3, 6, 7, 1, 4, 10, 16, 19, 1, 5, 15, 30, 45, 51, 1, 6, 21, 50, 90, 126, 141, 1, 7, 28, 77, 161, 266, 357, 393, 1, 8, 36, 112, 266, 504, 784, 1016, 1107, 1, 9, 45, 156, 414, 882, 1554, 2304, 2907, 3139, 1, 10, 55, 210, 615, 1452, 2850, 4740, 6765, 8350
Offset: 1

Views

Author

Reinhard Zumkeller, Aug 17 2005

Keywords

Comments

Consider a doubly infinite chessboard with squares labeled (n,k), ranks or rows n in Z, files or columns k in Z (Z denotes ...,-2,-1,0,1,2,... ); 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. A026300, A114929, A114972.
Triangle of numbers C^(2)(n-1,k), n>=1, of combinations with repetitions from elements {1,2,...,n} over k, such that every element i, i=1,...,n, appears in a k-combination either 0 or 1 or 2 times (cf. also A213742-A213745). - Vladimir Shevelev and Peter J. C. Moses, Jun 19 2012

References

  • Harrie Grondijs, Neverending Quest of Type C, Volume B - the endgame study-as-struggle.

Crossrefs

Row sums give A027914; central terms give A027908;
T(n, 0) = 0;
T(n, 1) = n for n>1;
T(n, 2) = A000217(n) for n>1;
T(n, 3) = A005581(n) for n>2;
T(n, 4) = A005712(n) for n>3;
T(n, 5) = A000574(n) for n>4;
T(n, 6) = A005714(n) for n>5;
T(n, 7) = A005715(n) for n>6;
T(n, 8) = A005716(n) for n>7;
T(n, 9) = A064054(n-5) for n>8;
T(n, n-5) = A098470(n) for n>4;
T(n, n-4) = A014533(n-3) for n>3;
T(n, n-3) = A014532(n-2) for n>2;
T(n, n-2) = A014531(n-1) for n>1;
T(n, n-1) = A005717(n) for n>0;
T(n, n) = central terms of A027907 = A002426(n).

Programs

  • Maple
    T := (n,k) -> simplify(GegenbauerC(k, -n, -1/2)):
    for n from 0 to 9 do seq(T(n,k), k=0..n) od; # Peter Luschny, May 09 2016
  • Mathematica
    Table[GegenbauerC[k, -n, -1/2], {n,0,10}, {k,0,n}] // Flatten (* G. C. Greubel, Feb 28 2017 *)

Formula

(1 + x + x^2)^n = Sum(T(n,k)*x^k: 0<=k<=n) + Sum(T(n,k)*x^(2*n-k): 0<=k
T(n, k) = A027907(n, k) = Sum_{i=0,..,(k/2)} binomial(n, n-k+2*i) * binomial(n-k+2*i, i), 0<=k<=n.
T(n, k) = GegenbauerC(k, -n, -1/2). - Peter Luschny, May 09 2016

Extensions

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

A095660 Pascal (1,3) triangle.

Original entry on oeis.org

3, 1, 3, 1, 4, 3, 1, 5, 7, 3, 1, 6, 12, 10, 3, 1, 7, 18, 22, 13, 3, 1, 8, 25, 40, 35, 16, 3, 1, 9, 33, 65, 75, 51, 19, 3, 1, 10, 42, 98, 140, 126, 70, 22, 3, 1, 11, 52, 140, 238, 266, 196, 92, 25, 3, 1, 12, 63, 192, 378, 504, 462, 288, 117, 28, 3, 1, 13, 75, 255, 570, 882, 966, 750, 405, 145, 31, 3
Offset: 0

Author

Wolfdieter Lang, May 21 2004

Keywords

Comments

This is the third member, q=3, in the family of (1,q) Pascal triangles: A007318 (Pascal (q=1)), A029635 (q=2) (but with T(0,0)=2, not 1).
This is an example of a Riordan triangle (see A053121 for a comment and the 1991 Shapiro et al. reference on the Riordan group) with o.g.f. of column no. m of the type g(x)*(x*f(x))^m with f(0)=1. Therefore the o.g.f. for the row polynomials p(n,x) = Sum_{m=0..n} T(n,m)*x^m is G(z,x) = g(z)/(1-x*z*f(z)). Here: g(x) = (3-2*x)/(1-x), f(x) = 1/(1-x), hence G(z,x) = (3-2*z)/(1-(1+x)*z).
The SW-NE diagonals give Sum_{k=0..ceiling((n-1)/2)} T(n-1-k,k) = A000285(n-2), n>=2, with n=1 value 3. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs.
Central terms: T(2*n,n) = A028329(n) = A100320(n) for n > 0, A028329 are the central terms of triangle A028326. - Reinhard Zumkeller, Apr 08 2012
Let P be Pascal's triangle, A007318 and R the Riordan array, A097805. Then Pascal triangle (1,q) = ((q-1) * R) + P. Example: Pascal triangle (1,3) = (2 * R) + P. - Gary W. Adamson, Sep 12 2015

Examples

			Triangle starts:
  3;
  1,  3;
  1,  4,  3;
  1,  5,  7,   3;
  1,  6, 12,  10,   3;
  1,  7, 18,  22,  13,   3;
  1,  8, 25,  40,  35,  16,   3;
  1,  9, 33,  65,  75,  51,  19,   3;
  1, 10, 42,  98, 140, 126,  70,  22,   3;
  1, 11, 52, 140, 238, 266, 196,  92,  25,   3;
  1, 12, 63, 192, 378, 504, 462, 288, 117,  28,  3;
  1, 13, 75, 255, 570, 882, 966, 750, 405, 145, 31, 3;
		

Crossrefs

Row sums: A000079(n+1), n>=1, 3 if n=0. Alternating row sums are [3, -2, followed by 0's].
Column sequences (without leading zeros) give for m=1..9 with n>=0: A000027(n+3), A055998(n+1), A006503(n+1), A095661, A000574, A095662, A095663, A095664, A095665.
Cf. A097805.

Programs

  • Haskell
    a095660 n k = a095660_tabl !! n !! k
    a095660_row n = a095660_tabl !! n
    a095660_tabl = [3] : iterate
       (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [1,3]
    -- Reinhard Zumkeller, Apr 08 2012
    
  • Magma
    A095660:= func< n,k | n eq 0 select 3 else (1+2*k/n)*Binomial(n,k) >;
    [A095660(n,k): k in [0..n], n in [1..12]]; // G. C. Greubel, May 02 2021
    
  • Maple
    T(n,k):=piecewise(n=0,3,0Mircea Merca, Apr 08 2012
  • Mathematica
    {3}~Join~Table[(1 + 2 k/n) Binomial[n, k], {n, 11}, {k, 0, n}] // Flatten (* Michael De Vlieger, Sep 14 2015 *)
  • Sage
    def A095660(n,k): return 3 if n==0 else (1+2*k/n)*binomial(n,k)
    flatten([[A095660(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, May 02 2021

Formula

Recursion: T(n, m)=0 if m>n, T(0, 0)= 3; T(n, 0)=1 if n>=1; T(n, m) = T(n-1, m) + T(n-1, m-1).
G.f. column m (without leading zeros): (3-2*x)/(1-x)^(m+1), m>=0.
T(n,k) = (1+2*k/n) * binomial(n,k), for n>0. - Mircea Merca, Apr 08 2012
Closed-form formula for arbitrary left and right borders of Pascal like triangle see A228196. - Boris Putievskiy, Aug 19 2013

A005716 Coefficient of x^8 in expansion of (1+x+x^2)^n.

Original entry on oeis.org

1, 15, 90, 357, 1107, 2907, 6765, 14355, 28314, 52624, 93093, 157950, 258570, 410346, 633726, 955434, 1409895, 2040885, 2903428, 4065963, 5612805, 7646925, 10293075, 13701285, 18050760, 23554206, 30462615, 39070540, 49721892, 62816292, 78816012, 98253540
Offset: 4

Keywords

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 78.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Programs

  • Magma
    I:=[1, 15, 90, 357, 1107, 2907, 6765, 14355, 28314]; [n le 9 select I[n] else 9*Self(n-1)-36*Self(n-2)+84*Self(n-3)-126*Self(n-4)+126*Self(n-5)-84*Self(n-6)+36*Self(n-7)-9*Self(n-8)+Self(n-9): n in [1..40]]; // Vincenzo Librandi, Jun 16 2012
    
  • Magma
    /* By definition: */ P:=PolynomialRing(Integers()); [ Coefficients((1+x+x^2)^n)[9]: n in [4..32] ]; // Bruno Berselli, Jun 17 2012
  • Maple
    A005716:=-(6*z-9*z**2+3*z**3+1)/(z-1)**9; # Conjectured by Simon Plouffe in his 1992 dissertation.
    A005716 := n -> GegenbauerC(`if`(8A005716(n)), n=4..20); # Peter Luschny, May 10 2016
  • Mathematica
    CoefficientList[Series[(1+6*x-9*x^2+3*x^3)/(1-x)^9,{x,0,40}],x] (* Vincenzo Librandi, Jun 16 2012 *)

Formula

a(n) = binomial(n+1, 5)*(n^2+23*n-84)*(n+10)/336, n >= 4.
G.f.: (x^4)*(1+6*x-9*x^2+3*x^3)/(1-x)^9. (Numerator polynomial is N3(8, x) from A063420).
a(n) = A027907(n, 8), n >= 4 (ninth column of trinomial coefficients).
a(n) = A111808(n,8) for n>7. - Reinhard Zumkeller, Aug 17 2005
a(n) = 9*a(n-1) - 36*a(n-2) + 84*a(n-3) - 126*a(n-4) + 126*a(n-5) - 84*a(n-6) + 36*a(n-7) - 9*a(n-8) + a(n-9). Vincenzo Librandi, Jun 16 2012
a(n) = binomial(n,4) + 10*binomial(n,5) + 15*binomial(n,6) + 7*binomial(n,7) + binomial(n,8) (see our comment in A026729). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
a(n) = GegenbauerC(N, -n, -1/2) where N = 8 if 8Peter Luschny, May 10 2016

Extensions

More terms from Vladeta Jovovic, Oct 02 2000

A005714 Coefficient of x^6 in expansion of (1+x+x^2)^n.

Original entry on oeis.org

1, 10, 45, 141, 357, 784, 1554, 2850, 4917, 8074, 12727, 19383, 28665, 41328, 58276, 80580, 109497, 146490, 193249, 251713, 324093, 412896, 520950, 651430, 807885, 994266, 1214955, 1474795, 1779121, 2133792, 2545224, 3020424, 3567025, 4193322, 4908309, 5721717
Offset: 3

Keywords

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 78.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Programs

  • Magma
    I:=[1, 10, 45, 141, 357, 784, 1554]; [n le 7 select I[n] else 7*Self(n-1)-21*Self(n-2)+35*Self(n-3)-35*Self(n-4)+21*Self(n-5)-7*Self(n-6)+Self(n-7): n in [1..40]]; // Vincenzo Librandi, Jun 16 2012
    
  • Magma
    /* By definition: */ P:=PolynomialRing(Integers()); [ Coefficients((1+x+x^2)^n)[7]: n in [3..35] ]; // Bruno Berselli, Jun 17 2012
  • Maple
    A005714:=-(1+3*z-4*z**2+z**3)/(z-1)**7; # Conjectured by Simon Plouffe in his 1992 dissertation.
    A005714 := n -> GegenbauerC(`if`(6A005714(n)), n=3..20); # Peter Luschny, May 10 2016
  • Mathematica
    a[n_] := Coefficient[(1 + x + x^2)^n, x, 6]; Table[a[n], {n, 3, 35}]
    CoefficientList[Series[(1+3*x-4*x^2+x^3)/(1-x)^7,{x,0,40}],x] (* Vincenzo Librandi, Jun 16 2012 *)

Formula

a(n) = binomial(n, 3)*(n^3+18*n^2+17*n-120) /120.
G.f.: (x^3)*(1+3*x-4*x^2+x^3)/(1-x)^7. (Numerator polynomial is N3(6, x) from A063420).
a(n) = A027907(n, 6), n >= 3 (seventh column of trinomial coefficients).
a(n) = A111808(n,6) for n>5. - Reinhard Zumkeller, Aug 17 2005
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). Vincenzo Librandi, Jun 16 2012
a(n) = binomial(n,3) + 6*binomial(n,4) + 5*binomial(n,5) + binomial(n,6) (see our comment in A026729). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
a(n) = GegenbauerC(N, -n, -1/2) where N = 6 if 6Peter Luschny, May 10 2016
E.g.f.: exp(x)*x^3*(120 + 180*x + 30*x^2 + x^3)/720. - Stefano Spezia, Mar 28 2023

Extensions

More terms from Vladeta Jovovic, Oct 02 2000

A278309 T(n,k)=Number of nXk 0..2 arrays with rows and columns in lexicographic nondecreasing order but with exactly one mistake.

Original entry on oeis.org

0, 3, 3, 16, 32, 16, 51, 294, 294, 51, 126, 2089, 4558, 2089, 126, 266, 11486, 70795, 70795, 11486, 266, 504, 51562, 986014, 2360544, 986014, 51562, 504, 882, 197981, 11557658, 79562696, 79562696, 11557658, 197981, 882, 1452, 672365, 114457714
Offset: 1

Author

R. H. Hardin, Nov 17 2016

Keywords

Comments

Table starts
...0......3........16............51..............126.................266
...3.....32.......294..........2089............11486...............51562
..16....294......4558.........70795...........986014............11557658
..51...2089.....70795.......2360544.........79562696..........2506281752
.126..11486....986014......79562696.......6345491150........507575149862
.266..51562..11557658....2506281752.....507575149862.....100825279690194
.504.197981.114457714...69684770828...38819080346585...20065923383306483
.882.672365.979384739.1689884963173.2710823731820118.3886257287342627627

Examples

			Some solutions for n=3 k=4
..1..2..1..2. .0..2..2..1. .1..1..1..2. .1..2..2..2. .0..1..1..2
..2..0..0..1. .2..0..1..2. .2..2..2..0. .0..0..1..2. .0..0..2..2
..2..1..1..1. .2..1..2..2. .0..1..1..1. .0..1..2..2. .1..1..2..1
		

Crossrefs

Column 1 is A000574(n+1).

Formula

Empirical for column k:
k=1: a(n) = (1/120)*n^5 + (1/8)*n^4 + (5/24)*n^3 - (1/8)*n^2 - (13/60)*n
k=2: [polynomial of degree 17]
k=3: [polynomial of degree 53]
k=4: [polynomial of degree 161]

A278627 T(n,k)=Number of nXk 0..2 arrays with rows in nondecreasing lexicographic order and columns in nonincreasing lexicographic order, but with exactly one mistake.

Original entry on oeis.org

0, 3, 3, 16, 46, 16, 51, 357, 357, 51, 126, 1952, 4754, 1952, 126, 266, 8518, 49503, 49503, 8518, 266, 504, 31605, 439446, 1069536, 439446, 31605, 504, 882, 103546, 3438414, 21121532, 21121532, 3438414, 103546, 882, 1452, 307087, 24103803
Offset: 1

Author

R. H. Hardin, Nov 24 2016

Keywords

Comments

Table starts
....0......3........16............51..............126.................266
....3.....46.......357..........1952.............8518...............31605
...16....357......4754.........49503...........439446.............3438414
...51...1952.....49503.......1069536.........21121532...........387542112
..126...8518....439446......21121532........978005050.........43853346948
..266..31605...3438414.....387542112......43853346948.......4902306226424
..504.103546..24103803....6594175430....1892563134910.....540194658701142
..882.307087.153073965..103536313036...77595353266488...58237100230743229
.1452.838936.888863183.1496475492375.2984253200734849.6049936223396297740

Examples

			Some solutions for n=3 k=4
..1..1..0..0. .2..2..1..1. .1..1..1..0. .2..1..0..0. .1..2..1..1
..1..0..1..0. .1..1..2..2. .2..2..2..1. .2..1..1..0. .1..2..2..0
..2..1..1..2. .2..0..2..2. .2..1..0..1. .0..0..2..2. .2..1..0..1
		

Crossrefs

Column 1 is A000574(n+1).

Formula

Empirical for column k:
k=1: a(n) = (1/120)*n^5 + (1/8)*n^4 + (5/24)*n^3 - (1/8)*n^2 - (13/60)*n
k=2: [polynomial of degree 15]
k=3: [polynomial of degree 43]
k=4: [polynomial of degree 125]

A005715 Coefficient of x^7 in expansion of (1+x+x^2)^n.

Original entry on oeis.org

4, 30, 126, 393, 1016, 2304, 4740, 9042, 16236, 27742, 45474, 71955, 110448, 165104, 241128, 344964, 484500, 669294, 910822, 1222749, 1621224, 2125200, 2756780, 3541590, 4509180, 5693454, 7133130, 8872231, 10960608, 13454496
Offset: 4

Keywords

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 78.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Programs

  • Magma
    I:=[4, 30, 126, 393, 1016, 2304, 4740, 9042]; [n le 8 select I[n] else 8*Self(n-1)-28*Self(n-2)+56*Self(n-3)-70*Self(n-4)+56*Self(n-5)-28*Self(n-6)+8*Self(n-7)-Self(n-8): n in [1..40]]; // Vincenzo Librandi, Jun 16 2012
    
  • Magma
    /* By definition: */ P:=PolynomialRing(Integers()); [ Coefficients((1+x+x^2)^n)[8]: n in [4..33] ]; // Bruno Berselli, Jun 17 2012
  • Maple
    A005715:=(z-2)*(z**2-2)/(z-1)**8; # Conjectured by Simon Plouffe in his 1992 dissertation.
    A005715 := n -> GegenbauerC(`if`(7A005715(n)), n=4..20); # Peter Luschny, May 10 2016
  • Mathematica
    CoefficientList[Series[(x-2)*(x^2-2)/(1-x)^8,{x,0,40}],x] (* Vincenzo Librandi, Jun 16 2012 *)

Formula

a(n) = binomial(n, 4)*(n^3+27*n^2+116*n-120)/210, n >= 4.
G.f.: (x^4)*(x-2)*(x^2-2)/(1-x)^8. (Numerator polynomial is N3(7, x) from A063420).
a(n) = A027907(n, 7), n >= 4 (eighth column of trinomial coefficients).
a(n) = A111808(n,7) for n>6. - Reinhard Zumkeller, Aug 17 2005
a(n) = 8*a(n-1) -28*a(n-2) +56*a(n-3) -70*a(n-4) +56*a(n-5) -28*a(n-6) +8*a(n-7) -a(n-8). Vincenzo Librandi, Jun 16 2012
a(n) = 4*binomial(n,4) + 10*binomial(n,5) + 6*binomial(n,6) + binomial(n,7) (see our comment in A026729). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
a(n) = GegenbauerC(N, -n, -1/2) where N = 7 if 7Peter Luschny, May 10 2016

Extensions

More terms from Vladeta Jovovic, Oct 02 2000
Showing 1-10 of 11 results. Next