結果

問題 No.96 圏外です。
ユーザー pekempeypekempey
提出日時 2016-02-18 07:38:15
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 181 ms / 5,000 ms
コード長 3,108 bytes
コンパイル時間 2,238 ms
コンパイル使用メモリ 175,208 KB
実行使用メモリ 112,108 KB
最終ジャッジ日時 2024-06-06 16:15:31
合計ジャッジ時間 5,202 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
99,328 KB
testcase_01 AC 37 ms
99,232 KB
testcase_02 AC 38 ms
99,328 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 40 ms
99,584 KB
testcase_05 AC 42 ms
99,696 KB
testcase_06 AC 43 ms
100,016 KB
testcase_07 AC 44 ms
100,272 KB
testcase_08 AC 48 ms
100,716 KB
testcase_09 AC 53 ms
101,248 KB
testcase_10 AC 56 ms
101,888 KB
testcase_11 AC 62 ms
102,784 KB
testcase_12 AC 69 ms
103,904 KB
testcase_13 AC 83 ms
104,528 KB
testcase_14 AC 91 ms
106,176 KB
testcase_15 AC 111 ms
107,008 KB
testcase_16 AC 119 ms
109,304 KB
testcase_17 AC 127 ms
111,700 KB
testcase_18 AC 181 ms
111,480 KB
testcase_19 AC 147 ms
111,488 KB
testcase_20 AC 121 ms
109,168 KB
testcase_21 AC 151 ms
112,108 KB
testcase_22 AC 159 ms
110,860 KB
testcase_23 AC 170 ms
110,860 KB
testcase_24 AC 37 ms
99,280 KB
testcase_25 AC 100 ms
108,032 KB
testcase_26 AC 118 ms
110,592 KB
testcase_27 AC 107 ms
109,056 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:97:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   97 |                 scanf("%lld %lld", &x[i], &y[i]);
      |                 ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define rep(i, a) for (int i = 0; i < (a); i++)
#define rep2(i, a, b) for (int i = (a); i < (b); i++)
#define repr(i, a) for (int i = (a) - 1; i >= 0; i--)
#define repr2(i, a, b) for (int i = (b) - 1; i >= (a); i--)
template<class T1, class T2> bool chmin(T1 &a, T2 b) { return b < a && (a = b, true); }
template<class T1, class T2> bool chmax(T1 &a, T2 b) { return a < b && (a = b, true); }
typedef long long ll;

typedef double D;
typedef complex<D> P;
const D eps = 1e-8;
D dot(P a, P b) { return real(conj(a) * b); }
D cross(P a, P b) { return imag(conj(a) * b); }

vector<P> convexfull(vector<P> &ps) {
	auto comp = [](P a, P b) {
		if (a.real() != b.real()) return a.real() < b.real();
		return a.imag() < b.imag();
	};
	int n = ps.size();
	sort(ps.begin(), ps.end(), comp);
	int k = 0;
	vector<P> qs(n * 2);
	for (int i = 0; i < n; i++) {
		while (k > 1 && cross(qs[k - 1] - qs[k - 2], ps[i] - qs[k - 1]) <= 0) k--;
		qs[k++] = ps[i];
	}
	for (int i = n - 2, t = k; i >= 0; i--) {
		while (k > t && cross(qs[k - 1] - qs[k - 2], ps[i] - qs[k - 1]) <= 0) k--;
		qs[k++] = ps[i];
	}
	qs.resize(k - 1);
	return qs;
}

double farthest(vector<P> &ps) {
	auto qs = convexfull(ps);
	int n = qs.size();
	if (n == 2) return abs(qs[0] - qs[1]);
	auto comp = [](P a, P b) {
		if (a.real() != b.real()) return a.real() < b.real();
		return a.imag() < b.imag();
	};
	int i = 0, j = 0;
	for (int k = 0; k < n; k++) {
		if (!comp(qs[i], qs[k])) i = k;
		if (comp(qs[j], qs[k])) j = k;
	}
	double ret = 0;
	int si = i, sj = j;
	while (i != sj || j != si) {
		ret = max(ret, abs(qs[i] - qs[j]));
		if (cross(qs[(i + 1) % n] - qs[i], qs[(j + 1) % n] - qs[j]) < 0) {
			i = (i + 1) % n;
		} else {
			j = (j + 1) % n;
		}
	}
	return ret;
}

struct UnionFind {
	vector<int> parent, size;
	UnionFind(int n) : parent(n), size(n, 1) { for (int i = 0; i < n; i++) parent[i] = i; }
	int operator[](int x) { return parent[x] == x ? x : (parent[x] = operator[](parent[x])); }
	bool merge(int x, int y) {
		if ((x = operator[](x)) == (y = operator[](y))) return false;
		parent[x] = parent[y] = size[x] > size[y] ? parent[x] : parent[y];
		size[x] = size[y] = size[x] + size[y];
		return true;
	}
};

ll hypot2(ll x, ll y) {
	return x * x + y * y;
}

int main() {
	int n;
	cin >> n;
	UnionFind uf(n);

	if (n == 0) {
		cout << 1 << endl;
		return 0;
	} else if (n == 1) {
		cout << 2 << endl;
		return 0;
	}

	vector<ll> x(n);
	vector<ll> y(n);
	static vector<int> bucket[2020][2020];
	rep(i, n) {
		scanf("%lld %lld", &x[i], &y[i]);
		x[i] += 10010;
		y[i] += 10010;
		bucket[x[i] / 10][y[i] / 10].push_back(i);
	}

	rep(i, n) {
		for (int dy = -1; dy <= 1; dy++) {
			for (int dx = -1; dx <= 1; dx++) {
				int bx = x[i] / 10;
				int by = y[i] / 10;
				for (int j : bucket[bx + dx][by + dy]) {
					if (hypot2(x[j] - x[i], y[j] - y[i]) <= 100) uf.merge(i, j);
				}
			}
		}
	}
	
	vector<vector<P>> ps(n);
	rep(i, n) ps[uf[i]].emplace_back(x[i], y[i]);

	double ans = 0;
	rep(i, n) if (!ps[i].empty()) chmax(ans, farthest(ps[i]));
	ans += 2;
	printf("%.20f\n", ans);

	return 0;
}
0