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=;
 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=; x=;
     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
     if(ch=='-') f=-,ch=getchar();
     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; 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[];
 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))<=) {
         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)<=) {
         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()<) {
         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))<) 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)<) 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))>=) 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,,) read(p[i].x),read(p[i].y);
     ans=abs(cross(p[]-p[],p[]-p[]));
     For(i,,) {
         ins(s1,p[i],);
         ins(s2,pt(-p[i].x,-p[i].y),);
     }
     read(n);
     For(i,,n) {
         pt P;
         read(P.x); read(P.y);
         ins(s1,P,);
         ins(s2,pt(-P.x,-P.y),);
         printf("%lld\n",ans);
     }
     Formylove;
 }