题目:http://acm.hdu.edu.cn/showproblem.php?pid=1521
因为要求解的问题是从n种物品中抽取m个物品,由于每种物品的数量为多个所以有一定难度。
假如这m个物品中包含的k种物品,且这k种物品中的每个物品的数量为ni(i=1,2,3.....k);
由基本的数学中的排列组合知识易知:结果为 ( m ! ) /(n1!*n2!*n3!........*nk!)
本题目采用指数型母函数解题,以x的次幂表示每种物品所取的数目,由上面的公式易知,x项的系数为x项指数的阶乘。
下面是本题目的代码:
#include<iostream>
#include<cstdio>
using namespace std;
const int Max = 11;
double fac[Max],ans[Max],tmp[Max];
void InitFac()
{
fac[0]=1;
for(int i=1;i<Max;i++)
fac[i]=fac[i-1]*i;
}
int main()
{
InitFac();
int n,m,i,j,k;
int Num[Max];
while(cin>>n>>m)
{
memset(Num,0,sizeof(Num));
for(i=0;i<n;i++)
cin>>Num[i];
memset(ans,0,sizeof(ans));
memset(tmp,0,sizeof(tmp));
for(i=0;i<=Num[0];i++)
ans[i]=1.0/fac[i];
for(i=1;i<n;i++)
{
for(j=0;j<Max;j++)
for(k=0;(k<=Num[i]&&k+j<Max);k++)
tmp[k+j]+=ans[j]/fac[k];
for(j=0;j<Max;j++)
ans[j]=tmp[j],tmp[j]=0;
}
printf("%.lf\n",ans[m]*fac[m]);
}
return 0;
}