Problem:
Problem D: Circle and Points

You are given N points in the xy-plane. You have a
circle of radius one and move it on the xy-plane, so as to
enclose as many of the points as possible. Find how many points can be
simultaneously enclosed at the maximum. A point is considered
enclosed by a circle when it is inside or on the circle.

Fig 1. Circle and Points

Input

The input consists of a series of data sets, followed by a single line
only containing a single character '0', which indicates the end of the
input. Each data set begins with a line containing an integer N,
which indicates the number of points in the data set. It is
followed by N lines describing the coordinates of the
points. Each of the N lines has two decimal fractions X
and Y, describing the x- and y-coordinates of a
point, respectively. They are given with five digits after the decimal
point.

You may assume 1 <= N <= 300, 0.0 <= X <= 10.0, and 0.0
<= Y <= 10.0. No two points are closer than 0.0001. No two
points in a data set are approximately at a distance of 2.0. More
precisely, for any two points in a data set, the distance d
between the two never satisfies 1.9999 <= d <= 2.0001. Finally,
no three points in a data set are simultaneously very close to a
single circle of radius one. More precisely,
let P1, P2, and P3 be any three points in a
data set, and d1, d2, and d3
the distances from an arbitrarily selected point in the xy-plane to each of them
respectively. Then it never simultaneously holds that 0.9999 <=
di <= 1.0001 (i = 1, 2, 3).

Output

For each data set, print a single line containing the maximum number
of points in the data set that can be simultaneously enclosed by a
circle of radius one. No other characters including leading and
trailing spaces should be printed.

Sample Input

3
6.47634 7.69628
5.16828 4.79915
6.69533 6.20378
6
7.15296 4.08328
6.50827 2.69466
5.91219 3.86661
5.29853 4.16097
6.10838 3.46039
6.34060 2.41599
8
7.90650 4.01746
4.10998 4.18354
4.67289 4.01887
6.33885 4.28388
4.98106 3.82728
5.12379 5.16473
7.84664 4.67693
4.02776 3.87990
20
6.65128 5.47490
6.42743 6.26189
6.35864 4.61611
6.59020 4.54228
4.43967 5.70059
4.38226 5.70536
5.50755 6.18163
7.41971 6.13668
6.71936 3.04496
5.61832 4.23857
5.99424 4.29328
5.60961 4.32998
6.82242 5.79683
5.44693 3.82724
6.70906 3.65736
7.89087 5.68000
6.23300 4.59530
5.92401 4.92329
6.24168 3.81389
6.22671 3.62210
0

Output for the Sample Input

2
5
5
11

Solution in JAVA (judge verdict: Time Limit Exceeded)

import java.io.IOException;
import java.util.InputMismatchException;

public class Main {
	double EPS = 1e-12;
	double INF = 1e12;

	class Point {
		double x, y;

		Point(double x, double y) {
			this.x = x;
			this.y = y;
		}
	}

	class Segment {
		Point a, b;

		Segment(Point a, Point b) {
			this.a = a;
			this.b = b;
		}
	}

	class Line {
		Point a, b;

		Line(Point a, Point b) {
			this.a = a;
			this.b = b;
		}
	}

	class Circle {
		Point c;
		double r;

		Circle(Point c, double r) {
			this.c = c;
			this.r = r;
		}
	}

	double norm(Point a) {
		return Math.hypot(a.x, a.y);
	}

	Point sub(Point a, Point b) {
		return new Point(a.x - b.x, a.y - b.y);
	}

	double dot(Point a, Point b) {
		return a.x * b.x + a.y * b.y;
	}

	double cross(Point a, Point b) {
		return a.x * b.y - a.y * b.x;
	}

	double abs(double a) {
		return 0 < a ? a : -a;
	}

