hdu 4962 Closed Paths 计算几何 2014 Multi-University Training Contest 9-1003

题意：

有n个点pi，表示一个环路的折点。一个女孩从p1走到p2，再走到p3，。。。，最后从pn走回p1，走过的路可以重叠或者部分重叠。现在要求在x轴上找的一个点作为男孩的位置，这个位置不能是女孩的路线上，男孩在女孩的散步过程中一直盯着女孩，所以会逆时针ap度，以及顺时针选择bp度，我们可以发现ap-bp是360度的倍数，现在要求找到x轴的一个点使得(ap-bp)/360最大，即逆时针转的最大圈数。输出最大圈数。

题解：

我们观察图形，可以发现一个点的旋转次数等于其左边的折线自上而下穿越x轴的次数减去自下而上穿越的次数（因为图像成环，所以上下穿越次数必定相同，那么计算前面的次数就可以知道女孩绕该点转了几圈，也就是男孩转的圈数）。

之后，只要找出图形与x轴的所有交点，然后判断交点是自上而下就标记+1，自下而上标记-1，点在x轴上未穿越记为0（如线在x轴上，则记两个端点），将所有点根据x轴排序。接着保存所有在x轴上的线段，根据x坐标排序；扫描所有的交点，找两个交点的中心，判断中心是否在路线上（用前面保存的在x轴上的线段判断），在则无效；不在则求该坐标前的所有交点标记总和即可。输出最大的标记总和。

注意：如果点不在图形内部，则必定是旋转0次的。

代码：

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <map>
#include <set>
#include <queue>
using namespace std;
//基础点和向量运算
struct Point{
    double x,y;
    Point(double x=0,double y=0):x(x),y(y){}
};
typedef Point Vector;
Vector operator + (Vector A,Vector B){return Vector(A.x+B.x,A.y+B.y);}
Vector operator - (Vector A,Vector B){return Vector(A.x-B.x,A.y-B.y);}
Vector operator * (Vector A,double p){return Vector(A.x*p,A.y*p);}
Vector operator / (Vector A,double p){return Vector(A.x/p,A.y/p);}
bool operator <(const Point& a, const Point& b)
{
    return a.x<b.x||(a.x==b.x&&a.y<b.y);
}
const double eps=1e-10;
int dcmp(double x)//判断正负，或者等于0
{
    if(fabs(x)<eps)return 0;else return x<0?-1:1;
}
bool operator==(const Point& a,const Point &b)
{
    return dcmp(a.x-b.x)==0&&dcmp(a.y-b.y)==0;
}
double Dot(Vector A, Vector B){return A.x*B.x+A.y*B.y;}//点积
double Length(Vector A){return sqrt(Dot(A,A));}//OA长
double Angle(Vector A,Vector B){return acos(Dot(A,B)/Length(A)/Length(B));}//OA和OB的夹角
double Cross(Vector A,Vector B){return A.x*B.y-A.y*B.x;}//叉积

//点和直线
//P+tv表示一条直线，P为点，tv为方向向量
Point GetLineIntersection(Point P,Vector v,Point Q,Vector w)//求直线交点，确保存在交点，即Cross(v,w)非0
{
    Vector u=P-Q;
    double t=Cross(w,u)/Cross(v,w);
    return P+v*t;
}



const int maxn=1e3+10;
struct node{
    double x;
    int dir;
}e[maxn];
Point p[maxn],f[maxn];
int t;
int cmp(node a,node b)
{
    return a.x<b.x;
}
int cmp2(Point a,Point b)
{
    if(a.x==b.x)return a.y<b.y;
    return a.x<b.x;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int i,j,k,n,tot=0;
        scanf("%d",&n);
        for(i=0;i<n;i++)
            scanf("%lf%lf",&p[i].x,&p[i].y);
        p[n]=p[0];
        p[n+1]=p[1];
        Point x(1,0),out,O(0,0),c;
        t=0;
        for(i=1;i<=n;i++)
        {
            if(dcmp(p[i].y)==0)
            {
                //printf("*%d\n",i);
                if(dcmp(p[i-1].y)*dcmp(p[i+1].y)<0)
                {
                    e[t].dir=dcmp(p[i-1].y);
                    e[t++].x=p[i].x;
                }
                else
                {
                    e[t].dir=0;
                    e[t++].x=p[i].x;
                }
            }
        }
        for(i=0;i<n;i++)
        {
            if(dcmp(p[i].y)==0&&dcmp(p[i+1].y)==0)
            {
                f[tot].x=min(p[i].x,p[i+1].x);
                f[tot++].y=max(p[i].x,p[i+1].x);
                continue;
            }
            if(dcmp(p[i].y)!=0&&dcmp(p[i+1].y)!=0)
            {
                if(dcmp(p[i].y*p[i+1].y)<0)
                {
                    out=GetLineIntersection(O,x,p[i],p[i+1]-p[i]);
                    //printf("out:%f %f %d %d\n",out.x,out.y,i,i+1);
                    e[t].x=out.x;
                    if(dcmp(p[i].y)>0)e[t++].dir=1;
                    else e[t++].dir=-1;
                }
            }
        }
        sort(e,e+t,cmp);
        sort(f,f+tot,cmp2);
        //for(i=0;i<t;i++)printf("**%d %f %d\n",i,e[i].x,e[i].dir);
        //for(i=0;i<tot;i++)printf("***%d %f %f\n",i,f[i].x,f[i].y);
        int ans=0,now=0;
        for(i=0,j=0;i<t-1;i++)
        {

            c.x=(e[i].x+e[i+1].x)/2;
            c.y=0;
            now+=e[i].dir;
            while(j<tot&&dcmp(e[i].x-f[j].y)>0)j++;
            if(j<tot&&dcmp(e[i].x-f[j].x)>=0)continue;
            ans=max(ans,now);
        }
        printf("%d\n",ans);
    }
    return 0;
}
/*
10
5
1 1
-1 0
1 -1
0 1
0 -1
6
-1 1
-1 -1
1 -1
1 1
2 0
-2 0
*/