Clifford algebras and approximating the permanent

(MATH) We study approximation algorithms for the permanent of an n x n (0,1) matrix A based on the following simple idea: obtain a random matrix B by replacing each 1-entry of A independently by ± e, where e is a random basis element of a suitable algebra; then output |det(B)|2. This estimator is always unbiased, but it may have exponentially large variance. In our first main result we show that, if we take the algebra to be a Clifford algebra of dimension polynomial in n, then we get an estimator with small variance. Hence only a constant number of trials suffices to estimate the permanent to good accuracy. The idea of using Clifford algebras is a natural extension of earlier work by Godsil and Gutman, Karmarkar et al., and Barvinok, who used the real numbers, complex numbers and quaternions respectively.(MATH) The above result implies that, in principle, this approach gives a fully-polynomial randomized approximation scheme for the permanent, provided |det(B)|2 can be efficiently computed in the Clifford algebras. Since these algebras are non-commutative it is not clear how to do this. However, our second main result shows how to compute in polynomial time an estimator with the same mean and variance over the 4-dimensional algebra (which is the quaternions, and is non-commutative); in addition to providing some hope that the computations can be performed in higher dimensions, this quaternion algorithm provides an exponential improvement in the variance over that of the 2-dimensional complex version studied by Karmarkar et al.

[1]  Richard J. Lipton,et al.  A Monte-Carlo Algorithm for Estimating the Permanent , 1993, SIAM J. Comput..

[2]  Lars Eilstrup Rasmussen,et al.  Approximating the Permanent: A Simple Approach , 1994, Random Struct. Algorithms.

[3]  Mark Jerrum,et al.  Approximating the Permanent , 1989, SIAM J. Comput..

[4]  Alex Samorodnitsky,et al.  A Deterministic Strongly Polynomial Algorithm for Matrix Scaling and Approximate Permanents , 2000, Comb..

[5]  Andrei Z. Broder,et al.  How hard is it to marry at random? (On the approximation of the permanent) , 1986, STOC '86.

[6]  Jeremy J. Gray,et al.  Years ago , 1996 .

[7]  Manfred Knebusch,et al.  The algebraic theory of quadratic forms , 1980 .

[8]  Helmer Aslaksen,et al.  Quaternionic determinants , 1996 .

[9]  Alex Samorodnitsky,et al.  A Deterministic Strongly Polynomial Algorithm for Matrix Scaling and Approximate Permanents , 1998, STOC '98.

[10]  Alexander Barvinok,et al.  New Permanent Estimators via Non-Commutative Determinants , 2000 .

[11]  A. Sinclair,et al.  On approximating the permanent and other #p-complete problems , 1998 .

[12]  Noam Nisan,et al.  Lower bounds for non-commutative computation , 1991, STOC '91.

[13]  Leslie G. Valiant,et al.  The Complexity of Computing the Permanent , 1979, Theor. Comput. Sci..

[14]  Eric Vigoda,et al.  A polynomial-time approximation algorithm for the permanent of a matrix with non-negative entries , 2001, STOC '01.

[15]  I. Gel'fand,et al.  A theory of noncommutative determinants and characteristic functions of graphs , 1992 .

[16]  A. Sinclair,et al.  Approximating the number of monomer-dimer coverings of a lattice , 1996 .

[17]  Mark Jerrum An analysis of a Monte Carlo algorithm for estimating the permanent , 1993, IPCO.

[18]  Alexander I. Barvinok,et al.  Polynomial Time Algorithms to Approximate Permanents and Mixed Discriminants Within a Simply Exponential Factor , 1999, Random Struct. Algorithms.