[atARC115D]Odd Degree

考虑对于一棵树$G$,这个问题的答案——

当$k$为奇数时答案显然为0,否则从$V$中任选$k$个点,以任意一点为根,从底往上不难发现子图数量唯一

换言之,当$k$为偶数时,每一个合法(恰有$k$个奇度数的点)子图恰好对应于一种选择方案,即${|V|\choose k}$

当$G$是一张连通图时,继续来分析答案——

首先$k$仍要是偶数,且仍然考虑任选$k$点,并求出其一棵生成树

对于生成树以外的边,任意选每一条边是否加入子图,之后同样可以通过这棵生成树构造出一组方案,换言之每一组选点方案恰对应于$2^{|E|-|V|+1}$个子图,答案即${|V|\choose k}2^{|E|-|V|+1}$

有多个连通块时,显然每一个连通块独立,用$f_{i,j}$表示前$i$个连通块中选$j$个点,枚举最后一个所选的点转移即可,复杂度为$o(n^{2})$,可以通过

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define N 5005
 4 #define mod 998244353
 5 int n,m,x,y,mi[N],fac[N],inv[N],f[N],szV[N],szE[N],dp[N][N];
 6 int c(int n,int m){
 7     return 1LL*fac[n]*inv[m]%mod*inv[n-m]%mod;
 8 }
 9 int find(int k){
10     if (k==f[k])return k;
11     return f[k]=find(f[k]);
12 }
13 void merge(int x,int y){
14     x=find(x),y=find(y);
15     if (x==y)szE[x]++;
16     else{
17         f[x]=y;
18         szV[y]+=szV[x];
19         szE[y]+=szE[x]+1;
20     }
21 }
22 int main(){
23     mi[0]=fac[0]=inv[0]=inv[1]=1;
24     for(int i=1;i<N;i++)mi[i]=2*mi[i-1]%mod;
25     for(int i=1;i<N;i++)fac[i]=1LL*fac[i-1]*i%mod;
26     for(int i=2;i<N;i++)inv[i]=1LL*(mod-mod/i)*inv[mod%i]%mod;
27     for(int i=1;i<N;i++)inv[i]=1LL*inv[i-1]*inv[i]%mod;
28     scanf("%d%d",&n,&m);
29     for(int i=1;i<=n;i++){
30         f[i]=i;
31         szV[i]=1,szE[i]=0;
32     }
33     for(int i=1;i<=m;i++){
34         scanf("%d%d",&x,&y);
35         merge(x,y);
36     }
37     dp[0][0]=1;
38     int scc=0;
39     for(int i=1;i<=n;i++)
40         if (f[i]==i){
41             scc++;
42             for(int j=0;j<=n;j++)
43                 for(int k=0;k<=min(j,szV[i]);k+=2)
44                     dp[scc][j]=(dp[scc][j]+1LL*mi[szE[i]-szV[i]+1]*c(szV[i],k)%mod*dp[scc-1][j-k])%mod;
45         }
46     for(int i=0;i<=n;i++)printf("%d\n",dp[scc][i]);
47 }
View Code

 

posted @ 2021-03-28 21:59  PYWBKTDA  阅读(222)  评论(0编辑  收藏  举报