The SK model — random Ising spins with all-to-all random couplings — is solved by Parisi's replica symmetry breaking ansatz. Below Tₒ, the system freezes into an ultrametric landscape of exponentially many metastable states. The overlap distribution P(q) transitions from a delta function (replica symmetric) to a broad distribution — a fingerprint of spin glass order.