2544.cpp

自己的ac代码 在acm.zju.edu.cn 上的题目

#include<stdio.h>
#include<math.h>
struct point 
{
	double x,y,z;
};
int n,m;
point a[105],b[105];
double sqr(double x)
{
	return x*x;
}

double dis(point a,point b)
{
	return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y)+sqr(a.z-b.z));
}

int main()
{
	int i,j,k;
	point cen;
	cen.x=0;cen.y=0;cen.z=0;
	while(scanf("%d%d",&n,&m))
	{
		if(n==0 && m==0)break;
		for(i=1;i<=n;i++)
			scanf("%lf%lf%lf",&a[i].x,&a[i].y,&a[i].z);
		for(j=1;j<=m;j++)
			scanf("%lf%lf%lf",&b[j].x,&b[j].y,&b[j].z);
		int ans=0;
		for(j=1;j<=m;j++)
		{
			for(i=1;i<=n;i++)
			{
				double o,p,q;
				o=dis(a[i],b[j]);
				p=dis(a[i],cen);
				q=40000.0/acos(-1.0)/2.0;
				if(p*p-q*q>=o*o){ans++;break;}
			}
		}
		printf("%d\n",ans);
	}
	return 0;
}