题意:
维护一个凸包,支持动态加点。
题解:
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));
}
}