Angle Beats HDU - 6731

做这道之前可以考虑做一下hdu5785

刚看到这一题想到的做法是(n+q)2log(n+q)离线操作，对每个点做一次极角排序，然后找到匹配的点，对于找90度我们用尺取法，定义两个指针k1，k2。我们让k1遍历每个点k2作为小于90度的第一个点的坐标，下一个点开始时，k1从k2这个点开始遍历。

然后我刚开始用的是atan2这个极角排序，后来wa了应该是被卡精度了

参考了https://blog.csdn.net/EIP_silly/article/details/101673375

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

using namespace std;

const int MAXN=4444;
const double eps=1e-8;
int sgn(double x)
{
    if(fabs(x)<eps) return 0;
    if(x<0) return -1;
    else return 1;
}
struct Point{
    double x,y,angle;
    int fg,id;
    Point(){}
    Point(double _x,double _y){
        x=_x;y=_y;
    }
    bool operator <(const Point &b)const {
        return angle<b.angle;
    }
    double operator ^(const Point &b)const {
        return x*b.y-y*b.x;
    }
    double operator *(const Point &b)const{
        return x*b.x+y*b.y;
    }
};
int qua(Point a)
{
    if(a.x>0 && a.y>=0)  return 1;
    if(a.x<=0 && a.y>0)  return 2;
    if(a.x<0 && a.y<=0)  return 3;
    if(a.x>=0 && a.y<0)  return 4;
}

bool cmp1(const Point &a,const Point &b)
{
    return (a^b) > 0;//叉积
}

bool cmp(const Point &a,const Point &b)  //先按象限从小到大排序 再按极角从小到大排序
{
    if(qua(a)==qua(b))//返回值就是象限
        return cmp1(a,b);
    else return qua(a)<qua(b);
}
Point p[MAXN],p1[MAXN],p2[MAXN];
int ans[MAXN];
int main()
{
    int n,q;
    while(scanf("%d%d",&n,&q)==2)
    {
        memset(ans,0,sizeof(ans));
        for (int i = 0; i < n; i++) scanf("%lf%lf", &p[i].x, &p[i].y), p[i].fg = 0;
        for (int i = 0; i < q; i++) scanf("%lf%lf", &p1[i].x, &p1[i].y), p1[i].fg = 1, p1[i].id = i;
        for (int i = 0; i < q; i++)
        {
            for(int j=0;j<n;j++)
            {
                p2[j].x=p[j].x-p1[i].x;
                p2[j].y=p[j].y-p1[i].y;
                p2[j].angle=atan2(p2[j].y,p2[j].x);
            }
            sort(p2,p2+n,cmp);
            int k1=0,k2=0;
            for(int j=0;j<n;j++)
            {
                k1=k2;
                if(k1==j) k2=k1=(j+1)%n;
                while(k1!=j&&(p2[k1]*p2[j]>=0)&&(p2[j]^p2[k1])>=0)
                {
                    if(p2[k1]*p2[j]==0)
                    {
                        ans[i]++;
                        k1=(k1+1)%n;
                    }
                    else k2=k1=(k1+1)%n;
                }

            }
        }
        for(int i=0;i<n;i++)
        {
            int cont=0;
            for(int j=0;j<n;j++)
            {
                if(j==i) continue;
                p2[cont].x=p[j].x-p[i].x;
                p2[cont].y=p[j].y-p[i].y;
                p2[cont].fg=p[j].fg;
                p2[cont].id=p[j].id;
                p2[cont].angle=atan2(p2[cont].y,p2[cont].x);
                cont++;
            }
            for(int j=0;j<q;j++)
            {
                p2[cont].x=p1[j].x-p[i].x;
                p2[cont].y=p1[j].y-p[i].y;
                p2[cont].id=p1[j].id;
                p2[cont].fg=p1[j].fg;
                p2[cont].angle=atan2(p2[cont].y,p2[cont].x);
                cont++;
            }
            sort(p2,p2+cont,cmp);
            int k1=0,k2=0;
            for(int j=0;j<cont;j++)
            {
                k1=k2;
                if(k1==j) k2=k1=(j+1)%cont;
                while(k1!=j&&(p2[k1]*p2[j]>=0)&&(p2[j]^p2[k1])>=0)
                {
                    if(p2[k1]*p2[j]==0)
                    {
                        int m=p2[k1].fg+p2[j].fg;
                        if(m==1)
                        {
                            if(p2[k1].fg) ans[p2[k1].id]++;
                            else if(p2[j].fg) ans[p2[j].id]++;
                        }
                        k1=(k1+1)%cont;
                    }
                    else k2=k1=(k1+1)%cont;
                }
            }
        }
        for(int i=0;i<q;i++)
            printf("%d\n",ans[i]);
    }
    return 0;
}