黑魔法师之门
(magician.pas/c/cpp)
题目描述
经过了16个工作日的紧张忙碌,未来的人类终于收集到了足够的能源。然而在与Violet星球的战争中,由于Z副官的愚蠢,地球的领袖applepi被邪恶的黑魔法师Vani囚禁在了Violet星球。为了重启Nescafé这一宏伟的科技工程,人类派出了一支由XLk、Poet_shy和lydrainbowcat三人组成的精英队伍,穿越时空隧道,去往Violet星球拯救领袖applepi。
applepi被囚禁的地点只有一扇门,当地人称它为“黑魔法师之门”。这扇门上画着一张无向无权图,而打开这扇门的密码就是图中每个点的度数大于零且都是偶数的子图的个数对1000000009取模的值。此处子图 (V, E) 定义为:点集V和边集E都是原图的任意子集,其中E中的边的端点都在V中。
但是Vani认为这样的密码过于简单,因此门上的图是动态的。起初图中只有N个顶点而没有边。Vani建造的门控系统共操作M次,每次往图中添加一条边。你必须在每次操作后都填写正确的密码,才能够打开黑魔法师的牢狱,去拯救伟大的领袖applepi。
输入格式
第一行包含两个整数N和M。
接下来M行,每行两个整数A和B,代表门控系统添加了一条无向边 (A, B)。
输出格式
输出一共M行,表示每次操作后的密码。
样例输入
4 8
3 1
3 2
2 1
2 1
1 3
1 4
2 4
2 3
样例输出
0
0
1
3
7
7
15
31
样例说明
第三次添加之后,存在一个满足条件的子图 {1, 2, 3}(其中1, 2, 3是数据中边的标号)。
第四次添加之后,存在三个子图 {1, 2, 3},{1, 2, 4},{3, 4}。
……
数据范围与约定
对于30% 的数据,N, M≤10。
对于100% 的数据,N≤200000,M≤300000。
题目中定义的子环,经过分析可以发现,它就是一个一个环。
我们可以求出元环个数n来,再根据排列组合求出元环组合(子图)总数为2^n-1.
我考试的时候卡在了后面那句话上,因此就没做出来。
即子图总数为2^n-1,因为我以为这些环必须是有交点的。实际上不用,再仔细看一下定义,可以明白,下面这种情况,其实他的子图个数是3而不是2。
当时我就卡在了这里,我想找到所有有交点的环,再用排列组合来算,但是这实际上很难处理。。。都是题意没有搞清楚的过。
多读题,仔细理解,总会有收获的。
另一个我没有突破的,实际上可能和以上问题有关,限制了自己的想法。其实,我们从它的边是一条一条加入的,能够看出一些端倪。另外貌似noip要求的在线算法不多,常常会用到并查集。而这道题,确实是要用到并查集。
和环有关的算法,了解到并查集的优势后,现在能够使用的就有三个了。Tarjan,floyd,并查集。(Tarjan放在这里其实是比较牵强的)
而并查集,我们只需要在加入新的边之前,判断是否已经在一个联通分量里,如果在,增加一条边,必定会增加一个元环。
题解里的思路和我的不同。顺便贴上题解。
黑魔法师之门
题目来源:XLk 提供
大致思路:
实际上每次操作后的答案就是2^(图中”元”环的个数)。
元环的意思如右图所示,(1-2-3-4-1)和(3-4-5-3) 是元环,
1-2-3-5-4-1 不是,因为它可以看做由上述的两个环合成。
因为一个环里每个点的度数都是大于零的偶数,我们可以这
样来构造答案:每个环有选和不选两种选择,如果选择了该
环,那么环上所有边的“选择次数”+1。最后取所有“选择次数”为奇数的边构成一个边集,
就是一个答案。可以证明这样构造出来的解不重复且涵盖了所有情况。因此答案就是2^(图
中”元”环的个数)。实现方法非常简单,只需要一个并查集即可。
具体实现方法:
并查集维护连通性,初始化ans=1。
加入一条边(x,y)时,如果x 和y 在同一集合内,ans*=2。
每次询问输出ans-1。
时间复杂度O(M α(N)),α(N)代表并查集的复杂度。
#include <cstdio>
#include <iostream>
using std::cout;
long fa[200010];
long getroot(long l)
{
if (fa[l] == l)
return fa[l];
return fa[l] = getroot(fa[l]);
}
inline void merge(long a,long b)
{
fa[getroot(a)] = getroot(b);
}
long getint()
{
long rs=0;bool sgn=1;char tmp;
do tmp=getchar();
while (!isdigit(tmp)&&tmp-'-');
if (tmp=='-'){sgn=0;tmp=getchar();}
do rs=(rs<<3)+(rs<<1)+tmp-'0';
while (isdigit(tmp=getchar()));
return sgn?rs:-rs;
}
int main()
{
freopen("magician.in","r",stdin);
freopen("magician.out","w",stdout);
long n = getint();
long m = getint();
long ans = 1;
for (long i=1;i<n+1;i++)
fa[i] = i;
for (long i=1;i<m+1;i++)
{
long a = getint();
long b = getint();
if (getroot(a) == getroot(b))
{
ans <<= 1;
if (ans > 1000000008)
ans -= 1000000009;
}
else
merge(a,b);
printf("%ld\n",ans-1);
}
return 0;
}