BZOJ 1137 半平面交

半平面交的板子

//By SiriusRen
#include <bits/stdc++.h>
#define double long double
using namespace std;
const int N=;
const double eps=1e-;
int n,m,xx,yy,tot;
double Ans;
vector<int>vec[N];
struct Point{double x,y;}point[N];
struct Line{Point a,b;double angle;}line[N],q[N];
void addline(Line &l,Point a,Point b){
    l.a=a,l.b=b,l.angle=atan2(b.y-a.y,b.x-a.x);
}
Point operator-(Point a,Point b){
    Point c;c.x=a.x-b.x,c.y=a.y-b.y;return c;
}
double operator*(Point a,Point b){
    return a.x*b.y-a.y*b.x;
}
bool operator<(Line a,Line b){
    if(a.angle==b.angle)return (b.b-a.a)*(b.a-a.a)>;
    return a.angle<b.angle;
}
Point inter(Line a,Line b){
    double k1,k2,t;
    k1=(a.b-b.a)*(b.b-b.a);
    k2=(b.b-b.a)*(a.a-b.a);
    t=k1/(k1+k2);
    Point ans;
    ans.x=a.b.x+(a.a.x-a.b.x)*t;
    ans.y=a.b.y+(a.a.y-a.b.y)*t;
    return ans;
}
double dis(Point x,Point y){
    x=y-x;
    return sqrt(x.x*x.x+x.y*x.y);
}
bool judge(Line a,Line b,Line t){
    Point p=inter(a,b);
    return (t.a-p)*(t.b-p)<;
}
void bpmj(){
    sort(line+,line++tot);
    n=;
    for(int i=;i<=tot;i++){
        if(abs(line[i].angle-line[i-].angle)>eps)n++;
        line[n]=line[i];
    }
    int r=,l=;
    q[]=line[],q[]=line[];
    for(int i=;i<=n;i++){
        while(l<r&&judge(q[r],q[r-],line[i]))r--;
        while(l<r&&judge(q[l],q[l+],line[i]))l++;
        q[++r]=line[i];
    }
    while(l<r&&judge(q[r],q[r-],q[l]))r--;
    while(l<r&&judge(q[l],q[l+],q[r]))l++;
    q[r+]=q[l],tot=;
    for(int i=l;i<=r;i++)point[++tot]=inter(q[i],q[i+]);
    point[++tot]=point[];
    for(int i=;i<tot;i++)Ans+=dis(point[i],point[i+]);
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;i++)
        scanf("%Lf%Lf",&point[i].x,&point[i].y);
    Ans=-dis(point[],point[n]);
    for(int i=;i<=m;i++){
        scanf("%d%d",&xx,&yy);
        if(xx>yy)swap(xx,yy);
        vec[xx].push_back(yy);
    }
    for(int i=,j,k;i<=n;i++){
        sort(vec[i].begin(),vec[i].end());
        for(j=n,k=vec[i].size()-;j>i&&~k;j--,k--)
            if(vec[i][k]!=j)break;
        if(i==&&j==n){printf("%.10Lf\n",dis(point[],point[n]));return ;}
        if(j>i)addline(line[++tot],point[j],point[i]);
    }addline(line[++tot],point[],point[n]);
    bpmj();
    printf("%.10Lf\n",Ans);
}