Codeforces 568B Symmetric and Transitive

http://codeforces.com/contest/568/problem/B

题意:题意还挺绕的，其实就是说:要你求出一个图，要求保证其中有至少一个点不连任何边，然后其他连边的点构成的每个联通块都必须构成完全连通图

思路:f[i][j]代表i个点，构成j个联通块的方案数

f[i][j]=f[i][j-1]*j(代表与其中一个联通块合并)+f[i-1][j-1](代表新开一个联通块)

然后答案是Σc[n][i]*f[i][j]

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<iostream>
#define ll long long
const ll Mod=1000000007;
ll f[4005][4005],c[4005][4005];
int n;
ll Pow(ll x,ll y){
    ll res=1;
    while (y){
        if (y&1) res=(res*x)%Mod;
        y/=2;
        x=(x*x)%Mod;
    }
    return res;
}
int read(){
    int t=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while ('0'<=ch&&ch<='9'){t=t*10+ch-'0';ch=getchar();}
    return t*f;
}
int main(){
    n=read();
    ll ans=0;
    for (int i=1;i<=n;i++) c[i][0]=c[i][i]=1;
    for (int i=2;i<=n;i++)
     for (int j=1;j<i;j++)
      c[i][j]=(c[i-1][j]+c[i-1][j-1])%Mod;
    f[0][0]=1;  
    for (int i=1;i<=n;i++)
     for (int j=1;j<=i;j++)
      f[i][j]=((1LL*f[i-1][j]*j)%Mod+f[i-1][j-1])%Mod;
    for (int i=0;i<n;i++)
     for (int j=0;j<=i;j++)
      ans=(ans+f[i][j]*c[n][i])%Mod;
    printf("%lld
",ans);  
    return 0;
}