結果

問題 No.2179 Planet Traveler
ユーザー kwm_tkwm_t
提出日時 2023-01-06 22:52:40
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 123 ms / 3,000 ms
コード長 2,274 bytes
コンパイル時間 4,500 ms
コンパイル使用メモリ 268,020 KB
実行使用メモリ 13,428 KB
最終ジャッジ日時 2023-08-20 16:16:26
合計ジャッジ時間 7,438 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,356 KB
testcase_01 AC 1 ms
4,388 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 1 ms
4,376 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 1 ms
4,384 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 2 ms
4,380 KB
testcase_08 AC 2 ms
4,380 KB
testcase_09 AC 1 ms
4,376 KB
testcase_10 AC 1 ms
4,384 KB
testcase_11 AC 40 ms
11,768 KB
testcase_12 AC 123 ms
13,428 KB
testcase_13 AC 118 ms
11,888 KB
testcase_14 AC 35 ms
11,784 KB
testcase_15 AC 35 ms
11,816 KB
testcase_16 AC 35 ms
11,980 KB
testcase_17 AC 118 ms
12,436 KB
testcase_18 AC 121 ms
11,448 KB
testcase_19 AC 118 ms
12,244 KB
testcase_20 AC 16 ms
4,380 KB
testcase_21 AC 110 ms
11,916 KB
testcase_22 AC 72 ms
7,816 KB
testcase_23 AC 62 ms
7,512 KB
testcase_24 AC 95 ms
8,208 KB
testcase_25 AC 74 ms
8,436 KB
testcase_26 AC 104 ms
12,412 KB
testcase_27 AC 17 ms
4,376 KB
testcase_28 AC 52 ms
8,524 KB
testcase_29 AC 95 ms
7,600 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
//using mint = modint1000000007;
//const int mod = 1000000007;
//using mint = modint998244353;
//const int mod = 998244353;
const int INF = 1e9;
const long long LINF = 1e18;
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rep2(i,l,r)for(int i=(l);i<(r);++i)
#define rrep(i, n) for (int i = (n-1); i >= 0; --i)
#define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i)
#define all(x) (x).begin(),(x).end()
#define allR(x) (x).rbegin(),(x).rend()
#define endl "\n"
#define P pair<int,int>
template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }
template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }
long long sqrt(long long n) {
	long long ok = 0;
	long long ng = 2 * INF + 10;
	while (1 < abs(ok - ng)) {
		long long c = (ng + ok) / 2;
		if ((c * c) <= n) ok = c;
		else ng = c;
	}
	//cout << n << " " << ok << endl;
	return ok;
}
int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	int n; cin >> n;
	vector<long long>x(n), y(n), t(n);
	rep(i, n) cin >> x[i] >> y[i] >> t[i];
	vector<pair<long long, P>>v;
	rep(i, n)rep2(j, i + 1, n) {
		if (t[i] == t[j]) {
			long long d2 = 0;
			d2 += (x[i] - x[j]) *  (x[i] - x[j]);
			d2 += (y[i] - y[j]) *  (y[i] - y[j]);
			v.push_back({ d2, {i,j} });
			//cout << i << j << "_" << d2 << endl;
		}
		else {
			long long d1, d0;
			d0 = x[i] * x[i] + y[i] * y[i];
			d1 = x[j] * x[j] + y[j] * y[j];
			if (d0 > d1)swap(d0, d1);
			long long ng = -1;
			long long ok = LINF + 1;
			long long sq = sqrt(4 * d0 * d1);
			while (1 < abs(ok - ng)) {
				long long mid = (ng + ok) / 2;
				auto check = [&]() {
					// (sqrt(d2)-sqrt(d1))^2<=mid
					// d1 + d2 - 2 * sqrt(d1*d2) <= mid * mid
					// d1 + d2 + mid * mid <= sqrt(4 * d1 * d2);
					return (d0 + d1 - mid <= sq);
				};
				if (check()) ok = mid;
				else ng = mid;
			}
			//cout << i << j << " " << ok << endl;
			v.push_back({ ok, {i,j} });
		}
	}
	sort(all(v));
	dsu uf(n);
	rep(i, v.size()) {
		uf.merge(v[i].second.first, v[i].second.second);
		if (uf.same(0, n - 1)) {
			cout << v[i].first << endl;
			break;
		}
	}
	return 0;
}
0