The smallest edge length of a square to include point <code>(x, y)</code> is <code>max(abs(x), abs(y)) * 2</code>.

Sort the points by <code>max(abs(x), abs(y))</code> and try each edge length, check the included point tags.