BZOJ1249：SGU277 HERO 动态凸包（set）

传送门

题意：
维护一个凸包，支持动态加点。

题解：
set维护上下凸包，代码稍微有点麻烦。

注意设置比较函数的时候上下凸包应该一样，才能保证正确的弹出。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
struct IO{
    streambuf *ib,*ob;
    inline void init(){
        ios::sync_with_stdio(false);
        cin.tie(NULL);cout.tie(NULL);
        ib=cin.rdbuf();ob=cout.rdbuf();
    }
    inline ll read(){
        char ch=ib->sbumpc();ll i=0,f=1;
        while(!isdigit(ch)){if(ch=='-')f=-1;ch=ib->sbumpc();}
        while(isdigit(ch)){i=(i<<1)+(i<<3)+ch-'0';ch=ib->sbumpc();}
        return i*f;
    }
    inline void W(ll x){
        static int buf[50];
        if(!x){ob->sputc('0');ob->sputc('\n');return;}
        if(x<0){ob->sputc('-');x=-x;}
        while(x){buf[++buf[0]]=x%10;x/=10;}
        while(buf[0]){ob->sputc(buf[buf[0]--]+'0');}
        ob->sputc('\n');
    }
}io;
ll area;
struct point{
    ll x,y;
    point(ll x=0,ll y=0):x(x),y(y){}
    friend inline point operator -(const point &a,const point &b){return point(a.x-b.x,a.y-b.y);}
    friend inline ll operator *(const point &a,const point &b){return a.x*b.y-a.y*b.x;}
};
struct cmp{
    inline bool operator()(const point &a,const point &b){
        return a.x<b.x||(a.x==b.x&&a.y<b.y);
    }
};
set<point,cmp>S_up;
set<point,cmp>S_down;
typedef set<point>::iterator it;
inline void maintain_down(point p){
    it r=S_down.lower_bound(p);
    if(r==S_down.begin()){
        it nxt=++r;r--;
        while(nxt!=S_down.end()&&(*r-p)*((*nxt-p))<=0)area-=((*nxt)*(*r)),S_down.erase(r),r=(nxt++);
        area+=((*r)*p);S_down.insert(p);
    }else if(r==S_down.end()){
        it nxt=(--r);r--;
        while(nxt!=S_down.begin()&&(*r-p)*(*nxt-p)<=0)area-=((*nxt)*(*r)),S_down.erase(nxt),nxt=(r--);
        area+=(p*(*nxt));S_down.insert(p);
    }else {
        it l=(--r);r++;
        if((*l-p)*(*r-p)>=0)return;
        area-=(*r*(*l));
        while(r!=(--S_down.end())){
            it nxt=(++r);r--;
            if((*r-p)*(*nxt-p)<=0)area-=(*nxt*(*r)),S_down.erase(r),r=(nxt++);
            else break;
        }
        area+=((*r)*p);
        while(l!=S_down.begin()){
            it nxt=--l;l++;
            if((*nxt-p)*(*l-p)<=0)area-=(*l*(*nxt)),S_down.erase(l),l=(nxt--);
            else break;
        }
        area+=(p*(*l));
        S_down.insert(p);
    }
}
inline void maintain_up(point p){
    it r=S_up.lower_bound(p);
    if(r==S_up.begin()){
        it nxt=++r;r--;
        while(nxt!=S_up.end()&&(*r-p)*((*nxt-p))>=0)area-=((*r)*(*nxt)),S_up.erase(r),r=(nxt++);
        area+=(p*(*r));S_up.insert(p);
    }else if(r==S_up.end()){
        it nxt=(--r);r--;
        while(nxt!=S_up.begin()&&(*r-p)*(*nxt-p)>=0)area-=((*r)*(*nxt)),S_up.erase(nxt),nxt=(r--);
        area+=(*nxt*p);S_up.insert(p);
    }else {
        it l=(--r);r++;
        if((*l-p)*(*r-p)<=0)return;
        area-=(*l*(*r));
        while(r!=(--S_up.end())){
            it nxt=(++r);r--;
            if((*r-p)*(*nxt-p)>=0)area-=(*r*(*nxt)),S_up.erase(r),r=(nxt++);
            else break;
        }
        area+=(p*(*r));
        while(l!=S_up.begin()){
            it nxt=--l;l++;
            if((*nxt-p)*(*l-p)>=0)area-=(*nxt*(*l)),S_up.erase(l),l=(nxt--);
            else break;
        }
        area+=(*l*p);
        S_up.insert(p);
    }
}
int main(){
    io.init();point p1,p2,p3;
    p1.x=io.read(),p1.y=io.read(),p2.x=io.read(),p2.y=io.read(),p3.x=io.read(),p3.y=io.read();
    if(p1.x>p2.x||(p1.x==p2.x&&p1.y>p2.y))swap(p1,p2);
    if(p1.x>p3.x||(p1.x==p3.x&&p1.y>p3.y))swap(p1,p3);
    if(p3.x<p2.x||(p3.x==p2.x&&p3.y<p2.y))swap(p2,p3);
    S_up.insert(p1);S_up.insert(p3);S_down.insert(p1);S_down.insert(p3);
    ((p1-p2)*(p3-p2)>0)?(S_up.insert(p2),area=p1*p2+p2*p3+p3*p1):(S_down.insert(p2),area=p1*p3+p3*p2+p2*p1);
    int q=io.read();
    while(q--){
        p1.x=io.read(),p1.y=io.read();
        maintain_up(p1);maintain_down(p1);
        io.W(abs(area));
    }
}