OFFSET
0,4
LINKS
FORMULA
Moebius transform of A005648. - Christian G. Bower
a(n) = 1/(4*n) * Sum_{d|n} mu(n/d) * (2*n*binomial(2*floor(d/2), floor(d/2)) + binomial(2*d, d)) for n > 0. - Andrew Howroyd, Sep 28 2017
MATHEMATICA
a[n_] := If[n == 0, 1, Sum[MoebiusMu[n/d] (2n Binomial[ 2 Quotient[d, 2], Quotient[d, 2]] + Binomial[2d, d]), {d, Divisors[n]}]/(4n)];
Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Aug 28 2019, from PARI *)
PROG
(PARI) a(n) = if(n<1, n==0, sumdiv(n, d, moebius(n/d)*(2*n*binomial(2*(d\2), d\2) + binomial(2*d, d)))/(4*n)); \\ Andrew Howroyd, Sep 28 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved