This work was produced for the Algorithms for Optimization, Inference and Learning course in Fall '19. It is based on M. Mézard, F. Ricci-Tersenghi, R. Zecchina, Two Solutions to Diluted p-Spin Models and XORSAT Problems, Journal of Statistical Physics, 2003.
Belief Propagation equations are implemented in their Min-Sum form to look for solutions to XOR-SAT instances drawn from a random ensemble in the thermodynamic limit of variables.
The problem is equivalent to that of finding the Ground State for a p-spin Ising model (an Ising model where spins interact not in pairs but in p-uples) defined over a factor graph.
A brief explanation of the main results follows
Recursively remove from the graph all leaves and the factors (also called hyperedges) they are attached to, until there is no leaf left. The subgraph resulting at the end, called core, is the frozen part of the graph, formed by those variable whose value is forced by the rigidness of the system.
The density 
is the ratio of number of constraints and number of variables. As the density varies from 0 to 1 and beyond, the system undergoes two phase transitions
- Dynamic transition: At some density
(~0.82 for p=3), the core size jumps from 0 to a value extensive with N, the number of variables - SAT/UNSAT transition: From
(~0.92 for p=3) up, the core contains more constraints than variables and the probability of an instance being satisfiable becomes exponentially small 
Energy in the p-spin Ising model corresponds to the number of unsatisfied constraints in the XOR-SAT formula. For < 
, all solutions are concentrated in one cluster, meaning that each one can be reached from any other with a finite number of spin flips. In this regime BP, which is a local algorithm, is able to find solutions. Instead, for < < 
, systems are still solvable, but solutions are scattered into an extensive (with respect to N) number of clusters, therefore more sophisticated techniques such as Survey Propagation are required in order
to approximate the actual behavior.
- Julia Language: Julia
- Phase transitions in combinatorial optimization: "Simplest random K-satisfiability problem"
- Survey Propagation: "Survey propagation: an algorithm for satisfiability"