hdu 4045 2011北京赛区网络赛F 组合数+斯特林数 ***

插板法基础知识

斯特林数见百科

 #include<iostream>
 #include<cmath>
 #include<cstdio>
 #include<cstring>
 #define LL long long
 #define eps 1e-7
 #define MOD 1000000007
 using namespace std;
 int c[][]={},stir2[][]={};
 int main(){
     for(int i=;i<=;i++){
         c[i][]=c[i][i]=;
         for(int j=;j<i;j++)
             c[i][j]=(c[i-][j-]+c[i-][j])%MOD;
     }
     for(int i=;i<=;i++){
         stir2[i][]=;
         stir2[i][i]=;
         for(int j=;j<i;j++)
             stir2[i][j]=((LL)j*stir2[i-][j]+stir2[i-][j-])%MOD;
     }
     int n,r,k,m;
     while(cin>>n>>r>>k>>m){
         int sum=;
         if(n-((r-)*k+)<){
             cout<<<<endl;
             continue;
         }
         for(int i=;i<=min(r,m);i++)
             sum=(sum+stir2[r][i])%MOD;
         cout<<((LL)c[n-((r-)*k+)+r+-][r]*sum)%MOD<<endl;
     }
     return ;
 }