标签:style string 小朋友 ring 方案 编程 使用 int 求助
@发源于 小朋友最近特别喜欢球。有一天他脑子抽了,从口袋里拿出了N个不同的球,想把它们放到M个相同的盒子里,并且要求每个盒子中至少要有一个球,他好奇有几种放法,于是尝试编程实现,但由于他天天不好好学习,只会上B站看游泳教练,于是他向你求助。
输入格式:
多组数据,每行两个数N,M。
输出格式:
每组数据一行,表示方案数。
4 2 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)种方法。
高精度使用重载运算符会很方便
1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 #include<cstring> 5 using namespace std; 6 struct node{ 7 long long a[500],len; 8 node() 9 {memset(a,0,sizeof(a));len=0;} 10 }; 11 node f[101][101]; 12 node operator + (node c,node d) 13 { 14 node h; 15 h.len=max(c.len,d.len); 16 int i; 17 for(i=1;i<=h.len;i++){ 18 h.a[i]+=c.a[i]+d.a[i]; 19 if(h.a[i]>=10){ 20 h.a[i+1]+=h.a[i]/10; 21 h.a[i]%=10; 22 if(i==h.len)h.len++; 23 } 24 } 25 return h; 26 } 27 node operator * (node c,long long u) 28 { 29 node h; 30 h.len=c.len; 31 int i; 32 for(i=1;i<=h.len;i++){ 33 h.a[i]+=c.a[i]*u; 34 if(h.a[i]>=10){ 35 h.a[i+1]+=h.a[i]/10; 36 h.a[i]%=10; 37 if(i==h.len)h.len++; 38 } 39 } 40 return h; 41 } 42 int main() 43 {int n,m,i,j; 44 while (cin>>n>>m) 45 { 46 if (n<m) 47 { 48 printf("0\n"); 49 } 50 else if (n==m) 51 { 52 printf("1\n"); 53 } 54 else 55 { 56 for (i=0;i<=n;i++) 57 f[i][i].a[1]=1,f[i][1].a[1]=1,f[i][i].len=f[i][1].len=1; 58 for (i=2;i<=n;i++) 59 { 60 for (j=2;j<=i-1;j++) 61 f[i][j]=f[i-1][j]*j+f[i-1][j-1]; 62 } 63 64 for (i=f[n][m].len;i>=1;i--) 65 printf("%lld",f[n][m].a[i]); 66 printf("\n"); 67 } 68 69 } 70 }
标签:style string 小朋友 ring 方案 编程 使用 int 求助
原文地址:http://www.cnblogs.com/Y-E-T-I/p/7112834.html