112 research outputs found
Polynomial-time Solvable #CSP Problems via Algebraic Models and Pfaffian Circuits
A Pfaffian circuit is a tensor contraction network where the edges are
labeled with changes of bases in such a way that a very specific set of
combinatorial properties are satisfied. By modeling the permissible changes of
bases as systems of polynomial equations, and then solving via computation, we
are able to identify classes of 0/1 planar #CSP problems solvable in
polynomial-time via the Pfaffian circuit evaluation theorem (a variant of L.
Valiant's Holant Theorem). We present two different models of 0/1 variables,
one that is possible under a homogeneous change of basis, and one that is
possible under a heterogeneous change of basis only. We enumerate a series of
1,2,3, and 4-arity gates/cogates that represent constraints, and define a class
of constraints that is possible under the assumption of a ``bridge" between two
particular changes of bases. We discuss the issue of planarity of Pfaffian
circuits, and demonstrate possible directions in algebraic computation for
designing a Pfaffian tensor contraction network fragment that can simulate a
swap gate/cogate. We conclude by developing the notion of a decomposable
gate/cogate, and discuss the computational benefits of this definition
Exact Algorithm for Sampling the 2D Ising Spin Glass
A sampling algorithm is presented that generates spin glass configurations of
the 2D Edwards-Anderson Ising spin glass at finite temperature, with
probabilities proportional to their Boltzmann weights. Such an algorithm
overcomes the slow dynamics of direct simulation and can be used to study
long-range correlation functions and coarse-grained dynamics. The algorithm
uses a correspondence between spin configurations on a regular lattice and
dimer (edge) coverings of a related graph: Wilson's algorithm [D. B. Wilson,
Proc. 8th Symp. Discrete Algorithms 258, (1997)] for sampling dimer coverings
on a planar lattice is adapted to generate samplings for the dimer problem
corresponding to both planar and toroidal spin glass samples. This algorithm is
recursive: it computes probabilities for spins along a "separator" that divides
the sample in half. Given the spins on the separator, sample configurations for
the two separated halves are generated by further division and assignment. The
algorithm is simplified by using Pfaffian elimination, rather than Gaussian
elimination, for sampling dimer configurations. For n spins and given floating
point precision, the algorithm has an asymptotic run-time of O(n^{3/2}); it is
found that the required precision scales as inverse temperature and grows only
slowly with system size. Sample applications and benchmarking results are
presented for samples of size up to n=128^2, with fixed and periodic boundary
conditions.Comment: 18 pages, 10 figures, 1 table; minor clarification
Counting Shortest Two Disjoint Paths in Cubic Planar Graphs with an NC Algorithm
Given an undirected graph and two disjoint vertex pairs and
, the Shortest two disjoint paths problem (S2DP) asks for the minimum
total length of two vertex disjoint paths connecting with , and
with , respectively.
We show that for cubic planar graphs there are NC algorithms, uniform
circuits of polynomial size and polylogarithmic depth, that compute the S2DP
and moreover also output the number of such minimum length path pairs.
Previously, to the best of our knowledge, no deterministic polynomial time
algorithm was known for S2DP in cubic planar graphs with arbitrary placement of
the terminals. In contrast, the randomized polynomial time algorithm by
Bj\"orklund and Husfeldt, ICALP 2014, for general graphs is much slower, is
serial in nature, and cannot count the solutions.
Our results are built on an approach by Hirai and Namba, Algorithmica 2017,
for a generalisation of S2DP, and fast algorithms for counting perfect
matchings in planar graphs
Faster Geometric Algorithms via Dynamic Determinant Computation
The computation of determinants or their signs is the core procedure in many
important geometric algorithms, such as convex hull, volume and point location.
As the dimension of the computation space grows, a higher percentage of the
total computation time is consumed by these computations. In this paper we
study the sequences of determinants that appear in geometric algorithms. The
computation of a single determinant is accelerated by using the information
from the previous computations in that sequence.
We propose two dynamic determinant algorithms with quadratic arithmetic
complexity when employed in convex hull and volume computations, and with
linear arithmetic complexity when used in point location problems. We implement
the proposed algorithms and perform an extensive experimental analysis. On one
hand, our analysis serves as a performance study of state-of-the-art
determinant algorithms and implementations. On the other hand, we demonstrate
the supremacy of our methods over state-of-the-art implementations of
determinant and geometric algorithms. Our experimental results include a 20 and
78 times speed-up in volume and point location computations in dimension 6 and
11 respectively.Comment: 29 pages, 8 figures, 3 table
- …