（整数划分）（母函数）

容斥原理一种常见的思维方法，有时情况比较多，相互有重叠。为了使重叠部分不被重复计算，可以这样处理：先不考虑重叠的情况，把所有对象的数目的计算出来，然后减去重复重复计算的数目。计数方法为容斥原理。

整数划分问题

整数划分，不同划分方法的总数叫作划份数。
这道题是不超过该划分数，还有的题是就是求该划分数的情况。
考虑整数和划分数的大小关系
将从整数n划分为m分：
当n和m都等于1时，返回1；
当n等于m时，就不断地返回1+（n，n-1）
当m大的时候就只返回（n，n）
当n比较大的时候，分为划分数中一个数为m和划分数中每个数都小于m
。。。。。。。。。。。。。。。。。。。。。。。。。
但递归求解hdu1028TLE
最后一行程序递归翻倍，复杂度O(n²)
才有了之后的Dp。。。。。。。。。。。。。。。
1.递归求解(2ⁿ)：

#include <iostream>

using namespace std;
int part(int n,int m)//将整数n划分成最大数不超过m的划分数
{
    if(n==1||m==1) return 1;
    else if(n<m) return part(n,n);
    else if(n==m) return 1+part(n,n-1);
    else return part(n-m,m)+part(n,m-1);//导致超时
}
int main()
{
    int n;
    while(cin>>n)
    {
        cout<<part(n,n)<<endl;
    }
    return 0;
}

以上第三种情况，return part(n-m,m)+part(n,m-1);
（划分有一个数为m和划分有一个数小于m)

(递归改成递推式)
分类的情况一样的。。。。。。。。。
2.Dp划分(n²)：

const int N=200;
int dp[N+1][N+1];//dp[n][m]将整数n划分成最大数不超过m的划分数
void part()
{
    for(int i=1;i<=N;i++)
        for(int j=1;j<=N;j++)
        if((i==1)||(j==1)) dp[i][j]=1;
    else if(i<j) dp[i][j]=dp[i][i];
    else if(i==j) dp[i][j]=dp[i][j-1]+1;
    else dp[i][j]=dp[i][j-1]+dp[i-j][j];
}

母函数：
原理：把组合问题的加法和幂级数的乘幂对应起来
用母函数解决整数划分问题：

const int N=200;
int c1[N+1];
int c2[N+1];
void part()
{
    int i,j,k;
    for(i=0;i<=N;i++)//第一部分的系数都是1
    {
        c1[i]=1;c2[i]=0;
    }
    for(k=2;k<=N;k++)//第二部分开始展开
    {
        for(i=0;i<=N;i++)
            for(j=0;j+i<=N;j+=k;)
            c2[i+j]+=c2[i];
        for(i=0;i<=N;i++)//更新本次结果
        {
            c1[i]=c2[i];
            c2[i]=0;
        }
    }
}

这个是"普通型"母函数用于求组合数方案，还会有一种“指数型”母函数，用于求排列数
对组合数方案的题目，如果能用普通的母函数求解，一般也能用Dp求解，但众所周知Dp的难点在于递推关系，别人都说想不到就做不出来，母函数相对直观一点。
hdu1521（指数型母函数）
当同时选多个物品时，把公式相乘，其展开项为排列情况。分母阶乘的形式，是表示它体现了排列和组合的关系K!倍的关系。

#include <iostream>
#include <cstring>
using namespace std;
int f(int n)
{
    int i,s=1;
    for(i=1;i<=n;i++)
    {
        s=s*i;
    }
    return s;
}

int main()
{
    int n,m,i,j,k,num[12];
    double c1[12],c2[12];
    while(cin>>n>>m)
    {
        memset(c1,0,sizeof(c1));
        memset(c2,0,sizeof(c2));
        for(i=0;i<n;i++)
        {
            cin>>num[i];
        }
        for(i=0;i<=num[0];i++)
        {
            c1[i]=1.0/f(i);
        }
        for(i=1;i<n;i++)
        {
            for(j=0;j<=n;j++)
                for(k=0;k<=num[i]&&k+j<=n;k++)
                c2[k+j]+=c1[j]/f(k);
            for(j=0;j<=n;j++)
            {
                c1[j]=c2[j];
                c2[j]=0;
            }
        }
        cout<<(c1[m]*f(m))<<endl;
    }
    return 0;
}