#795880 (C++17) No.2179 Planet Traveler

• 問題: No.2179 Planet Traveler
• ユーザー: MasKoaTS
• 提出日時: 2022-09-06 20:36:00
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 556 ms / 3,000 ms
• コード長: 1,220 bytes
• コンパイル時間: 4,177 ms
• コンパイル使用メモリ: 268,432 KB
• 実行使用メモリ: 8,448 KB
• 最終ジャッジ日時: 2024-05-07 20:28:33
• 合計ジャッジ時間: 12,639 ms
• ジャッジサーバーID （参考情報）: judge4 / judge1

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	2 ms; 5,248 KB
testcase_02	AC	1 ms; 5,376 KB
testcase_03	AC	1 ms; 5,376 KB
testcase_04	AC	2 ms; 5,376 KB
testcase_05	AC	2 ms; 5,376 KB
testcase_06	AC	2 ms; 5,376 KB
testcase_07	AC	2 ms; 5,376 KB
testcase_08	AC	1 ms; 5,376 KB
testcase_09	AC	2 ms; 5,376 KB
testcase_10	AC	2 ms; 5,376 KB
testcase_11	AC	514 ms; 8,320 KB
testcase_12	AC	556 ms; 8,320 KB
testcase_13	AC	552 ms; 8,320 KB
testcase_14	AC	541 ms; 8,448 KB
testcase_15	AC	526 ms; 8,448 KB
testcase_16	AC	511 ms; 8,320 KB
testcase_17	AC	547 ms; 8,448 KB
testcase_18	AC	548 ms; 8,320 KB
testcase_19	AC	547 ms; 8,448 KB
testcase_20	AC	30 ms; 5,376 KB
testcase_21	AC	494 ms; 8,064 KB
testcase_22	AC	264 ms; 6,272 KB
testcase_23	AC	219 ms; 5,888 KB
testcase_24	AC	396 ms; 7,424 KB
testcase_25	AC	286 ms; 6,528 KB
testcase_26	AC	465 ms; 7,808 KB
testcase_27	AC	34 ms; 5,376 KB
testcase_28	AC	168 ms; 5,376 KB
testcase_29	AC	404 ms; 7,296 KB

ソースコード

#include <math.h>
#include <bits/stdc++.h>
#include <atcoder/all>
#define rep(i, l, n) for (int i = (l); i < (n); i++)
#define squ(n) (n) * (n)
using namespace std;
using namespace atcoder;
using ll = long long;
using u64 = unsigned long long;
template <class T>	using V = vector<T>;

u64 isqrt(u64 n) {
	u64 ok = 0ll;
	u64 ng = 4000000000ll;
	while (ng - ok > 1ll) {
		u64 mid = (ok + ng) >> 1ll;
		if (mid * mid <= n) {
			ok = mid;
		}
		else {
			ng = mid;
		}
	}
	return ok;
}


int main(void) {
	int n;	cin >> n;
	V<V<ll> > p(n, V<ll>(3));
	rep(i, 0, n) {
		ll x, y, t;	cin >> x >> y >> t;
		p[i] = { x,y,t };
	}

	V<V<u64> > dist(n, V<u64>(n, 1000000000000000000ll));
	rep(i, 0, n) {
		dist[i][i] = 0ll;
	}
	rep(i, 0, n - 1) {
		u64 r1 = (u64)(squ(p[i][0]) + squ(p[i][1]));
		rep(j, i + 1, n) {
			u64 r2 = (u64)(squ(p[j][0]) + squ(p[j][1]));
			u64 d = (u64)(squ(p[i][0] - p[j][0]) + squ(p[i][1] - p[j][1]));
			if (p[i][2] != p[j][2]) {
				d = r1 + r2 - isqrt(4ll * r1 * r2);
			}
			dist[i][j] = dist[j][i] = d;
		}
	}

	rep(k, 0, n) {
		rep(i, 0, n) {
			rep(j, 0, n) {
				dist[i][j] = min(dist[i][j], max(dist[i][k], dist[k][j]));
			}
		}
	}
	u64 ans = dist[0][n - 1];
	cout << ans << endl;

	return 0;
}