Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1212.0025 · Electronic colloquium on computational complexity

Generalizing and Derandomizing Gurvits's Approximation Algorithm for the Permanent

Scott Aaronson, Travis Hance

Published Nov 30, 2012• 4 citations
arXivHTMLPDF

Abstract

Around 2002, Leonid Gurvits gave a striking randomized algorithm to approximate the permanent of an n*n matrix A. The algorithm runs in O(n^2/eps^2) time, and approximates Per(A) to within eps*||A||^n additive error. A major advantage of Gurvits's algorithm is that it works for arbitrary matrices, not just for nonnegative matrices. This makes it highly relevant to quantum optics, where the permanents of bounded-norm complex matrices play a central role. Indeed, the existence of Gurvits's algorithm is why, in their recent work on the hardness of quantum optics, Aaronson and Arkhipov (AA) had to talk about sampling problems rather than estimation problems. In this paper, we improve Gurvits's algorithm in two ways. First, using an idea from quantum optics, we generalize the algorithm so that it yields a better approximation when the matrix A has either repeated rows or repeated columns. Translating back to quantum optics, this lets us classically estimate the probability of any outcome of an AA-type experiment---even an outcome involving multiple photons "bunched" in the same mode---at least as well as that probability can be estimated by the experiment itself. (This does not, of course, let us solve the AA sampling problem.) It also yields a general upper bound on the probabilities of "bunched" outcomes, which resolves a conjecture of Gurvits and might be of independent physical interest. Second, we use eps-biased sets to derandomize Gurvits's algorithm, in the special case where the matrix A is nonnegative. More interestingly, we generalize the notion of eps-biased sets to the complex numbers, construct "complex eps-biased sets," then use those sets to derandomize even our generalization of Gurvits's algorithm to the multirow/multicolumn case (again for nonnegative A). Whether Gurvits's algorithm can be derandomized for general A remains an outstanding problem.

DOI
10.48550/arxiv.1212.0025

Related papers

Generalizing and derandomizing Gurvits's approximation algorithm for the permanent

Scott Aaronson, Travis Hance

Around 2002, Leonid Gurvits gave a striking randomized algorithm to approximate the permanent of an $n\times n$ matrix A. The algorithm runs in $O( n^{2}/\varepsilon^{2})$ time, and approximates $\operatorname*{Per}( A)$ to within $\pm\varepsilon\left\Vert A\right\Vert ^{n}$ additive error. A major advantage of Gurvits's algorithm is that it works for arbitrary matrices, not just for nonnegative matrices. \ This makes it highly relevant to quantum optics, where the permanents of bounded-norm complex matrices play a central role. (In particular, $n\times n$ permanents arise as the transition amplitudes for $n$ identical photons.) Indeed, the existence of Gurvits's algorithm is why, in their recent work on the \textit{hardness} of quantum optics, Aaronson and Arkhipov (AA) had to talk about sampling problems rather than ordinary decision problems. In this paper, we improve Gurvits's algorithm in two ways. First, using an idea from quantum optics, we generalize the algorithm so that it yields a better approximation when the matrix A has either repeated rows or repeated columns. Translating back to quantum optics, this lets us classically estimate the probability of any outcome of an AA-type experiment -- even an outcome involving multiple photons "bunched" in the same mode -- at least as well as that probability can be estimated by the experiment itself. It also yields a general upper bound on the probabilities of "bunched" outcomes, which resolves a conjecture of Gurvits and might be of independent physical interest. Second, we use $\varepsilon$-biased sets to derandomize Gurvits's algorithm, in the special case where the matrix A is nonnegative. More interestingly, we generalize the notion of $\varepsilon$-biased sets to the complex numbers, construct "complex $\varepsilon$-biased sets", then use those sets to derandomize even our generalization of Gurvits's algorithm to the case (again for nonnegative A) where some rows or columns are identical. Whether Gurvits's algorithm can be derandomized for general A remains an outstanding problem.

