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.

Previous Showing 11-20 of 36 results. Next

A110813 A triangle of pyramidal numbers.

Original entry on oeis.org

1, 3, 1, 5, 4, 1, 7, 9, 5, 1, 9, 16, 14, 6, 1, 11, 25, 30, 20, 7, 1, 13, 36, 55, 50, 27, 8, 1, 15, 49, 91, 105, 77, 35, 9, 1, 17, 64, 140, 196, 182, 112, 44, 10, 1, 19, 81, 204, 336, 378, 294, 156, 54, 11, 1, 21, 100, 285, 540, 714, 672, 450, 210, 65, 12, 1, 23, 121, 385, 825
Offset: 0

Views

Author

Paul Barry, Aug 05 2005

Keywords

Comments

Triangle A029653 less first column. In general, the product (1/(1-x),x/(1-x))*(1+m*x,x) yields the Riordan array ((1+(m-1)x)/(1-x)^2,x/(1-x)) with general term T(n,k)=(m*n-(m-1)*k+1)*C(n+1,k+1)/(n+1). This is the reversal of the (1,m)-Pascal triangle, less its first column. - Paul Barry, Mar 01 2006
The column sequences give, for k=0..10: A005408 (odd numbers), A000290 (squares), A000330, A002415, A005585, A040977, A050486, A053347, A054333, A054334, A057788.
Linked to Chebyshev polynomials by the fact that this triangle with interpolated zeros in the rows and columns is a scaled version of A053120.
Row sums are A033484. Diagonal sums are A001911(n+1) or F(n+4)-2. Factors as (1/(1-x),x/(1-x))*(1+2x,x). Inverse is A110814 or (-1)^(n-k)*A104709.
This triangle is a subtriangle of the [2,1] Pascal triangle A029653 (omit there the first column).
Subtriangle of triangles in A029653, A131084, A208510. - Philippe Deléham, Mar 02 2012
This is the iterated partial sums triangle of A005408 (odd numbers). Such iterated partial sums of arithmetic progression sequences have been considered by Narayana Pandit (see the Mar 20 2015 comment on A000580 where the MacTutor History of Mathematics archive link and the Gottwald et al. reference, p. 338, are given). - Wolfdieter Lang, Mar 23 2015

Examples

			The number triangle T(n, k) begins
n\k  0   1   2   3    4    5    6   7   8  9 10 11
0:   1
1:   3   1
2:   5   4   1
3:   7   9   5   1
4:   9  16  14   6    1
5:  11  25  30  20    7    1
6:  13  36  55  50   27    8    1
7:  15  49  91 105   77   35    9   1
8:  17  64 140 196  182  112   44  10   1
9:  19  81 204 336  378  294  156  54  11  1
10: 21 100 285 540  714  672  450 210  65 12  1
11: 23 121 385 825 1254 1386 1122 660 275 77 13  1
... reformatted by _Wolfdieter Lang_, Mar 23 2015
As a number square S(n, k) = T(n+k, k), rows begin
  1,   1,   1,   1,   1,   1, ...
  3,   4,   5,   6,   7,   8, ...
  5,   9,  14,  20,  27,  35, ...
  7,  16,  30,  50,  77, 112, ...
  9,  25,  55, 105, 182, 294, ...
		

Crossrefs

