3069: [Pa2011]Hard Choice 艰难的选择

Description

Byteasar是一个很纠结的人。每次他经过Bytetown的时候都知道有至少2条不同的路径可以选择，这导致他必须花很长时间来决定走哪条路。Byteasar最近听说了Bytetown的修路计划，他可能是唯一一个为此感到高兴的人——他有机会消除他的烦恼。

在Byteasar一共有n个岔口，连接着m条双向道路。两条路径完全不同当且仅当他们没有公共的道路（但是允许经过相同的岔口）。

Byteasar想知道：对于两个岔口x y，是否存在一对完全不同的路径。

Input

第一行3个整数：n, m, z (2<=n<=100000, 1<=m,z<=100000)，分别代表：n个岔口，m条边，事件数z。岔口编号为1~n。

下面m行：ai, bi (1<=ai,bi<= n, ai!=bi)，描述一条边

然后下面z行描述事件：ti, ci, di (t='Z' or 'P', 1<=ci,di<=n, ci!=di)。事件按照时间排序。

• 当t='Z'，表示删除一条边(ci, di)，保证这条边之前没有被删除。注意，边可以被全部删除！
• 当t='P'，询问是否存在从ci到di的一对完全不同的路径。

Output

对于每组询问，如果存在，输出TAK，否则输出NIE。

逆着操作顺序加边，同时并查集维护连通性和边双连通分量
预处理出按加边顺序得到的生成森林，在上面把每个连通块定向为有根树并标上深度
正式加边时，若两侧不连通则标为联通，否则将边的两端在生成树上的路径并成同一个边双连通分量
查询时可以直接判断两点是否在同个边双连通分量内

#include<bits/stdc++.h>
char buf[],*ptr=buf-;
int _(){
    int x=,c=*++ptr;
    while(c<)c=*++ptr;
    while(c>)x=x*+c-,c=*++ptr;
    return x;
}
int _o(){
    int c=*++ptr;
    while(c<'A'||c>'Z')c=*++ptr;
    return c;
}
const int N=;
int es[N*],enx[N*],ev[N*],e0[N],ep=,del[N];
int n,m,q,qs[N][],f1[N],f2[N],fa[N],ap=,dep[N];
bool ed[N],as[N];
int get(int*f,int x){
    int a=x,c;
    while(x!=f[x])x=f[x];
    while(x!=f[a])c=f[a],f[a]=x,a=c;
    return x;
}
struct edge{
    int a,b,id;
    bool operator<(edge e)const{return a!=e.a?a<e.a:b<e.b;}
    void chk(){
        int x=get(f1,a),y=get(f1,b);
        if(x!=y){
            f1[x]=y;
            ev[id<<]=ev[id<<|]=;
        }
    }
    void ins(){
        int x=get(f1,a),y=get(f1,b);
        if(x!=y){
            f1[x]=y;
            return;
        }
        a=get(f2,a);b=get(f2,b);
        while(a!=b){
            if(dep[a]<dep[b])b=f2[b]=get(f2,fa[b]);
            else a=f2[a]=get(f2,fa[a]);
        }
    }
}e[N];
void dfs(int w){
    ed[w]=;
    for(int i=e0[w];i;i=enx[i]){
        int u=es[i];
        if(!ed[u]&&ev[i])fa[u]=w,dep[u]=dep[w]+,dfs(u);
    }
}
void ae(){
    int a=_(),b=_();
    if(a>b)std::swap(a,b);
    e[ep>>]=(edge){a,b,ep>>};
    es[ep]=b;enx[ep]=e0[a];e0[a]=ep++;
    es[ep]=a;enx[ep]=e0[b];e0[b]=ep++;
}
void query(int a,int b){
    as[ap++]=(get(f2,a)==get(f2,b));
}
int main(){
    fread(buf,,sizeof(buf),stdin);
    n=_();m=_();q=_();
    for(int i=;i<=m;++i)ae();
    std::sort(e+,e+m+);
    for(int i=,a,b;i<=q;++i){
        if(qs[i][]=(_o()=='Z')){
            a=_();b=_();
            if(a>b)std::swap(a,b);
            a=std::lower_bound(e+,e+m+,(edge){a,b})-e;
            ++del[a];
            qs[i][]=a;
        }else{
            qs[i][]=_();
            qs[i][]=_();
        }
    }
    for(int i=;i<=n;++i)f1[i]=i;
    for(int i=;i<=m;++i){
        while(del[i])i+=del[i];
        if(i<=m)e[i].chk();
    }
    for(int i=q;i;--i)if(qs[i][])e[qs[i][]].chk();
    for(int i=;i<=n;++i)if(!ed[i])dfs(i);
    for(int i=;i<=n;++i)f1[i]=f2[i]=i;
    for(int i=;i<=m;++i){
        while(del[i])i+=del[i];
        if(i<=m)e[i].ins();
    }
    for(int i=q;i;--i)if(qs[i][])e[qs[i][]].ins();else query(qs[i][],qs[i][]);
    while(ap)puts(as[--ap]?"TAK":"NIE");
    return ;
}