A Fast Fourier Transform for the Johnson graph

Abstract

The set X of k-subsets of an n-set has a natural graph structure where two k-subsets are connected if and only if the size of their intersection is k-1. This is known as the Johnson graph. The symmetric group Sn acts on the space of complex functions on X and this space has a multiplicity-free decomposition as sum of irreducible representations of Sn, so it has a well-defined Gelfand-Tsetlin basis up to scalars. The Fourier transform on the Johnson graph is defined as the change of basis matrix from the delta function basis to the Gelfand-Tsetlin basis. The direct application of this matrix to a generic vector requires nk2 arithmetic operations. We show that this matrix can be factorized as a product of n-1 orthogonal matrices, each one with at most two nonzero elements in each column. The factorization is based on the construction of n-1 intermediate bases which are parametrized via the Robinson-Schensted insertion algorithm. This factorization shows that the number of arithmetic operations required to apply this matrix to a generic vector is bounded above by 2(n-1) nk. We show that each one of these sparse matrices can be constructed using O(nk) arithmetic operations. Our construction does not depend on numerical methods. Instead, they are obtained by solving small linear systems with integer coefficients derived from the Jucys-Murphy operators. Then both the construction and the succesive application of all these n-1 matrices can be performed using O(n nk) operations. As a consequence, we show that the problem of computing all the weights of the isotypic components of a given function can be solved in O(n nk) operations, improving the previous bound O(k2 nk) when k asymptotically dominates n.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…