Programs

  • Mathematica
    Table[2*Binomial[n + 1, k + 1] - Binomial[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, Oct 19 2017 *)
  • PARI
    for(n=0,10, for(k=0,n, print1(2*binomial(n+1, k+1) - binomial(n,k), ", "))) \\ G. C. Greubel, Oct 19 2017

Formula

Number triangle T(n, k) = C(n, k)*(2n-k+1)/(k+1) = 2*C(n+1, k+1) - C(n, k); Riordan array ((1+x)/(1-x)^2, x/(1-x)); As a number square read by antidiagonals, T(n, k)=C(n+k, k)(2n+k+1)/(k+1).
Equals A007318 * an infinite bidiagonal matrix with 1's in the main diagonal and 2's in the subdiagonal. - Gary W. Adamson, Dec 01 2007
Binomial transform of an infinite lower triangular matrix with all 1's in the main diagonal, all 2's in the subdiagonal and the rest zeros. - Gary W. Adamson, Dec 12 2007
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k) - T(n-2,k-1), T(0,0)=T(1,1)=1, T(1,0)=3, T(n,k)=0 if k<0 or if k>n. - Philippe Deléham, Nov 30 2013
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(7 + 9*x + 5*x^2/2! + x^3/3!) = 7 + 16*x + 30*x^2/2! + 50*x^3/3! + 77*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - Peter Bala, Dec 21 2014
T(n, k) = ps(1, 2; k, n-k) with ps(a, d; k, n) = sum(ps(a, d; k-1, j), j=0..n) and input ps(a, d; 0, j) = a + d*j. See the iterated partial sums comment from Mar 23 2015 above. - Wolfdieter Lang, Mar 23 2015
From Franck Maminirina Ramaharo, May 21 2018: (Start)
T(n,k) = coefficients in the expansion of ((x + 2)*(x + 1)^n - 2)/x.
T(n,k) = A135278(n,k) + A135278(n-1,k).
T(n,k) = A097207(n,n-k).
G.f.: (y + 1)/((y - 1)*(x*y + y - 1)).
E.g.f.: ((x + 2)*exp(x*y + y) - 2*exp(y))/x.
(End)

A057788 Expansion of (1+x)/(1-x)^12.

Original entry on oeis.org

1, 13, 90, 442, 1729, 5733, 16744, 44200, 107406, 243542, 520676, 1058148, 2057510, 3848222, 6953544, 12183560, 20764055, 34512075, 56071470, 89224590, 139299615, 213696795, 322561200, 479634480, 703323660, 1018031196, 1455797448, 2058314440, 2879378332
Offset: 0

Views

Author

N. J. A. Sloane, Nov 04 2000

Keywords

Comments

1/2^10 of twelfth unsigned column of triangle A053120 (T-Chebyshev, rising powers, zeros omitted).
If a 2-set Y and an (n-3)-set Z are disjoint subsets of an n-set X then a(n-12) is the number of 12-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 08 2007
11-dimensional square numbers, tenth partial sums of binomial transform of [1,2,0,0,0,...]. a(n) = sum_{i=0..n} C(n+10,i+10)*b(i), where b(i)=[1,2,0,0,0,...]. - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
2*a(n) is number of ways to place 10 queens on an (n+10) X (n+10) chessboard so that they diagonally attack each other exactly 45 times. The maximal possible attack number, p=binomial(k,2) =45 for k=10 queens, is achievable only when all queens are on the same diagonal. In graph-theory representation they thus form the corresponding complete graph. - Antal Pinter, Dec 27 2015

Crossrefs

Partial sums of A054334.
Sixth column of A111125 (s=5, without leading zeros). - Wolfdieter Lang, Oct 18 2012

Programs

  • GAP
    List([0..30], n -> (2*n+11)*Binomial(n+10, 10)/11); # G. C. Greubel, Dec 02 2018
  • Magma
    [Binomial(n+10,10)*(2*n+11)/11: n in [0..40]]; // Vincenzo Librandi, Feb 14 2016
    
  • Maple
    A057788 := proc(n)
            1/39916800*(2*n+11) *(n+10) *(n+9) *(n+8) *(n+7) *(n+6) *(n+5) *(n+4) *(n+3) *(n+2) *(n+ 1) ; end proc: # R. J. Mathar, Mar 22 2011
  • Mathematica
    Table[(2*n+11)*Binomial[n+10, 10]/11, {n,0,40}] (* G. C. Greubel, Dec 02 2018 *)
    CoefficientList[Series[(1 + x) / (1 - x)^12, {x, 0, 40}], x] (* Vincenzo Librandi, Feb 14 2016 *)
    LinearRecurrence[{12,-66,220,-495,792,-924,792,-495,220,-66,12,-1},{1,13,90,442,1729,5733,16744,44200,107406,243542,520676,1058148},30] (* Harvey P. Dale, Sep 07 2022 *)
  • PARI
    Vec((1+x)/(1-x)^12+O(x^99)) \\ Charles R Greathouse IV, Sep 23 2012
    
  • Sage
    [(2*n+11)*binomial(n+10, 10)/11 for n in range(40)] # G. C. Greubel, Dec 02 2018
    

