POJ 1106 极角排序

题意：

给定半圆的圆心和半径和n个点，求能在半圆内的最多的点数。

题解：

先把所有的在圆外的点删除，然后求出所有的点关于圆心的极角，然后排序，因为是环形的，极角复制一份加在原来的后面，维护双指针即可~

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>

#define PI 3.141592653589
#define EPS 1e-7
#define N 600

using namespace std;

struct PO
{
    double x,y,angle;
}p[N];

int n,ans;
double sr;

inline bool cmp(const PO &a,const PO &b)
{
    return a.angle<b.angle;
}

inline int doublecmp(double x)
{
    if(x>EPS) return 1;
    else if(x<-EPS) return -1;
    return 0;
}

inline double get_dis2(PO &a,PO &b)
{
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}

inline void read()
{
    n=0;
    int tmp;scanf("%d",&tmp);
    for(int i=1;i<=tmp;i++)
    {
        ++n;
        scanf("%lf%lf",&p[n].x,&p[n].y);
        if(doublecmp(get_dis2(p[n],p[0])-sr*sr)==1) n--;
        else p[n].angle=atan2(p[n].y-p[0].y,p[n].x-p[0].x);
    }
}

inline void go()
{
    ans=0;
    sort(p+1,p+1+n,cmp);
    for(int i=1;i<=n;i++) p[i+n].angle=p[i].angle+2*PI;
    for(int g1=1,g2=1;g1<=n&&g2<=n+n;g2++)
        if(doublecmp(p[g2].angle-p[g1].angle-PI)==1)
        {
            ans=max(g2-g1,ans);
            g1++;
        }
    printf("%d\n",ans);
}

int main()
{
    while(scanf("%lf%lf%lf",&p[0].x,&p[0].y,&sr),doublecmp(sr)>=0) read(),go();
    return 0;
}