Problem:
問題 ４

　平面上の n 個の点 P1, ..., Pn が与えられたとき，
距離が最小の２点を求めたい．

　入力ファイルの１行目には整数 n が書いてある．
２行目から n+1 行目のそれぞれには，
２つの正整数 x, y が１つの半角空白文字を区切りとして書いてある．
i+1 行目の x, y はそれぞれ Pi の x 座標，
Pi の y 座標である．
これら n 点の中から最も近い２点を選んだとき，
この２点間の距離の２乗を出力せよ．

　ただし，
2≤n≤500000 かつ -10000≤x≤10000, -10000≤y≤10000 とし，
５つの入力ファイルのうち３つでは n≤100 である．
また，どの入力ファイルにおいても，全ての座標は異なるものとする．

出力ファイルにおいては，
出力の最後の行にも改行コードを入れること．

入出力例

入力例１

2
0 0
1 1

出力例１

2

入力例２

3
5 5
0 0
-3 -4

出力例２

25

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

Solution in C++11 (judge verdict: Accepted)

#include <iostream>
#include <ccomplex>
#include <vector>
#include <cassert>
#include <utility>
#include <algorithm>
#include <string>
#include <deque>
#include <queue>
#include <functional>
#include <cmath>
#include <iomanip>
#include <map>
#include <numeric>
#include <list>
#include <assert.h>
#include <math.h>
#include <valarray>
#include <stdio.h>
#include <algorithm>
#include <set>
#include <complex>
#include <list>
#include <time.h>
#include <stack>
#include <locale>
#include <clocale>
#include <ctype.h>
#include <wchar.h>

using namespace std;
typedef long long int LL;
typedef pair<long long int, long long int> pii;
typedef pair<double, double> pdd;

#define SORT(c) sort((c).begin(),(c).end())
#define BACKSORT(c) sort((c).begin(),(c).end(),std::greater<LL>())
#define FOR(i,a,b) for(LL i=(a);i<(b);++i)
#define REP(i,n)  FOR(i,0,n)

int xcomp(vector<int>& a, vector<int>& b) {
	return a[0] < b[0];
}

int ycomp(vector<int>& a, vector<int>& b) {
	return a[1] < b[1];
}

#define XSORT(c) sort((c).begin(),(c).end(),xcomp)
#define YSORT(c) sort((c).begin(),(c).end(),ycomp)

#define ALL(x) (x).begin(),(x).end()

int henkan(int a) {
	if (a == 10000) {
		return 199;
	}
	else {
		return (a + 10000) / 100;
	}
}

LL mod = 1000000007;
int main() {

	int n;
	cin >> n;

	if (n <= 100) {
		vector<vector<int>> v(n, vector<int>(2));
		REP(i, n) {
			cin >> v[i][0] >> v[i][1];

		}
		int min = 1000000000;
		REP(i, n) {
			FOR(j, i+1, n) {
				int r = (v[i][0] - v[j][0])*(v[i][0] - v[j][0]) + (v[i][1] - v[j][1])*(v[i][1] - v[j][1]);
				if (r < min) { min = r; }
			}
		}
		cout << min << endl;
	}
	else {

		vector<vector<vector<vector<int>>>> v(200, vector<vector<vector<int>>>(200, vector<vector<int>>(0)));
		REP(i, n) {
			int a, b;
			cin >> a >> b;
			v[henkan(a)][henkan(b)].push_back({ a,b });
		}
		int min = 100000;
		REP(i, 199) {

			REP(j, 199) {
				vector<vector<int>> vv(0);

				vv.insert(vv.end(), v[i][j].begin(), v[i][j].end());
				vv.insert(vv.end(), v[i][j+1].begin(), v[i][j + 1].end());
				vv.insert(vv.end(), v[i + 1][j].begin(), v[i + 1][j].end());
				vv.insert(vv.end(), v[i + 1][j + 1].begin(), v[i + 1][j + 1].end());

				REP(ii, vv.size()) {
					FOR(jj, ii + 1, vv.size()) {
						int r = (vv[ii][0] - vv[jj][0])*(vv[ii][0] - vv[jj][0]) + (vv[ii][1] - vv[jj][1])*(vv[ii][1] - vv[jj][1]);
						if (r < min) { min = r; }
					}
				}
			}
		}
		cout << min << endl;
	}
}