母函数浅析

母函数

母函数可以用于数列,概率,组合数学中,是个极其强大的工具。

定义

对于一个数列 \(a_0,a_1,a_2,a_3\dots\) ,定义 \(G(x)=a_0+a_1x+a_2x^2+a_3x^3\dots\) 为其母函数(这里对是否无穷项并没有什么要求)。

看起来没什么用。

但是我们看一个例子:

已知母函数 \(G(x)=C^{0}_{n}+C^{1}_{n}x+C^{2}_{n}x^2+C^{3}_{n}x^3\dots\) 。

很明显,这是个二项式展开式,那我们不妨将其写成 \(G(x)=(x+1)^n\) 。这里母函数的作用就体现出来了。

我们可以将一个无穷式变成一个收敛形式,便于计算处理。

普通型母函数

首先来看看几个常见普通母函数:

  1. \(G(x)=\frac{1}{1-x}=1+x+x^2+x^3+\dots\)

    由等比数列求和公式可知 \(G(x)=\frac{a_1(1-q^n)}{1-q}\) ,这里的话 \(q=x,a_1=1\) 。由于 \(x\) 可以是任何数,我们不妨设其小于 \(1\) 。那么 \(q^n\) 趋近于 \(0\) ,就得到收敛形式。

  2. \(G(x)=\frac{1}{1-kx}=1+kx+k^2x^2+k^3x^3\dots\)

    证明同上。

  3. \(G(x)=\frac{1}{(1-x)^{n}}=C^{0}_{n+0-1}+C^{1}_{n+1-1}x+C^{2}_{n+2-1}x^2+C^{3}_{n+3-1}x^3\dots +C^{k}_{n+k-1}x^k\dots\)

    证明:

    数学归纳法。

    假设 \(n\) 的情况成立,那么对于 \(n+1\) :
    \[
    G(x)=\frac{1}{(1-x)^{n+1}}=\frac{1}{(1-x)^{n}}\frac{1}{(1-x)}\\bf 此时对于第\mit \;k\;\bf 项有:\\ \mit
    a_k=\sum^{k}_{i=1}C^{i}_{n+i-1}=C^{k}_{n+1+k-1}
    \]
    这里用到了 \(C^{n}_{m}=C^{n-1}_{m-1}+C^{n}_{m-1}\)

    最后的递归边界由于 \(C^{0}_{n-1}=C^{0}_{n+1-1}\) 所以就可以滚雪球上去。


虽然好像母函数主要是被用在组合上,但是用其来解数列递推式也非常香(初赛的递推数列题可以用母函数暴力推)。

这里以斐波那契数列为例:

? 已知 \(f_{n+1}=f_n+f_{n-1}\) 。

? 那么对于其母函数 \(T\) ,有:
\[
\begin{align}
T=f_0+f_1x+f_2x^2+f_3x^3+\dots\xT=f_0x+f_1x^2+f_2x^3+\dots\x^2T=f_0x^2+f_1x^3+\dots\\end{align}
\\therefore
T-f_0-f_1x=xT-f_0x+x^2T\\bf 解得: \mit \qquad T=\frac{1}{1-x-x^2}\\]
? 那么我们有了斐波那契数列母函数的收敛形式,只要将其向几个基本的母函数上靠,期望就能得到通项:

\[
\begin{align}
T=\frac{1}{-(\frac{\sqrt{5}-1}{2}-x)(\frac{-\sqrt{5}-1}{2}-x)}\\end{align}
\]

? 可以看出其中有黄金分割率,看起来离正解不远了,我们设其为 \(\phi_1\) 和 \(\phi_2\) 然后继续:

\[
T=\frac{1}{-\phi_1\phi_2(1-\frac{1}{\phi_1}x)(1-\frac{1}{\phi_2}x)}\=\frac{1}{(1-\frac{1}{\phi_1}x)(1-\frac{1}{\phi_2}x)}\\]
? 然后再裂个项:
\[
T=\frac{1}{\sqrt{5}} (-\frac{\phi_2}{1-\frac{1}{\phi_1}x}+\frac{\phi_1}{1-\frac{1}{\phi_2}x})
\]
? 根据几个基本母函数,斐波那契的通项为:
\[
F_n=\frac{1}{\sqrt{5}} [(\frac{1+\sqrt{5}}{2})^{n+1}-(\frac{1-\sqrt{5}}{2})^{n+1}]
\]