Quantum Information and Computation•May 1, 2014•20 citations

Approximating the Permanent of a Random Matrix with Polynomially Small Mean: Zeros and Universality

Frederic Koehler, Pui Kuen Leung

We study algorithms for approximating the permanent of a random matrix when the entries are slightly biased away from zero. This question is motivated by the goal of understanding the classical complexity of linear optics and \emph{boson sampling} (Aaronson and Arkhipov '11; Eldar and Mehraban '17). Barvinok's interpolation method enables efficient approximation of the permanent, provided one can establish a sufficiently large zero-free region for the polynomial $\mathrm{per}(zJ + W)$, where $J$ is the all-ones matrix and $W$ is a random matrix with independent mean-zero entries. We show that when the entries of $W$ are standard complex Gaussians, all zeros of the random polynomial $\mathrm{per}(zJ + W)$ lie within a disk of radius $\tilde{O}(n^{-1/3})$, which yields an approximation algorithm when the bias of the entries is $\tildeΩ(n^{-1/3})$. Previously, there were no efficient algorithms at biases smaller than $1/\mathrm{polylog}(n)$, and it was unknown whether there typically exist zeros $z$ with $|z| \ge 1$. As a complementary result, we show that the bulk of the zeros, namely $(1 - ε)n$ of them, have magnitude $Θ(n^{-1/2})$. This prevents our interpolation method from contradicting the conjectured average-case hardness of approximating the permanent. We also establish analogous zero-free regions for the hardcore model on general graphs with complex vertex fugacities. In addition, we prove universality results establishing zero-free regions for random matrices $W$ with i.i.d. subexponential entries.

cs.DSCross-listedarXiv:2604.01367•Apr 3, 2026•0 citations

Computing matrix permanent with collective boson operators

Joonsuk Huh

Computing permanents of matrices are known to be a classically hard problem that the computational cost grows exponentially with the size of the matrix increases. So far, there exist a few classical algorithms to compute the matrix permanents in deterministic and in randomized ways. By exploiting the series expansion of products of boson operators regarding collective boson operators, a generalized algorithm for computing permanents is developed that the algorithm can handle the arbitrary matrices with repeated columns and rows. In a particular case, the formula is reduced to Glynn's form. Not only the algorithm can be used for a deterministic direct calculation of the matrix permanent but also can be expressed as a sampling problem like Gurvits's randomized algorithm.

•May 27, 2016•0 citations

Quantum estimation bound of the Gaussian matrix permanent

Joonsuk Huh

Exact calculation and even multiplicative error estimation of matrix permanent are challenging for both classical and quantum computers. Regarding the permanents of random Gaussian matrices, the additive error estimation is closely linked to boson sampling and achieving multiplicative error estimation requires exponentially many samplings. Our formula for matrix permanents and its corresponding quantum expression have enabled better estimation of the average additive error for random Gaussian matrices compared to Gurvits' classical sampling algorithm. The well-known Ryser formula has been converted into a quantum permanent estimator. When dealing with real random Gaussian square matrices of size $N$, the quantum estimator can approximate the matrix permanent with an additive error smaller than $\ensuremath{\epsilon}{(\sqrt{\mathrm{e}N})}^{N}$, where $\ensuremath{\epsilon}$ is the estimation precision. In contrast, Gurvits' classical sampling algorithm has an estimation error of $\ensuremath{\epsilon}{(2\sqrt{N})}^{N}$, which is exponentially larger ($1.{2}^{N}$) than the quantum method. As expected, the quantum additive error bound fails to reach the multiplicative error bound of ${(2\ensuremath{\pi}N)}^{1/4}\ensuremath{\epsilon}{(\sqrt{N/\mathrm{e}})}^{N}$. Additionally, the quantum permanent estimator can be up to quadratically faster than the classical estimator when using quantum phase estimation-based amplitude estimation.

arXiv:2205.01328•Physical Review A•Jan 13, 2025•6 citations

Weak Permanent Anti-Concentration for Random Gaussian Matrices in Boson Sampling

Fei Meng, Bin Cheng, Jianan Li +1 more

Recent demonstrations of quantum computational advantage have been driven largely by sampling problems. A prominent model, boson sampling, involves sampling from the output distribution of a linear optical network. However, its classical hardness hinges on two plausible yet less-studied conjectures: the average-case hardness of approximating Gaussian permanents, and the permanent anti-concentration conjecture (PACC). The PACC is a purely mathematical assertion regarding the distributional properties of random Gaussian matrices. While the typical magnitude of the permanent has been established for discrete random matrices, the complex Gaussian case, which governs transition amplitudes in linear optical networks, has remained open. Here, we establish a weak anti-concentration bound by upper-bounding the probability that a random Gaussian permanent is superexponentially smaller than its standard deviation. Tightening this bound to an inverse-polynomial fraction would prove the original PACC. As a corollary, we establish the typical magnitude of Gaussian permanents, on par with Tao and Vu's seminal result for Bernoulli matrices. Combined with the Aaronson-Arkhipov framework, our result implies that classically simulating boson sampling to within a superexponentially small total variation distance would collapse the polynomial hierarchy, assuming the remaining conjectures hold.

quant-pharXiv:2607.22088•Jul 27, 2026•0 citations

Approximating Permanent of Random Matrices with Vanishing Mean: Made Better and Simpler

Zhengfeng Ji, Zhihan Jin, Pinyan Lu

The algorithm and complexity of approximating the permanent of a matrix is an extensively studied topic. Recently, its connection with quantum supremacy and more specifically BosonSampling draws a special attention to the average-case approximation problem of the permanent of random matrices with zero or small mean value for each entry. Eldar and Mehraban (FOCS 2018) gave a quasi-polynomial time algorithm for random matrices with mean at least 1/polyloglog(n). In this paper, we improve the result by designing a deterministic quasi-polynomial time algorithm and a PTAS for random matrices whose module of mean is at least 1/ polylog(n). We note that if the algorithm can be further improved to work with a mean value that is a sufficiently small 1/poly(n), it will disprove a central conjecture for quantum supremacy. Our algorithm is also much simpler and has a better and flexible trade-off for running time. The running time can be quasi-polynomial in both n and 1/∊, or PTAS (polynomial in n but exponential in 1/∊), where ∊ is the approximation parameter.

arXiv:1911.11962•Society for Industrial and Applied Mathematics eBooks•Jan 1, 2021•2 citations
Around 2002, Leonid Gurvits gave a striking randomized algorithm to approximate the permanent of an nn matrix A. The algorithm runs in O(n^2/eps^2) time, and approximates Per(A) to within eps||A||^n additive error. A major advantage of Gurvits's algorithm is that it works for arbitrary matrices, not just for nonnegative matrices. This makes it highly relevant to quantum optics, where the permanents of bounded-norm complex matrices play a central role. Indeed, the existence of Gurvits's algorithm is why, in their recent work on the hardness of quantum optics, Aaronson and Arkhipov (AA) had to talk about sampling problems rather than estimation problems. In this paper, we improve Gurvits's algorithm in two ways. First, using an idea from quantum optics, we generalize the algorithm so that it yields a better approximation when the matrix A has either repeated rows or repeated columns. Translating back to quantum optics, this lets us classically estimate the probability of any outcome of an AA-type experiment---even an outcome involving multiple photons "bunched" in the same mode---at least as well as that probability can be estimated by the experiment itself. (This does not, of course, let us solve the AA sampling problem.) It also yields a general upper bound on the probabilities of "bunched" outcomes, which resolves a conjecture of Gurvits and might be of independent physical interest. Second, we use eps-biased sets to derandomize Gurvits's algorithm, in the special case where the matrix A is nonnegative. More interestingly, we generalize the notion of eps-biased sets to the complex numbers, construct "complex eps-biased sets," then use those sets to derandomize even our generalization of Gurvits's algorithm to the multirow/multicolumn case (again for nonnegative A). Whether Gurvits's algorithm can be derandomized for general A remains an outstanding problem.
Around 2002, Leonid Gurvits gave a striking randomized algorithm to approximate the permanent of an n×nn\times nn×n matrix A. The algorithm runs in O(n2/ε2)O( n^{2}/\varepsilon^{2})O(n2/ε2) time, and approximates Per⁡(A)\operatorname*{Per}( A)Per(A) to within ±ε∥A∥n\pm\varepsilon\left\Vert A\right\Vert ^{n}±ε∥A∥n additive error. A major advantage of Gurvits's algorithm is that it works for arbitrary matrices, not just for nonnegative matrices. \ This makes it highly relevant to quantum optics, where the permanents of bounded-norm complex matrices play a central role. (In particular, n×nn\times nn×n permanents arise as the transition amplitudes for nnn identical photons.) Indeed, the existence of Gurvits's algorithm is why, in their recent work on the \textit{hardness} of quantum optics, Aaronson and Arkhipov (AA) had to talk about sampling problems rather than ordinary decision problems. In this paper, we improve Gurvits's algorithm in two ways. First, using an idea from quantum optics, we generalize the algorithm so that it yields a better approximation when the matrix A has either repeated rows or repeated columns. Translating back to quantum optics, this lets us classically estimate the probability of any outcome of an AA-type experiment -- even an outcome involving multiple photons "bunched" in the same mode -- at least as well as that probability can be estimated by the experiment itself. It also yields a general upper bound on the probabilities of "bunched" outcomes, which resolves a conjecture of Gurvits and might be of independent physical interest. Second, we use ε\varepsilonε-biased sets to derandomize Gurvits's algorithm, in the special case where the matrix A is nonnegative. More interestingly, we generalize the notion of ε\varepsilonε-biased sets to the complex numbers, construct "complex ε\varepsilonε-biased sets", then use those sets to derandomize even our generalization of Gurvits's algorithm to the case (again for nonnegative A) where some rows or columns are identical. Whether Gurvits's algorithm can be derandomized for general A remains an outstanding problem.
We study algorithms for approximating the permanent of a random matrix when the entries are slightly biased away from zero. This question is motivated by the goal of understanding the classical complexity of linear optics and \emph{boson sampling} (Aaronson and Arkhipov '11; Eldar and Mehraban '17). Barvinok's interpolation method enables efficient approximation of the permanent, provided one can establish a sufficiently large zero-free region for the polynomial per(zJ+W)\mathrm{per}(zJ + W)per(zJ+W), where JJJ is the all-ones matrix and WWW is a random matrix with independent mean-zero entries. We show that when the entries of WWW are standard complex Gaussians, all zeros of the random polynomial per(zJ+W)\mathrm{per}(zJ + W)per(zJ+W) lie within a disk of radius O~(n−1/3)\tilde{O}(n^{-1/3})O~(n−1/3), which yields an approximation algorithm when the bias of the entries is Ω~(n−1/3)\tildeΩ(n^{-1/3})Ω~(n−1/3). Previously, there were no efficient algorithms at biases smaller than 1/polylog(n)1/\mathrm{polylog}(n)1/polylog(n), and it was unknown whether there typically exist zeros zzz with ∣z∣≥1|z| \ge 1∣z∣≥1. As a complementary result, we show that the bulk of the zeros, namely (1−ε)n(1 - ε)n(1−ε)n of them, have magnitude Θ(n−1/2)Θ(n^{-1/2})Θ(n−1/2). This prevents our interpolation method from contradicting the conjectured average-case hardness of approximating the permanent. We also establish analogous zero-free regions for the hardcore model on general graphs with complex vertex fugacities. In addition, we prove universality results establishing zero-free regions for random matrices WWW with i.i.d. subexponential entries.
Computing permanents of matrices are known to be a classically hard problem that the computational cost grows exponentially with the size of the matrix increases. So far, there exist a few classical algorithms to compute the matrix permanents in deterministic and in randomized ways. By exploiting the series expansion of products of boson operators regarding collective boson operators, a generalized algorithm for computing permanents is developed that the algorithm can handle the arbitrary matrices with repeated columns and rows. In a particular case, the formula is reduced to Glynn's form. Not only the algorithm can be used for a deterministic direct calculation of the matrix permanent but also can be expressed as a sampling problem like Gurvits's randomized algorithm.
Exact calculation and even multiplicative error estimation of matrix permanent are challenging for both classical and quantum computers. Regarding the permanents of random Gaussian matrices, the additive error estimation is closely linked to boson sampling and achieving multiplicative error estimation requires exponentially many samplings. Our formula for matrix permanents and its corresponding quantum expression have enabled better estimation of the average additive error for random Gaussian matrices compared to Gurvits' classical sampling algorithm. The well-known Ryser formula has been converted into a quantum permanent estimator. When dealing with real random Gaussian square matrices of size NNN, the quantum estimator can approximate the matrix permanent with an additive error smaller than $\ensuremath{\epsilon}{(\sqrt{\mathrm{e}N})}^{N}$, where $\ensuremath{\epsilon}$ is the estimation precision. In contrast, Gurvits' classical sampling algorithm has an estimation error of $\ensuremath{\epsilon}{(2\sqrt{N})}^{N}$, which is exponentially larger (1.2N1.{2}^{N}1.2N) than the quantum method. As expected, the quantum additive error bound fails to reach the multiplicative error bound of ${(2\ensuremath{\pi}N)}^{1/4}\ensuremath{\epsilon}{(\sqrt{N/\mathrm{e}})}^{N}$. Additionally, the quantum permanent estimator can be up to quadratically faster than the classical estimator when using quantum phase estimation-based amplitude estimation.
Recent demonstrations of quantum computational advantage have been driven largely by sampling problems. A prominent model, boson sampling, involves sampling from the output distribution of a linear optical network. However, its classical hardness hinges on two plausible yet less-studied conjectures: the average-case hardness of approximating Gaussian permanents, and the permanent anti-concentration conjecture (PACC). The PACC is a purely mathematical assertion regarding the distributional properties of random Gaussian matrices. While the typical magnitude of the permanent has been established for discrete random matrices, the complex Gaussian case, which governs transition amplitudes in linear optical networks, has remained open. Here, we establish a weak anti-concentration bound by upper-bounding the probability that a random Gaussian permanent is superexponentially smaller than its standard deviation. Tightening this bound to an inverse-polynomial fraction would prove the original PACC. As a corollary, we establish the typical magnitude of Gaussian permanents, on par with Tao and Vu's seminal result for Bernoulli matrices. Combined with the Aaronson-Arkhipov framework, our result implies that classically simulating boson sampling to within a superexponentially small total variation distance would collapse the polynomial hierarchy, assuming the remaining conjectures hold.
The algorithm and complexity of approximating the permanent of a matrix is an extensively studied topic. Recently, its connection with quantum supremacy and more specifically BosonSampling draws a special attention to the average-case approximation problem of the permanent of random matrices with zero or small mean value for each entry. Eldar and Mehraban (FOCS 2018) gave a quasi-polynomial time algorithm for random matrices with mean at least 1/polyloglog(n). In this paper, we improve the result by designing a deterministic quasi-polynomial time algorithm and a PTAS for random matrices whose module of mean is at least 1/ polylog(n). We note that if the algorithm can be further improved to work with a mean value that is a sufficiently small 1/poly(n), it will disprove a central conjecture for quantum supremacy. Our algorithm is also much simpler and has a better and flexible trade-off for running time. The running time can be quasi-polynomial in both n and 1/∊, or PTAS (polynomial in n but exponential in 1/∊), where ∊ is the approximation parameter.