差别

这里会显示出您选择的修订版和当前版本之间的差别。

--- 2020-2021:teams:legal_string:jxm2001:contest:arc_115 [2021/05/02 21:20]
jxm2001
+++ 2020-2021:teams:legal_string:jxm2001:contest:arc_115 [2021/05/02 21:36] (当前版本)
jxm2001
@@ 行 7: / 行 7: @@
 ==== 题意 ====
-给定一个图，询问满足度数为奇数的点的个数恰好为 $k(1\sim k\sim n)$ 的生成图数量。
+给定一个图，询问满足度数为奇数的点的个数恰好为 $k(1\sim n)$ 的生成图数量。
 其中生成图的定义为原图的点集 $+$ 原图的边集的子集。
@@ 行 29: / 行 29: @@
 然后等价于找 $n$ 个点中有 $k’$ 个点的生成图，根据数学归纳法，这是唯一的。
-接下来考虑连通图的情况，假设图中有 $n$ 个点和 $m$ 条边。
+接下来考虑连通图的情况，假设图中有 $n$ 个点和 $m$ 条边，先确定一棵生成树。
+对于非树边，有选于不选两种方案，同时非树边选与不选同时改变两个点的奇偶性。
+当 $k$ 为偶数时，任意指定 $k$ 个点为度数为奇数的点，先考虑所有非树边选与不选对所有点奇偶性的影响，共 $2^{m-n+1}$ 种情况。
+最后等价于构造图有 $k'$ 个点是度数为奇数的点的生成树的子图，此时方案数唯一。于是总方案数为 $2^{m-n+1}{n\choose k}$。
+对于不连通图，考虑分治 $\text{NTT}$ 处理背包的合并，时间复杂度 $O(n\log^2 n)$。
+<hidden 查看代码>
+<code cpp>
+const int Mod=998244353,MAXN=5005;
+int quick_pow(int a,int k){
+	int ans=1;
+	while(k){
+		if(k&1)ans=1LL*ans*a%Mod;
+		a=1LL*a*a%Mod;
+		k>>=1;
+	}
+	return ans;
+}
+namespace Poly{
+	const int G=3;
+	int rev[MAXN<<2];
+	int build(int k){
+		int n,pos=0;
+		while((1<<pos)<=k)pos++;
+		n=1<<pos;
+		_for(i,0,n)rev[i]=(rev[i>>1]>>1)|((i&1)<<(pos-1));
+		return n;
+	}
+	void NTT(int *f,int n,bool type){
+		_for(i,0,n)if(i<rev[i])
+		swap(f[i],f[rev[i]]);
+		int t1,t2;
+		for(int i=1;i<n;i<<=1){
+			int w=quick_pow(G,(Mod-1)/(i<<1));
+			for(int j=0;j<n;j+=(i<<1)){
+				int cur=1;
+				_for(k,j,j+i){
+					t1=f[k],t2=1LL*cur*f[k+i]%Mod;
+					f[k]=(t1+t2)%Mod,f[k+i]=(t1-t2)%Mod;
+					cur=1LL*cur*w%Mod;
+				}
+			}
+		}
+		if(!type){
+			reverse(f+1,f+n);
+			int div=quick_pow(n,Mod-2);
+			_for(i,0,n)f[i]=(1LL*f[i]*div%Mod+Mod)%Mod;
+		}
+	}
+	void Mul(int *f,int _n,int *g,int _m){
+		static int temp1[MAXN<<2],temp2[MAXN<<2];
+		int n=build(_n+_m-2);
+		_for(i,0,_n)temp1[i]=f[i];
+		_for(i,_n,n)temp1[i]=0;
+		_for(i,0,_m)temp2[i]=g[i];
+		_for(i,_m,n)temp2[i]=0;
+		NTT(temp1,n,true);NTT(temp2,n,true);
+		_for(i,0,n)temp1[i]=1LL*temp1[i]*temp2[i]%Mod;
+		NTT(temp1,n,false);
+		_for(i,0,_n+_m-1)f[i]=temp1[i];
+	}
+}
+int frac[MAXN],inv[MAXN],pow2[MAXN];
+int C(int n,int m){
+	return 1LL*frac[n]*inv[m]%Mod*inv[n-m]%Mod;
+}
+vector<int> g[MAXN];
+int blk_id[MAXN],blk_cnt;
+int pool[MAXN<<1],blk_sz[MAXN],blk_edge[MAXN],*blk_s[MAXN];
+void dfs(int u){
+	blk_id[u]=blk_cnt;
+	_for(i,0,g[u].size()){
+		int v=g[u][i];
+		if(blk_id[v])continue;
+		dfs(v);
+	}
+}
+void solve(int lef,int rig){
+	if(lef==rig)return;
+	int mid=lef+rig>>1;
+	solve(lef,mid);
+	solve(mid+1,rig);
+	int len1=1,len2=1;
+	_rep(i,lef,mid)len1+=blk_sz[i];
+	_rep(i,mid+1,rig)len2+=blk_sz[i];
+	Poly::Mul(blk_s[lef],len1,blk_s[mid+1],len2);
+}
+int main()
+{
+	frac[0]=pow2[0]=1;
+	_for(i,1,MAXN)frac[i]=1LL*i*frac[i-1]%Mod,pow2[i]=2LL*pow2[i-1]%Mod;
+	inv[MAXN-1]=quick_pow(frac[MAXN-1],Mod-2);
+	for(int i=MAXN-1;i;i--)inv[i-1]=1LL*inv[i]*i%Mod;
+	int n=read_int(),m=read_int();
+	_for(i,0,m){
+		int u=read_int(),v=read_int();
+		g[u].push_back(v);
+		g[v].push_back(u);
+	}
+	_rep(i,1,n){
+		if(!blk_id[i]){
+			blk_cnt++;
+			dfs(i);
+		}
+		blk_sz[blk_id[i]]++;
+		blk_edge[blk_id[i]]+=g[i].size();
+	}
+	int *pos=pool;
+	_rep(i,1,blk_cnt){
+		blk_s[i]=pos;
+		pos+=blk_sz[i]+1;
+		int base=pow2[blk_edge[i]/2-blk_sz[i]+1];
+		_rep(j,0,blk_sz[i]){
+			if(j%2==0)
+			blk_s[i][j]=1LL*C(blk_sz[i],j)*base%Mod;
+			else
+			blk_s[i][j]=0;
+		}
+	}
+	solve(1,blk_cnt);
+	_rep(i,0,n)
+	enter(blk_s[1][i]);
+	return 0;
+}
+</code>
+</hidden>
 ===== E - LEQ and NEQ =====

CVBB ACM Team

用户工具

站点工具

差别

页面工具