Formula

a(n) = 2*C(n+11, 11) - C(n+10, 10). - Paul Barry, Mar 04 2003
a(n) = C(n+10,10) + 2*C(n+10,11). - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
a(n) = C(n+10,10)*(2n+11)/11. - Antal Pinter, Dec 27 2015
a(n) = 12*a(n-1)-66*a(n-2)+220*a(n-3)-495*a(n-4)+792*a(n-5)-924*a(n-6)+792*a(n-7)-495*a(n-8)+220*a(n-9)-66*a(n-10)+12*a(n-11)-a(n-12) for n >11. - Vincenzo Librandi, Feb 14 2016
a(n) = (2*n+11)*binomial(n+10, 10)/11. - G. C. Greubel, Dec 02 2018
From Amiram Eldar, Jan 26 2022: (Start)
Sum_{n>=0} 1/a(n) = 419751541/13230 - 2883584*log(2)/63.
Sum_{n>=0} (-1)^n/a(n) = 720896*Pi/63 - 237793798/6615. (End)

A156308 Inverse of triangle S(n,m) defined by sequence A156290, n >= 1, 1 <= m <= n.

Original entry on oeis.org

1, 4, 1, 9, 6, 1, 16, 20, 8, 1, 25, 50, 35, 10, 1, 36, 105, 112, 54, 12, 1, 49, 196, 294, 210, 77, 14, 1, 64, 336, 672, 660, 352, 104, 16, 1, 81, 540, 1386, 1782, 1287, 546, 135, 18, 1, 100, 825, 2640, 4290, 4004, 2275, 800, 170, 20, 1
Offset: 1

Views

Author

Hartmut F. W. Hoft, Feb 07 2009

Keywords

Comments

From Wolfdieter Lang, Jun 26 2011: (Start)
This triangle S(n,m) appears as U_m(n) in the Knuth reference on p. 285. It is related to the Riordan triangle T_m(n) = A111125(n,m) by S(n,m) = A111125(n,m) - A111125(n-1,m), n >= m >= 1 (identity on p. 286).
Also, S(n,m)-S(n-1,m) = A111125(n-1,m-1), n >= 2, m >= 1 (identity on p. 286). (End)
These polynomials may be expressed in terms of the Faber polynomials of A263916 and are embedded in A127677 and A208513. - Tom Copeland, Nov 06 2015

Examples

			Triangle starts:
  n=1:  1;
  n=2:  4,  1;
  n=3:  9,  6,  1;
  n=4: 16, 20,  8,  1;
  ...
		

Crossrefs

Same as triangle A208513 with the first column truncated.
Columns: A000290 (m=1), A002415 (m=2), A040977 (m=3), A053347 (m=4), A054334 (m=5).

Programs

  • Magma
    [(n/k)*Binomial(n+k-1, 2*k-1): k in [1..n], n in [1..12]]; // G. C. Greubel, Feb 01 2022
  • Mathematica
    S[m_] := Flatten[Table[k/j Binomial[k + j - 1, 2 j - 1], {k, 1, m}, {j, 1, k}]]
  • Sage
    flatten([[(n/k)*binomial(n+k-1, 2*k-1) for k in (1..n)] for n in (1..12)]) # G. C. Greubel, Feb 01 2022
    

Formula

S(n, m) = (n/m) * binomial(n + m - 1, 2*m - 1).
From Peter Bala, May 01 2014: (Start)
The n-th row o.g.f. is polynomial R(n,x) = 2/x*( T(n,(x + 2)/2) - 1 ), where T(n,x) is Chebyshev polynomial of the first kind. They form a divisibility sequence: if n divides m then R(n,x) divides R(m,x) in the ring Z[x].
R(2*n,x) = (x + 4)*U(n-1,(x + 2)/2)^2;
R(2*n + 1,x) = ( U(n,(x + 2)/2) + U(n-1,(x + 2)/2) )^2.
O.g.f.: Sum_{n >= 0} R(n,x)*z^n = z*(1 + z)/( (1 - z)*(1 - (x + 2)*z + z^2) ). (End)
The polynomial R(n,x) defined above satisfies (x + 1/x - 2) * R(n, x + 1/x - 2) = x^n + 1/x^n - 2. - Alexander Burstein, May 23 2021

