HDU 1140 War on Weather

http://acm.hdu.edu.cn/showproblem.php?pid=1140

根据球的切点做简单的判断，一开始找到答案没有跳出循环，各种错。。。

#include <iostream>
#include <cmath>
using namespace std ;
const double PI=acos(-1.0) ;//π的表示方法 
struct point 
{
    double x,y,z ;
} ;
double dis(point p1,point p2)
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y)+(p1.z-p2.z)*(p1.z-p2.z)) ;
}
int main()
{
    int k,m ;
    while(scanf("%d%d",&k,&m),(k||m))
    {
        point O ;
        O.x=O.y=O.z=0 ;
        point p[101] ;
        for(int i=0;i<k;i++)
            scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].z) ;
        double r=20000.0/PI ;
        int ans=0 ;
        while(m--)
        {
            point target ;
            scanf("%lf%lf%lf",&target.x,&target.y,&target.z) ;
            for(int i=0;i<k;i++)
                if(dis(p[i],O)*dis(p[i],O)>=r*r+dis(p[i],target)*dis(p[i],target))
                { 
                    ans++ ;
                    break ;
                }
        }
        printf("%d\n",ans) ;
    }
    return 0 ;
}