点在直线那一端

 poj 1106   给定的点中，求半圆的覆盖最多的点数。

 

Point sec,temp,all[320];
double r;
int multiply(Point p2,Point p0)//叉积，其中有一点是原点。
{
    return((p2.x-p0.x)*p0.y-p0.x*(p2.y-p0.y));
}

int fun(int n)
{
    int coun,i,j,Max = 0;
    for(i = 0;i < n;i++)
    {
        coun = 0;
        for(j = 0;j < n;j++)
        {
            if(i!=j)
            {
                if(multiply(all[j],all[i])>=0)//判断两点是否在直线的同一端，直线过原点。
                    coun++;
            }
        }
        Max = max(coun,Max);
    }
    return Max;
}
int main()
{
    int n,m,i,j,k;
    double dis;
    while (scanf("%d%d%lf",&sec.x,&sec.y,&r)!=EOF,r>-eps)
    {
        scanf("%d",&n);
        k = 0;
        for(i = 0;i < n;i++)
        {
            scanf("%d%d",&temp.x,&temp.y);
            temp.x -= sec.x;//将圆心移至原点处。
            temp.y -= sec.y;
            dis = (double)(temp.x*temp.x+temp.y*temp.y);
            if(r*r-dis>-eps)
            {
                all[k++] = temp;
            }
        }
        m = fun(k);
        if(k==0)
            printf("0\n");
        else
            printf("%d\n",m+1);
    }
    return 0;
}