HDU 1521 排列组合（指数型母函数）

HDU - 1521

指数型母函数也超级神奇，看到这样一句话：普通型母函数主要是来求组合的方案数，而指数型母函数是求多重排列数

运用指数型母函数我们可以将求解多重排列数，即排列组合混合运算

这个博客讲得很好，例子很生动： 指数型母函数博客

关于指数型母函数，想说一下自己的一些见解

和普通型母函数相似，选择的数的个数用x的幂来表示，只不过我们在进行运算的时候不一定是单纯的排列数，可能有重复的，这个时候我们需要把重复的情况除出去：

假设有n个元素，其中a1，a2，····，an互不相同，进行全排列，可得n!个不同的排列。若其中某一元素a1重复了n1次，全排列出来必有重复元素，其中真正不同的排列数应为n!/n1!，即其重复度为n1!

同样理由a1重复了n1次，a2重复了n2次，····，ak重复了nk次，n1+n2+····+nk=n。对于这样的n个元素进行全排列，可得不同排列的个数实际上是

（摘自百度百科）

假设是这样，需要挑9个数，假设其中一种情况是有两类数，分别有4个相同的和5个相同的，此时结果应该是9！/（4！*5！），我们先把所有类似该类符合情况的分母（1 / （4！*5！））加起来，最后乘以9！即为最终总的情况个数，非常神奇也很好理解，原理就是把所有组成9个数的情况中重复的情况除去，再把各种符合条件的情况相加，代码实现也和普通型母函数十分相似

直接挂代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<set>
#include<map>
#include<queue>
#include<cmath>
using namespace std;
//普通型母函数主要是来求组合的方案数，而指数型母函数是求多重排列数
#define maxn 15
int fac[15];
void init()
{
    fac[0] = fac[1] = 1;
    for(int i = 2; i <= 10; i ++)
        fac[i] = fac[i - 1] * i;
}
double n1[maxn * 10], n2[maxn * 10];
int num[maxn];
int main()
{
    int n, m;
    init();
    while(scanf("%d%d",&n,&m) != EOF)
    {
        //init();
        memset(n1, 0, sizeof(n1));
        memset(n2, 0, sizeof(n2));
        for(int i = 1; i <= n; i ++)
            scanf("%d",&num[i]);
        for(int i = 0; i <= num[1]; i ++)
            n1[i] = 1.0 / fac[i];//注意是1.0,否则结果默认为0
        for(int i = 2; i <= n; i ++)
        {
            for(int j = 0; j <= m; j ++)
            {
                for(int k = 0; k <= num[i]; k ++)
                    if(k + j <= m)
                    n2[j + k] += n1[j] / fac[k];
            }
            //strcpy(n1, n2, sizeof(n1));
            for(int i = 0; i <= m; i ++)
                n1[i] = n2[i];
            memset(n2, 0, sizeof(n2));
        }
        printf("%.0lf\n",n1[m]*fac[m]);
    }
    return 0;
}