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.
%I A144401 #5 Mar 17 2013 05:01:01 %S A144401 1,1,1,1,2,1,1,3,3,2,1,4,6,6,3,1,5,10,13,11,4,1,6,15,24,27,18,6,1,7, %T A144401 21,40,55,51,30,9,1,8,28,62,100,116,94,50,13,1,9,36,91,168,231,234, %U A144401 171,81,19,1,10,45,128,266,420,505,460,303,130,28,1,11,55,174,402,714,987,1065 %N A144401 Padovan ( A000931) version of A038137: expansion of polynomials as antidiagonal: p(x,n)=1/(1-x-x^3)^n. %C A144401 Row sums are: 1, 2, 4, 9, 20, 44, 97, 214, 472, 1041, 2296, 5064, 11169, 24634, 54332 (cf. A008998). %C A144401 These polynomials are sort of pseudo-combinations with the last element Padovan instead of one. %C A144401 If you subtract the binomial triangle sequence you get: %C A144401 {0}, %C A144401 {0, 0}, %C A144401 {0, 0, 0}, %C A144401 {0, 0, 0, 1}, %C A144401 {0, 0, 0, 2, 2}, %C A144401 {0, 0, 0, 3, 6, 3}, %C A144401 {0, 0, 0, 4, 12, 12, 5}, %C A144401 {0, 0, 0, 5, 20, 30, 23, 8}, %C A144401 {0, 0, 0, 6, 30, 60, 66, 42, 12} %F A144401 p(x,n)=1/(1-x-x^3)^n; t(n,m)=anti_diagonal_expansion(p(x,n)). %e A144401 {1}, %e A144401 {1, 1}, %e A144401 {1, 2, 1}, %e A144401 {1, 3, 3, 2}, %e A144401 {1, 4, 6, 6, 3}, %e A144401 {1, 5, 10, 13, 11, 4}, %e A144401 {1, 6, 15, 24, 27, 18, 6}, %e A144401 {1, 7, 21, 40, 55, 51, 30, 9}, %e A144401 {1, 8, 28, 62, 100, 116, 94, 50, 13}, %e A144401 {1, 9, 36, 91, 168, 231, 234, 171, 81, 19}, %e A144401 {1, 10, 45, 128, 266, 420, 505, 460, 303, 130, 28}, %e A144401 {1, 11, 55, 174, 402, 714, 987, 1065, 879, 527, 208, 41}, %e A144401 {1, 12, 66, 230, 585, 1152, 1792, 2220, 2175, 1640, 906, 330, 60}, %e A144401 {1, 13, 78, 297, 825, 1782, 3072, 4278, 4815, 4320, 3006, 1539, 520, 88}, %e A144401 {1, 14, 91, 376, 1133, 2662, 5028, 7752, 9807, 10122, 8391, 5424, 2586, 816, 129} %t A144401 Clear[f, b, a, g, h, n, t]; f[t_, n_] = 1/(1 - t - t^3)^n; a = Table[Table[SeriesCoefficient[Series[f[t, m], {t, 0, 30}], n], {n, 0, 30}], {m, 1, 31}]; b = Table[Table[a[[n - m + 1]][[m]], {m, 1, n }], {n, 1, 15}]; Flatten[b] %Y A144401 Cf. A000931, A038137. %K A144401 nonn,tabl,uned %O A144401 1,5 %A A144401 _Roger L. Bagula_ and _Gary W. Adamson_, Oct 03 2008