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-2 of 2 results.

A081581 Pascal-(1,6,1) array.

Original entry on oeis.org

1, 1, 1, 1, 8, 1, 1, 15, 15, 1, 1, 22, 78, 22, 1, 1, 29, 190, 190, 29, 1, 1, 36, 351, 848, 351, 36, 1, 1, 43, 561, 2339, 2339, 561, 43, 1, 1, 50, 820, 5006, 9766, 5006, 820, 50, 1, 1, 57, 1128, 9192, 28806, 28806, 9192, 1128, 57, 1, 1, 64, 1485, 15240, 68034, 116208, 68034, 15240, 1485, 64, 1
Offset: 0

Views

Author

Paul Barry, Mar 23 2003

Keywords

Comments

One of a family of Pascal-like arrays. A007318 is equivalent to the (1,0,1)-array. A008288 is equivalent to the (1,1,1)-array. Rows include A016993, A081591, A081592. Coefficients of the row polynomials in the Newton basis are given by A013614.

Examples

			Rows start as:
  1,  1,   1,    1,    1, ... A000012;
  1,  8,  15,   22,   29, ... A016993;
  1, 15,  78,  190,  351, ... A081591;
  1, 22, 190,  848, 2339, ...
  1, 29, 351, 2339, 9766, ...
The triangle starts as:
  1;
  1,  1;
  1,  8,   1;
  1, 15,  15,    1;
  1, 22,  78,   22,    1;
  1, 29, 190,  190,   29,   1;
  1, 36, 351,  848,  351,  36,  1;
  1, 43, 561, 2339, 2339, 561, 43, 1;
		

Crossrefs

Cf. Pascal (1,m,1) array: A123562 (m = -3), A098593 (m = -2), A000012 (m = -1), A007318 (m = 0), A008288 (m = 1), A081577 (m = 2), A081578 (m = 3), A081579 (m = 4), A081580 (m = 5), A081582 (m = 7), A143683 (m = 8).

Programs

  • Magma
    A081581:= func< n,k,q | (&+[Binomial(k, j)*Binomial(n-j, k)*q^j: j in [0..n-k]]) >;
    [A081581(n,k,6): k in [0..n], n in [0..12]]; // G. C. Greubel, May 26 2021
    
  • Mathematica
    Table[Hypergeometric2F1[-k, k-n, 1, 7], {n,0,10}, {k,0,n}]//Flatten (* Jean-François Alcover, May 24 2013 *)
  • PARI
    t(n, k) = sum(j=0, n-k, binomial(n-k, j)*binomial(k, j)*7^j) \\ Michel Marcus, May 24 2013
    
  • Sage
    flatten([[hypergeometric([-k, k-n], [1], 7).simplify() for k in (0..n)] for n in (0..12)]) # G. C. Greubel, May 26 2021

Formula

Square array T(n, k) defined by T(n, 0) = T(0, k) = 1, T(n, k) = T(n, k-1) + 6*T(n-1, k-1) + T(n-1, k).
Rows are the expansions of (1+6*x)^k/(1-x)^(k+1).
T(n, k) = Hypergeometric2F1([-k, k-n], [1], 7). - Jean-François Alcover, May 24 2013
E.g.f. for the n-th subdiagonal, n = 0,1,2,..., equals exp(x)*P(n,x), where P(n,x) is the polynomial Sum_{k = 0..n} binomial(n,k)*(7*x)^k/k!. For example, the e.g.f. for the second subdiagonal is exp(x)*(1 + 14*x + 49*x^2/2) = 1 + 15*x + 78*x^2/2! + 190*x^3/3! + 351*x^4/4! + 561*x^5/5! + .... - Peter Bala, Mar 05 2017
From G. C. Greubel, May 26 2021: (Start)
T(n, k, m) = Sum_{j=0..n-k} binomial(k,j)*binomial(n-j,k)*m^j, for m = 6.
Sum_{k=0..n} T(n, k, 6) = A083099(n+1). (End)

A081591 Third row of Pascal-(1,6,1) array A081581.

Original entry on oeis.org

1, 15, 78, 190, 351, 561, 820, 1128, 1485, 1891, 2346, 2850, 3403, 4005, 4656, 5356, 6105, 6903, 7750, 8646, 9591, 10585, 11628, 12720, 13861, 15051, 16290, 17578, 18915, 20301, 21736, 23220, 24753, 26335, 27966, 29646, 31375, 33153, 34980, 36856, 38781, 40755
Offset: 0

Views

Author

Paul Barry, Mar 23 2003

Keywords

Comments

1. Smallest triangular number T(k) (other than the trivial adjacent ones) such that T(n) + T(k) is a square. T(n-1) and T(n+1) are trivial triangular numbers such that T(n) + T(n-1) and T(n) + T(n+1) both are squares. 0+1 = 1, 1+15 = 16, 3+78 = 81, 6+190 = 196 etc. 2. (7n+5)-th triangular number. - Amarnath Murthy, Jun 20 2003

Crossrefs

Programs

  • Magma
    [(2-21*n+49*n^2)/2: n in [0..50]]; // Vincenzo Librandi, Jun 18 2011
    
  • Mathematica
    Table[(2-21n+49n^2)/2,{n,0,40}] (* or *) LinearRecurrence[{3,-3,1},{1,15,78},40] (* Harvey P. Dale, Aug 03 2012 *)
  • PARI
    a(n)=(2-21*n+49*n^2)/2 \\ Charles R Greathouse IV, Jun 17 2017

Formula

a(n) = (2 - 21*n + 49*n^2)/2.
G.f.: (1+6*x)^2/(1-x)^3.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3); a(0)=1, a(1)=15, a(2)=78. - Harvey P. Dale, Aug 03 2012
E.g.f.: exp(x)*(2 + 28*x + 49*x^2)/2. - Elmo R. Oliveira, Jun 09 2025
Showing 1-2 of 2 results.