hdu 4045 Machine scheduling

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4045

题目大意:从n个机器(编号1~n)中选出r个两两编号差不小于k的机器分成x(1<=x<=m)组,求共有几种分法.

题目思路:
把题拆成两部分:
第一部分:选出r个两两编号差不小于k的机器.
第二部分:把r个机器分成x(1<=x<=m)组,很明显的斯特林第二类数.

第一部分比较难,个人的数学比较差,想了超久...
可以把第一部分转换为n个人(一模一样的人)排队,其中有r个人两两间至少隔k个人.
第一步:把r个人先取出排好,然后再相邻两人之间插入k个人,这样就满足了条件,然后还剩下n-r-(r-1)*k个人.
第二步:可以知道已排好的队伍的头,尾和k个人所在的那些位置,无论怎么插人,都满足题意,所以我们可以把剩下的人分成r组(可以插的位置只有r个,k个人在一起算一个可插位置),方法就是最普遍的隔板法.
通过前两步,我们就把问题转换成了在n-r-(r-1)*k+1个位置(剩下n-r-(r-1)*k个人,头尾皆可插板)插入r块板,有几种插法.得到公式:C(n-(r-1)*(k-1),r).

代码:

#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <ctype.h>
#include <math.h>
#include <time.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;
#define ull unsigned __int64
//#define ll __int64
//#define ull unsigned long long
#define ll long long
#define son1 New(p.xl,xm,p.yl,ym),(rt<<2)-2
#define son2 New(p.xl,xm,min(ym+1,p.yr),p.yr),(rt<<2)-1
#define son3 New(min(xm+1,p.xr),p.xr,p.yl,ym),rt<<2
#define son4 New(min(xm+1,p.xr),p.xr,min(ym+1,p.yr),p.yr),rt<<2|1
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define middle (l+r)>>1
#define MOD 1000000007
#define esp (1e-8)
const int INF=0x3F3F3F3F;
const double DINF=10000.00;
//const double pi=acos(-1.0);
const int M=1001;
ll min(ll x,ll y){return x<y? x:y;}
int max(int x,int y){return x>y? x:y;}
void swap(ll& x,ll& y){ll t=x;x=y;y=t;}
int T,cas;

int n,r,k,m;
ll s[M][M],C[M][M];

void run(){
	//scanf("%d%d%d%d",&n,&r,&k,&m);
	ll sum=0,nn=n-(r-1)*(k-1);
	for(int i=1;i<=m;i++) sum=(sum+s[r][i])%MOD;
	//printf("%I64d %I64d\n",sum1,sum2);
	printf("%I64d\n", nn>0? C[nn][r]*sum%MOD:0);
}


int main(){
	//freopen("1.in","r",stdin);
	//freopen("1.out","w",stdout);
	//run();
	//scanf("%d",&T);for(cas=1;cas<=T;cas++) run();
	for(int i=1;i<M;i++){
		s[i][i]=1,s[i][0]=0;
		for(int j=1;j<i;j++)
			s[i][j]=(s[i-1][j]*(ll)j+s[i-1][j-1])%MOD;
	}
	for(int i=1;i<M;i++){
		C[i][0]=C[i][i]=1;
		for(int j=1;j<=(i>>1);j++){
			C[i][j]=C[i][i-j]=(C[i-1][j]+C[i-1][j-1])%MOD;
		}
	}
	while(~scanf("%d%d%d%d",&n,&r,&k,&m))  run();
	return 0;
}