BNU校赛总决赛J 小白兔小灰兔 相交计算几何模板

J 小白兔小灰兔

时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 32768K，其他语言65536K
Special Judge, 64bit IO Format: %lld

题目描述

老山羊伯伯在地里收白菜，小白兔和小灰兔看见了就一起来帮忙。

他们干了半天，终于干完了。

羊伯伯：小灰兔，这车白菜送给你！

小灰兔：谢谢羊伯伯！

羊伯伯：小白兔，我也送你一车白菜！

小白兔：我不要白菜！给我一包白菜种子吧！

羊伯伯：好！给你~

小白兔：谢谢羊伯伯~

小灰兔把白菜拉到家里之后，就跟大家梦想中的生活一样，躺着啥都不干，吃吃吃，吃了玩~（好想一辈子都这样啊~小灰兔心想。）

小白兔把种子拿回去，打算开始勤劳地种白菜。然而他发现不是所有土地都能用来种白菜，只有被阳光照到的地方可以种白菜。

小白兔生活的星球可以看作二维平面中的一个简单多边形，太阳可以看作一个点。小白兔想知道，这个星球上一共有长度多少的土地可以用来种白菜。

输入描述:

第一行是一个正整数T(≤ 20)，表示测试数据组数，

对于每组测试数据，

第一行是一个整数n(3≤ n≤ 10)，表示简单多边形的顶点数，

接下来n行，每行是两个整数x_i,y_i(-10≤ x_i,y_i≤10)，按照逆时针绕向给出简单多边形的n个顶点的坐标，

最后一行是两个整数x,y(-10≤ x,y≤10)，表示太阳的坐标，保证太阳在多边形外且不在多边形任意一条边所在直线上。

输出描述:

对于每组测试数据，输出一行，包含一个实数，表示可以用来种白菜的土地的总长度，要求相对误差不超过，

也就是说，令输出结果为a，标准答案为b，若满足，则输出结果会被认为是正确答案。

示例1

输入

2
4
0 0
2 0
2 2
0 2
4 1
4
0 0
2 0
2 2
0 2
4 4

输出

2.000000000000
4.000000000000

---

#include<bits/stdc++.h>
#define clr(x) memset(x,0,sizeof(x))
#define clr_1(x) memset(x,-1,sizeof(x))
#define mod 1000000007
#define INF 0x3f3f3f3f
#define LL long long
#define pb push_back
#define pbk pop_back
#define ls(i) (i<<1)
#define rs(i) (i<<1|1)
#define mp make_pair
using namespace std;
typedef double db;
const int N=5e5+10;
const db eps=1e-9;
int equzero(db t)//带精度大小判断
{
    if(t>eps) return 1;
    if(t<-eps) return -1;
    return 0;
}
struct Point
{
    db x,y;
    Point(){}
    Point(db _x,db _y):x(_x),y(_y) {}
    Point operator +(const Point &b) const //+
    {
        return Point(x+b.x,y+b.y);
    }
    Point operator -(const Point &b) const//-
    {
        return Point(x-b.x,y-b.y);
    }
    Point operator *(const db &b) const// 放大
    {
        return Point(x*b,y*b);
    }
    Point operator /(const db &b) const//缩小
    {
        return Point (x/b,y/b);
    }
    db dot(const Point &b) const//点积
    {
        return x*b.x+y*b.y;
    }
    db det(const Point &b) const//叉积
    {
        return x*b.y-b.x*y;
    }
    bool operator ==(const Point &b) const
    {
        return  equzero(x-b.x)==0 && equzero(y-b.y)==0;
    }
    bool operator !=(const Point &b) const
    {
        return !(Point(x,y)==b);
    }
    db mo() //模
    {
        return sqrt(x*x+y*y);
    }
}pt[N],u,v,mid;
struct Line
{
    Point a,b;
    Line (){}
    Line (Point _a,Point _b):a(_a),b(_b) {};
    bool quickrej(const Line &t) const//快速排斥判断
    {
        return equzero(min(a.x,b.x)-max(t.a.x,t.b.x))<=0 && equzero(min(t.a.x,t.b.x)-max(a.x,b.x))<=0 && equzero(min(a.y,b.y)-max(t.a.y,t.b.y))<=0 && equzero(min(t.a.y,t.b.y)-max(a.y,b.y))<=0;
    }
    db len() const
    {
        return (b-a).mo();
    }
    bool stra(const Line &t) const //跨立判断
    {
        return equzero((a-t.a).det(t.b-t.a)*(b-t.a).det(t.b-t.a))<0 && equzero((t.a-a).det(b-a)*(t.b-a).det(b-a))<0 ;
    }
    bool gfxj(const Line &t) const //规范相交
    {
        return quickrej(t) && stra(t);
    }
    bool on(const Point &t) const //点在线段上
    {
        Line p=Line(a,b);
        if(p.a.x>p.b.x)
            swap(p.a,p.b);
        if(equzero(t.x-p.a.x)<0 || equzero(t.x-p.b.x)>0)
            return false;
        if(equzero((p.a-t).det(p.b-t))==0)
            return true;
        else
            return false;
    }
    bool kddxj(const Line &t) const //可多点相交
    {
        return quickrej(t) && (stra(t) || on(t.a) || on(t.b) || t.on(a) || t.on(b));
    }
    bool para(const Line &t) const // 平行或在一条直线上
    {
        return equzero((a.y-b.y)*(t.a.x-t.b.x)-(t.a.y-t.b.y)*(a.x-b.x))==0 ;
    }
    db operator & (const Line &t) const //返回切点占直线比例,距a比例
    {
        if(equzero((a-b).det(t.a-t.b))==0)return -INF;
        db tk=((a-t.a).det(t.a-t.b))/((a-b).det(t.a-t.b));
        return equzero(tk)>=0 && equzero(tk-1)<=0 ? tk : -INF;
    }
}le[N],lu,lv;
int n,m,T;
vector<db> cp;
db ans,ins,rate;
bool flag;
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=0;i<=n;i++)
            scanf("%lf%lf",&pt[i].x,&pt[i].y);
        for(int i=0;i<n;i++)
            le[i]=Line(pt[i],pt[(i+1)%n]);
        ans=0;
        for(int i=0;i<n;i++)
        {
            cp.clear();
            for(int j=0;j<n;j++)
            {
                ins=le[i]&Line(pt[n],pt[j]);
                if(ins>=-1) cp.pb(ins);
            }
            sort(cp.begin(),cp.end());
            rate=0;
            for(int j=0;j<cp.size()-1;j++)
            {
                lu=Line(pt[n],le[i].a+(le[i].b-le[i].a)*((cp[j]+cp[j+1])/2));
                flag=0;
                for(int k=0;k<n;k++)
                {
                    if(lu.gfxj(le[k]))
                    {
                        flag=1;
                        break;
                    }
                }
                if(!flag)
                {
                    rate+=cp[j+1]-cp[j];
                }
            }
            ans+=rate*le[i].len();
        }
        printf("%.12f
",ans);
    }
    return 0;
}