3090613_ce.c

北大大牛代码 1240道题的原代码 超级权威

#include <stdio.h>
#include <math.h>

int cnt;
double x[101], y[101], z[101], d[101];
double r = 20000.00*20000.00/(acos(-1)*acos(-1));

int main()
{
	int k, m, i;
	double a, b, c, d1;

	while(scanf("%d%d",&k,&m)==2)
	{
		if(m==0&&k==0)
			break;
		for(i = 0; i < k; i++)
		{
			scanf("%lf%lf%lf",&x[i],&y[i],&z[i]);
			d[i] = sqrt(x[i]*x[i]+y[i]*y[i]+z[i]*z[i]-r);
		}
		cnt = 0;
		while(m-- > 0)
		{
			scanf("%lf%lf%lf",&a,&b,&c);
			for(i = 0; i < k; i++)
			{
				d1 = sqrt((a-x[i])*(a-x[i])+(b-y[i])*(b-y[i])+(c-z[i])*(c-z[i]));
				d1 = d1-d[i];
				if(d1 <= 1e-8)
				{
					cnt++;
					break;
				}
			}
		}
		printf("%d\n",cnt);
	}
	return 0;
}