OFFSET
7,2
LINKS
Marius A. Burtea, Table of n, a(n) for n = 7..5000
N. J. A. Sloane, Transforms
FORMULA
Moebius transform of C(n-1,6).
G.f.: Sum_{k>=1} mu(k) * x^(7*k) / (1 - x^k)^7. - Ilya Gutkovskiy, Feb 05 2020
MAPLE
with(numtheory):
a:= n-> add(mobius(n/d)*binomial(d-1, 6), d=divisors(n)):
seq(a(n), n=7..50); # Alois P. Heinz, Feb 05 2020
MATHEMATICA
a[n_]:=DivisorSum[n, Binomial[#-1, 6] MoebiusMu[n/#]&]; Array[a, 37, 7] (* or *) a[n_]:=Sum[Boole[Divisible[n, k]] MoebiusMu[n/k] Binomial[k-1, 6], {k, 1, n}]; Table[a[n], {n, 7, 45}] (* Vincenzo Librandi, Feb 07 2020 *)
PROG
(Magma) [&+[MoebiusMu(n div d)*Binomial(d-1, 6):d in Divisors(n)]:n in[7..41]]; // Marius A. Burtea, Feb 07 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved