258 research outputs found
Bisection of Bounded Treewidth Graphs by Convolutions
In the Bisection problem, we are given as input an edge-weighted graph G. The task is to find a partition of V(G) into two parts A and B such that ||A| - |B|| <= 1 and the sum of the weights of the edges with one endpoint in A and the other in B is minimized. We show that the complexity of the Bisection problem on trees, and more generally on graphs of bounded treewidth, is intimately linked to the (min, +)-Convolution problem. Here the input consists of two sequences (a[i])^{n-1}_{i = 0} and (b[i])^{n-1}_{i = 0}, the task is to compute the sequence (c[i])^{n-1}_{i = 0}, where c[k] = min_{i=0,...,k}(a[i] + b[k - i]).
In particular, we prove that if (min, +)-Convolution can be solved in O(tau(n)) time, then Bisection of graphs of treewidth t can be solved in time O(8^t t^{O(1)} log n * tau(n)), assuming a tree decomposition of width t is provided as input. Plugging in the naive O(n^2) time algorithm for (min, +)-Convolution yields a O(8^t t^{O(1)} n^2 log n) time algorithm for Bisection. This improves over the (dependence on n of the) O(2^t n^3) time algorithm of Jansen et al. [SICOMP 2005] at the cost of a worse dependence on t. "Conversely", we show that if Bisection can be solved in time O(beta(n)) on edge weighted trees, then (min, +)-Convolution can be solved in O(beta(n)) time as well. Thus, obtaining a sub-quadratic algorithm for Bisection on trees is extremely challenging, and could even be impossible. On the other hand, for unweighted graphs of treewidth t, by making use of a recent algorithm for Bounded Difference (min, +)-Convolution of Chan and Lewenstein [STOC 2015], we obtain a sub-quadratic algorithm for Bisection with running time O(8^t t^{O(1)} n^{1.864} log n)
qTorch: The Quantum Tensor Contraction Handler
Classical simulation of quantum computation is necessary for studying the
numerical behavior of quantum algorithms, as there does not yet exist a large
viable quantum computer on which to perform numerical tests. Tensor network
(TN) contraction is an algorithmic method that can efficiently simulate some
quantum circuits, often greatly reducing the computational cost over methods
that simulate the full Hilbert space. In this study we implement a tensor
network contraction program for simulating quantum circuits using multi-core
compute nodes. We show simulation results for the Max-Cut problem on 3- through
7-regular graphs using the quantum approximate optimization algorithm (QAOA),
successfully simulating up to 100 qubits. We test two different methods for
generating the ordering of tensor index contractions: one is based on the tree
decomposition of the line graph, while the other generates ordering using a
straight-forward stochastic scheme. Through studying instances of QAOA
circuits, we show the expected result that as the treewidth of the quantum
circuit's line graph decreases, TN contraction becomes significantly more
efficient than simulating the whole Hilbert space. The results in this work
suggest that tensor contraction methods are superior only when simulating
Max-Cut/QAOA with graphs of regularities approximately five and below. Insight
into this point of equal computational cost helps one determine which
simulation method will be more efficient for a given quantum circuit. The
stochastic contraction method outperforms the line graph based method only when
the time to calculate a reasonable tree decomposition is prohibitively
expensive. Finally, we release our software package, qTorch (Quantum TensOR
Contraction Handler), intended for general quantum circuit simulation.Comment: 21 pages, 8 figure
On the expressive power of planar perfect matching and permanents of bounded treewidth matrices
Valiant introduced some 25 years ago an algebraic model of computation along
with the complexity classes VP and VNP, which can be viewed as analogues of the
classical classes P and NP. They are defined using non-uniform sequences of
arithmetic circuits and provides a framework to study the complexity for
sequences of polynomials. Prominent examples of difficult (that is,
VNP-complete) problems in this model includes the permanent and hamiltonian
polynomials. While the permanent and hamiltonian polynomials in general are
difficult to evaluate, there have been research on which special cases of these
polynomials admits efficient evaluation. For instance, Barvinok has shown that
if the underlying matrix has bounded rank, both the permanent and the
hamiltonian polynomials can be evaluated in polynomial time, and thus are in
VP. Courcelle, Makowsky and Rotics have shown that for matrices of bounded
treewidth several difficult problems (including evaluating the permanent and
hamiltonian polynomials) can be solved efficiently. An earlier result of this
flavour is Kasteleyn's theorem which states that the sum of weights of perfect
matchings of a planar graph can be computed in polynomial time, and thus is in
VP also. For general graphs this problem is VNP-complete. In this paper we
investigate the expressive power of the above results. We show that the
permanent and hamiltonian polynomials for matrices of bounded treewidth both
are equivalent to arithmetic formulas. Also, arithmetic weakly skew circuits
are shown to be equivalent to the sum of weights of perfect matchings of planar
graphs.Comment: 14 page
- …