OFFSET
1,2
FORMULA
A(n,k) = -Sum_{j=1..n} (-k)^gcd(n,j).
G.f. of column k: k * Sum_{j>=1} phi(j) * x^j / (1 + k*x^j).
EXAMPLE
Square array begins:
1, 2, 3, 4, 5, 6, 7, ...
0, -2, -6, -12, -20, -30, -42, ...
3, 12, 33, 72, 135, 228, 357, ...
0, -16, -84, -264, -640, -1320, -2436, ...
5, 40, 255, 1040, 3145, 7800, 16835, ...
0, -60, -714, -4056, -15540, -46500, -117390, ...
7, 140, 2205, 16408, 78155, 279972, 823585, ...
PROG
(PARI) a(n, k) = -sumdiv(n, d, eulerphi(n/d)*(-k)^d);
CROSSREFS
KEYWORD
sign,tabl
AUTHOR
Seiichi Manyama, Apr 12 2025
STATUS
approved
