A322264 Square array A(n,k), n >= 1, k >= 0, read by antidiagonals: A(n,k) = denominator of Sum_{d|n} 1/d^k.
1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 8, 9, 4, 1, 1, 16, 27, 16, 5, 1, 1, 32, 81, 64, 25, 1, 1, 1, 64, 243, 256, 125, 18, 7, 1, 1, 128, 729, 1024, 625, 6, 49, 8, 1, 1, 256, 2187, 4096, 3125, 648, 343, 64, 9, 1, 1, 512, 6561, 16384, 15625, 648, 2401, 512, 81, 5, 1, 1, 1024, 19683, 65536, 78125, 23328, 16807, 4096, 729, 10, 11, 1
Offset: 1
Examples
Square array begins: 1, 1, 1, 1, 1, 1, ... 2, 3/2, 5/4, 9/8, 17/16, 33/32, ... 2, 4/3, 10/9, 28/27, 82/81, 244/243, ... 3, 7/4, 21/16, 73/64, 273/256, 1057/1024, ... 2, 6/5, 26/25, 126/125, 626/625, 3126/3125, ... 4, 2, 25/18, 7/6, 697/648, 671/648, ...
Crossrefs
Programs
-
Mathematica
Table[Function[k, Denominator[DivisorSigma[-k, n]]][i - n], {i, 0, 12}, {n, 1, i}] // Flatten Table[Function[k, Denominator[DivisorSigma[k, n]/n^k]][i - n], {i, 0, 12}, {n, 1, i}] // Flatten Table[Function[k, Denominator[SeriesCoefficient[Sum[x^j/(j^k (1 - x^j)), {j, 1, n}], {x, 0, n}]]][i - n], {i, 0, 12}, {n, 1, i}] // Flatten
Formula
G.f. of column k: Sum_{j>=1} x^j/(j^k*(1 - x^j)) (for rationals Sum_{d|n} 1/d^k).
Dirichlet g.f. of column k: zeta(s)*zeta(s+k) (for rationals Sum_{d|n} 1/d^k).
A(n,k) = denominator of sigma_k(n)/n^k.