OFFSET
0,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1047 (terms 0..350 from Alois P. Heinz)
FORMULA
a(n) = Sum_{d|n} 9^d * mu(n/d) for n>0, a(0) = 1.
G.f.: 1 + 9 * Sum_{k>=1} mu(k) * x^k / (1 - 9*x^k). - Ilya Gutkovskiy, Apr 15 2021
MAPLE
with(numtheory):
a:= n-> `if`(n=0, 1, add(9^d*mobius(n/d), d=divisors(n))):
seq(a(n), n=0..30);
PROG
(PARI) a(n) = if (n==0, 1, sumdiv(n, d, 9^d*moebius(n/d))); \\ Michel Marcus, Apr 15 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 21 2012
STATUS
approved