bzoj1249: SGU277 HERO 动态凸包

动态维护凸包面积。

//Achen
#include<bits/stdc++.h>
#define For(i,a,b) for(int i=(a);i<=(b);i++)
#define Rep(i,a,b) for(int i=(a);i>=(b);i--)
#define Formylove return 0
const int N=200007;
typedef long long LL;
typedef double db;
using namespace std;
int n;
LL ans;

template<typename T> void read(T &x) {
    char ch=getchar(); T f=1; x=0;
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') f=-1,ch=getchar();
    for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0'; x*=f;
}

struct pt {
    LL x,y;
    pt(){}
    pt(LL x,LL y):x(x),y(y){}
    friend bool operator <(const pt&A,const pt&B) {
        return A.x<B.x||(A.x==B.x&&A.y<B.y);
    }
}p[10];
pt operator -(const pt&A,const pt&B) { return pt(A.x-B.x,A.y-B.y); }
pt operator +(const pt&A,const pt&B) { return pt(A.x+B.x,A.y+B.y); }
LL dot(pt A,pt B) { return A.x*B.x+A.y*B.y; }
LL cross(pt A,pt B) { return A.x*B.y-A.y*B.x; }

set<pt>s1,s2;
#define IT set<pt>::iterator 
void delr(set<pt>&s,pt P,IT l,IT r,int f) {
    while(cross(P-(*l),(*r)-(*l))<=0) {
        if(f) ans+=cross((*r)-(*l),P-(*l));
        s.erase(r); r=l; 
        if(l!=s.begin()) l--;
        else break;
    }
}

void dell(set<pt>&s,pt P,IT l,IT r,int f) {
    while(cross((*r)-P,(*l)-P)<=0) {
        if(f) ans+=cross((*l)-P,(*r)-P);
        s.erase(*l); l=r; r++;
        if(r==s.end()) break;
    }
}

void ins(set<pt>&s,pt P,int f) {
    IT l,r,it;
    if(s.size()<2) {
        s.insert(P); return ;
    }
    it=s.lower_bound(P);
    if(it==s.end()) {
        pt st=*s.begin(); --it;
        if(f) ans-=cross((*it)-st,P-st);
        r=it; l=r; --l;
        if(cross(P-(*l),(*r)-(*l))<0) delr(s,P,l,r,f);
    }
    else {
        if(it==s.begin()) {
            pt ed=*s.rbegin();
            if(f) ans-=cross(P-ed,(*it)-ed);
            l=it; r=l; ++r;
            if(cross((*r)-P,(*l)-P)<0) dell(s,P,l,r,f);
        }
        else {
            r=it; l=r; --l;
            if((*r).x==P.x&&(*r).y==P.y) return ;
            if(cross(P-(*l),(*r)-(*l))>=0) return ;
            if(f) ans+=cross((*l)-P,(*r)-P);
            r=l; 
            if(l!=s.begin()) {
                --l;
                delr(s,P,l,r,f);
            }
            l=it; r=l; ++r;
            if(r!=s.end()) dell(s,P,l,r,f);
        }
    }
    s.insert(P);
}

int main() {
    //freopen("1.in","r",stdin);
    //freopen("1.out","w",stdout);
    For(i,1,3) read(p[i].x),read(p[i].y);
    ans=abs(cross(p[2]-p[1],p[3]-p[1]));
    For(i,1,3) {
        ins(s1,p[i],0);
        ins(s2,pt(-p[i].x,-p[i].y),0);
    }
    read(n);
    For(i,1,n) {
        pt P;
        read(P.x); read(P.y);
        ins(s1,P,1);
        ins(s2,pt(-P.x,-P.y),1);
        printf("%lld
",ans);
    }
    Formylove;
}