Title
Adaptive simulated annealing: A near-optimal connection between sampling and counting
Abstract
We present a near-optimal reduction from approximately counting the cardinality of a discrete set to approximately sampling elements of the set. An important application of our work is to approximating the partition function Z of a discrete system, such as the Ising model, matchings or colorings of a graph. The standard approach to estimating the partition function Z(\beta *) at some desired inverse temperature \beta * is to define a sequence, which we call a cooling schedule, \beta 0 = 0 \le \beta 1 \le\cdots\le \beta \ell= \beta * where Z(0) is trivial to compute and the ratios Z(\beta i + 1)/Z(\beta i) are easy to estimate by sampling from the distribution corresponding to Z(\beta i). Previous approaches required a cooling schedule of length {\rm O}*(1nA) where A = Z(0), thereby ensuring that each ratio Z(\beta i + 1)/Z(\beta i) is bounded. We present a cooling schedule of length \ell= {\rm O}*\left( {\sqrt {1nA} } \right). For well-studied problems such as estimating the partition function of the Ising model, or approximating the number of colorings or matchings of a graph, our cooling schedule is of length {\rm O}*\left( {\sqrt n } \right) and the total number of samples required is {\rm O}*\left( n \right). This implies an overall savings of a factor of roughly n in the running time of the approximate counting algorithm compared to the previous best approach. A similar improvement in the length of the cooling schedule was recently obtained by Lovász and Vempala in the context of estimating the volume of convex bodies. While our reduction is inspired by theirs, the discrete analogue of their result turns out to be significantly more difficult. Whereas a fixed schedule suffices in their setting, we prove that in the discrete setting we need an adaptive schedule, i. e., the schedule depends on Z. More precisely, we prove any non-adaptive cooling schedule has length at least {\rm O}*\left( {1nA} \right), and we present an algorithm to find an adaptive schedule of length {\rm O}*\left( {\sqrt {1nA} } \right) and a nearly matching lower bound.
Year
DOI
Venue
2009
10.1145/1516512.1516520
IEEE Symposium on Foundations of Computer Science
Keywords
DocType
Volume
fixed schedule suffices,ratio z,discrete set,markov chain monte carlo,near-optimal connection,length o,adaptive schedule,partition function z,ising model,cooling schedule,counting,simulated annealing,discrete analogue,adaptive simulated annealing,ratios z,discrete system,discrete mathematics,data structure,partition function
Journal
56
Issue
ISSN
ISBN
3
0272-5428
0-7695-3010-9
Citations 
PageRank 
References 
26
1.67
16
Authors
3
Name
Order
Citations
PageRank
Daniel Štefankovič120116.34
Santosh Vempala23546523.21
Eric Vigoda374776.55