题意:
给出
n
,
m
(
0
≤
n
≤
50
,
0
≤
m
≤
50
)
n,m(0\leq n\leq 50,0\leq m\leq50)
n,m(0≤n≤50,0≤m≤50)
将n只带编号的猪放入m个不带编号的房子里面,要求每个房子至少放一只猪,求方案数
题解:
第二类斯特林数
S(n,m)=S(n-1,m)*m+S(n-1,m-1);
要用高精,这个数据是会爆long long的,另外记忆化搜索下,否则时间也会爆炸
code:
#include<bits/stdc++.h>
#define ll long long
using namespace std;
inline int read(){
int x=0;char ch=getchar();
while(!isdigit(ch))ch=getchar();
while(isdigit(ch))x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
return x;
}
struct Bigint{
bool iszero;
int p[505],len;
Bigint(){memset(p,0,sizeof p);len=0;iszero=true;}
Bigint(int x){
memset(p,0,sizeof p);len=0;iszero=true;
if(x){
iszero=false;
while(x){
p[len++]=x%10;
x/=10;
}
}
}
}dp[55][55];
void print(Bigint p){
if(p.iszero){
printf("0");
return;
}
for(int i=p.len-1;i>=0;i--)
printf("%d",p.p[i]);
}
Bigint operator + (const Bigint &a,const Bigint &o){
Bigint re;
int carry=0;re.iszero=false;
for(int i=0;i<max(a.len,o.len);i++,re.len++){
re.p[i]=a.p[i]+o.p[i]+carry;
carry=re.p[i]/10;
re.p[i]%=10;
}
while(carry){
re.p[re.len++]=carry%10;
carry/=10;
}
return re;
}
Bigint operator * (Bigint a,int o){
Bigint re;
re.iszero=false;
int carry=0;
for(int i=0;i<a.len;i++,re.len++){
re.p[i]=a.p[i]*o+carry;
carry=re.p[i]/10;
re.p[i]%=10;
}
while(carry){
re.p[re.len++]=carry%10;
carry/=10;
}
return re;
}
Bigint S(int n,int m){
if(n<m)return Bigint(0);
if(n==m||m==1)return Bigint(1);
if(!dp[n][m].iszero)return dp[n][m];
dp[n][m].iszero=false;
return dp[n][m]=S(n-1,m)*m+S(n-1,m-1);
}
int main(){
int n=read(),m=read();
print(S(n,m));
return 0;
}