LA 2963 Hypertransmission

题意：需要在n个星球上各安装一个广播装置，作用范围均为R。每个星球广播A节目或者B节目，若一个星球收听到的和自己节目不同的广播数量大于和自己广播相同的数量则说明这个星球是不稳定的，求在不稳定星球数量最多的情况下R的最小值。

分析:先求出所有星球间两两的距离后排序，遍历一遍维护最大值，O（n*n）。

#include<iostream>
#include<string>
#include<algorithm>
#include<cstdlib>
#include<cstdio>
#include<set>
#include<map>
#include<vector>
#include<cstring>
#include<stack>
#include<cmath>
#include<queue>
#define INF 0x3f3f3f3f
#define MAXN 1005
#define sqr(x) ((x)*(x)) 
using namespace std;
int n,tot,x[MAXN],y[MAXN],z[MAXN],p[MAXN],f1[MAXN],f2[MAXN],jud[MAXN];
struct Edge
{
	int u,v;
	double dis;
	Edge() { }
	Edge(int x,int y,double val)
	{
		u = x;
		v = y;
		dis = val;
	}
	friend bool operator < (Edge a,Edge b)
	{
		return a.dis < b.dis;
	}
} s[MAXN*MAXN];
int main()
{
	while(scanf("%d",&n) != EOF)
	{
		tot = 0;
		memset(f1,0,sizeof(f1));
		memset(f2,0,sizeof(f2));
		memset(jud,0,sizeof(jud));
		for(int i = 1;i <= n;i++) scanf("%d%d%d%d",x+i,y+i,z+i,p+i);
		for(int i = 1;i < n;i++)
		 for(int j = i + 1;j <= n;j++)
		  s[++tot] = Edge(i,j,sqrt(sqr(*(x+i)-*(x+j))+sqr(*(y+i)-*(y+j))+sqr(*(z+i)-*(z+j))));
		sort(s+1,s+1+tot);
		int now = 0,ans = 0;
		double ans2 = 0;
		for(int i = 1;i <= tot;i++)
		{
			int u = s[i].u,v = s[i].v;
			if(jud[u]) 
			{
				now--;
				jud[u] = false;
			}
			if(jud[v]) 
			{
				now--;
				jud[v] = false;
			} 
			if(p[u] == p[v])
			{
				f1[u]++;
				f1[v]++;
			}
			else 
			{
				f2[u]++;
				f2[v]++;
			}
			if(f1[u] + 1 < f2[u])
			{
				jud[u] = true;
				now++;
			}
			if(f1[v] + 1 < f2[v])
			{
				jud[v] = true;
				now++;
			}
			if(i != tot && s[i].dis == s[i+1].dis) continue;
			if(now > ans)
			{
				ans = now;
				ans2 = s[i].dis; 
			}
		}
		printf("%d\n",ans);
		printf("%.6lf\n",ans2);
	}	
}