Extensions

Edited by Max Alekseyev, Mar 05 2018

A101104 a(1)=1, a(2)=12, a(3)=23, and a(n)=24 for n>=4.

Original entry on oeis.org

1, 12, 23, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24
Offset: 1

Views

Author

Cecilia Rossiter, Dec 15 2004

Keywords

Comments

Original name: The first summation of row 4 of Euler's triangle - a row that will recursively accumulate to the power of 4.

Crossrefs

For other sequences based upon MagicNKZ(n,k,z):
..... | n = 1 | n = 2 | n = 3 | n = 4 | n = 5 | n = 6 | n = 7
---------------------------------------------------------------------------
z = 0 | A000007 | A019590 | .......MagicNKZ(n,k,0) = A008292(n,k+1) .......
z = 1 | A000012 | A040000 | A101101 | thisSeq | A101100 | ....... | .......
z = 2 | A000027 | A005408 | A008458 | A101103 | A101095 | ....... | .......
z = 3 | A000217 | A000290 | A003215 | A005914 | A101096 | ....... | .......
z = 4 | A000292 | A000330 | A000578 | A005917 | A101098 | ....... | .......
z = 5 | A000332 | A002415 | A000537 | A000583 | A022521 | ....... | A255181
Cf. A101095 for an expanded table and more about MagicNKZ.

Programs

  • Mathematica
    MagicNKZ = Sum[(-1)^j*Binomial[n+1-z, j]*(k-j+1)^n, {j, 0, k+1}];Table[MagicNKZ, {n, 4, 4}, {z, 1, 1}, {k, 0, 34}]
    Join[{1, 12, 23},LinearRecurrence[{1},{24},56]] (* Ray Chandler, Sep 23 2015 *)

Formula

a(k) = MagicNKZ(4,k,1) where MagicNKZ(n,k,z) = Sum_{j=0..k+1} (-1)^j*binomial(n+1-z,j)*(k-j+1)^n (cf. A101095). That is, a(k) = Sum_{j=0..k+1} (-1)^j*binomial(4, j)*(k-j+1)^4.
a(1)=1, a(2)=12, a(3)=23, and a(n)=24 for n>=4. - Joerg Arndt, Nov 30 2014
G.f.: x*(1+11*x+11*x^2+x^3)/(1-x). - Colin Barker, Apr 16 2012

Extensions

New name from Joerg Arndt, Nov 30 2014
Original Formula edited and Crossrefs table added by Danny Rorabaugh, Apr 22 2015

A101095 Fourth difference of fifth powers (A000584).

Original entry on oeis.org

1, 28, 121, 240, 360, 480, 600, 720, 840, 960, 1080, 1200, 1320, 1440, 1560, 1680, 1800, 1920, 2040, 2160, 2280, 2400, 2520, 2640, 2760, 2880, 3000, 3120, 3240, 3360, 3480, 3600, 3720, 3840, 3960, 4080, 4200, 4320, 4440, 4560, 4680, 4800, 4920, 5040, 5160, 5280
Offset: 1

Views

Author

Cecilia Rossiter, Dec 15 2004

Keywords

Comments

Original Name: Shells (nexus numbers) of shells of shells of shells of the power of 5.
The (Worpitzky/Euler/Pascal Cube) "MagicNKZ" algorithm is: MagicNKZ(n,k,z) = Sum_{j=0..k+1} (-1)^j*binomial(n + 1 - z, j)*(k - j + 1)^n, with k>=0, n>=1, z>=0. MagicNKZ is used to generate the n-th accumulation sequence of the z-th row of the Euler Triangle (A008292). For example, MagicNKZ(3,k,0) is the 3rd row of the Euler Triangle (followed by zeros) and MagicNKZ(10,k,1) is the partial sums of the 10th row of the Euler Triangle. This sequence is MagicNKZ(5,k-1,2).

Crossrefs

