OFFSET
1,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..1275 (first 50 antidiagonals)
FORMULA
G.f. of column k: Sum_{j>=1} j*sigma_k(j)*x^j/(1 - x^j).
L.g.f. of column k: -log(Product_{j>=1} (1 - x^j)^sigma_k(j)).
A(n,k) = Sum_{d|n} d^(k+1)*sigma_1(n/d).
EXAMPLE
Square array begins:
1, 1, 1, 1, 1, 1, ...
5, 7, 11, 19, 35, 67, ...
7, 13, 31, 85, 247, 733, ...
17, 35, 95, 311, 1127, 4295, ...
11, 31, 131, 631, 3131, 15631, ...
35, 91, 341, 1615, 8645, 49111, ...
MATHEMATICA
Table[Function[k, Sum[d DivisorSigma[k, d], {d, Divisors[n]}]][i - n], {i, 0, 11}, {n, 1, i}] // Flatten
Table[Function[k, SeriesCoefficient[Sum[j DivisorSigma[k, j] x^j/(1 - x^j), {j, 1, n}], {x, 0, n}]][i - n], {i, 0, 11}, {n, 1, i}] // Flatten
PROG
(PARI) T(n, k)={sumdiv(n, d, d^(k+1)*sigma(n/d))}
for(n=1, 10, for(k=0, 8, print1(T(n, k), ", ")); print); \\ Andrew Howroyd, Nov 26 2018
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Ilya Gutkovskiy, Nov 26 2018
STATUS
approved