小朋友的球

题目描述

@发源于 小朋友最近特别喜欢球。有一天他脑子抽了，从口袋里拿出了N个不同的球，想把它们放到M个相同的盒子里，并且要求每个盒子中至少要有一个球，他好奇有几种放法，于是尝试编程实现，但由于他天天不好好学习，只会上B站看游泳教练，于是他向你求助。

输入输出格式

输入格式：

 

多组数据，每行两个数N，M。

输出格式：

每组数据一行，表示方案数。

输入输出样例

输入样例#1：

4 2
1 1

输出样例#1：

7
1

说明

【样例解释】

N=4,M=2

1,2 3 4

2,1 3 4

3,1 2 4

4,1 2 3

1 2,3 4

1 3,2 4

1 4,2 3

对于20%的数据，满足1≤N，M≤10；

对于100%的数据，满足1≤N，M≤100，数据组数≤10。

stirling数，递推公式s[i][j]=s[i-1][j]*j+s[i-1][j-1]

 

S(p,k)的一个组合学解释是：将p个物体划分成k个非空的不可辨别的（可以理解为盒子没有编号）集合的方法数。

 

k!S(p,k)是把p个人分进k间有差别(如：被标有房号）的房间(无空房）的方法数。

 

S(p,k)的递推公式是：S(p,k)=k*S(p-1,k)+S(p-1,k-1) ,1<= k<=p-1

 

边界条件：S(p,p)=1 ,p>=0 S(p,0)=0 ,p>=1

 

递推关系的说明：

 

考虑第p个物品，p可以单独构成一个非空集合，此时前p-1个物品构成k-1个非空的不可辨别的集合，方法数为S(p-1,k-1)；

 

也可以前p-1种物品构成k个非空的不可辨别的集合，第p个物品放入任意一个中，这样有k*S(p-1,k)种方法。

高精度使用重载运算符会很方便

 

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
struct node{
    long long a[500],len;
    node()
    {memset(a,0,sizeof(a));len=0;}
};
node f[101][101];
node operator + (node c,node d)
{
    node h;
    h.len=max(c.len,d.len);
    int i;
    for(i=1;i<=h.len;i++){
        h.a[i]+=c.a[i]+d.a[i];
        if(h.a[i]>=10){
            h.a[i+1]+=h.a[i]/10;
            h.a[i]%=10;
            if(i==h.len)h.len++;
        }
    }
    return h;
}
node operator * (node c,long long u)
{
    node h;
    h.len=c.len;
    int i;
    for(i=1;i<=h.len;i++){
        h.a[i]+=c.a[i]*u;
        if(h.a[i]>=10){
            h.a[i+1]+=h.a[i]/10;
            h.a[i]%=10;
            if(i==h.len)h.len++;
        }
    }
    return h;
}
int main()
{int n,m,i,j;
    while (cin>>n>>m)
    {
        if (n<m)
        {
            printf("0\n");
        }
        else if (n==m)
        {
            printf("1\n");
        }
        else 
        {
            for (i=0;i<=n;i++)
            f[i][i].a[1]=1,f[i][1].a[1]=1,f[i][i].len=f[i][1].len=1;
            for (i=2;i<=n;i++)
            {
                for (j=2;j<=i-1;j++)
                 f[i][j]=f[i-1][j]*j+f[i-1][j-1];
            }

      for (i=f[n][m].len;i>=1;i--)
      printf("%lld",f[n][m].a[i]);
      printf("\n");
        }
        
    }
}