Fourth differences of A000584, third differences of A022521, second differences of A101098, and first differences of A101096.
For other sequences based upon MagicNKZ(n,k,z):
...... | n = 1 | n = 2 | n = 3 | n = 4 | n = 5 | n = 6 | n = 7 | n = 8
--------------------------------------------------------------------------------------
z = 0 | A000007 | A019590 | ....... MagicNKZ(n,k,0) = T(n,k+1) from A008292 .......
z = 1 | A000012 | A040000 | A101101 | A101104 | A101100 | ....... | ....... | .......
z = 2 | A000027 | A005408 | A008458 | A101103 | thisSeq | ....... | ....... | .......
z = 3 | A000217 | A000290 | A003215 | A005914 | A101096 | ....... | ....... | .......
z = 4 | A000292 | A000330 | A000578 | A005917 | A101098 | ....... | ....... | .......
z = 5 | A000332 | A002415 | A000537 | A000583 | A022521 | ....... | A255181 | .......
z = 12 | A001288 | A057788 | ....... | A254870 | A254471 | A254683 | A254646 | A254642
z = 13 | A010965 | ....... | ....... | ....... | A254871 | A254472 | A254684 | A254647
z = 14 | A010966 | ....... | ....... | ....... | ....... | A254872 | ....... | .......
--------------------------------------------------------------------------------------
Cf. A047969.

Programs

  • Magma
    I:=[1,28,121,240,360]; [n le 5 select I[n] else 2*Self(n-1)-Self(n-2): n in [1..50]]; // Vincenzo Librandi, May 07 2015
    
  • Mathematica
    MagicNKZ=Sum[(-1)^j*Binomial[n+1-z, j]*(k-j+1)^n, {j, 0, k+1}];Table[MagicNKZ, {n, 5, 5}, {z, 2, 2}, {k, 0, 34}]
    CoefficientList[Series[(1 + 26 x + 66 x^2 + 26 x^3 + x^4)/(1 - x)^2, {x, 0, 50}], x] (* Vincenzo Librandi, May 07 2015 *)
    Join[{1,28,121,240},Differences[Range[50]^5,4]] (* or *) LinearRecurrence[{2,-1},{1,28,121,240,360},50] (* Harvey P. Dale, Jun 11 2016 *)
  • PARI
    a(n)=if(n>3, 120*n-240, 33*n^2-72*n+40) \\ Charles R Greathouse IV, Oct 11 2015
  • Sage
    [1,28,121]+[120*(k-2) for k in range(4,36)] # Danny Rorabaugh, Apr 23 2015
    

Formula

a(k+1) = Sum_{j=0..k+1} (-1)^j*binomial(n + 1 - z, j)*(k - j + 1)^n; n = 5, z = 2.
For k>3, a(k) = Sum_{j=0..4} (-1)^j*binomial(4, j)*(k - j)^5 = 120*(k - 2).
a(n) = 2*a(n-1) - a(n-2), n>5. G.f.: x*(1+26*x+66*x^2+26*x^3+x^4) / (1-x)^2. - Colin Barker, Mar 01 2012

Extensions

MagicNKZ material edited, Crossrefs table added, SeriesAtLevelR material removed by Danny Rorabaugh, Apr 23 2015
Name changed and keyword 'uned' removed by Danny Rorabaugh, May 06 2015

A157703 G.f.s of the z^p coefficients of the polynomials in the GF2 denominators of A156925.

Original entry on oeis.org

1, 1, 5, 5, 2, 62, 152, 62, 2, 91, 1652, 5957, 5957, 1652, 91, 52, 5240, 77630, 342188, 551180, 342188, 77630, 5240, 52, 12, 8549, 424921, 5629615, 28123559, 61108544, 61108544, 28123559, 5629615, 424921, 8549, 12
Offset: 0

Views

Author

Johannes W. Meijer, Mar 07 2009

Keywords

Comments

