OFFSET
2,1
LINKS
G. C. Greubel, Rows n=2..100 of triangle, flattened
FORMULA
a(n) = (R*(R+3)-2*(n-3))^2*(R*(1-R)+2*(n+1))^2/16 where R = floor((sqrt(8*n-15)-1)/2). - Luce ETIENNE, Jun 04 2017
G.f.: (x*(1 + 11*x + 11*x^2 + x^3)*(-1 + y)^4 - 2*(-1 + x)^2*x*(1 + x)*(-1 + y)^2*y*(1 + y) + (-1 + x)^4*y*(1 + 11*y + 11*y^2 + y^3))/((-1 + x)^5*(-1 + y)^5). - Stefano Spezia, Oct 21 2018
EXAMPLE
Triangle begins:
9;
64, 25;
225, 144, 49;
576, 441, 256, 81;
1225, 1024, 729, 400, 121;
2304, 2025, 1600, 1089, 576, 169;
...
MAPLE
a:=(n, k)->(n^2-k^2)^2: seq(seq(a(n, k), k=1..n-1), n=2..12); # Muniru A Asiru, Oct 21 2018
MATHEMATICA
Table[n^2 - k^2, {n, 2, 12}, {k, n - 1}]^2 // Flatten (* Michael De Vlieger, Jun 06 2017, after Alonso del Arte at A120070 *)
f[n_] := Table[SeriesCoefficient[(x (1 + 11 x + 11 x^2 + x^3) (-1 + y)^4 - 2 (-1 + x)^2 x (1 + x) (-1 + y)^2 y (1 + y) + (-1 + x)^4 y (1 + 11 y + 11 y^2 + y^3))/((-1 + x)^5 (-1 + y)^5) , {x, 0, i}, {y, 0, j}], {i, n, n}, {j, 1, n-1}]; Flatten[Array[f, 11]] (* Stefano Spezia, Oct 21 2018 *)
PROG
(PARI) for (n=2, 10, for(k=1, n-1, print1((n^2-k^2)^2, ", ")); print()); \\ Michel Marcus, Jun 04 2017
(GAP) Flat(List([2..12], n->List([1..n-1], k->n^4-2*k^2*n^2+k^4))); # Muniru A Asiru, Oct 21 2018
(Magma) [(Floor((Sqrt(8*n-15)-1)/2)*(Floor((Sqrt(8*n-15)-1)/2)+3)-2*(n-3))^2*(Floor((Sqrt(8*n-15)-1)/2)*(1-Floor((Sqrt(8*n-15)-1)/2))+2*(n+1))^2/16: n in [2..30]]; // G. C. Greubel, Oct 20 2018
CROSSREFS
KEYWORD
AUTHOR
Paul Curtz, Sep 20 2009
EXTENSIONS
More terms from Michel Marcus, Jun 04 2017
STATUS
approved