OFFSET
0,8
COMMENTS
The name expresses the 'row view' of the array. The 'column view' regards the array as the collection of the inverse Möbius transforms of the power sequences k^n = 0^n, 1^n, 2^n, .... (n >= 0). Viewed this way, the array is a generalization of the number of divisors sequence tau (A000005), to which it reduces in the case k = 1.
The array has offset (0, 0). It uses the usual definition of 'k divides n' as described in Apostol, rather than the shortened version, which restricts to values k > 0 as some programs do (but not SageMath). Such a restriction makes sense in the context of rational numbers but not in the case of natural numbers.
REFERENCES
Tom M. Apostol, Introduction to Analytic Number Theory, Springer 1976, p. 14.
FORMULA
A(n, k) = Sum_{j=0..n} divides(j, n) * k^j, where divides(k, n) <-> [k = n or (k > 0 and n mod k = 0)], and '[ ]' denotes the Iverson bracket.
The columns are the inverse Möbius transforms of the powers x^n, x >= 0.
EXAMPLE
Array A(n, k) starts:
[0] 1, 1, 1, 1, 1, 1, 1, 1, 1, ... A000012
[1] 0, 1, 2, 3, 4, 5, 6, 7, 8, ... A001477
[2] 0, 2, 6, 12, 20, 30, 42, 56, 72, ... A002378
[3] 0, 2, 10, 30, 68, 130, 222, 350, 520, ... A034262
[4] 0, 3, 22, 93, 276, 655, 1338, 2457, 4168, ...
[5] 0, 2, 34, 246, 1028, 3130, 7782, 16814, 32776, ...
[6] 0, 4, 78, 768, 4180, 15780, 46914, 118048, 262728, ...
[7] 0, 2, 130, 2190, 16388, 78130, 279942, 823550, 2097160, ...
[8] 0, 4, 278, 6654, 65812, 391280, 1680954, 5767258, 16781384, ...
.
Triangle T(n, k) starts:
[0] 1;
[1] 0, 1;
[2] 0, 1, 1;
[3] 0, 2, 2, 1;
[4] 0, 2, 6, 3, 1;
[5] 0, 3, 10, 12, 4, 1;
[6] 0, 2, 22, 30, 20, 5, 1;
[7] 0, 4, 34, 93, 68, 30, 6, 1;
[8] 0, 2, 78, 246, 276, 130, 42, 7, 1;
[9] 0, 4, 130, 768, 1028, 655, 222, 56, 8, 1;
MAPLE
divides := (k, n) -> ifelse(k = n or (k > 0 and irem(n, k) = 0), 1, 0):
A := (n, k) -> local j; add(divides(j, n) * k^j, j = 0 ..n):
for n from 0 to 8 do seq(A(n, k), k = 0..8) od;
# If we introduce the 'inverse Möbius transform' InvMoebius acting on s ...
InvMoebius := (s, n) -> local j; add(divides(j, n) * s(j), j = 0 ..n):
# ... the transposed array is given by applying InvMoebius to the powers r^m:
seq(lprint(seq(InvMoebius(m -> r^m, n), n = 0..8)), r = 0..8);
# For instance we see that the number of divisors is the inverse
# Moebius transform of the constant sequence s = 1.
PROG
(SageMath)
def A(n, k): return sum(j.divides(n) * k^j for j in (0..n))
for n in srange(9): print([A(n, k) for k in (0..8)])
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Peter Luschny, Jun 27 2023
STATUS
approved