BZOJ 2300 防线修建

http://www.lydsy.com/JudgeOnline/problem.php?id=2300

题意：给点，有以下操作:删去一个点，询问这些点构成凸包的周长。

思路:用splay维护上凸壳，splay排序的关键字是X，同时还要记录每个点与左右点的斜率，当加入一个点时，我们找到它应该在的X的位置，一个个减掉左边，一个个减掉右边，符合条件就是当左边的斜率大于右边的斜率，就表明它在凸包上。

PS:如果当前加入点的左边斜率小于右边斜率，说明它在凸包内，可以直接删掉这个点。

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
struct node{
     int id,val,opt;
}q[];
struct point{
    double x,y;
}p[];
const double eps=1e-;
int tmp,root;
double c[],xx[],yy[],ans,rk[],lk[];
int n,ch[][],fa[],vis[];
int read(){
    int t=,f=;char ch=getchar();
    while (ch<''||ch>''){if (ch=='-') f=-;ch=getchar();}
    while (''<=ch&&ch<=''){t=t*+ch-'';ch=getchar();}
    return t*f;
}
void rotate(int x,int &rt){
     int y=fa[x],z=fa[y],l,r;
     if (ch[y][]==x) l=;else l=;r=l^;
     if (y!=rt){
        if (ch[z][]==y) ch[z][]=x;else ch[z][]=x;
     }else rt=x;
     fa[x]=z;fa[y]=x;fa[ch[x][r]]=y;
     ch[y][l]=ch[x][r];ch[x][r]=y;
}
void splay(int x,int &rt){
     while (x!=rt){
        int y=fa[x],z=fa[y];
        if (y!=rt){
           if (ch[z][]==y^ch[y][]==x) rotate(x,rt);
           else rotate(y,rt);
        }
        rotate(x,rt);
     }
}
void find(int k,double X){
    if (!k) return;
    if (xx[k]<X+eps) tmp=k,find(ch[k][],X);
    else find(ch[k][],X);
}
int suc(int x){
    x=ch[x][];
    while (ch[x][]) x=ch[x][];
    return x;
}
int pre(int x){
    x=ch[x][];
    while (ch[x][]) x=ch[x][];
    return x;
}
double sqr(double x){
    return x*x;
}
double dis(int x,int y){
    if (!x||!y) return 0.0;
    return sqrt(sqr(xx[x]-xx[y])+sqr(yy[x]-yy[y]));
}
void getslope(int x,int y){
    if (!x){lk[y]=1e9;return;}
    if (!y){rk[x]=-1e9;return;}
    if (fabs(xx[x]-xx[y])<1e-){
        if (yy[x]<yy[y]) rk[x]=lk[y]=1e9;
        else rk[x]=lk[y]=-1e9;
        return;
    }
    rk[x]=lk[y]=(yy[x]-yy[y])/(xx[x]-xx[y]);
}
void insert(double X,double Y,int id){
    tmp=;
    find(root,X);int x=tmp,y=;
    if (!x){
       x=root;
       while (ch[x][]) x=ch[x][];
       splay(x,root);
       y=x;x=;
    }else splay(x,root),splay(y=suc(x),ch[x][]);
    xx[id]=X;yy[id]=Y;
    if (y) fa[id]=y,ch[y][]=id;
    else fa[id]=x,ch[x][]=id;
    getslope(x,id);getslope(id,y);
    ans-=dis(x,y);
    if (rk[id]>=lk[id]){getslope(x,y);fa[id]=;ch[y][]=;ans+=dis(x,y);return;}
    splay(id,root);
    root=id;
    x=pre(id);
    ans+=dis(x,id);
    while (lk[x]<=rk[x]&&x){
        y=pre(x);
        splay(y,ch[x][]);
        ans-=dis(x,id);
        ans-=dis(x,y);
        ans+=dis(y,id);
        fa[x]=ch[x][]=ch[x][]=;
        ch[id][]=y;fa[y]=id;
        getslope(y,id);
        x=y;
    }
    x=suc(id);
    ans+=dis(x,id);
    while (lk[x]<=rk[x]&&x){
       y=suc(x);
       splay(y,ch[x][]);
       ans-=dis(x,id);
       ans-=dis(x,y);
       ans+=dis(y,id);
       fa[x]=ch[x][]=ch[x][]=;
       ch[id][]=y;fa[y]=id;
       getslope(id,y);
       x=y;
    }
}
int main(){
    int m=read(),x=read(),y=read();
    n=read();
    for (int i=;i<=n;i++)
      scanf("%lf%lf",&p[i].x,&p[i].y);
    int T=read();
    for (int i=;i<=T;i++){
      q[i].opt=read();
      if (q[i].opt==) continue;
      q[i].id=i;q[i].val=read();vis[q[i].val]=;
    }
    insert(,,);
    insert(m,,);
    insert(x,y,);
    for (int i=;i<=n;i++)
        if (!vis[i]){
              insert(p[i].x,p[i].y,i+);
        }
    int top=;
    for (int i=T;i>=;i--){
        if (q[i].opt==){
            c[++top]=ans;
        }else{
            insert(p[q[i].val].x,p[q[i].val].y,q[i].val+);
        }
    }
    for (int i=top;i>=;i--)
        printf("%.2f\n",c[i]);
}