Description
有N个点(编号1到N)组成的无向图,已经为你连了M条边。请你再连K条边,使得所有的点的度数都是偶数。求有多少种连的方法。要求你连的K条边中不能有重边,但和已经连好的边可以重。不允许自环的存在。求连边的方法数。我们只关心它模10007的余数。
Input
输入的第一行有三个自然数,分别表示点数N,已经连好的边数M,和你要连的边数K。保证K≤N(N-1)/2
接下来M行每行两个整数x,y,描述了一条连接x和y的边。
30%的数据满足:N≤200
100%的数据满足:N≤1000,M≤N,K≤1000,K≤N(N-1)/2
Output
输出一个整数,表示连边的方法数模10007的余数
Sample Input
5 1 4
1 2
Sample Output
13
Hint
【样例说明】
以下是13种连边的方法(只显示你连的边):
{(1,2),(1,3),(1,4),(3,4)}
{(1,2),(1,3),(1,5),(3,5)}
{(1,2),(1,4),(1,5),(4,5)}
{(1,2),(2,3),(2,4),(3,4)}
{(1,2),(2,3),(2,5),(3,5)}
{(1,2),(2,4),(2,5),(4,5)}
{(1,2),(3,4),(3,5),(4,5)}
{(1,3),(2,4),(3,5),(4,5)}
{(1,3),(2,5),(3,4),(4,5)}
{(1,4),(2,3),(3,5),(4,5)}
{(1,4),(2,5),(3,4),(3,5)}
{(1,5),(2,3),(3,4),(4,5)}
{(1,5),(2,4),(3,4),(3,5)}
思路
DP。
令f[i][j]表示加i条边,j个奇点的方案
- 当两端点的度数都为奇数,则奇数点的个数-2,so f[i][j]+=f[i-1][j-2]*C(n-j+2,2)
- 当两端点的度数都为偶数,则奇数点的个数+2,so f[i][j]+=f[i-1][j+2]*C(j+2,2)
- 当两端点的度数一奇一偶,则奇数点的个数不变 , so f[i][j]+=f[i-1][j]*(n-j)*j
减去重复:
f[i][j]-=f[i-2][j]*(C(n,2)-i+2)
有序变无序
ans/=i
代码
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int b[1001],cnt,n,m,k;
int mo=10007;
long long f[1001][1001],a[1001];
long long C(int x)
{
if (x<2) return 0;
return (x-1)*x/2;
}
int main()
{
int i,j,u,v;
cin>>n>>m>>k;
a[0]=1;a[1]=1;
for (i=2;i<=1000;i++)
a[i]=(mo-mo/i)*a[mo%i]%mo;
for (i=1;i<=m;i++)
{
scanf("%d%d",&u,&v);
b[u]++;b[v]++;
}
for (i=1;i<=n;i++) if (b[i]%2==1) cnt++;
f[0][cnt]=1;
for (i=1;i<=k;i++)
{
for (j=0;j<=n;j++)
{
if (j>=2)
f[i][j]+=f[i-1][j-2]*C(n-j+2)%mo;
f[i][j]%=mo;
if (j+2<=n)
f[i][j]+=f[i-1][j+2]*C(j+2)%mo;
f[i][j]%=mo;
f[i][j]+=f[i-1][j]*j*(n-j)%mo;
f[i][j]%=mo;
if (i>=2) f[i][j]-=f[i-2][j]*(C(n)-i+2)%mo;
f[i][j]=(f[i][j]+mo)%mo;
f[i][j]*=a[i];f[i][j]%=mo;
}
}
cout<<f[k][0];
}