UVALive 2963 Hypertransmission

题目链接：https://vjudge.net/problem/UVALive-2963

题目大意：有n个星球，每个星球坐标为（xi，yi，zi），可以看成一个点。每个星球广播A类节目或B类节目，广播范围为R（以该星球为中心半径为R的球体）。令N+(i)表示星球i听到的和自己广播相同节目的星球数（包括自己），N-(i)表示星球i听到的和自己广播不同节目的星球数。如果N+(i)<N-(i)，则称星球i是不稳定的。求不稳定星球的最大数目以及在此前提下R的最小值。

思路：容易知道，R的值肯定为某两个星球之间的距离。因此将所有距离从小到大排序，并存储每个距离所对应的两个点。按顺序枚举每个距离，根据这些距离所对应的点更新N+(i)和N-(i)即可。又因为半径大的球会包含半径小的球，所以前面的结果可以适用于后面。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cctype>
#include <cstdlib>
#include <map>
#include <iostream>
#include <queue>
using namespace std;
typedef long long LL;
const int dr[] = {0,0,-1,1};
const int dc[] = {-1,1,0,0};
const int maxn = 1000 + 10;

struct Point
{
   int x, y, z, t;
}p[maxn];

struct Triple
{
	int d, pa, pb;
	bool operator < (const Triple& rhs) const {
		return d < rhs.d;
	}
	bool operator == (const Triple &rhs) const {
		return d == rhs.d;
	}
	Triple(int d = 0, int pa = 0, int pb = 0) : d(d),pa(pa),pb(pb){}
}tri[maxn*maxn];

int dist(const Point& A, const Point& B)
{
	return (A.x-B.x)*(A.x-B.x) + (A.y-B.y)*(A.y-B.y) + (A.z-B.z)*(A.z-B.z);
}

int vis[maxn];

int main()
{
	int n;
	while(scanf("%d", &n) == 1)
	{
		for(int i = 0; i < n; i++) scanf("%d%d%d%d", &p[i].x, &p[i].y, &p[i].z, &p[i].t);
		int cnt = 0;
	    for(int i = 0; i < n; i++)
	       for(int j = i+1; j < n; j++)
	       	  tri[cnt++] = Triple(dist(p[i],p[j]), i, j);
	    sort(tri, tri+cnt);
        fill(vis, vis+n+2, 1);
        int temp = 0, ans = 0, R = 0, j;
        for(int i = 0; i < cnt;) {
           for(j = i; j < cnt && tri[j].d == tri[i].d; j++)
           {
           	  int ta = tri[j].pa, tb = tri[j].pb;
           	  if(p[ta].t != p[tb].t)
           	  {
           	  	 if(--vis[ta] == -1) temp++;
           	  	 if(--vis[tb] == -1) temp++;
           	  }
           	  else {
           	  	 if(++vis[ta] == 0) temp--;
           	  	 if(++vis[tb] == 0) temp--;
           	  }
           }
           if(temp > ans) {
           	  ans = temp;
           	  R = tri[i].d;
           }
           i = j;
        }
        printf("%d\n%.4f\n", ans, sqrt(R*1.0));
	}
	return 0;
}