OFFSET
1,7
COMMENTS
This should be related to the Coxeter transformation for the Tamari lattices.
The source sequence is 1, -1, -2, 3, 6, -10, -20, 35, 70, -126, ... (A001405). Its Mobius transform is 1, -2, -3, 4, 5, -6, -21, 32, 72, -130, -253, 468, 923, ... and division of each term through n generates a(n). - R. J. Mathar, Jul 23 2012
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
F. Chapoton, Le module dendriforme sur le groupe cyclique, Ann. Inst. Fourier (Grenoble) 58 (2008), no. 7, 2333-2350. In French.
EXAMPLE
G.f. = x - x^2 - x^3 + x^4 + x^5 - x^6 - 3*x^7 + 4*x^8 + 8*x^9 - 13*x^10 + ...
MAPLE
with(numtheory):
a:= n-> add(mobius(n/d)*[1$2, -1$2][1+irem(d, 4)]*
binomial(d-1, iquo(d-1, 2)), d=divisors(n))/n:
seq(a(n), n=1..50); # Alois P. Heinz, Apr 05 2013
MATHEMATICA
a[ n_] := If[ n < 1, 0, DivisorSum[ n, MoebiusMu[ n/#] (-1)^Quotient[ #, 2] Binomial[ # - 1, Quotient[ # - 1, 2]] &] / n]; (* Michael Somos, Sep 14 2015 *)
PROG
(Sage)
def lam(n):
return (-1)**binomial(n, 2) * binomial(n - 1, (n - 1) // 2)
def a(n):
return sum(moebius(n // d) * lam(d) for d in divisors(n)) // n
[a(n) for n in range(1, 20)]
(PARI) {a(n) = if( n<1, 0, sumdiv( n, d, moebius(n/d) * (-1)^(d\2) * binomial(d-1, (d-1)\2)) / n)}; /* Michael Somos, Dec 23 2014 */
CROSSREFS
KEYWORD
sign
AUTHOR
F. Chapoton, Jun 09 2010
STATUS
approved