|
|
A165447
|
|
T(n,k) = n^4 - 2*k^2*n^2 + k^4 = A120070(n, k)^2.
|
|
1
|
|
|
9, 64, 25, 225, 144, 49, 576, 441, 256, 81, 1225, 1024, 729, 400, 121, 2304, 2025, 1600, 1089, 576, 169, 3969, 3600, 3025, 2304, 1521, 784, 225, 6400, 5929, 5184, 4225, 3136, 2025, 1024, 289, 9801, 9216, 8281, 7056, 5625, 4096, 2601, 1296, 361, 14400, 13689, 12544, 11025, 9216, 7225, 5184, 3249, 1600, 441, 20449, 19600, 18225, 16384, 14161, 11664, 9025, 6400, 3969, 1936, 529
(list;
table;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
2,1
|
|
LINKS
|
|
|
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
|
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
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|