The formula for the PDGF2(z;n) polynomials in the GF2 denominators of A156925 can be found below.
The general structure of the GFKT2(z;p) that generate the z^p coefficients of the PDGF2(z; n) polynomials can also be found below. The KT2(z;p) polynomials in the numerators of the GFKT2(z;p) have a nice symmetrical structure.
The sequence of the number of terms of the first few KT2(z;p) polynomials is: 1, 1, 2, 5, 6, 9, 12, 13, 16, 19, 22, 23, 26. The first differences follow a simple pattern. The positions of the 1's follow the Lazy Caterer's sequence A000124 with one exception, here a(0) = 0.
A Maple algorithm that generates relevant GFKT2(z;p) information can be found below.

Examples

			Some PDGF2 (z;n) are:
  PDGF2(z;n=3) = (1-z)^3*(1-2*z)^2*(1-3*z)
  PDGF2(z;n=4) = (1-z)^4*(1-2*z)^3*(1-3*z)^2*(1-4*z)
The first few GFKT2's are:
  GFKT2(z;p=0) = 1/(1-z)
  GFKT2(z;p=1) = -z/(z-1)^4
  GFKT2(z;p=2) = z^2*(5+5*z)/(1-z)^7
Some KT2(z,p) polynomials are:
  KT2(z;p=2) = 5+5*z
  KT2(z;p=3) = 2+62*z+152*z^2+62*z^3+2*z^4
  KT2(z;p=4) = 91+1652*z+5957*z^2+5957*z^3+1652*z^4+91*z^5
		

Crossrefs

Originator sequence A156925.
See A000292 for the z^1 coefficients and A040977 for the z^2 coefficients divided by 5.
Row sums equal A025035.

Programs

  • Maple
    p:=2; fn:=sum((-1)^(n1+1)*binomial(3*p+1,n1) *a(n-n1),n1=1..3*p+1): fk:=rsolve(a(n) = fn,a(k)): for n2 from 0 to 3*p+1 do fz(n2):=product((1-m*z)^(n2+1-m),m=1..n2): a(n2):= coeff(fz(n2),z,p): end do: b:=n-> a(n): seq(b(n), n=0..3*p+1); a(n)=fn; a(k)=sort(simplify(fk)); GFKT2(p):=sum((fk)*z^k,k=0..infinity); q2:=ldegree((numer(GFKT2(p)))): KT2(p):=sort((-1)^p*simplify((GFKT2(p)*(1-z)^(3*p+1))/z^q2),z, ascending);

Formula

PDGF2(z;n) = Product_{m=1..n} (1-m*z)^(n+1-m) with n = 1, 2, 3, ...
GFKT2(z;p) = (-1)^(p)*(z^q2)*KT2(z, p)/(1-z)^(3*p+1) with p = 0, 1, 2, ...
The recurrence relation for the z^p coefficients a(n) is a(n) = Sum_{k=1..3*p+1} (-1)^(k+1)*binomial(3*p + 1, k)*a(n-k) with p = 0, 1, 2, ... .

A027790 a(n) = 10*(n+1)*binomial(n+3,5)/3.

Original entry on oeis.org

10, 80, 350, 1120, 2940, 6720, 13860, 26400, 47190, 80080, 130130, 203840, 309400, 456960, 658920, 930240, 1288770, 1755600, 2355430, 3116960, 4073300, 5262400, 6727500, 8517600, 10687950, 13300560, 16424730, 20137600
Offset: 2

Views

Author

thi ngoc dinh (via R. K. Guy)

Keywords

Comments

Number of 9-subsequences of [ 1, n ] with just 3 contiguous pairs.

Crossrefs

Programs

  • Magma
    A027790:= func< n | 10*(n+1)*Binomial(n+3,5)/3 >;
    [A027790(n): n in [2..40]]; // G. C. Greubel, Feb 21 2025
    
  • Mathematica
    Table[10(n+1) Binomial[n+3,5]/3,{n,2,30}] (* or *) LinearRecurrence[{7,-21,35,-35,21,-7,1},{10,80,350,1120,2940,6720,13860},30] (* Harvey P. Dale, Jan 15 2015 *)
  • SageMath
    def A027790(n): return 10*(n+1)*binomial(n+3,5)//3
    print([A027790(n) for n in range(2,41)]) # G. C. Greubel, Feb 21 2025

Formula

