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 ;

 }