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, j/d odd} d^k - (d - 1)^k.
EXAMPLE
Square array begins:
1, 1, 1, 1, 1, 1, 1, ...
2, 4, 8, 16, 32, 64, 128, ...
4, 10, 28, 82, 244, 730, 2188, ...
5, 17, 65, 257, 1025, 4097, 16385, ...
7, 27, 127, 627, 3127, 15627, 78127, ...
9, 41, 225, 1313, 7809, 46721, 280065, ...
11, 55, 353, 2419, 16841, 117715, 823673, ...
MATHEMATICA
T[n_, k_] := Sum[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 17 2021 *)
PROG
(PARI) T(n, k) = sum(j=1, n, (n\(2*j-1))^k);
(PARI) T(n, k) = sum(j=1, n, sumdiv(j, d, j/d%2*(d^k-(d-1)^k)));
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Seiichi Manyama, Dec 16 2021
STATUS
approved