n*a(n) is the number of n-member subsets of {1,2,3,...,2*n-1} that sum to 1 mod n, cf. A145855. - Vladeta Jovovic, Oct 28 2008
a(n) is the number of orbits under the S_n action on a set closely related to the set of parking functions. See Konvalinka-Tewari reference below. - Vasu Tewari, Mar 17 2020
A131868 := proc(n) local a, d ; a := 0 ; for d in numtheory[divisors](n) do a := a+(-1)^(n+d)*numtheory[mobius](n/d)*binomial(2*d, d) ; od: a/2/n^2 ; end: seq(A131868(n), n=1..30) ; # R. J. Mathar, Oct 24 2007
MATHEMATICA
a = {}; For[n = 1, n < 30, n++, b = Divisors[n]; s = 0; For[j = 1, j < Length[b] + 1, j++, s = s + (-1)^(n + b[[j]])*MoebiusMu[n/b[[j]]]* Binomial[2*b[[j]], b[[j]]]]; AppendTo[a, s/(2*n^2)]]; a (* Stefan Steinerberger, Oct 26 2007 *)