bzoj5303[HAOI2018]反色游戏（tarjan，线性基）

题目描述
小C和小G经常在一起研究搏弈论问题，有一天他们想到了这样一个游戏． 有一个n个点m条边的无向图，初始时每个节点有一个颜色，要么是黑色，要 么是白色．现在他们对于每条边做出一次抉择：要么将这条边连接的两个节点都 反色（黑变白，白变黑），要么不作处理．他们想把所有节点都变为白色，他们 想知道在2^m
种决策中，有多少种方案能达成这个目标． 小GG认为这个问题太水了，于是他还想知道，对于第ii个点，在删去这个点及 与它相连的边后，新的答案是多少． 由于答案可能很大，你只需要输出答案对10^9+7取模后的结果．
输入输出格式
输入格式：
从文件game.in中读入数据. 第一行一个整数T，表示数据组数． 每组数据第一行两个整数n,m表示点数和边数. 接下来m行，每行两个整数u，v,描述无向图的一条边． 接下来一行一个长度为n的0/1串，如果第i个字符为0表示第i个点为白色，否 则为黑色．
输出格式：
输出到文件game.out中. 每组数据输出一行n+1个整数，第一个整数表示不删去任何点时的答案．接 下来n个整数，第i个表示删去第ii个点时的答案．
输入输出样例
输入样例#1：
2
5 5
1 2
2 3
3 4
2 4
3 5
00000
5 4
1 2
2 3
2 4
2 5
11111
输出样例#1：
2 2 1 1 1 2
0 1 0 1 1 1
说明
对于所有数据，有$1\le T\le 5,1\le n,m\le 10^5$
$1\le u,v\le n$，没有重边和自 环．

---

$solution$：
考虑无解的情况。

可以发现，每个联通块里面只要有偶数个黑点就是合法的，只要两两路径取反即可。如果是奇数个黑点就是不合法的。

考虑怎么计算方案数

对于一张图，我们把边看成未知数，可以列出$n$个异或方程，那么最后的答案就是$2^{自由元}$

分治消元好像不是很行。

可以发现对于一个大小为$n$的联通块，取$n-1$个主元后最后一行肯定可以被消完，并且必定能选取$n-1$条边，这个可以放到生成树上考虑

那么答案就是$2^{m-n+联通块个数}$

现在要做的就是统计删去每个点后是否合法以及新的联通块个数

上$tarjan$就行了

#include<bits/stdc++.h>
using namespace std;
#define rep(i,j,k) for(int i = j;i <= k;++i)
#define repp(i,j,k) for(int i = j;i >= k;--i)
#define rept(i,x) for(int i = linkk[x],y = e[i].y;i;i = e[i].n,y = e[i].y)
#define P pair<int,int>
#define Pil pair<int,ll>
#define Pli pair<ll,int>
#define Pll pair<ll,ll>
#define pb push_back 
#define pc putchar
#define mp make_pair
#define file(k) memset(k,0,sizeof(k))
#define ll long long
int rd()
{
	int num = 0;char c = getchar();bool flag = true;
	while(c < '0'||c > '9') {if(c == '-') flag = false;c = getchar();}
	while(c >= '0' && c <= '9') num = num*10+c-48,c = getchar();
	if(flag) return num;else return -num;
}
const int p = 1e9+7;
int n,m,tot;
int linkk[101000],t,a[101000];
int dfn[101000],low[101000],root;
int st[101000],top;
int mi[201000],mx;
int v[101000],bel[101000];//表示i子树内的权值 
int du[101000],out[101000];
bool cut[101000],ck[101000],only[101000];
char s[101000];
struct node{int n,y;}e[201000];
inline int max(int a,int b){return a>b?a:b;}
inline int min(int a,int b){return a<b?a:b;}
inline int mul(int a,int b){return (a*=b)>=p?a-=p:a;}
void insert(int x,int y)
{
	e[++t].y = y;e[t].n = linkk[x];linkk[x] = t;du[y]++;
	e[++t].y = x;e[t].n = linkk[y];linkk[y] = t;du[x]++;
}
void tarjan(int x)
{
	dfn[x] = low[x] = ++tot;bel[x] = root;
	v[x] = a[x];ck[x] = false;cut[x] = false;only[x] = false;out[x] = 0;
	if(x == root && linkk[x] == 0)
	{
		only[x] = true;
		return;
	}
	rept(i,x)
		if(!dfn[y])
		{
			tarjan(y);
			low[x] = min(low[x],low[y]);
			v[x] += v[y];
		    if(low[y] >= dfn[x])
		    {
		    	out[x]++;
		    	if(x != root || out[x] > 1) cut[x] = true;
		        ck[x] |= (v[y]&1); 
			}
		}
		else low[x] = min(low[x],dfn[y]);
	out[x]++;
	if(x == root) out[x]--;
}
void init()
{
	n = rd();m = rd();
	
	rep(i,1,n) linkk[i] = dfn[i] = du[i] = 0;
	t = 0;tot = 0;
	
	rep(i,1,m)
	{
		int x = rd(),y = rd();
		insert(x,y);
	}
	int flag = 0;
	int num = 0;
	scanf("%s",s+1);
	rep(i,1,n) a[i] = s[i] == '1';
	rep(i,1,n) if(!dfn[i]) root = i,tarjan(i),flag += v[root]&1,num++;
	rep(i,1,n) ck[i] |= ( v[bel[i]] - a[i]) & 1;
	if(flag) printf("0");else printf("%d",mi[m-n+num]);
	rep(x,1,n) if(!cut[x])
	{
		if(!a[x]) {if(flag) printf(" 0");else printf(" %d",mi[m-du[x]-n+1+num-only[x]]);}
		else if( flag+ ((v[bel[x]] - a[x])&1?1:-1)) printf(" 0");
		     else printf(" %d",mi[m-du[x] - n+1 +num]);
	}
	else
	{
		if(ck[x]) printf(" 0");
		else if(flag-(v[bel[x]]&1)) printf(" 0"); 
		else printf(" %d",mi[m-du[x]-n+num+out[x]]);
	}
	printf("\n");
}
int main()
{
    int T = rd();mi[0] = 1;
    rep(i,1,200000) mi[i] = mul(mi[i-1],2);
    while(T--) init();
	return 0;
}