Description 對Samuel星球的探險已經取得了非常巨大的成就，於是科學家們將目光投向了Samuel星球所在的星系——一個巨大的由千百萬星球構成的Samuel星系。 星際空間站的Samuel II巨型電腦經過長期探測，已經鎖定了Samuel星系中許多星球的空間坐標，並對這些星球從1開始編 ...
Submit: 587 Solved: 259
[Submit][Status][Discuss]
Description
對Samuel星球的探險已經取得了非常巨大的成就，於是科學家們將目光投向了Samuel星球所在的星系——一個巨大的由千百萬星球構成的Samuel星系。 星際空間站的Samuel II巨型電腦經過長期探測，已經鎖定了Samuel星系中許多星球的空間坐標，並對這些星球從1開始編號1、2、3……。 一些先遣飛船已經出發，在星球之間開闢探險航線。 探險航線是雙向的，例如從1號星球到3號星球開闢探險航線，那麼從3號星球到1號星球也可以使用這條航線。 例如下圖所示：
在5個星球之間，有5條探險航線。 A、B兩星球之間，如果某條航線不存在，就無法從A星球抵達B星球，我們則稱這條航線為關鍵航線。 顯然上圖中，1號與5號星球之間的關鍵航線有1條：即為4-5航線。 然而，在宇宙中一些未知的磁暴和行星的衝撞，使得已有的某些航線被破壞，隨著越來越多的航線被破壞，探險飛船又不能及時回覆這些航線，可見兩個星球之間的關鍵航線會越來越多。 假設在上圖中，航線4-2（從4號星球到2號星球）被破壞。此時，1號與5號星球之間的關鍵航線就有3條：1-3，3-4，4-5。 小聯的任務是，不斷關註航線被破壞的情況，並隨時給出兩個星球之間的關鍵航線數目。現在請你幫助完成。
Input
第一行有兩個整數N，M。表示有N個星球（1< N < 30000），初始時已經有M條航線（1 < M < 100000）。隨後有M行，每行有兩個不相同的整數A、B表示在星球A與B之間存在一條航線。接下來每行有三個整數C、A、B。C為1表示詢問當前星球A和星球B之間有多少條關鍵航線；C為0表示在星球A和星球B之間的航線被破壞，當後面再遇到C為1的情況時，表示詢問航線被破壞後，關鍵路徑的情況，且航線破壞後不可恢復； C為-1表示輸入文件結束，這時該行沒有A,B的值。被破壞的航線數目與詢問的次數總和不超過40000。Output
對每個C為1的詢問，輸出一行一個整數表示關鍵航線數目。 註意：我們保證無論航線如何被破壞，任意時刻任意兩個星球都能夠相互到達。在整個數據中，任意兩個星球之間最多只可能存在一條直接的航線。Sample Input
5 51 2
1 3
3 4
4 5
4 2
1 1 5
0 4 2
1 5 1
-1
Sample Output
13
HINT
Source
時間倒流
首先把所有的操作全都讀進來
然後把最終形態的樹建出來，
這樣刪邊操作就變成了加邊操作
維護聯通性用LCT，縮點的話用並查集維護祖先，然後暴力改祖先即可
access操作稍微有一些改動，具體看代碼吧
// luogu-judger-enable-o2 #include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int MAXN = 1e5+10; inline int read() { char c = getchar(); int x = 0, f = 1; while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();} while(c >= '0' && c <= '9') {x = x * 10 + c - '0'; c = getchar();} return x * f; } int N, M; int father[MAXN]; int find(int x) { if(father[x] == x) return father[x]; else return father[x] = find(father[x]); } int unionn(int x, int y) { int fx = find(x), fy = find(y); father[fx] = fy; } #define ls(x) T[x].ch[0] #define rs(x) T[x].ch[1] #define fa(x) T[x].f struct node { int f, ch[2], r, siz; }T[MAXN]; int ident(int x) { return T[fa(x)].ch[0] == x ? 0 : 1; } void update(int x) { T[x].siz = T[ls(x)].siz + T[rs(x)].siz + 1; } void pushdown(int x) { if(T[x].r) { T[ls(x)].r ^= 1, T[rs(x)].r ^= 1; swap(ls(x), rs(x)); T[x].r = 0; } } void connect(int x, int fa, int how) { T[x].f = fa; T[fa].ch[how] = x; } bool IsRoot(int x) { return T[fa(x)].ch[0] != x && T[fa(x)].ch[1] != x; } void rotate(int x) { int Y = fa(x), R = fa(Y), Yson = ident(x), Rson = ident(Y); int B = T[x].ch[Yson ^ 1]; T[x].f = R; if(!IsRoot(Y)) connect(x, R, Rson); connect(B, Y, Yson); connect(Y, x, Yson ^ 1); update(Y);update(x); } int st[MAXN]; void splay(int x) { int y = x, top = 0; st[++top] = y; while(!IsRoot(y)) st[++top] = y = fa(y); while(top) pushdown(st[top--]); for(int y = fa(x); !IsRoot(x); rotate(x), y = fa(x)) if(!IsRoot(y)) rotate( ident(x) == ident(y) ? y : x); } void access(int x) { for(int y = 0; x;y = x, x = fa(y) = find(fa(x)))//這裡要改一下 splay(x), rs(x) = y, update(x); } void makeroot(int x) { access(x); splay(x); T[x].r ^= 1; } void link(int x, int y) { //makeroot(x); T[x].f = y; } int findroot(int x) { //makeroot(x); access(x);splay(x); pushdown(x); while(ls(x)) pushdown(x = ls(x)); splay(x); return x; } void delet(int x, int to) { if(x) father[x] = to, delet(ls(x), to), delet(rs(x), to); } void merge(int x, int y) { if(x == y) return ; makeroot(x); if(findroot(y) != x) { link(x, y); return ; } delet(rs(x), x); rs(x) = 0; update(x); } int split(int x,int y) { makeroot(x); access(y); splay(y); } struct Edge { int u,v; bool operator < (const Edge &rhs) const{ return u < rhs.u || (u == rhs.u && v < rhs.v); } }E[MAXN]; int vis[MAXN]; struct Query { int opt, x, y; }Q[MAXN]; int Qnum = 0; int ans[MAXN]; int main() { #ifdef WIN32 freopen("a.in", "r", stdin); #else //freopen("lane.in","r",stdin); //freopen("lane.out","w",stdout); #endif N = read(), M = read(); for(int i = 1; i <= N; i++) father[i] = i; for(int i = 1; i <= M; i++) { E[i].u = read(), E[i].v = read(); if(E[i].u > E[i].v) swap(E[i].u, E[i].v); } sort(E + 1, E + M + 1); while(scanf("%d", &Q[++Qnum].opt) && Q[Qnum].opt != -1) { Q[Qnum].x = read(), Q[Qnum].y = read(); if(Q[Qnum].x > Q[Qnum].y) swap(Q[Qnum].x, Q[Qnum].y); if(Q[Qnum].opt == 0) vis[lower_bound(E + 1, E + M + 1, (Edge){Q[Qnum].x, Q[Qnum].y} ) - E] = 1; } Qnum--; for(int i = 1; i <= M; i++) if(!vis[i]) merge( find(E[i].u), find(E[i].v) ); int ansnum = 0; for(int i = Qnum; i >= 1; i--) { Q[i].x = find(Q[i].x); Q[i].y = find(Q[i].y); if(Q[i].opt == 0) merge(Q[i].x, Q[i].y); else split(Q[i].x, Q[i].y), ans[++ansnum] = T[Q[i].y].siz - 1; } while(ansnum) printf("%d\n", ans[ansnum--]); return 0; }
