【题解】P2521 [HAOI2011]防线修建(动态凸包)

【题解】P2521 [HAOI2011]防线修建(动态凸包)

凸包是易插入不好删除的东西，按照剧情所以我们时光倒流

然后问题就是维护凸包的周长，支持加入

本来很简单，但是计算几何就是一些小地方经验不足容易WA和RE

然后代码注释里有一些经验

//@winlere
#include
#include
#include
#include
#include
#include
#include
#define pf(x) ((x)*(x))
#define DEBUG(s,a) cerr<<#s" = "<<(s)<<" \n"[(a)==1]
#define getchar() (__c==__ed?(__ed=__buf+fread(__c=__buf,1,1<<18,stdin),*__c++):*__c++)

using namespace std;  typedef long long ll;   char __buf[1<<18],*__c=__buf,*__ed=__buf;
inline int qr(){
    int ret=0,f=0,c=getchar();
    while(!isdigit(c))f|=c==45,c=getchar();
    while(isdigit(c)) ret=ret*10+c-48,c=getchar();
    return f?-ret:ret;
}

const int maxn=1e5+5;
struct NODE{//向量 坐标用一个结构体代替就行了
    ll x,y;//用ll吧 别小气
    NODE operator - (NODE b)const{return {x-b.x,y-b.y};}
    ll operator % (NODE b)const{return x*b.y-y*b.x;}
    bool operator < (NODE b)const{return x > ve;
double len,ans[maxn];
bool in[maxn];
set s;

void insert(NODE x){
    vector ve;
    if((*s.upper_bound(x)-x)%(x-*--s.upper_bound(x))<0) return;//先判断是否能够插入，这样以后好写一些
    for(set::iterator t=--s.upper_bound(x),temp;t!=s.begin();t=temp){
        temp=prev(t);
        if((x-*temp)%(*t-*temp)<0) ve.push_back(*t);//避免一边用iterator一边erase
        else break;
    }
    reverse(ve.begin(),ve.end());
    for(set::iterator t=s.upper_bound(x),temp;t!=--s.end();t=temp){
        temp=next(t);
        if((*temp-x)%(*temp-*t)>0) ve.push_back(*t);//同上
        else break;
    }
    if(ve.size()==0)//多写几个if比仔细比巧妙统一快得多，对于=0,=1都特判一下，这里只要=0特判
        len+=(*s.upper_bound(x)^x)+(*--s.upper_bound(x)^x)-(*--s.upper_bound(x)^*s.upper_bound(x));
    else {
        auto l=*--s.find(ve.front()),r=*++s.find(ve.back());
        len+=(l^x)+(r^x)-(l^ve.front())-(r^ve.back());
        for(int t=1,ed=ve.size();t){op,qr()}:(pair){op,t});
        if(op==1) in[ve.back().second]=1;
    }
    for(int t=1;t<=m;++t)
        if(!in[t])
            insert(data[t]);
    reverse(ve.begin(),ve.end());
    for(auto t:ve)
        if(t.first==1) insert(data[t.second]);
        else ans[t.second]=len;
    reverse(ve.begin(),ve.end());
    for(auto t:ve)
        if(t.first==2) printf("%.2lf\n",ans[t.second]);
    return 0;
}