G.f.: 10*(1+x)*x^2/(1-x)^7.
a(n) = binomial(n+1, 3)*binomial(n+3, 3) = A000292(n-1)*A000292(n+1). - Zerinvary Lajos, May 13 2005
a(n) = 10*A040977(n). - R. J. Mathar, May 22 2013
From Amiram Eldar, Jan 06 2021: (Start)
Sum_{n>=2} 1/a(n) = 3*Pi^2/2 - 235/16.
Sum_{n>=2} (-1)^n/a(n) = 3*Pi^2/4 - 117/16. (End)
E.g.f.: (1/36)*x^2*(180 + 300*x + 135*x^2 + 21*x^3 + x^4)*exp(x). - G. C. Greubel, Feb 21 2025

A060556 Bisection of triangle A060098: odd-indexed members of column sequences of A060098 (not counting leading zeros).

Original entry on oeis.org

1, 1, 2, 1, 6, 3, 1, 12, 16, 4, 1, 20, 50, 32, 5, 1, 30, 120, 140, 55, 6, 1, 42, 245, 448, 316, 86, 7, 1, 56, 448, 1176, 1284, 622, 126, 8, 1, 72, 756, 2688, 4170, 3102, 1113, 176, 9, 1, 90, 1200, 5544, 11550, 12122
Offset: 0

Views

Author

Wolfdieter Lang, Apr 06 2001

Keywords

Comments

Row sums give A060557. Column sequences without leading zeros give for m=0..5: A000012 (powers of 1), A002378 = 2*A000217, A004320, 4*A040977, A060558, 2*A060559.
Companion triangle (even-indexed members) A060102.
With offset 1 for n and k, T(n,k) is the number of (1-2-3)-avoiding trapezoidal words of length n that contain n+1-k 1s. A trapezoidal word (following Riordan) is a sequence (a_1,a_2,...,a_n) of integers with 1 <= a_i <= 2i-1. For example, T(3,3)=3 counts 122, 132, 133 and T(4,2)=12 counts 1112, 1113, 1114, 1115, 1116, 1117, 1121, 1131, 1141, 1151, 1211, 1311. - David Callan, Aug 25 2009

Examples

			{1}; {1,2}; {1,6,3}; {1,12,16,4}; ...; Po(3,x) = 3 + x.
		

Formula

a(n, m)= A060098(2*n+1-m, m).
G.f. for column m: (x^m)*Po(m+1, x)/(1-x)^(2*m+1), with Po(n, x) = Sum_{j=0..floor(n/2)} binomial(n, 2*j+1)*x^j (odd members of row n of Pascal triangle A007318).
a(n, m) = Sum_{j=0..floor((m+1)/2)} binomial(n-j+m, 2*m)*binomial(m+1, 2*j+1), n >= m >= 0, otherwise zero.

A133112 Triangular array, read by rows, associated with sums of certain Vandermonde determinants.

Original entry on oeis.org

1, 2, 1, 3, 4, 1, 4, 10, 8, 1, 5, 20, 35, 16, 1, 6, 35, 112, 126, 32, 1, 7, 56, 294, 672, 462, 64, 1, 8, 84, 672, 2772, 4224, 1716, 128, 1, 9, 120, 1386, 9504, 28314, 27456, 6435, 256, 1, 10, 165, 2640, 28314, 151008, 306735, 183040, 24310, 512, 1
Offset: 1

Views

Author

Peter Bala, Sep 18 2007

Keywords

Comments

Appears to be the same as A073165 read as a triangular array (excluding the first column).

Examples

			Triangle starts:
  1;
  2 1;
  3 4 1;
  4 10 8 1;
  5 20 35 16 1;
		

Crossrefs

A000292 (column 2), A040977 (column 3), A133111 (column 4). Cf. A103905.

Formula

T(n,k)=1/(1!*2! ... *k!)*Sum_{1 <= x_1, ...,x_k <= n} |det V(x_1, ...,x_k)|, where V(x_1, ...,x_k) is the Vandermonde matrix of order k. For example, T(n,2) = 1/2*Sum_{1<=i,j<= n} |i-j|. T(n,k) = 1/(1!*2! ... *k!)*Sum_{1 <= x_1, ...,x_k <= n} |(Product_{1 <= i < j <= k} (x_j - x_i) )|.

