母函数

母函数又称生成函数。定义是给出序列:a0,a1,a2,.......ak,......,那么函数G(x)=a0+a1*x+a2*x²+......ak*x^k称为序列a0,a1,a2,.......ak,......的母函数(即生成函数)。

例如：序列1,2,3.......n的生成函数为：G(x)=x+2x²+3x³+........nxⁿ。点此链接:百度百科

特别的当序列为:1，1，1，1，.......1，这个生成函数为:G(x)=x+x²+x³+.......+xⁿ=(1-xⁿ)/(1-x),当-1<x<1时G(x)=1/(1-x)

1/(1-x)ⁿ=1+C(n,1)x+C(n+1,2)x²+C(n+2,3)x³+...+C(n+k-1,k)x^k+...可以把生成函数还原为数列。

==========================================================================================================================

例 1：使用母函数求出斐波那契数列的通项公式。Fib(n)=Fib(n-1)+Fib(n-2),这里假设Fib(1)=1,Fib(2)=1;

求解这种递推关系的方法是：①、将递推关系变成母函数方程；②、求解母函数方程；③、将母函数变成幂级数形式。

所以斐波那契数列的生成函数为:G(x)=x+x²+2x³+3x⁴+5x⁵+8x⁶..........。

等式两边同时*x有：xG(x)=x²+x³+2x⁴+3x⁵+5x⁶+8x⁷+.......。

相加有：G(x)+xG(x)=x+2x²+3x³+5x⁴+8x⁵+13x⁶+........。

我们对比G(x)可以得到:G(x)+xG(x)=G(x)/x-1;所以我们可以得到:G(x)=x/(1-x-x²)。

可以令:1-x-x²=0,得到两根为:a=(1-√5)/2,b=(1+√5)/2,所以我们可以知道:1-x-x²=(x-x1)(x-x2)=(1-ax)(1-bx);

假设x/(1-x-x²)=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[bⁿ⁺¹-aⁿ⁺¹]。

==========================================================================================================================

例题2：若有1克、2克、3克、4克的砝码各一 枚，能称出哪几种重量？各有几种可能方案？

构造母函数，如果用x的指数表示称出的重量，则：
1个1克的砝码可以用函数1+x表示，(前面的这个1表示1克的砝码个数为0)
1个2克的砝码可以用函数1+x²表示，
1个3克的砝码可以用函数1+x³表示，
1个4克的砝码可以用函数1+x⁴表示，

那么几种砝码的组合情况的用乘积表示有:(1+x)(1+x²)(1+x³)(1+x⁴)=1+x+x²+2x³+2x⁴+2x⁵+2x⁶+2x⁷+x⁸+x⁹+x¹⁰ ,系数即为方案数。

例称出重量为6的物品：①、1，2，3；②、2，4两种方案。

==========================================================================================================================

例题3：求用1分、2分、3分的邮票贴出不同数值的方案数？

这个相对于上面的那个例子是：这个邮票可以重复。可知其生成函数为：G(x)=(1+x+x²+....)(1+x²+x⁴+....)(1+x³+x⁶+...),同理展开后其系数即为方案数。

==========================================================================================================================

例题4：德.梅其里亚克称重问题

(1)重为a1,a2,a3.....ak的砝码，如何放在天平的两端，记可称重量为n的物体的不同方式为Cn,则Cn的母函数为:

G(x)=(x^-a1+1+x^a1)(x^-a2+1+x^a2).........(x^-ak+1+x^ak) ------ x^-a1表示砝码a1和物体放在同一个托盘内,x^a1表示砝码和物体放在不同的托盘内,1则为不用这个砝码。

(2)重为a1,a2,a3....ak的砝码，如只可以放在天平的一端，记可称重量为n的物体的不同方式为Cn，则Cn的母函数为：

G(x)=(1+x^a1)(1+x^a2).........(1+x^ak)

==========================================================================================================================

例题5：数的划分，将整数分解为若干个整数(相当于将n个苹果放在n个无区别的盘子里，每个盘子可以放多个，也可以不放)，上一篇博文中有提到。

假设1出现的次数为记为a1,2出现的次数记为a2.........k出现的次数记为ak,那么生成函数为：

G(x)=(1+x+x²+x³+x⁴+.....)(1+x²+x⁴+x⁶+x⁸+......)(1+x³+x⁶+x⁹+....)........(1+xⁿ)

前面的1+x²+x⁴+x⁶+x⁸+......意思是当出现一个2时为x²，当出现两个2时为x⁴.....,为什么当出现n时，只有两项(1+xⁿ)，因为是将数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等等~有空再做下，还有很多东西不懂，应该先归纳下其它的知识~~