ZOJ Problem Set - 1041

这个传说中的计算几何，让我搞的话真有点困难，因为神马叉乘判断神马的，以前只是知道名字，具体怎么用从来没有做过，嘿嘿！

这算一次尝试吧，自己想了想，感觉应该直接暴力，怕超时，看了看别人的报告，发现居然可以搞定！

因为要求在一个半圆里面，最多的点的个数，这就要求我们怎么样才知道它是在我们现在所画的半圆内，这样的话，需要我们首先把不在这个圆外的点去除，针对在这个圆内的点，对每个点和圆心构成直 径 的顺时针半圆，用叉乘来判断是否在里面，这样，暴解就搞定了！！！（PS:怎么样通过叉乘判断是在顺时针的半圆内呢，这个是通过那个叉乘的数学性质了，如果两个向量的角度小于180，大于0，则那个乘积为正，反之为负的说！嘿嘿！）

#include<stdio.h> #define N 155 struct node{ double x,y; int sign; }s[N]; double x,y,r; double is(int i) { return (s[i].x-x)*(s[i].x-x)+(s[i].y-y)*(s[i].y-y); } double xmult(int i,int j) { return (s[i].x-x)*(s[j].y-y)-(s[j].x-x)*(s[i].y-y); } int main() { int i,j,n,count,max; while(scanf("%lf%lf%lf",&x,&y,&r),r>0) { scanf("%d",&n); r*=r; for(i=0;i<n;i++) s[i].sign=0; for(i=0;i<n;i++) { scanf("%lf%lf",&s[i].x,&s[i].y); if(is(i)<=r)s[i].sign=1; } max=0; for(i=0;i<n;i++) { if(s[i].sign==1) { count=0; for(j=0;j<n;j++) { if(s[j].sign==1&&xmult(i,j)>=0) count++; } if(count>max)max=count; } } printf("%d/n",max); } return 0; }