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.

A136218 Triangle, read by rows of A136219(n) terms, where row n+1 is generated from row n by first inserting zeros in row n at positions {[m*(m+7)/6], m>=0} and then taking partial sums, starting with a '1' in row 0.

Original entry on oeis.org

1, 1, 1, 1, 3, 3, 3, 2, 2, 1, 1, 15, 15, 15, 12, 12, 9, 9, 6, 6, 4, 2, 2, 1, 108, 108, 108, 93, 93, 78, 78, 63, 63, 51, 39, 39, 30, 21, 21, 15, 9, 9, 5, 3, 1, 1, 1036, 1036, 1036, 928, 928, 820, 820, 712, 712, 619, 526, 526, 448, 370, 370, 307, 244, 244, 193, 154, 115, 115, 85
Offset: 0

Views

Author

Paul D. Hanna, Dec 23 2007

Keywords

Comments

A variant of the triple factorial triangle A136213. Compare to square array A136217, which is generated by a complementary process.

Examples

			Triangle begins:
1;
1,1,1;
3,3,3,2,2,1,1;
15,15,15,12,12,9,9,6,6,4,2,2,1;
108,108,108,93,93,78,78,63,63,51,39,39,30,21,21,15,9,9,5,3,1,1;
1036,1036,1036,928,928,820,820,712,712,619,526,526,448,370,370,307,244,244,193,154,115,115,85,64,43,43,28,19,10,10,5,2,1;
12569,12569,12569,11533,11533,10497,10497,9461,9461,8533,7605,7605,6785,5965,5965,5253,4541,4541,3922,3396,2870,2870,2422,2052,1682,1682,1375,1131,887,887,694,540,425,310,310,225,161,118,75,75,47,28,18,8,8,3,1;
...
Number of terms in rows is given by A136219, which starts:
[1,3,7,13,22,33,47,64,84,106,131,159,190,224,261,301,343,388,...].
To generate row 3, start with row 2:
[3,3,3,2,2,1,1];
insert zeros at positions [0,1,3,5,7,10] to get:
[0,0,3,0,3,0,3,0,2,2,0,1,1],
then take reverse partial sums (from right to left) to obtain row 3:
[15,15,15,12,12,9,9,6,6,4,2,2,1].
For row 4, insert zeros in row 3 at positions [0,1,3,5,7,10,13,16,20]:
[0,0,15,0,15,0,15,0,12,12,0,9,9,0,6,6,0,4,2,2,0,1]
then take reverse partial sums to obtain row 4:
[108,108,108,93,93,78,78,63,63,51,39,39,30,21,21,15,9,9,5,3,1,1].
Continuing in this way will generate all the rows of this triangle.
Amazingly, column 0 of this triangle = column 0 of triangle P=A136220:
1;
1, 1;
3, 2, 1;
15, 10, 3, 1;
108, 75, 21, 4, 1;
1036, 753, 208, 36, 5, 1;
12569, 9534, 2637, 442, 55, 6, 1;
185704, 146353, 40731, 6742, 805, 78, 7, 1; ...
where column k of P^3 = column 0 of P^(3k+3) such that
column 0 of P^3 = column 0 of P shift one place left.
		

Crossrefs

Cf. A136221 (column 0), A136219; A136213, A136220.

Programs

  • PARI
    {T(n,k)=local(A=[1],B);if(n>0,for(i=1,n,m=1;B=[0]; for(j=1,#A,if(j+m-1==(m*(m+7))\6,m+=1;B=concat(B,0));B=concat(B,A[j])); A=Vec(Polrev(Vec(Pol(B)/(1-x+O(x^#B)))))));if(k+1>#A,0,A[k+1])} /* for(n=0,6,for(k=0,2*n^2,if(T(n,k)==0,break,print1(T(n,k),",")));print("")) */

A136169 a(n) = 2*a(n-1) - [(n+1)/3] for n>0 with a(0) = 1.

Original entry on oeis.org

1, 2, 3, 5, 9, 16, 30, 58, 113, 223, 443, 882, 1760, 3516, 7027, 14049, 28093, 56180, 112354, 224702, 449397, 898787, 1797567, 3595126, 7190244, 14380480, 28760951, 57521893, 115043777, 230087544, 460175078, 920350146, 1840700281, 3681400551
Offset: 0

Views

Author

Paul D. Hanna, Dec 23 2007

Keywords

Crossrefs

Cf. A136219.

Programs

  • Mathematica
    nxt[{n_,a_}]:={n+1, 2a-Floor[(n+2)/3]}; NestList[nxt,{0,1},50][[All,2]] (* Harvey P. Dale, Aug 09 2021 *)
  • PARI
    {a(n)=if(n==0,1,2*a(n-1)-((n+1)\3))}

Formula

G.f.: (1 - x^2*(1+x+x^2)/(1-x^3)^2 )/(1-2*x).
Showing 1-2 of 2 results.