HDU1521排列组合

题目：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;
}