OFFSET
1,3
FORMULA
G.f. of column k: (1/(1 - x)) * Sum_{j>=1} (j^k - (j - 1)^k) * x^j/(1 + x^(2*j)).
T(n,k) = Sum_{j=1..n} Sum_{d|j} A101455(j/d) * (d^k - (d - 1)^k).
EXAMPLE
Square array begins:
1, 1, 1, 1, 1, 1, 1, ...
2, 4, 8, 16, 32, 64, 128, ...
2, 8, 26, 80, 242, 728, 2186, ...
3, 15, 63, 255, 1023, 4095, 16383, ...
5, 25, 125, 625, 3125, 15625, 78125, ...
5, 33, 209, 1281, 7745, 46593, 279809, ...
5, 45, 335, 2385, 16775, 117585, 823415, ...
MATHEMATICA
T[n_, k_] := Sum[(-1)^(j + 1) * Floor[n/(2*j - 1)]^k, {j, 1, n}]; Table[T[k, n - k + 1], {n, 1, 11}, {k, 1, n}] // Flatten (* Amiram Eldar, Dec 18 2021 *)
PROG
(PARI) T(n, k) = sum(j=1, n, (-1)^(j+1)*(n\(2*j-1))^k);
(PARI) T(n, k) = sum(j=1, n, sumdiv(j, d, kronecker(-4, j/d)*(d^k-(d-1)^k)));
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Seiichi Manyama, Dec 18 2021
STATUS
approved