Diophantine problems arise frequently in physics, in for example anomaly cancellation conditions, string consistency conditions and so forth. We present methods to solve such problems to high order on annealers that are based on the quadratic Ising Model. This is the intrinsic framework for both quantum annealing and for common forms of classical simulated annealing. We demonstrate the method on so-called Taxicab numbers (discovering some apparently new ones), and on the realistic problem of anomaly cancellation in U(1) extensions of the Standard Model.