A217476 Coefficient triangle for the square of the monic integer Chebyshev T-polynomials A127672.

Original entry on oeis.org

4, 0, 1, 4, -4, 1, 0, 9, -6, 1, 4, -16, 20, -8, 1, 0, 25, -50, 35, -10, 1, 4, -36, 105, -112, 54, -12, 1, 0, 49, -196, 294, -210, 77, -14, 1, 4, -64, 336, -672, 660, -352, 104, -16, 1, 0, 81, -540, 1386, -1782, 1287, -546, 135, -18, 1, 4, -100, 825, -2640, 4290, -4004, 2275, -800, 170, -20, 1
Offset: 0

Views

Author

Wolfdieter Lang, Oct 17 2012

Keywords

Comments

The monic integer T-polynomials, called R(n,x) (in Abramowitz-Stegun C(n,x)), with their coefficient triangle given in A127672, when squared, become polynomials in y=x^2:
R(n,x)^2 = sum(T(n,k)*y^k,m=0..n).
R(n,x)^2 = 2 + R(2*n,x). From the bisection of the R-(or T-)polynomials, the even part. Directly from the R(m*n,x)=R(m,R(n,x)) property for m=2.
The o.g.f. is G(z,y) := sum((R(n,sqrt(y))^2)*z^n ,n=0..infinity) = (4 + (4 - 3*y)*z + y*z^2)/((1 +(2-y)*z + z^2)*(1-z)). From the bisection.
The o.g.f.s of the columns k>=1 are x^k*(1-x)/(1+x)^(2*k+1),
and for k=0 the o.g.f. is 4/(1-x^2).
Hetmaniok et al. (2015) refer to these as "modified Chebyshev" polynomials. - N. J. A. Sloane, Sep 13 2016

Examples

			The triangle begins:
n\k 0    1    2      3     4      5     6     7    8   9  10
0:  4
1:  0    1
2:  4   -4    1
3:  0    9   -6      1
4:  4  -16   20     -8     1
5:  0   25  -50     35   -10      1
6:  4  -36  105   -112    54    -12     1
7:  0   49 -196    294  -210     77   -14     1
8:  4  -64  336   -672   660   -352   104   -16    1
9:  0   81 -540   1386 -1782   1287  -546   135  -18   1
10: 4 -100  825  -2640  4290  -4004  2275  -800  170 -20   1
...
n=2:  R(2,x) = -2 + y, R(2,x)^2 = 4 -4*y + y^2, with y=x^2.
n=3:  R(3,x) = 3*x - x^3, R(3,x)^2 = 9*y - 6*y^2 +y^3, with y=x^2.
T(4,1) = 8*(-1)^3*binomial(5,3)/5 = -16.
T(4,0) = 2 + 8*(-1)^4*binomial(4,4)/4 = 4.
T(n,1) = (-1)^(n-1)*2*n*(n+1)!/((n-1)!*2!*(n+1)) = -((-1)^n)*n^2 = A162395(n), n >= 1.
T(n,2) = (-1)^n*A002415(n), n >= 0.
T(n,3) = -(-1)^n*A040977(n-3), n >= 3.
T(n,4) = (-1)^n*A053347(n-4), n >= 4.
T(n,5) = -(-1)^n*A054334(n-5), n >= 5.
		

References

  • E Hetmaniok, P Lorenc, S Damian, et al., Periodic orbits of boundary logistic map and new kind of modified Chebyshev polynomials in R. Witula, D. Slota, W. Holubowski (eds.), Monograph on the Occasion of 100th Birthday Anniversary of Zygmunt Zahorski. Wydawnictwo Politechniki Slaskiej, Gliwice 2015, pp. 325-343.

Crossrefs

Cf. A127672, A158454 (square of S-polynomials), A128495 (sum of square of S-polynomials).

Formula

T(n,k) = [x^(2*k)]R(n,x)^2, with R(n,x) the monic integer version of the Chebyshev T(n,x) polynomial.
T(n,k) = 0 if n=1. ([k=0] means 1 if k=0 else 0).
Previous Showing 11-20 of 36 results. Next