Problem:
Problem D: Graph Construction

n 匹のうさぎがいて, 0 番からn − 1 番の番号がついた小屋に1 匹ずつ住んでいる.

あるとき, 秘密の組織によって地下通路の建設工事が行われる, という情報がうさぎたちのもとへ入った. 地下通路を使えば, うさぎたちは他のうさぎたちの小屋へ遊びに行けるようになって嬉しい.

通路は両方向に進むことができ, また通路同士は交わらない. 諸事情により, 一度建設された通路が破壊されてしまうこともある. 通路の建設や破壊の工事は1 本ずつ行われ, 各工事中はうさぎは自分の小屋に留まっているものとする.

うさぎたちは, 工事のいくつかの段階において, あるうさぎとあるうさぎが遊べるかどうかを事前に知りたい. うさぎたちは仲良しなので, 遊びに行くときに他のうさぎの小屋を自由に通ることができる. プログラミングが好きなうさぎたちは, 昔似たような問題を解いたので簡単だろうと思って挑戦し出したが, なかなか効率の良いプログラムが書けない. うさぎの代わりにこの問題を解け.

Input

入力の一行目にはn とk がスペース区切りで与えられる。2 ≤ n ≤ 40 000, 1 ≤ k ≤ 40 000

つづくk 行には工事情報と質問が合わせて, 時間順に与えられる.

“1 u v” — 小屋u とv を結ぶ通路が建設される. 小屋u とv を結ぶ通路がないときにのみ現れる.

“2 u v” — 小屋u とv を結ぶ通路が破壊される. 小屋u とv を結ぶ通路があるときにのみ現れる.

“3 u v” — 小屋u とv のうさぎが遊べるかどうか判定せよ.

0 ≤ u < v < n

Output

入力に現れる各質問について, 遊べるなら”YES”を, そうでないなら”NO”を1 行ずつ出力せよ.

Sample Input 1

4 10
1 0 1
1 0 2
3 1 2
2 0 1
1 2 3
3 0 1
1 0 1
2 0 2
1 1 3
3 0 2

Sample Output 1

YES
NO
YES

Solution in C++11 (judge verdict: Accepted)

#include <bits/stdc++.h>
#define ll long long
#define INF 1000000005
#define MOD 1000000007
#define EPS 1e-10
#define rep(i,n) for(int i=0;i<(int)(n);++i)
#define rrep(i,n) for(int i=(int)(n)-1;i>=0;--i)
#define srep(i,s,t) for(int i=(int)(s);i<(int)(t);++i)
#define each(a,b) for(auto (a): (b))
#define all(v) (v).begin(),(v).end()
#define len(v) (int)(v).size()
#define zip(v) sort(all(v)),v.erase(unique(all(v)),v.end())
#define cmx(x,y) x=max(x,y)
#define cmn(x,y) x=min(x,y)
#define fi first
#define se second
#define pb push_back
#define show(x) cout<<#x<<" = "<<(x)<<endl
#define spair(p) cout<<#p<<": "<<p.fi<<" "<<p.se<<endl
#define svec(v) cout<<#v<<":";rep(kbrni,v.size())cout<<" "<<v[kbrni];cout<<endl
#define sset(s) cout<<#s<<":";each(kbrni,s)cout<<" "<<kbrni;cout<<endl
#define smap(m) cout<<#m<<":";each(kbrni,m)cout<<" {"<<kbrni.first<<":"<<kbrni.second<<"}";cout<<endl

using namespace std;

typedef pair<int,int> P;
typedef pair<ll,ll> pll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<double> vd;
typedef vector<P> vp;
typedef vector<string> vs;

const int MAX_N = 40005;

unordered_map<int,int> mp;

class UF_pop {
public:
	vector<int> par,nrank; vector<vector<int> > hist;
    int sz;
	UF_pop(int node_size){
        sz = node_size;
        par = nrank = vector<int>(sz,0);
        rep(i,sz){
            par[i] = i;
        }
    }
	void reinit(){
        rep(i,sz){
            nrank[i] = 0; par[i] = i;
        }
    }
	int find(int x){
        return (par[x] == x)?(x):find(par[x]);
    }
	void pop(){
		if(len(hist)) {
			auto a = hist.back();
			hist.pop_back();
			par[a[0]] = a[0]; nrank[a[0]] = a[1];
			par[a[2]] = a[2]; nrank[a[2]] = a[3];
		}
	}
	void unite(int x,int y) {
		x = find(x); y = find(y);
		hist.pb({x,nrank[x],y,nrank[y]});
		if(x == y) return;
		if(nrank[x] < nrank[y]){
            par[x] = y;
        }else{
            nrank[x] += (nrank[x] == nrank[y]);
			par[y] = x;
        }
	}
    bool same(int x,int y){
        return find(x) == find(y);
    }
};

void dfs(int l, int r, UF_pop& uf, vi& vid, vi& qid, const vvi& vec, const vvi& que, vi& res) {
	vi lvid,rvid;
	int num=0;
	int mid = (l+r) / 2;
    int st = lower_bound(all(qid),l) - qid.begin();
    int ed = lower_bound(all(qid),r) - qid.begin();
    if(ed <= st){
        return;
    }
	if(vid.empty()){
        srep(i,st,ed){
            if(uf.same(que[i][0],que[i][1])){
                res[i] = 1;
            }
        }
        return;
    }
	each(i,vid){
		const vi& v = vec[i];
		if(v[2] <= l && r <= v[3]) {
			uf.unite(v[0],v[1]);
			++num;
		}else{
			if(v[2] < mid){
                lvid.pb(i);
            }
			if(v[3] > mid){
                rvid.pb(i);
            }
		}
	}
	dfs(l,mid,uf,lvid,qid,vec,que,res);
	dfs(mid,r,uf,rvid,qid,vec,que,res);
	while(num){
        uf.pop();
        --num;
    }
}

int main(){
    int n,q;
    scanf("%d%d",&n,&q);
    vi vid,qid;
    vvi vec,que;
	rep(i,q){
		int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        if(x == 1){
            mp[y*MAX_N+z] = i;
        }else if(x == 2){
            vec.pb({y,z,mp[y*MAX_N+z],i});
            mp.erase(y*MAX_N+z);
        }else{
            qid.pb(i);
            que.pb({y,z});
        }
	}
	each(it,mp){
        vec.pb({it.fi/MAX_N,it.fi%MAX_N,it.se,q});
    }
	vid.resize(len(vec));
	iota(all(vid),0);
	UF_pop uf(n);
    vi res(len(que),0);
	dfs(0,q,uf,vid,qid,vec,que,res);
    rep(i,len(res)){
        if(res[i]){
            cout << "YES\n";
        }else{
            cout << "NO\n";
        }
    }
    return 0;
}