CONCENTRATION OF THE RUNNING TIME 211
Theorem 11.2. Let X be a d raw from the distribution with parameter value
β
, and
η
≤ 4e
ε
−2
/(2 + ln(2E[H(X)])). Then if for all i, Z
s(i+1)
/Z
s(i)
≤ e
η
,then
V
rel
(R)=V
rel
(S) ≤
ε
. (11.27)
The proof is relatively complex, see [58] for details.
For many problems, it is possible to obtain an upper bound on E[X ] based on the
problem input, which can then be used to get a usable bound on
η
. For instance, in
the Ising model, E[X ] is upper bounded by the number of nodes of the graph.
11.6 Concentration of the running time
Many of the results shown in this text have been of the form of bounds on the ex-
pected running time of the procedure. While it is also usually possible to bound the
standard deviation of the running time using similar procedures, such arguments tend
to be far more complex than the analysis of the expected value.
Therefore, it is useful to understand how well the running time concentrates
around the mean for various forms of perfect simulation. In general, the goal is to
obtain an exponential tail on the running time.
Definition 11.3. Let X be a nonnegative random variable with finite mean where
there exist nonnegative constants
α
1
and
α
2
such that for all sufficien tly large k,
P(X ≥ kE[X]) ≤
α
1
k
−
α
2
. Then say that X has a polynomial tail. Suppose there ex-
ist positive constants
α
1
and
α
2
such that for sufficiently large k, P(X ≥ kE[X ]) ≤
α
1
exp(−
α
2
k). Then say that X has an exponential tail.
In general it is not possible to guarantee an exponential tail on the running time
of a particular perfect simulation procedure. That being said, for interruptible algo-
rithms, or for CFTP, it is possible to guarantee a tail that goes down faster than any
polynomial in k. In order to describe these algorithms, it is helpful to have the notion
of quantile.
Definition 11.4. For a random variable T , an
α
-quantile of T is any number a such
that P(T ≤ a) ≥
α
and P(T ≥ a) ≤ 1 −
α
. Denote the smallest such
α
-quantile by
q
α
(T ).
In this section, the quantile q
3/4
(T ) will be the one considered, so P(T ≥
q
3/4
(T )) ≤ 1/4. For random variables that are finite with probability 1, q
3/4
(T ) is
finite whether or not T has a finite mean. In the situation th at T does have a finite
mean E[T ], then by Markov’s inequality q
3/4
(T ) ≤ 4E[T ].
Theorem 11.3. Consider an interruptible perfect simulation algorithm A .LetT
A
denote the random number of steps taken by A . Then there exists an interruptible
algorithm B that takes T
B
steps where E[T
B
] ≤ 6q
3/4
(T ), and for k ≥ 1, P(T
B
≥
kE[T
B
]) ≤ exp(−k).
To prove the theorem, it is necessary to build an algorithm with the desired con-
centration properties. Before getting to that algorithm, first a subroutine is needed
that takes input t, and then runs algorithm A forward as most t steps. If algorithm A
returns an answer X,thenBounded
time A returns X as well. Otherwise it returns
the symbol ⊥, indicating that it failed to complete in time.