#include #define rep(i, l, n) for (int i = (l); i < (n); i++) #define squ(n) (n) * (n) using namespace std; using ll = long long; using u64 = unsigned long long; using Pair = pair; template using V = vector; template using PQ = priority_queue, greater >; u64 isqrt(u64 n) { u64 ok = 0ull; u64 ng = 4000000000ull; while (ng - ok > 1ull) { u64 mid = (ok + ng) >> 1ull; if (mid * mid <= n) { ok = mid; } else { ng = mid; } } return ok; } int main(void) { int n; cin >> n; V > p(n, V(3)); rep(i, 0, n) { ll x, y, t; cin >> x >> y >> t; p[i] = { x,y,t }; } V > route(n, V({})); 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(4ull * r1 * r2); } route[i].push_back({ (u64)j,d }); route[j].push_back({ (u64)i,d }); } } V path(n, 1000000000000000000ull); PQ que = {}; que.push({ 0ull,0ull }); while (que.empty() == false) { Pair pa = que.top(); que.pop(); u64 v = pa.second, d = pa.first; if (path[v] <= d) { continue; } path[v] = d; if (v == n - 1) { break; } for (Pair npa : route[v]) { u64 nv = npa.first, nd = npa.second; que.push({ max(d,nd),nv }); } } u64 ans = path[n - 1]; cout << ans << endl; return 0; }