#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using namespace atcoder; typedef long long ll; #define rep(i, n) for (int i = 0; i < (int)(n); i++) #define repr(i, n) for (int i = (int)(n) - 1; i >= 0; i--) #define repk(i, k, n) for (int i = k; i < (int)(n); i++) #define all(v) v.begin(), v.end() #define mod1 1000000007 #define mod2 998244353 #define mod3 100000007 #define vi vector #define vs vector #define vc vector #define vl vector #define vb vector #define vvi vector> #define vvc vector> #define vvl vector> #define vvb vector> #define vvvi vector>> #define vvvl vector>> #define pii pair #define pil pair #define pli pair #define pll pair #define vpii vector> #define vpll vector> #define vvpii vector>> #define vvpll vector>> template void debug(T e) { cerr << e << endl; } template void debug(vector &v) { rep(i, v.size()) { cerr << v[i] << " "; } cerr << endl; } template void debug(vector> &v) { rep(i, v.size()) { rep(j, v[i].size()) { cerr << v[i][j] << " "; } cerr << endl; } } template void debug(vector> &v) { rep(i, v.size()) { cerr << v[i].first << " " << v[i].second << endl; } } template void debug(set &st) { for (auto itr = st.begin(); itr != st.end(); itr++) { cerr << *itr << " "; } cerr << endl; } template void debug(multiset &ms) { for (auto itr = ms.begin(); itr != ms.end(); itr++) { cerr << *itr << " "; } cerr << endl; } template void debug(map &mp) { for (auto itr = mp.begin(); itr != mp.end(); itr++) { cerr << itr->first << " " << itr->second << endl; } } void debug_out() { cerr << endl; } template void debug_out(Head H, Tail... T) { cerr << H << " "; debug_out(T...); } using mint = modint998244353; bool s1(ll x1, ll y1, ll x2, ll y2, ll s) { return (x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2) <= s; } bool s2(ll x1, ll y1, ll x2, ll y2, ll s) { ll left = x1 * x1 + y1 * y1 + x2 * x2 + y2 * y2 - s; if (x1 * x1 + y1 * y1 + x2 * x2 + y2 * y2 - s < 0) { return true; } return left * left <= 4 * (x1 * x1 + y1 * y1) * (x2 * x2 + y2 * y2); } int main() { ll N; cin >> N; vector X(N); vector Y(N); vector T(N); for (ll i = 0; i < N; i++) { cin >> X[i] >> Y[i] >> T[i]; } ll ng = -1; ll ok = 4000000000; while (ok - ng > 1) { ll s = (ok + ng) / 2; queue que; que.push(0); vector visited(N, false); visited[0] = true; while (que.size()) { ll p = que.front(); que.pop(); for (ll i = 0; i < N; i++) { if (visited[i]) { continue; } if (T[i] == T[p]) { if (s1(X[i], Y[i], X[p], Y[p], s)) { visited[i] = true; que.push(i); } } else { if (s2(X[i], Y[i], X[p], Y[p], s)) { visited[i] = true; que.push(i); } } } } if (visited[N - 1]) { ok = s; } else { ng = s; } } cout << ok << endl; }