Abstract
Despite significant advances in characterizing the highly nonconvex landscapes of constraint satisfaction problems, the good performance of certain algorithms in solving hard combinatorial optimization tasks remains poorly understood. This gap in understanding stems largely from the lack of theoretical tools for analyzing their out-of-equilibrium dynamics. To address this challenge, we develop a system of approximate master equations that capture the behavior of local search algorithms in constraint satisfaction problems. Our framework shows excellent qualitative agreement with the phase diagrams of two paradigmatic algorithms: Focused Metropolis Search (FMS) and greedy-WalkSAT (G-WalkSAT) for random 3-SAT. The equations not only confirm the numerical observation that G-WalkSAT's algorithmic threshold is nearly parameter-independent but also successfully predict FMS's threshold beyond the clustering transition. We also exploit these equations in a decimation scheme, demonstrating that the computed marginals encode valuable information about the local structure of the solution space explored by stochastic algorithms. Notably, our decimation approach achieves a threshold that surpasses the clustering transition, outperforming conventional methods like Belief Propagation-guided decimation. These results challenge the prevailing assumption that long-range correlations are always necessary to describe efficient local search dynamics and open a path to designing efficient algorithms to solve combinatorial optimization problems.