hdu4045 第二类斯特林，至少k个间隔选数

分两步，第一步n个数选r个，r个数的间隔数大于等于k-1，这个有公式！搞不懂，只能写进模板

r个不同数分成相同m组，可以为空，箱子装球问题7，第二类斯特林数的部分和

#include<iostream>
#include<stdio.h>
using namespace std;
typedef long long  ll;
const int mod=1e9+7;
ll s2[1005][1005];//第二类斯特灵
void initS()
{
	for(int i=1; i<=1000; i++)
	{
		s2[i][i]=s2[i][1]=1;//边界
	}
	for(int i=2; i<=1000; i++)
	{
		for(int j=1; j<i; j++)
		{
			s2[i][j]=s2[i-1][j-1]+1ll*j*s2[i-1][j];//递推式
			s2[i][j]%=mod;
		}
	}
}
ll B[1005][1005];
void initB()
{
	for(int i=1; i<=1000; i++)
	{
		B[i][1]=1;
		for(int j=2; j<=i; j++)
		{
			B[i][j]=B[i][j-1]%mod+s2[i][j]%mod;
			B[i][j]%=mod;
		}

	}
}
ll c[2005][1005];
void initC()
{
	for(int i=0; i<=2000; i++)
	{
		c[i][0]=c[i][i]=1;
	}
	for(int i=1; i<=2000; i++)
	{
		for(int j=1; j<=i&&j<=1000; j++)
		{
			c[i][j]=c[i-1][j-1]+c[i-1][j];//递推公式，1000比较小还能用 mod可以不是质数
			c[i][j]%=mod;
		}
	}
}

int main()
{
	ll n,r,k,m;
	initS();
	initB();
	initC();
	while(~scanf("%I64d%I64d%I64d%I64d",&n,&r,&k,&m))
	{
		ll sum=0;
		ll rem=n-((r-1)*k+1);
		if(rem<0)
		{
			printf("0\n");
			continue;
		}
//		for(int i=1;i<=min(r,m);i++)
//		{
//			sum+=s2[r][i]%mod;
//			sum%=mod;
//		}
		if(r<m)
			m=r;
		sum=c[rem+r][r]*B[r][m]%mod;
		cout<<sum%mod<<endl;
	}
}