P1197 [JSOI2008]星球大战(并查集判断连通块+正难则反)
并查集本来就是连一对不同父亲的节点就的话连通块就少一个。
题目描述
很久以前,在一个遥远的星系,一个黑暗的帝国靠着它的超级武器统治者整个星系。某一天,凭着一个偶然的机遇,一支反抗军摧毁了帝国的超级武器,并攻下了星系中几乎所有的星球。这些星球通过特殊的以太隧道互相直接或间接地连接。
但好景不长,很快帝国又重新造出了他的超级武器。凭借这超级武器的力量,帝国开始有计划地摧毁反抗军占领的星球。由于星球的不断被摧毁,两个星球之间的通讯通道也开始不可靠起来。现在,反抗军首领交给你一个任务:给出原来两个星球之间的以太隧道连通情况以及帝国打击的星球顺序,以尽量快的速度求出每一次打击之后反抗军占据的星球的连通快的个数。(如果两个星球可以通过现存的以太通道直接或间接地连通,则这两个星球在同一个连通块中)。
输入输出格式
输入格式:
输入文件第一行包含两个整数,N (1 <= N <= 2M) 和M (1 <= M <= 200,000),分别表示星球的数目和以太隧道的数目。星球用0~N-1的整数编号。
接下来的M行,每行包括两个整数X, Y,其中(0<=X<>Y<N),表示星球X和星球Y之间有以太隧道。注意所有的以太隧道都是双向的。
接下来一行是一个整数K,表示帝国计划打击的星球个数。
接下来的K行每行一个整数X,满足0<=X<N,表示帝国计划打击的星球编号。帝国总是按输入的顺序依次摧毁星球的。
输出格式:
输出文件的第一行是开始时星球的连通块个数。
接下来的K行,每行一个整数,表示经过该次打击后现存星球的连通块个数。
输入输出样例
8 13 0 1 1 6 6 5 5 0 0 6 1 2 2 3 3 4 4 5 7 1 7 2 7 6 3 6 5 1 6 3 5 7
1 1 1 2 3 3
说明
[JSOI2008]
解答
参考洛谷题解
读完题目,我们首先会想到每次删去一个点,然后重新建图,算出联通块的个数。然而,这根本不信,时间上过不去,这时就得让我们将头旋转180度,从后往前处理,也就是我们经常用的接边了。当然,再接边前需要预处理,当前这条边应当在什么时候接上,给边标上序号,然后依次处理。最后逆序输出,就可以啦啦啦。如果还不懂,请看代码。
PS:帝国攻占的地盘不算是反抗军的地盘。
比较有趣的一道并查集。
我们在做组合数学或者集合的时候会说一句话,叫做正难则反。这里也是同样的。
首先,既然是对联通块的操作,那自然而然的想到用并查集解决。
先来思考一个问题:如果有n个块,然后我们要进行合并,合并之后有多少块?只需要对合并函数稍稍修改:
inline void uni(int x, int y) {
int xx = find(x), yy = find(y);
if(xx == yy) return;
--cnt;
fa[xx] = yy;
}
每一次合并,如果两个在同一集合中就不做处理,否则相当于连通块-1.
返回这道题,摧毁的过程是相当难以维护的,由于这道题是离线的,所以我们可以反过来重建。以样例为例,在摧毁的最后,节点0 2 4分别构成了3个连通块,注意这里答案并不是8,因为不存在剩下的5个节点!
然后每每恢复一个节点,就可以把和它相关联的连通块都合并一遍,但是要注意,这个时候的合并的前提,是对面重建过。实现见下。
1 #include<cstdio> 2 #include<cstring> 3 #include<iostream> 4 #include<algorithm> 5 #define MAXN 200005 6 using namespace std; 7 int n,m,fa[MAXN*2],K,D[MAXN*2],vis[MAXN*2],ans[2*MAXN],num; 8 //num表示总联通块的个数,ans[]记住答案 ,vis[]给点标号 ,D[]被攻占的点 9 struct xcw{ 10 int x,y,c;//x,y表示这条边的两端点,c表示这条边的序号 11 bool operator <(const xcw b)const{return c<b.c;}//按序号从小到大排序 12 }a[MAXN]; 13 int read(){ 14 int ret=0,f=1;char ch=getchar(); 15 while(ch<'0'||ch>'9'){if(ch=='-') f=-f;ch=getchar();} 16 while(ch>='0'&&ch<='9') ret=ret*10+ch-48,ch=getchar(); 17 return ret*f; 18 } 19 int get(int x){return fa[x]==x?x:fa[x]=get(fa[x]);} 20 void mer(int x,int y){ 21 int fx=get(x),fy=get(y); 22 if(fx!=fy) fa[fx]=fy,num--;//两联通块合并,总联通块num-- 23 } 24 int main(){ 25 // freopen("a.in","r",stdin); 26 // freopen("a.out","w",stdout); 27 n=read();m=read();num=n; 28 for(int i=0;i<n;i++) fa[i]=i; 29 for(int i=1;i<=m;i++) a[i]=(xcw){read(),read(),0}; 30 K=read(); 31 for(int i=1;i<=K;i++) vis[D[i]=read()]=K-i+1;//给点标号,从大到小 32 for(int i=1;i<=m;i++) a[i].c=max(vis[a[i].x],vis[a[i].y]);//当前边的序号是两个点中最大的那个 33 sort(a+1,a+1+m);//给点按序号大小排序,方便处理,优化效率 34 for(int i=0,j=1;i<=K;i++){//枚举被攻占的星球序号 35 for(;a[j].c==i;j++) mer(a[j].x,a[j].y);//将当前序号全部合并 36 ans[i]=num-(K-i);//K-i是帝国攻占的星球个数(因为是逆序处理的) 37 } 38 for(int i=K;i>=0;i--) printf("%d ",ans[i]);//逆序输出 39 return 0; 40 }
代码具体过程:
D[]数组
i D[i] vis[D[i]]
1 1 5
2 6 4
3 3 3
4 5 2
5 7 1
vis[]数组:
0 1 2 3 4 5 6 7
0 5 0 3 0 2 4 1
排序前:
i a[i].x a[i].y a[i].c
1 0 1 5
2 1 6 5
3 6 5 4
4 5 0 2
5 0 6 4
6 1 2 5
7 2 3 3
8 3 4 3
9 4 5 2
10 7 1 5
11 7 2 1
12 7 6 4
13 3 6 4
排序后:
i a[i].x a[i].y a[i].c
1 7 2 1
2 5 0 2
3 4 5 2
4 2 3 3
5 3 4 3
6 6 5 4
7 0 6 4
8 7 6 4
9 3 6 4
10 0 1 5
11 1 6 5
12 1 2 5
13 7 1 5
结果:
i num ans[i]
0 8 3
1 7 3
2 5 2
3 3 1
4 2 1
5 1 1
1
1
1
2
3
3