?



除了解递推数列通项,普通型母函数还可以解一些无重复项的组合问题,质数分解问题:

? 将一个数 \(n\) 分解成质数和的形式,有多少种本质不同的方法。

? 根据题意构建方程:
\[
G(x)=(1+x^2+x^4+x^6+\cdots)\times(1+x^3+x^6+\cdots)\times(1+x^5+x^{10}\cdots)\times\cdots
\]
? 化简后 \(n\) 次项的系数即为方法数。

用这种方法甚至可以解决如凑零钱之类的 \(dp\) 问题,根据上例可以发现,系数是要手动爆枚的,所以复杂度上来说并没有什么实际意义。

指数型母函数

基本形式如下:
\[
G(x)=\sum^{\infty}_{i=0}\frac{x^i}{i!}=1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+\frac{x^4}{4!}+\cdots
\]
可以注意到阶乘的形式很特别,易联想到 \(Taylor\) 展开 ,因为 \(x\) 取值并没有什么关系,我们就取 \(e^x\) 作为其收敛形式(与 \(e^x\) 在 \(0\) 处的展开式相同)。

那么同样,常见的指数型母函数有:

  1. \(G(x)=\frac{e^x+e^{-x}}{2}=\sum^{\infty}_{i=0}\frac{x^{2i}}{(2i)!}=1+\frac{x^2}{2!}+\frac{x^4}{4!}+\cdots\)
  2. \(G(x)=\frac{e^x-e^{-x}}{2}=\sum^{\infty}_{i=0}\frac{x^{2i+1}}{(2i+1)!}=x+\frac{x^3}{3!}+\cdots\)

皆为 \(Taylor\) 展开式。



指数型母函数主要是为求解有重复项的组合问题,比如 \(1,1,2,3\) 可以组合成的四位数个数,明显是 \(\frac{4!}{2!}\) 个,这就是为何下面要除阶乘。

注意,最终我们需要的是 \(\frac{x^n}{n!}\) 前的系数,可以发现,下面没有阶乘的项乘出来会自动加个阶乘,而有阶乘的会把重复部分除掉,刚好达到预计效果。

如果还不明白可以找个例子手动操作一遍。

指数型母函数大多可以通过 \(Taylor\) 展开找到通项,这是其优点。

如:

? 求 \(1,3,5,7,9\) 可以组成的 \(n\) 位数个数。其中 \(3,7\) 分别出现偶数次。

? 依照题意构建母函数:
\[
T=(1 + \frac{x^2}{2!} + \frac{x^4}{4!} + \cdots )^2\times(1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+\frac{x^4}{4!}+\cdots)^3\=(\frac{e^x+e^{-x}}{2})^2\times e^{3x}= \frac{1}{4}(e^{5x}+2e^{3x}+e^x)\=\frac{1}{4} \sum^{\infty}_{i=0} (5^i+2\times3^i+1)\frac{x^i}{i!}
\]
? 那么第 \(n\) 项的系数是: \(\frac{5^n+2\times3^n+1}{4}\) 。



\(\frak by\;\;thorn\_\)

原文地址:https://www.cnblogs.com/thornblog/p/12207208.html

时间: 2024-10-17 22:47:17

母函数浅析的相关文章

浅析母函数

母函数,又称生成函数. 母函数就是一列用来展示一串数字的挂衣架. --赫伯特·唯尔夫   . 母函数问题一般是,对于一个问题,抽象出一个生成函数. 以换硬币问题为例,有一元的,四元的,九元的........硬币无数个,问y元钱有几种组合方式. 设一个数的指数为y(组成y元钱),系数为有几种组合方式,开始抽象为数学算式: 对于1元钱可组合的钱为,S1= x^0+x^1+x^2+x^3+....+x^300(上限为300是因为题目说y不会大于300) 对于2元钱可组合的钱为, S2=x^0+x^2+

HDU 2189 悼念512汶川大地震遇难同胞——来生一起走(母函数或完全背包)

悼念512汶川大地震遇难同胞--来生一起走 Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 3773    Accepted Submission(s): 1913 Problem Description 妈妈你别哭泪光照亮不了我们的路让我们自己慢慢的走 妈妈我会记住你和爸爸的模样记住我们的约定来生一起走 上面这首诗节选自一位诗人纪念遇难

