OFFSET
0,3
COMMENTS
a(n) is also the number of subspaces of R^n given by coordinate equalities of the form x_i = x_j, x_i = -x_j and x_i = 0, that are orthogonal to the vector of all 1's.
These subspaces play an important role in the field of network dynamical systems, where they correspond to anti-synchronization. That is, they capture the phenomenon where different cells in a network show the same or opposite behavior.
LINKS
James W. Swift, Table of n, a(n) for n = 0..500
Eddie Nijholt, Nándor Sieben, and James W. Swift, Invariant Synchrony and Anti-Synchrony Subspaces of Weighted Networks, arXiv:2206.00094 [math.DS], 2022.
FORMULA
E.g.f.: exp(I_0(2z)/2 - 1/2 + z), where I_0 is the modified Bessel function of the first kind.
a(0) = 1, a(n+1) = a(n) + Sum_{k>=0, m>=0, k + 2m + 1 = n} (n!/(k!*m!*(n-k-m)!)*a(k).
EXAMPLE
The set {1,2} has the trivial partition consisting only of the class {1,2}, together with the involution that sends {1,2} to itself. There is also the singleton partition with classes {1} and {2}, together with the involution that maps {1} to {2} and vice versa. The other involution on the two classes {1} and {2} has two fixed points and is therefore not counted, Hence, we find a(2) = 2.
Alternatively, in R^2 we have two subspaces given by coordinate equalities of the form x_i = x_j, x_i = -x_j and x_i = 0, that are orthogonal to the vector (1,1). These are the zero-space and the subspace given by x_1 = -x_2. In R^3 we find the zero-space and the three subspaces given by x_i = -x_j, x_k = 0 for {i,j,k} = {1,2,3}. This shows that a(3) = 4.
MATHEMATICA
nMax = 20; CoefficientList[ Series[Exp[BesselI[0, 2 x]/2 - 1/2 + x], {x, 0, nMax}], x] * Range[0, nMax]!
CROSSREFS
KEYWORD
nonn
AUTHOR
Eddie Nijholt, Jun 23 2022
STATUS
approved