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.
Square array begins:
(1),(1),(1),1,(1),1,(1),1,(1),1,1,(1),1,1,(1),1,1,(1),1,1,1,(1),1,1,1,...;
(1),(2),(3),4,(5),6,(7),8,(9),10,11,(12),13,14,(15),16,17,(18),19,20,21,..;
(4),(10),(18),28,(39),52,(66),82,(99),118,138,(159),182,206,(231),258,286,..;
(28),(80),(162),280,(418),600,(806),1064,(1350),1696,2074,(2485),2966,3484,..;
(280),(880),(1944),3640,(5714),8680,(12164),16840,(22194),29080,36824,(45474),.;
(3640),(12320),(29160),58240,(95064),151200,(219108),315440,(428652),581680,...;
(58240),(209440),(524880),1106560,(1864456),3082240,...;
where terms in parenthesis are at positions {[m*(m+5)/6], m>=0}
and are removed before taking partial sums to obtain the next row.
To generate the array, start with all 1's in row 0; from then on,
obtain row n+1 from row n by first removing terms in row n at
positions {[m*(m+5)/6], m>=0} and then taking partial sums.
For example, to generate row 2 from row 1:
[(1),(2),(3),4,(5),6,(7),8,(9),10,11,(12),13,14,(15),16,17,(18),...],
remove terms at positions [0,1,2,4,6,8,11,14,17,...] to get:
[4, 6, 8, 10,11, 13,14, 16,17, 19,20,21, 23,24,25, 27,28,29, ...]
then take partial sums to obtain row 2:
[4, 10, 18, 28,39, 52,66, 82,99, 118,138,159, 182,206,231, ...].
Continuing in this way will generate all the rows of this array.
t[n_, k_] := t[n, k] = Module[{a = 0, m = 0, c = 0, d = 0}, If[n == 0, a = 1, While[d <= k, If[c == Quotient[(m*(m + 5)), 6], m += 1, a += t[n - 1, c]; d += 1]; c += 1]]; a]; Table[t[n - k, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 06 2013, translated from Pari *)
{T(n, k)=local(A=0, m=0, c=0, d=0); if(n==0, A=1, until(d>k, if(c==(m*(m+5))\6, m+=1, A+=T(n-1, c); d+=1); c+=1)); A}
Table begins:
(1),(1),1,(1),1,1,(1),1,1,1,(1),1,1,1,1,(1),1,...;
(1),(2),3,(4),5,6,(7),8,9,10,(11),12,13,14,15,(16),...;
(3),(8),14,(22),31,41,(53),66,80,95,(112),130,149,169,190,...;
(14),(45),86,(152),232,327,(457),606,775,965,(1202),1464,1752,2067,...;
(86),(318),645,(1251),2026,2991,(4455),6207,8274,10684,(13934),17653,...;
(645),(2671),5662,(11869),20143,30827,(48480),70355,96990,128959,...;
(5662),(25805),56632,(126987),223977,352936,(582183),874664,1240239,...;
(56632),(280609),633545,(1508209),2748448,4438122,(7641111),11831184,...;
(633545),(3381993),7820115,(19651299),36837937,60743909,...; ...
where row n equals the partial sums of row n-1 after removing terms
at positions {m*(m+1)/2, m>=0} (marked by parenthesis in above table).
For example, to generate row 3 from row 2:
[3,8, 14, 22, 31,41, 53, 66,80,95, 112, 130,...]
remove terms at positions {0,1,3,6,10,...}, yielding:
[14, 31,41, 66,80,95, 130,149,169,190, ...]
then take partial sums to obtain row 3:
[14, 45,86, 152,232,327, 457,606,775,965, ...].
Continuing in this way generates all rows of this table.
RELATION TO POWERS OF A SPECIAL TRIANGULAR MATRIX.
Columns 0 and 1 are found in triangle T=A152400, which begins:
1;
1, 1;
3, 2, 1;
14, 8, 3, 1;
86, 45, 15, 4, 1;
645, 318, 99, 24, 5, 1;
5662, 2671, 794, 182, 35, 6, 1;
56632, 25805, 7414, 1636, 300, 48, 7, 1; ...
where column k of T = column 0 of matrix power T^(k+1) for k>=0.
Furthermore, matrix powers of triangle T=A152400 satisfy:
column k of T^(j+1) = column j of T^(k+1) for all j>=0, k>=0.
Column 3 of this square array = column 1 of T^2:
1;
2, 1;
8, 4, 1;
45, 22, 6, 1;
318, 152, 42, 8, 1;
2671, 1251, 345, 68, 10, 1;
25805, 11869, 3253, 648, 100, 12, 1; ...
RELATED TRIANGLE A127714 begins:
1;
1, 1, 1;
1, 2, 2, 3, 3, 3;
1, 3, 5, 5, 8, 11, 11, 14, 14, 14;
1, 4, 9, 14, 14, 22, 33, 44, 44, 58, 72, 72, 86, 86, 86;...
where right border = column 0 of this square array.
{T(n, k)=local(A=0, m=0, c=0, d=0); if(n==0, A=1, until(d>k, if(c==m*(m+1)/2, m+=1, A+=T(n-1, c); d+=1); c+=1)); A}
To generate the array, start with all 1's in row 0; from then on,
obtain row n+1 from row n by first removing terms in row n at
positions 0 and {(m+1)*(m+2)/2-2,m>0} and then taking partial sums.
This square array A begins:
(1), (1), 1, 1, (1), 1, 1, 1, (1), 1, 1, 1, 1, (1), 1, 1, 1, 1, 1, ...;
(1), (2), 3, 4, (5), 6, 7, 8, (9), 10, 11, 12, 13, (14), 15, 16, ...;
(3), (7), 13, 20, (28), 38, 49, 61, (74), 89, 105, 122, 140, (159),...;
(13), (33), 71, 120, (181), 270, 375, 497, (637), 817, 1019, 1244, ...;
(71), (191), 461, 836, (1333), 2150, 3169, 4413, (5906), 8001, ...;
(461), (1297), 3447, 6616, (11029), 19030, 29483, 42775, (59324),...;
(3447), (10063), 29093, 58576, (101351), 185674, 300329, 451277, ...;
(29093), (87669), 273343, 573672, (1024949), 1982310, 3330651, ...;
(273343), (847015), 2829325, 6159976, (11320359), 23009602, 39998897, ...;
where terms in parenthesis at positions {0,1,4,8,13,..} in a row
are removed before taking partial sums to obtain the next row.
...
RELATION TO SPECIAL TRIANGLE.
Triangle A104980 begins:
1;
1, 1;
3, 2, 1;
13, 7, 3, 1;
71, 33, 13, 4, 1;
461, 191, 71, 21, 5, 1;
3447, 1297, 461, 133, 31, 6, 1;
29093, 10063, 3447, 977, 225, 43, 7, 1; ...
in which column 0 and column 1 are found in square array A.
...
Matrix square of A104980 = triangle A104988 which begins:
1;
2, 1;
8, 4, 1;
42, 20, 6, 1;
266, 120, 38, 8, 1;
1954, 836, 270, 62, 10, 1;
16270, 6616, 2150, 516, 92, 12, 1;
151218, 58576, 19030, 4688, 882, 128, 14, 1; ...
where column 1 and column 2 are also found in square array A.
{T (n, k)=local (A=0, b=2, c=1, d=0); if (n==0, A=1, until (d>k, if (c==b* (b+1)/2-2, b+=1, A+=T (n-1, c); d+=1); c+=1)); A}
Comments