洛谷 P2147 [SDOI2008]洞穴勘测

以下这个做法应该是叫线段树分治。。。

根据修改操作预处理出每条边存在的时间区间[l,r]（以操作序号为时间），然后把所有形式化后的修改挂到线段树节点上。

处理完修改后，dfs一遍线段树，进入某个节点时把那个点上所有的修改操作做一遍连边（用按秩合并并查集），出来时再撤销那些连边；那么到达叶节点时，刚好就是完成了这个节点代表的时间所需要的一切修改操作

复杂度O(nlog^2n)，比lct要大

这个东西跟cdq分治一样是时间分治，但是好像不能降维的样子(?反正想了很久没想通)，不过能在有些有加入、删除、查询操作的题中，把删除操作规避掉

以下代码把线段树避免掉了（其实是处理过程中把树建出来。。。）

 #include<cstdio>
 #include<algorithm>
 #include<map>
 #include<vector>
 #define pb push_back
 using namespace std;
 struct Info
 {
     int fx,fy,dx,dy;
 }ttt[];
 int tlen;
 struct Q
 {
     bool type;
     int l,r,a,b;
 };
 bool ans[];
 int qnum,n,m;
 char tmp[];
 map<int,int> ma[];
 map<int,int>::iterator it;
 int fa[],dep[];
 int find(int x)
 {
     for(;x!=fa[x];x=fa[x]);
     return x;
 }
 void add(const Q &x)
 {
     int fx=find(x.a),fy=find(x.b);++tlen;
     ttt[tlen].fx=fx;ttt[tlen].fy=fy;
     ttt[tlen].dx=dep[fx];ttt[tlen].dy=dep[fy];
     if(dep[fx]>dep[fy])    fa[fy]=fx;
     else
     {
         fa[fx]=fy;
         if(dep[fx]==dep[fy])    dep[fy]++;
     }
 }
 void del(int tm)
 {
     for(int i=,fx,fy;i<=tm;i++)
     {
         fx=ttt[tlen].fx;fy=ttt[tlen].fy;
         fa[fx]=fx;fa[fy]=fy;
         dep[fx]=ttt[tlen].dx;dep[fy]=ttt[tlen].dy;
         --tlen;
     }
 }
 void solve(const vector<Q> &q,int l,int r)
 //当前线段树上区间为[l,r]，要处理的询问为q
 {
     int i,nn=,sz=q.size();
     if(l==r)
     {
         for(i=;i<sz;i++)
             if(q[i].type==)
                 add(q[i]),++nn;
         for(i=;i<sz;i++)
             if(q[i].type==)
                 ans[l]=find(q[i].a)==find(q[i].b);
         del(nn);
         return;
     }
     int mid=l+(r-l)/;vector<Q> qq1,qq2;
     for(i=;i<sz;i++)
     {
         if(q[i].type==)
         {
             if(q[i].l<=l&&r<=q[i].r)    add(q[i]),++nn;
             else
             {
                 if(q[i].l<=mid)    qq1.push_back(q[i]);
                 if(mid<q[i].r)    qq2.push_back(q[i]);
             }
         }
         else
         {
             if(q[i].l<=mid)    qq1.pb(q[i]);
             else    qq2.pb(q[i]);
         }
     }
     solve(qq1,l,mid);
     solve(qq2,mid+,r);
     del(nn);
 }
 vector<Q> q;
 bool type[];
 int main()
 {
     int i,a,b;
     scanf("%d%d",&n,&m);
     for(i=;i<=n;i++)    fa[i]=i;
     for(i=;i<=m;i++)
     {
         scanf("%s%d%d",tmp,&a,&b);if(a>b)    swap(a,b);
         if(tmp[]=='Q')    q.pb((Q){,i,,a,b}),type[i]=;
         else if(tmp[]=='C')    ma[a][b]=i;
         else if(tmp[]=='D')    q.pb((Q){,ma[a][b],i,a,b}),ma[a].erase(b);
     }
     for(i=;i<=n;i++)
         for(it=ma[i].begin();it!=ma[i].end();it++)
             q.pb((Q){,it->second,m,i,it->first});
     solve(q,,m);
     for(i=;i<=m;i++)
         if(type[i]==)
             puts(ans[i]?"Yes":"No");
     return ;
 }