OFFSET
1,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1048
FORMULA
a(n) = Sum_{j=1..n} Sum_{d|j} 9^(d-1) * mu(j/d).
a(n) = A143327(n,9).
a(n) = Sum_{j=1..n} A143325(j,9).
a(n) = A143326(n,9) / 9.
G.f.: (1/(1 - x)) * Sum_{k>=1} mu(k) * x^k / (1 - 9*x^k). - Ilya Gutkovskiy, Dec 11 2020
MAPLE
b:= n-> add(`if`(d=n, 9^(n-1), -b(d)), d=numtheory[divisors](n)):
a:= proc(n) option remember; b(n)+`if`(n<2, 0, a(n-1)) end:
seq(a(n), n=1..30);
PROG
(PARI) a(n) = sum(j=1, n, sumdiv(j, d, 9^(d-1)*moebius(j/d))); \\ Michel Marcus, Dec 11 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 05 2018
STATUS
approved