	// 交点の個数が0または無限の時はreturn null;
	Point[] crossPoint(Circle c1, Circle c2) {
		double d = norm(sub(c1.c, c2.c));
		if (d < EPS && abs(c1.r - c2.r) < EPS) {
			return null;
		} else if (c1.r + c2.r < d - EPS || d + EPS < abs(c1.r - c2.r)) {
			return null;
		} else {
			double a = (c1.r * c1.r - c2.r * c2.r + d * d) / (2 * d);
			double h = Math.sqrt(Math.max(c1.r * c1.r - a * a, 0));
			Point p1 = new Point(c1.c.x + a / d * (c2.c.x - c1.c.x), c1.c.y + a / d * (c2.c.y - c1.c.y));
			Point p2 = new Point(h / d * (c2.c.x - c1.c.x), h / d * (c2.c.y - c1.c.y));
			if (norm(p2) < EPS) {
				return new Point[] { p1 };
			} else {
				return new Point[] { new Point(p1.x - p2.y, p1.y + p2.x), new Point(p1.x + p2.y, p1.y - p2.x) };
			}
		}
	}

	boolean isCover(Circle c, Point p) {
		return norm(sub(c.c, p)) <= c.r + EPS;
	}

	void run() {
		MyScanner sc = new MyScanner();

		while (true) {
			int n = sc.nextInt();
			if (n == 0) {
				break;
			}
			Point[] p = new Point[n];
			Circle[] c = new Circle[n];
			for (int i = 0; i < n; i++) {
				p[i] = new Point(sc.nextDouble(), sc.nextDouble());
				c[i] = new Circle(p[i], 1.);
			}

			int max = 1;
			for (int i = 0; i < n; i++) {
				for (int j = i + 1; j < n; j++) {
					Point[] cp = crossPoint(c[i], c[j]);
					if (cp == null) {
						continue;
					}
					for (int k = 0; k < cp.length; k++) {
						Circle nc = new Circle(cp[k], 1.);
						int cnt = 0;
						for (int l = 0; l < n; l++) {
							if (isCover(nc, p[l])) {
								cnt++;
							}
						}
						max = Math.max(max, cnt);
					}
				}
			}
			System.out.println(max);
		}
	}

	public static void main(String[] args) {
		new Main().run();
	}

	public void mapDebug(int[][] a) {
		System.out.println("--------map display---------");
		for (int i = 0; i < a.length; i++) {
			for (int j = 0; j < a[i].length; j++) {
				System.out.printf("%3d ", a[i][j]);
			}
			System.out.println();
		}
		System.out.println("----------------------------" + '\n');
	}

	class MyScanner {
		int read() {
			try {
				return System.in.read();
			} catch (IOException e) {
				throw new InputMismatchException();
			}
		}

		boolean isSpaceChar(int c) {
			return c == ' ' || c == '\n' || c == '\r' || c == '\t' || c == -1;
		}

		boolean isEndline(int c) {
			return c == '\n' || c == '\r' || c == -1;
		}

		int nextInt() {
			return Integer.parseInt(next());
		}

		int[] nextIntArray(int n) {
			int[] array = new int[n];
			for (int i = 0; i < n; i++)
				array[i] = nextInt();
			return array;
		}

		long nextLong() {
			return Long.parseLong(next());
		}

		long[] nextLongArray(int n) {
			long[] array = new long[n];
			for (int i = 0; i < n; i++)
				array[i] = nextLong();
			return array;
		}

		double nextDouble() {
			return Double.parseDouble(next());
		}

		double[] nextDoubleArray(int n) {
			double[] array = new double[n];
			for (int i = 0; i < n; i++)
				array[i] = nextDouble();
			return array;
		}

		String next() {
			int c = read();
			while (isSpaceChar(c))
				c = read();
			StringBuilder res = new StringBuilder();
			do {
				res.appendCodePoint(c);
				c = read();
			} while (!isSpaceChar(c));
			return res.toString();
		}

		String[] nextStringArray(int n) {
			String[] array = new String[n];
			for (int i = 0; i < n; i++)
				array[i] = next();

			return array;
		}

		String nextLine() {
			int c = read();
			while (isEndline(c))
				c = read();
			StringBuilder res = new StringBuilder();
			do {
				res.appendCodePoint(c);
				c = read();
			} while (!isEndline(c));
			return res.toString();
		}
	}
}