hdu 4045 Machine scheduling 【第二类斯特林数】

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

---

 

题目大意：有N个机器，每天选出R个机器，而且每两个机器的编号差要大于等于K，每天将这R个机器最多分为M组工作，问最多                      有多少种方案。

思路：第一步先求出N个机器选R个满足条件的方案数；第二部将R个机器放入M个组的方案数，然后两者相乘；

先满足每两个机器之间至少有K-1个间隔，也就是还剩下X=n-((r-1)*k+1)个机器可以随意安排，把这些多余的插入到R个机器之间(加上两端共R+1个位置)。问题也就变为X个相同的球分到R+1个不同的组可以为空这种模型，不难推出是C(X+R，R)，网上有说直接用插板法公式。

第二部分：R个元素分为i个非空集合是第二类斯特林数，对i为1至m求和即可。
（第二类斯特林数定义：将N个不同的元素划分到M个集合的方案书）

 

---

AC代码：

#include<cstdio>
#include<iomanip>
#include<cstdlib>
#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#include<math.h>
#include<map>
#include<set>
#include<queue>
#include<vector>
#include<stack>
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
const double pi=acos(-1.0);
const ll mod=1e9+7;
const int M=100000;
const int N=2e3+10;

ll stl[N][N];//第二类斯特林数
ll c[N][N];//组合数C(i,j)

void work()
{  //递推
    for(int i=1;i<=2000;i++)
    {
        c[i][0]=c[i][i]=1;
        for(int j=1;j<i;j++)
            c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
    }
    stl[1][1]=1;
    for(int i=2;i<=1000;i++)
    {
        for(int j=1;j<=i;j++)
            stl[i][j]=(stl[i-1][j-1]+j*stl[i-1][j])%mod;
    }
}

int main()
{
    work();
    int n,r,k,m;
    while(~scanf("%d %d %d %d",&n,&r,&k,&m))
    {
        int x=n-((r-1)*k+1);
        if(x<0) printf("0\n");
        else
        {
            ll ans=c[x+r][r],sum=0;
            for(int i=1;i<=m;i++)
                sum=(sum+stl[r][i])%mod;
            ans=(sum*ans)%mod;
            printf("%lld\n",ans);
        }
    }
    return 0;
}