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 A157181 #2 Oct 12 2012 14:54:56 %S A157181 1,1,1,1,5,1,1,17,17,1,1,45,106,45,1,1,105,524,524,105,1,1,229,2231, %T A157181 4258,2231,229,1,1,481,8547,28771,28771,8547,481,1,1,989,30424,171283, %U A157181 290126,171283,30424,989,1,1,2009,102926,928070,2505074,2505074,928070 %N A157181 A new general triangle sequence based on the Eulerian form in three parts ( subtraction):m=3; t0(n,k)=If[n*k == 0, 1, Sum[(-1)^j Binomial[n + 1, j](k + 1 - j)^n, {j, 0, k + 1}]] t(n,k,m)=If[n == 0, 1, ( m*(n - k) + 1)*t0(n - 1 + 1, k - 1) + (m*k + 1)*t0(n - 1 + 1, k) - m*k*(n - k)*t0(n - 2 + 1, k - 1)]. %C A157181 Row sums are: %C A157181 {1, 2, 7, 36, 198, 1260, 9180, 75600, 695520, 7076160, 78926400,...}. %C A157181 The m=0 of the general sequence is A008518. %F A157181 m=3; %F A157181 t0(n,k)=If[n*k == 0, 1, Sum[(-1)^j Binomial[n + 1, j](k + 1 - j)^n, {j, 0, k + 1}]]; %F A157181 t(n,k,m)=If[n == 0, 1, ( m*(n - k) + 1)*t0(n - 1 + 1, k - 1) + %F A157181 (m*k + 1)*t0(n - 1 + 1, k) + %F A157181 m*k*(n - k)*t0(n - 2 + 1, k - 1)]. %e A157181 {1}, %e A157181 {1, 1}, %e A157181 {1, 5, 1}, %e A157181 {1, 17, 17, 1}, %e A157181 {1, 45, 106, 45, 1}, %e A157181 {1, 105, 524, 524, 105, 1}, %e A157181 {1, 229, 2231, 4258, 2231, 229, 1}, %e A157181 {1, 481, 8547, 28771, 28771, 8547, 481, 1}, %e A157181 {1, 989, 30424, 171283, 290126, 171283, 30424, 989, 1}, %e A157181 {1, 2009, 102926, 928070, 2505074, 2505074, 928070, 102926, 2009, 1}, %e A157181 {1, 4053, 336109, 4684096, 19330402, 30217078, 19330402, 4684096, 336109, 4053, 1} %t A157181 Clear[t, n, k, m]; %t A157181 t[n_, k_, m_] = (m*(n - k) + 1)*Binomial[n - 1, k - 1] + (m*k + 1)*Binomial[n - 1, k] - m*k*(n - k)*Binomial[n - 2, k - 1]; %t A157181 Table[t[n, k, m], {m, 0, 10}, {n, 0, 10}, {k, 0, n}]; %t A157181 Table[Flatten[Table[Table[t[n, k, m], {k, 0, n}], {n, 0, 10}]], {m, 0, 10}] %t A157181 Table[Table[Sum[t[n, k, m], {k, 0, n}], {n, 0, 10}], {m, 0, 10}]; %Y A157181 A008518 %K A157181 nonn,tabl %O A157181 0,5 %A A157181 _Roger L. Bagula_ and _Gary W. Adamson_, Feb 24 2009