Python之encode与decode浅析

 Python之encode与decode浅析 在 python 源代码文件中,如果你有用到非ASCII字符,则需要在文件头部进行字符编码的声明,声明如下: # code: UTF-8 因为python 只检查 #.coding 和编码字符串,为了美观等原因可以如下写法: #-*-coding:utf-8-*- 常见编码介绍: GB2312编码:适用于汉字处理.汉字通信等系统之间的信息交换. GBK编码:是汉字编码标准之一,是在 GB2312-80 标准基础上的内码扩展规范,使用了双字节编码.

浅析PHP的开源产品二次开发的基本要求

浅析PHP的开源产品二次开发的基本要求 第一, 基本要求:HTML(必须要非常熟悉),PHP(能看懂代码,能写一些小系统,如:留言板,小型CMS),Mysql(至少会一种数据库),Javascript(能看懂,能改现成的一些代码),Div+Css(能进行界面的调整,明白CSS是怎么使用的) 第二, 熟悉开源产品的使用,比如 Dedecms,你要知道怎么登录,怎么新建栏目,怎么添加文章,模板标签的使用方法,模型的概念和使用方法等等一些功能 第三, 要熟悉这个开源产品的数据库结构,还要理解里面核心文

word-break|overflow-wrap|word-wrap——CSS英文断句浅析

---恢复内容开始--- word-break|overflow-wrap|word-wrap--CSS英文断句浅析 一 问题引入 今天在再次学习 overflow 属性的时候,查看效果时,看到如下结果,内容在 div 中国换行了,可是两个 P 元素的内容并没有换行,搜索一番没有找到系统的答案,截图到群里请教大神,才知道是英文断句的问题,但是还是不太明白.之前没有遇到这种情况,为了彻底搞清楚,英文断句,又开始学习英文断句到底是怎么回事. 二 换行 每种语言里都有换行,就中文而言,我们最小语言单位

浅析vanish

浅析 VANISH --一种cache 第一部分:理解vanish的准备工作 1.对CDN的小剖析 CDN  content  delivery  network  内容分发(推送)网络,是在现有的Internet中增加一层新的网络架构,将网络内容发布到最接近用户的网络边缘(边缘服务器),使用户最近取得所需内容,解决网络拥挤状态,提高用户访问网站的速度. CDN网络架构主要有两部分组成,中心和边缘两部分,中心指CDN网管中心和DNS重定向解析中心,负责全局负载均衡.边缘主要指异地节点,CDN分发

健康,home? [java的内存浅析]

健康,home? [java的内存浅析] 摘要: 原创出处: http://www.cnblogs.com/Alandre/ 泥沙砖瓦浆木匠 希望转载,保留摘要,谢谢! 乐观上上,how can other kno u,u r yourself!I must be strong and carry on. -泥沙砖瓦浆木匠 一.闲谈下 201407月记着那时候身体垮了下来,呵呵.想说,对自己的说,也是对大家的负责吧.那时候胸疼胸闷,然后几乎累垮了,我还坚持了一星期,那一星期真的迷迷糊糊.完全不能

Mysql查询优化器浅析

--Mysql查询优化器浅析 -----------------------------2014/06/11 1 定义 Mysql查询优化器的工作是为查询语句选择合适的执行路径.查询优化器的代码一般是经常变动的,这和存储引擎不太一样.因此,需要理解最新版本的查询优化器是如何组织的,请参考相应的源代码.整体而言,优化器有很多相同性,对mysql一个版本的优化器做到整体掌握,理解起mysql新版本以及其他数据库的优化器都是类似的. 优化器会对查询语句进行转化,转化等价的查询语句.举个例子,优化器会将

Volley框架源码浅析(一)

尊重原创http://blog.csdn.net/yuanzeyao/article/details/25837897 从今天开始,我打算为大家呈现关于Volley框架的源码分析的文章,Volley框架是Google在2013年发布的,主要用于实现频繁而且粒度比较细小的Http请求,在此之前Android中进行Http请求通常是使用HttpUrlConnection和HttpClient进行,但是使用起来非常麻烦,而且效率比较地下,我想谷歌正式基于此种原因发布了Volley框架,其实出了Voll