2019.02.28【HAOI2018】【BZOJ5303】【洛谷P4494】反色游戏（割点） （高斯消元推结论）

BZOJ传送门

洛谷传送门

---

解析：

出题人的部分分给的太多了吧。。。光是$O(n^4)$高斯消元解异或方程组就有$50pts$。

不过还是可以理解的，毕竟谁会想到考高斯消元啊。

一个很显然的做法，高斯消元解异或方程组，设最后自由元的数量为$t$，显然答案就是$2^t$。

注意一点，如果最后某一个联通分量中有奇数个黑点就$gg$了。显然我们每次操作对黑点数量的改变都是偶数。

随便推几组样例会发现，答案就是$2^{m-n+c}$，其中$m,n$分别表示边数，点数，$C$表示联通分量的数量。

这个可以很轻松的发现。

考虑新加一条边对方程组的影响。

如果当前边连接的是两个不同的连通分量，那么这条边在方程组解出来就肯定不是自由元了。对应的$m+=1,C-=1$，答案没有变化。

不然的话，这条边相当于在原来联通分量中突然多出，解出来方程它会剩下，因为最后方程中没有该列系数为$1$的行了。它就是自由元。

现在已经可以用并查集在$O(n^2\alpha (n))$时间内解决这个问题拿到70分了。

我们考虑怎样能有更加优秀的复杂度？

计算出每个点删除后对联通分量和奇数黑点联通分量的影响！

边和点的个数的影响可以非常方便的计算。

连通分量数量的改变只需要看一看是不是割点或独立点就行了。

而我们同样可以在$tarjan$求割点的时候计算出每个点的$dfs$子树内部点权的异或和，可以知道$dfs$子树内部会增加多少个黑联通块，然后用树根的异或和去异或就可以得到子树外部。

实现和细节都相当简单。

---

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc get_char
#define gets get_s
#define cs const

namespace IO{
	inline char get_char(){
		static cs int Rlen=1<<20|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	
	inline void get_s(char *s){
		re char c;
		while(!isspace(c=gc()))*++s=c;
	}
	
	inline int getint(){
		re char c;
		while(!isdigit(c=gc()));re int num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return num;
	}
}
using namespace IO;

cs int mod=1e9+7;
inline int add(int a,int b){return a+b>=mod?a+b-mod:a+b;}
inline int dec(int a,int b){return a<b?a-b+mod:a-b;}
inline int mul(int a,int b){return (ll)a*b%mod;}

cs int N=1e5+5;
int n,m;
int pow2[N<<1];
int last[N],nxt[N<<1],to[N<<1],ecnt;
int du[N];
inline void addedge(int u,int v){
	++du[u],++du[v];
	nxt[++ecnt]=last[u],last[u]=ecnt,to[ecnt]=v;
	nxt[++ecnt]=last[v],last[v]=ecnt,to[ecnt]=u;
}

int black_bcc;
bool cut[N];
int around[N],black[N];
int dfn[N],low[N],dfs_clock;
int top,sta[N];
bool tag[N],sub[N];

void tarjan(int u,int fa){
	sta[++top]=u;
	dfn[u]=low[u]=++dfs_clock;
	int son=0;
	sub[u]=tag[u];
	for(int re e=last[u],v=to[e];e;v=to[e=nxt[e]])
	if(v^fa){
		if(!dfn[v]){
			tarjan(v,u);++son;
			low[u]=min(low[u],low[v]);
			tag[u]^=tag[v];
			if(low[v]>=dfn[u]){
				sub[u]^=tag[v];
				cut[u]=true;
				++around[u];
				if(tag[v])++black[u];
			}
		}
		else low[u]=min(low[u],dfn[v]);
	}
	if(son==1&&fa==0)cut[u]=around[u]=black[u]=0;
	else if(fa==0)around[u]--;
}

inline void init(){
	ecnt=0;black_bcc=0;
	memset(last,0,sizeof last);
	memset(dfn,0,sizeof dfn);
	memset(low,0,sizeof low);
	memset(cut,0,sizeof cut);
	memset(around,0,sizeof around);
	memset(black,0,sizeof black);
	memset(du,0,sizeof du);
}

int ans;
char s[N];
inline void solve(){
	n=getint();
	m=getint();
	for(int re i=1;i<=m;++i)addedge(getint(),getint());
	gets(s);
	for(int re i=1;i<=n;++i)tag[i]=s[i]^48;
	ans=m-n;
	for(int re i=1;i<=n;++i)if(!dfn[i]){
		++ans;
		tarjan(i,0);
		if(tag[i])++black_bcc;
		if(!du[i])around[i]=-1;
		for(;top;top--)if(cut[sta[top]]){
			int v=sta[top];
			if(tag[i]^sub[v])++black[v];
			if(tag[i])--black[v];
		}
		else black[sta[top]]=s[sta[top]]=='1'?(tag[i]?-1:1):0;
	}
	cout<<(black_bcc?0:pow2[ans])<<" ";
	for(int re i=1;i<=n;++i){
		black_bcc+=black[i],ans+=around[i]-du[i]+1;
		cout<<(black_bcc?0:pow2[ans])<<" ";
		black_bcc-=black[i],ans-=around[i]-du[i]+1;
	}
	cout<<"\n";
}

signed main(){
	pow2[0]=1;
	for(int re i=1;i<=2e5;++i)pow2[i]=mul(2,pow2[i-1]);
	for(int re T=getint();T--;){
		init();
		solve();
	}
	return 0;
}