SDNU 1011（第二类stirling数）

1011.盒子与球

Time Limit: 1000 MS    Memory Limit: 32768 KB
Total Submission(s): 667    Accepted Submission(s): 165

Description

现有r个互不相同的盒子和n个互不相同的球，要将这n个球放入r个盒子中，且不允许有空盒子。则有多少种放法？

Input

n, r(0 <= n, r <= 10)。

Output

有多少种放法。

Sample Input

3 2

Sample Output

6

Source

SDNU ACM-ICPC 2010复赛(2010级)

把从1到n标号的n个球放到k个无区别的盒子里，要求每个盒子里至少有一个小球，问不同的放法数量。例如，如果用A、B、C、D分别表示4个球，要分成两组（即放入无区别的盒子里），其方法有7种：
{A,B},{C,D}  {A},{B,C,D}
{A,C},{B,D}  {B},{A,C,D}
{A,D},{B,C}  {C},{A,B,D}
                    {D},{A,B,C}

这个数量可以用第二类斯特林 (Stirling) 数来计算，表示为S(n,k)，S(4,2)=7。第二类斯特林 (Stirling) 数也是计算机科学应用中很常见的公式。它有如下的递推公式：

整数参数n≥k≥0，且初始条件满足

式中是Knuth推荐的第二类斯特林数的表

S（n,k） =  S(n-1,k-1)  +  S(n-1,k)*k,  n>=k;

#include <stdio.h>
#include <math.h>
#include <string.h>

#define N 10
#define R 10

int dp[N][R];
int i,j;

int Stirling()
{
     memset(dp,0,sizeof(dp));
     for(i=1;i<=N;i++)
          dp[i][1]=1;
     for(i=2;i<=N;i++)
         for(j=1;j<=i;j++)
              dp[i][j]=dp[i-1][j-1]+dp[i-1][j]*j;
     return 0;
}

int main()
{
     int n,r;
     Stirling();
     while(scanf("%d%d",&n,&r)!=EOF)
     {
          int arr=1;
          for(i=1;i<=r;i++)
               arr*=i;
          printf("%d\n",dp[n][r]*arr);
     }
     return 0;
}