【poor 几何】UVA1606 转啊转啊转

题目

可怕，可怕
atan2教做人，之前根本不知道有这个函数
看了刘汝佳的代码，又看了这篇题解，这孩子在评论区放的图很亮
向量差积教做人，这个性质早忘了

#include<cmath>
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=1100;
struct point{
    int x,y,col;double rad;
    point(int x=0,int y=0):x(x),y(y){}
    bool  operator < (const point &a)const
        {return rad<a.rad;}
    bool  operator ^ (const point &a)const
        {return x*a.y-y*a.x>=0;}
    point operator - (const point &a)const
        {return point(x-a.x,y-a.y);}
}a[N],p[N],zero;
int n,L,R,cnt;
int main(){
    freopen("fuck.in","r",stdin);
    zero.x=zero.y=0;
    for (;~scanf("%d",&n)&&n;){
        for (int i=0;i<n;i++)
          scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].col);
        int ans=0;
        for (int i=0;i<n;i++){
            int P=0;

            for (int j=0;j<n;j++){
                if (i==j)continue;
                p[P]=a[j]-a[i];
                if (a[j].col)p[P]=zero-p[P];
                p[P].rad=atan2(p[P].y,p[P].x);
                P++;
            }
            sort(p,p+P);
            for (int j=0;j<P;j++)
            //printf("x=%d y=%d rad=%.2lf\n",p[j].x,p[j].y,p[j].rad); 
            L=0;R=0;cnt=2;
            for (;L<P;){
                if (L==R){R=(R+1)%P;cnt++;}
                while(L!=R&&p[L]^p[R]){R=(R+1)%P;cnt++;}
                L++; cnt--; ans=max(ans,cnt);
                //printf("L=%d R=%d ans=%d\n",L-1,R,ans);
            }
            //printf("base=%d ans=%d\n",i,ans);
        }
        printf("%d\n",ans);
    }
    return 0;
}