Fermat factorization method

Fermat factorization tests successive $r\ge \lceil \sqrt{N}\rceil$ until $r^2-N=s^2$, then returns $N=(r-s)(r+s)$.