母函数
母函数又称生成函数。定义是给出序列:a0,a1,a2,.......ak,......,那么函数G(x)=a0+a1*x+a2*x2+......ak*xk称为序列a0,a1,a2,.......ak,......的母函数(即生成函数)。
例如：序列1,2,3.......n的生成函数为：G(x)=x+2x2+3x3+........nxn。点此链接:百度百科
特别的当序列为:1，1，1，1，.......1，这个生成函数为:G(x)=x+x2+x3+.......+xn=(1-xn)/(1-x),当-1<x<1时G(x)=1/(1-x)
1/(1-x)n=1+C(n,1)x+C(n+1,2)x2+C(n+2,3)x3+...+C(n+k-1,k)xk+...可以把生成函数还原为数列。
==========================================================================================================================
例 1：使用母函数求出斐波那契数列的通项公式。Fib(n)=Fib(n-1)+Fib(n-2),这里假设Fib(1)=1,Fib(2)=1;
求解这种递推关系的方法是：①、将递推关系变成母函数方程；②、求解母函数方程；③、将母函数变成幂级数形式。
所以斐波那契数列的生成函数为:G(x)=x+x2+2x3+3x4+5x5+8x6..........。
等式两边同时*x有：xG(x)=x2+x3+2x4+3x5+5x6+8x7+.......。
相加有：G(x)+xG(x)=x+2x2+3x3+5x4+8x5+13x6+........。
我们对比G(x)可以得到:G(x)+xG(x)=G(x)/x-1;所以我们可以得到:G(x)=x/(1-x-x2)。
可以令:1-x-x2=0,得到两根为:a=(1-√5)/2,b=(1+√5)/2,所以我们可以知道:1-x-x2=(x-x1)(x-x2)=(1-ax)(1-bx);
假设x/(1-x-x2)=m/(1-ax)+n/(1-bx),通分有：x=m(1-bx)+n(1-ax).由系数关系可得m=-1/√5,n=1/√5,所以G(x)=-1/√5(1-bx)+1/√5(1-ax)。
我们可知：1/(1-bx)=1/[1-(1+√5)/2x]是以公比为(1+√5)/2的等比数列，1/(1-ax)是以公比为(1-√5)/2的等比数列，所以其通项公式为：Fib(n)=1/√5[bn+1-an+1]。
==========================================================================================================================
例题2：若有1克、2克、3克、4克的砝码各一 枚，能称出哪几种重量？各有几种可能方案？
构造母函数，如果用x的指数表示称出的重量，则：
1个1克的砝码可以用函数1+x表示，(前面的这个1表示1克的砝码个数为0)
1个2克的砝码可以用函数1+x2表示，
1个3克的砝码可以用函数1+x3表示，
1个4克的砝码可以用函数1+x4表示，
那么几种砝码的组合情况的用乘积表示有:(1+x)(1+x2)(1+x3)(1+x4)=1+x+x2+2x3+2x4+2x5+2x6+2x7+x8+x9+x10 ,系数即为方案数。
例称出重量为6的物品：①、1，2，3；②、2，4两种方案。
==========================================================================================================================
例题3：求用1分、2分、3分的邮票贴出不同数值的方案数？
这个相对于上面的那个例子是：这个邮票可以重复。可知其生成函数为：G(x)=(1+x+x2+....)(1+x2+x4+....)(1+x3+x6+...),同理展开后其系数即为方案数。
==========================================================================================================================
例题4：德.梅其里亚克称重问题
(1)重为a1,a2,a3.....ak的砝码，如何放在天平的两端，记可称重量为n的物体的不同方式为Cn,则Cn的母函数为:
G(x)=(x-a1+1+xa1)(x-a2+1+xa2).........(x-ak+1+xak) ------ x-a1表示砝码a1和物体放在同一个托盘内,xa1表示砝码和物体放在不同的托盘内,1则为不用这个砝码。
(2)重为a1,a2,a3....ak的砝码，如只可以放在天平的一端，记可称重量为n的物体的不同方式为Cn，则Cn的母函数为：
G(x)=(1+xa1)(1+xa2).........(1+xak)
==========================================================================================================================
例题5：数的划分，将整数分解为若干个整数(相当于将n个苹果放在n个无区别的盘子里，每个盘子可以放多个，也可以不放)，上一篇博文中有提到。
假设1出现的次数为记为a1,2出现的次数记为a2.........k出现的次数记为ak,那么生成函数为：
G(x)=(1+x+x2+x3+x4+.....)(1+x2+x4+x6+x8+......)(1+x3+x6+x9+....)........(1+xn)
前面的1+x2+x4+x6+x8+......意思是当出现一个2时为x2，当出现两个2时为x4.....,为什么当出现n时，只有两项(1+xn)，因为是将数n划分为若干项，所以不能超过该数，且由数1到n项数依次要<=n/k(k=1.2,3,4...n)。
还是以nyist 90(数的划分)为例：这里就直接套用网上的模板了
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAX=50;
#define CLR(arr,val) memset(arr,val,sizeof(arr))
int n,m,value[MAX],temp[MAX];
int main()
{ cin>>m;
while(m--)
{ cin>>n;
fill(value,value+MAX,1);//value用来存储系数
CLR(temp,0);//temp用来保存每一次的情况
for(int i=2;i<=n;i++)
{ for(int j=0;j<=n;j++)
for(int k=0;k+j<=n;k+=i) //控制每次系数的变化和每个数出现的最大项数
temp[k+j]+=value[j];
for(int j=0;j<=n;j++)
value[j]=temp[j],temp[j]=0;
}
cout<<value[n]<<endl;
}
return 0;
}
==========================================================================================================================
例题2：hdu 1085(硬币问题)
//有3种面额是1、2、5的硬币,输入3个数字代表每种硬币的枚数,求最小的不能由这些硬币组成的面额是多少？
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAX=8010;
#define CLR(arr,val) memset(arr,val,sizeof(arr))
int value[MAX],temp[MAX],num[3],coin[3]={1,2,5};
int main()
{ while(cin>>num[0]>>num[1]>>num[2])
{ if(num[0]+num[1]+num[2]==0) break;
int max=num[0]+2*num[1]+5*num[2];
CLR(value,0);
CLR(temp,0);
fill(value,value+num[0]+1,1);
for(int i=1;i<3;i++)
{ for(int j=0;j<=max;j++)
for(int k=0;k+j<=max&&k/coin[i]<=num[i];k+=coin[i])//注意不能超出个数
temp[k+j]+=value[j];
for(j=0;j<=max;j++)
value[j]=temp[j],temp[j]=0;
}
for(i=0;i<=max+1;i++)//遍历即可
if(value[i]==0) {cout<<i<<endl;break;}
}
return 0;
}
涉及到母函数的题目有：HDU 1171，1398，1709，2065，2069，2082，2152；POJ 3046,3716,3734等等~有空再做下，还有很多东西不懂，应该先归纳下其它的知识~~