#include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; const ll INF=1LL<<40; typedef pair P; typedef pair PP; const ll MOD=998244353; const double PI=acos(-1); /* // abs(√a - √b) <= √cかどうかの判定 bool f(ll a,ll b,ll c){ //ll t=(a*a + b*b + c*c)-(2*a*b + 2*b*c + 2*c*a); //t<=0 ll t=(a*a + b*b + c*c)/(2*a*b +2*b*c + 2*c*a); return t<=1; } */ /* bool f(ll a,ll b,ll c){ double t=abs(sqrt(a) - sqrt(b)); return t<=sqrt(c); } */ /* bool f(ll a,ll b,ll c){ //ll t=(a+b-c)/(4*a*b)*(a+b-c); ll t=(a+b-c)*(a+b-c)/(4*a*b); ll u=abs(a+b-c); if(u>(INF/u)) return false; return t<=1; } */ int main(){ int N; cin>>N; vector X(N),Y(N),T(N); for(int i=0;i>X[i]>>Y[i]>>T[i]; } auto f=[&](int i,int j){ if(T[i]==T[j]){ return (X[i]-X[j])*(X[i]-X[j])+(Y[i]-Y[j])*(Y[i]-Y[j]); }else{ //T[i]!=T[j] ll Ri=X[i]*X[i]+Y[i]*Y[i]; ll Rj=X[j]*X[j]+Y[j]*Y[j]; ll ub=1LL<<32; ll lb=-1; while(ub-lb>1){ ll mid=(ub+lb)/2; ll c=mid*mid; if(4*Ri*Rj=c lb=mid; } } ll r=lb; //return R1s+R2s - (ll)sqrt(4*R1s*R2s); return Ri+Rj - r; } }; ll ub=INF,lb=-1; //ubのときは必ずたどり着ける while(ub>lb+1){ ll mid=(ub+lb)/2; vector> G(N); for(int i=0;i