高效算法设计专项：LA 2963

刚开始完全没思路，看了训练指南上的提示说是扫描法，然后往这方面想。仔细分析题目发现可以用O(n^2)的时间解决。不难发现可能的结果只可能是n个点两两距离的，然后可以从小到大扫描这n*(n-1)/2条边，然后维护当前符合条件的planet的个数。需要注意处理边长相等的情况。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
#define M 1010
const double eps=1e-6;
int dcmp(double x)
{
if(fabs(x)<eps) return 0;
else return x<0?-1:1;
}
struct point
{
int x,y,z,p;
} a[M];
struct edge
{
int u,v;
double l;
edge(int u=0,int v=0,double l=0):u(u),v(v),l(l){}
bool operator<(const edge& tmp) const
{return dcmp(l-tmp.l)<0;}
} e[M*M];
double length(int i,int j)
{
int x=a[i].x-a[j].x,y=a[i].y-a[j].y,z=a[i].z-a[j].z;
return sqrt(x*x+y*y+z*z);
}
int n,m;
int c[M];
int main()
{
freopen("in.txt","r",stdin);
while(scanf("%d",&n)!=EOF)
{
for(int i=0;i<n;i++)
{
scanf("%d%d%d%d",&a[i].x,&a[i].y,&a[i].z,&a[i].p);
c[i]=1;
}
m=0;
for(int i=0;i<n;i++)
{
for(int j=i+1;j<n;j++)
{
e[m++]=edge(i,j,length(i,j));
}
}
sort(e,e+m);
int ans=0,cur=0;
double len=0;
for(int i=0;i<m;i++)
{
int u=e[i].u,v=e[i].v;
double l=e[i].l;
if(a[u].p!=a[v].p)
{
c[u]--;
c[v]--;
if(c[u]==-1) cur++;
if(c[v]==-1) cur++;
}
else
{
c[u]++;
c[v]++;
if(c[u]==0) cur--;
if(c[v]==0) cur--;
}
if(i!=m-1&&dcmp(l-e[i+1].l)==0) continue;
if(cur>ans)
{
ans=cur;
len=l;
}
}
printf("%d\n",ans);
printf("%.4f\n",len);
}
return 0;
}