[bzoj2169]连边

题目描述

有N个点（编号1到N）组成的无向图，已经为你连了M条边。请你再连K条边，使得所有的点的度数都是偶数。求有多少种连的方法。要求你连的K条边中不能有重边，但和已经连好的边可以重。不允许自环的存在。求连边的方法数。我们只关心它模10007的余数。

DP

设f[i,j]表示无顺序的连了i条边，有j个度数为奇数的点。
$f[i,j]=\frac{1}{i}(f[i-1,j-2]*C_{n-j+2}^2+f[i-1,j]*j*(n-j)+f[i-1,j+2]*C_{j+2}^2-f[i-2,j]*(C_n^2-(i-2)))$
是不是很巧妙呢？

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
const int maxn=1000+10,mo=10007;
int f[maxn][maxn],d[maxn];
int i,j,k,l,t,n,m,ans;
int qsm(int x,int y){
    if (!y) return 1;
    int t=qsm(x,y/2);
    t=t*t%mo;
    if (y%2) t=t*x%mo;
    return t;
}
int C(int n){
    return (n*(n-1)/2)%mo;
}
int main(){
    scanf("%d%d%d",&n,&m,&k);
    fo(i,1,m){
        scanf("%d%d",&j,&t);
        d[j]^=1;d[t]^=1;
    }
    t=0;
    fo(i,1,n) t+=d[i];
    f[0][t]=1;
    fo(i,1,k)
        fo(j,0,n){
            if (j>=2) f[i][j]=(f[i][j]+f[i-1][j-2]*C(n-j+2)%mo)%mo;
            f[i][j]=(f[i][j]+f[i-1][j]*j%mo*(n-j)%mo)%mo;
            if (j<=n-2) f[i][j]=(f[i][j]+f[i-1][j+2]*C(j+2)%mo)%mo;
            if (i>=2) f[i][j]=(f[i][j]-f[i-2][j]*(C(n)-(i-2))%mo)%mo;
            f[i][j]=f[i][j]*qsm(i,mo-2)%mo;
        }
    ans=f[k][0];
    //fo(i,1,k) ans=ans*qsm(i,mo-2)%mo;
    (ans+=mo)%=mo;
    printf("%d\n",ans);
}