OFFSET
1,3
COMMENTS
FORMULA
Inverse Moebius transform of an infinite lower triangular matrix with A063659, (1, 2, 3, 3, 5, 6, 7, 6, 8, 10, ...) in the main diagonal and the rest zeros.
EXAMPLE
First few rows of the triangle:
1;
1, 2;
1, 0, 3;
1, 2, 0, 3;
1, 0, 0, 0, 5;
1, 2, 3, 0, 0, 6;
1, 0, 0, 0, 0, 0, 7,
1, 2, 0, 3, 0, 0, 0, 6;
1, 0, 3, 0, 0, 0, 0, 0, 8;
...
MATHEMATICA
m = 14;
A051731 = Table[If[Mod[n, k] == 0, 1, 0], {n, m}, {k, m}];
A063659 = Table[Sum[MoebiusMu[GCD[n, k]]^2, {k, n}], {n, m}] // DiagonalMatrix;
Table[M[[n, k]], {n, m}, {k, n}] // Flatten (* Jean-François Alcover, Jan 18 2020 *)
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Gary W. Adamson, May 07 2007
EXTENSIONS
More terms from Jean-François Alcover, Jan 18 2020
STATUS
approved