OFFSET
9,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 9..1000
MAPLE
with(numtheory):
b:= proc(n, k) option remember; `if`(n=0, 1,
add(mobius(n/d)*k^d, d=divisors(n))/n)
end:
a:= n-> add(b(n, 9-j)*binomial(9, j)*(-1)^j, j=0..9):
seq(a(n), n=9..30);
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jan 25 2015
STATUS
approved