OFFSET
0
LINKS
August Ferdinand Möbius, Über eine besondere Art von Umkehrung der Reihen. Journal für die reine und angewandte Mathematik 9 (1832), 105-123.
FORMULA
T(n, k) = M(n, k)^(-1), where M(n, k) = [k <= n and k divides n], '(-1)' denotes matrix inversion, and '[ ]' denotes the Iverson bracket.
T(n, k) = k^n if k = 0, otherwise Moebius(n/k) if k divides n, otherwise 0.
Sum_{k=0..n} k*T(n, k) = phi(n) = A000010(n).
EXAMPLE
Triangle T(n, k) starts:
[0] 1;
[1] 0, 1;
[2] 0, -1, 1;
[3] 0, -1, 0, 1;
[4] 0, 0, -1, 0, 1;
[5] 0, -1, 0, 0, 0, 1;
[6] 0, 1, -1, -1, 0, 0, 1;
[7] 0, -1, 0, 0, 0, 0, 0, 1;
[8] 0, 0, 0, 0, -1, 0, 0, 0, 1;
[9] 0, 0, 0, -1, 0, 0, 0, 0, 0, 1;
MAPLE
A363914 := (n, k) -> ifelse(k = 0, k^n, ifelse(irem(n, k) = 0, NumberTheory:-mu(n/k), 0)): for n from 0 to 9 do seq(A363914(n, k), k = 0..n) od;
# By inverting the Moebius matrix:
divides := (k, n) -> ifelse(k > n, 0, ifelse(k = n or (k > 0 and irem(n, k) = 0), 1, 0)): M := Matrix(10, (n, k) -> divides(k - 1, n - 1)):
# The shift in k, n is necessary because Maple's 'Matrix' is (1, 1)-based.
LinearAlgebra:-MatrixInverse(M);
PROG
CROSSREFS
KEYWORD
sign,tabl
AUTHOR
Peter Luschny, Jul 01 2023
EXTENSIONS
Name edited by Peter Luschny, Jul 29 2023
STATUS
approved