For random k-SAT with prespecified solution, Stirling's asymptotic expansion in Eq. (8) shows that for y=O(n) when m grows as , which is much less that the maximum . In this case, the asymptotic behavior of the bound B is determined by the values near the maximum of the binomial coefficient, i.e., near y=n/2. Thus if , we have . From Eq. (29), at least when B< 1. This is the case for where
where . For instance, is 1.01 and 17.3 for k=2 and 3, respectively. Thus, the algorithm presented here is simple enough to allow an analytic bound on its behavior for highly constrained problems, thus demonstrating its asymptotic effectiveness in these cases. Other, more complex, structured quantum algorithms have only been evaluated empirically [29, 30], which is limited to small problems.
The algorithm's behavior with fewer constraints, i.e., , is not easily evaluted analytically since the bound provided by B is no longer useful. Instead, the behavior can be examined empirically using a classical simulation [34], which is however limited to problems with a relatively small number of variables. These empirical studies may eventually be extendable to larger problems using approximate evaluation techniques [9].
An example is given in Fig. 3. This shows good performance for highly constrained problems, as expected from the behavior of the lower bound. Performance is also good with few constraints; not because the algorithm is capturing the problem structure particularly well but rather because there are many solutions to weakly constrained problems. As with other classical and quantum search methods that use problem structure, the hardest cases are for problems with an intermediate number of constraints [32].
Figure 3: Probability to find a solution for random 3-SAT for n=10 (solid) and 20 (dashed) vs. m/n, on a log-log scale. Each point is an average over at least 100 problem instances, and includes error bars for the standard deviation in this estimate of the averge. The error bars are smaller than the size of the plotted points. The gray lines show the corresponding lower bounds .
From Fig. 4 we see that the nonzero asymptotic limit for the probability of a solution appears to continue for somewhat fewer constraints than expected from the value of . Below this value, the probability for finding a solution appears to decrease as a power of n, indicated by linear scaling on the log-log plot of Fig. 4 for . Similar empirical evaluations of the scaling with an intermediate number of constraints where the hard problems are concentrated, e.g., m = 4 n for 3-SAT, shows linear scaling on a log-plot, indicating exponential decrease in the probability to find a solution. Moreover, the resulting search costs in these cases are larger than those of other structured quantum search algorithms [29, 30]. Thus the structure of these harder cases differs enough from the simple 1-SAT problems that this algorithm is not effective for them.
Figure 4: Probability to find a solution for random soluble 3-SAT vs. n with, from top to bottom, , and , respectively, on a log-log plot. For each value of n, at least 100 problem instances were used. Error bars showing the expected error in the estimate are included but are smaller than the size of the plotted points.
In summary, the algorithm solves highly constrained problems with in O(1) steps for , and possibly for somewhat smaller values of as well. As the number of clauses is further reduced, the required number of steps appears to grow polynomially when and exponentially when m = O(n).