51NOD1835 完全图

传送门

分析

令f(i,j)表示i点完全图有j个联通块的方案数。

讨论有i-1个点已经固定了，我们拉出一个代表元素然后讨论它的集合大小然后组合数算一下就可以了。

$$ dp(i,j) = \sum_{k=1}^{i-1} C_{i-1}^{k-1} dp(i-k,j-1) dp(k,1) $$

$$ dp(i,1) = 2^{\frac{i(i-1)}{2}} - \sum{j=2}^i dp(i,j) $$

因为至少要删一条边，所以m=1时候最终答案要减1。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cctype>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<ctime>
#include<vector>
#include<set>
#include<map>
#include<stack>
using namespace std;
const long long mod = ;
#define add(x,y) x=(x+y)%mod
long long c[][],dp[][];
inline long long pw(long long x,long long p){
    long long res=;
    while(p){
      if(p&)res=res*x%mod;
      x=x*x%mod;
      p>>=;
    }
    return res;
}
int main(){
    long long n,m,i,j,k;
    scanf("%lld%lld",&n,&m);
    c[][]=;
    for(i=;i<=;i++)
      c[i][]=c[i][i]=;
    for(i=;i<=;i++)
      for(j=;j<i;j++)
        c[i][j]=(c[i-][j]+c[i-][j-])%mod;
    dp[][]=;
    for(i=;i<=n;i++){
      for(j=;j<=i;j++)
        for(k=;k<i;k++)
          add(dp[i][j],c[i-][k-]*dp[i-k][j-]%mod*dp[k][]%mod);
      dp[i][]=pw(,i*(i-)/);
      for(j=;j<=i;j++)
        dp[i][]=(dp[i][]-dp[i][j]+mod)%mod;
    }
    printf("%lld\n",m==?dp[n][m]-:dp[n][m]);
    return ;
}