bzoj 2169 连边 —— DP+容斥

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=2169

就和这篇博客说的一样：https://blog.csdn.net/WerKeyTom_FTD/article/details/70274470

注意每次是 /i 而不是 /(i!)，因为 i-1 时也已经去了重，现在就是对于新加一条边的多种方式带来一种局面去重，从每一种局面看，新加的边可以是任意一条，所以 /i。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
int const xn=,mod=;
int n,m,k,deg[xn];
ll f[xn][xn];
int rd()
{
  int ret=,f=; char ch=getchar();
  while(ch<''||ch>''){if(ch=='-')f=; ch=getchar();}
  while(ch>=''&&ch<='')ret=(ret<<)+(ret<<)+ch-'',ch=getchar();
  return f?ret:-ret;
}
ll pw(ll a,int b)
{
  ll ret=;
  for(;b;b>>=,a=(a*a)%mod)
    if(b&)ret=(ret*a)%mod;
  return ret;
}
ll calc(int x){return (ll)x*(x-)/;}
int main()
{
  n=rd(); m=rd(); k=rd(); int num=;
  for(int i=,x,y;i<=m;i++)x=rd(),y=rd(),deg[x]++,deg[y]++;
  for(int i=;i<=n;i++)if(deg[i]&)num++;
  f[][num]=;
  for(int i=;i<=k;i++)
    for(int j=;j<=n;j++)
      {
    f[i][j]=(f[i][j]+f[i-][j]*((ll)j*(n-j)%mod)%mod)%mod;
    if(j>=)f[i][j]=(f[i][j]+f[i-][j-]*calc(n-j+)%mod)%mod;//+2!
    if(j<=n-)f[i][j]=(f[i][j]+f[i-][j+]*calc(j+)%mod)%mod;
    if(i>=)f[i][j]=(f[i][j]-(f[i-][j]*(calc(n)-i+))%mod+mod)%mod;
    f[i][j]=(f[i][j]*pw(i,mod-))%mod;//
      }
  printf("%lld\n",f[k][]);
  return ;
}