Problem:
最古の遺跡

問題

昔, そこには集落があり, 多くの人が暮らしていた. 人々は形も大きさも様々な建物を建てた. だが, それらの建造物は既に失われ, 文献と, 遺跡から見つかった柱だけが建造物の位置を知る手がかりだった.

文献には神殿の記述がある. 神殿は上から見ると正確に正方形になっており, その四隅には柱があった. 神殿がどの向きを向いていたかはわからない. また, 辺上や内部に柱があったかどうかもわからない. 考古学者たちは, 遺跡から見つかった柱の中で, 正方形になっているもののうち, 面積が最大のものが神殿に違いないと考えた.

柱の位置の座標が与えられるので, 4 本の柱でできる正方形のうち面積が最大のものを探し, その面積を出力するプログラムを書け. なお, 正方形の辺は座標軸に平行とは限らないことに注意せよ.

例

下の図の例では, 10 本の柱があり, 座標 (4, 2), (5, 2), (5, 3), (4, 3) にある 4 本と座標 (1, 1), (4, 0), (5, 3), (2, 4) にある 4 本が正方形をなしている. 面積が最大の正方形は後者で, その面積は 10 である.

入力

入力は複数のデータセットからなる．各データセットは以下の形式で与えられる．入力はゼロ１つを含む行で終了する．

1 行目には, 遺跡から見つかった柱の本数 n が書かれている.

2 行目から n + 1 行目までの n 行の各々には, 柱の x 座標と y 座標が空白区切りで書かれている.

1 本の柱が 2 度以上現れることはない.

n は 1 ≤ n ≤ 3000 を満たす整数であり, 柱の x 座標と y 座標は 0 以上 5000 以下の整数である.

採点用データのうち, 配点の 30% 分は 1 ≤ n ≤ 100 を満たし, 配点の 60% 分は1 ≤ n ≤ 500 を満たす.

データセットの数は 10 を超えない．

出力

データセットごとに 1 個の整数を出力する. 4 本の柱からなる正方形が存在する場合は, そのような正方形のうち面積が最大のものの面積を出力し, そのような正方形が存在しない場合は 0 を出力せよ.

入出力例

入力例

10
9 4
4 3
1 1
4 2
2 4
5 8
4 0
5 3
0 5
5 2
10
9 4
4 3
1 1
4 2
2 4
5 8
4 0
5 3
0 5
5 2
0

出力例

10
10

上記問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C++ (judge verdict: Wrong Answer)

#include <vector>
#include <map>
#include <set>
#include <stack>
#include <queue>
#include <algorithm>
#include <utility>
#include <functional>
#include <sstream>
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <string>
#include <cstring>
#include <ctime>
#include <climits>
#include <fstream>
#include <list>
using namespace std;
inline int toInt(string s) { int v; istringstream sin(s); sin >> v; return v; }
template<class T> inline string toStr(T x) { ostringstream sout; sout << x; return sout.str(); }
typedef vector<int> vi;
typedef vector<vi>  vvi;
typedef vector<string> vs;
typedef pair<int, int> pii;
typedef long long ll;
#define ALL(a) (a).begin(),(a).end()
#define RALL(a) (a).rbegin(),(a).rend()
#define FOR(i,a,b) for(int i=(a);i<=(b);++i)
#define REP(i,n) FOR(i,0,(n)-1)
const double EPS = 1e-10;
const double PI = acos(-1.0);
const int INF = INT_MAX / 10;

#include <complex>
typedef complex<double> P;
#define EQ(x,y) (abs((x)-(y))<EPS)

P rotate(P t, P p, double r) {
	double ta = cos(r)*(t.real() - p.real()) - sin(r)*(t.imag() - p.imag()) + p.real();
	double tb = sin(r)*(t.real() - p.real()) + cos(r)*(t.imag() - p.imag()) + p.imag();
	return P(ta, tb);
}

P mid(P a, P b) {
	return P((a.real()+ b.real()) / 2, (a.imag() + b.imag()) / 2);
}

bool existPt[5001][5001];

bool isPillar(P a) {
	int far = floor(a.real()), fai = floor(a.imag());
	if (!EQ(a.real(), far) || !EQ(a.imag(), fai)) {
		return false;
	}
	if (far < 0 || 5000 < far || fai < 0 || 5000 < fai) {
		return false;
	}
	return existPt[far][fai];
}

bool isSquare(P a, P b) {
	P midpoint = mid(a, b);
	P pa = rotate(a, midpoint, PI / 2), pb = rotate(b, midpoint, PI / 2);
	return (isPillar(pa) && isPillar(pb));
}

int area(P a, P b) {
	return abs(a - b)*abs(a - b) / 2;
}

int main() {
	int n;
	while (cin >> n, n) {
		memset(existPt, false, sizeof(existPt));
		int x, y;
		vector<P> pts(n);
		REP(i, n) {
			cin >> x >> y;
			pts[i] = P(x, y);
			existPt[x][y] = true;
		}

		int ans = 0;
		REP(i, n) {
			FOR(j, i + 1, n - 1) {
				if (isSquare(pts[i], pts[j])) {
					ans = max(ans, area(pts[i], pts[j]));
				}
			}
		}

		cout << ans << endl;
	}
}