母函数[通俗易懂]

2022-07-08 20:08:21 浏览数 (1)

大家好,又见面了,我是全栈君。

在数学中,某个序列的母函数(Generating function,又称生成函数)是一种形式幂级数。其每一项的系数能够提供关于这个序列的信息。使用母函数解决这个问题的方法称为母函数方法。

我们首先来看下这个多项式乘法:

由此能够看出:

1.x的系数是a1,a2,…an 的单个组合的全体。 2. x^2的系数是a1,a2,…a2的两个组合的全体。 ……… n. x^n的系数是a1,a2,….an的n个组合的全体(仅仅有1个)。

令a1,a2…an都等于1,由此可得

这里先给出两句话。不懂的能够等看完这篇文章再回过头来看:

1.“把组合问题的加法法则和幂级数的乘幂相应起来” 2.“母函数的思想非常easy — 就是把离散数列和幂级数一 一相应起来,把离散数列间的相互结合关系相应成为幂级数间的运算关系。最后由幂级数形式来确定离散数列的构造. ”

母函数的定义

对于序列a0。a1,a2,…构造一函数:

第一种:

有1克、2克、3克、4克的砝码各一枚,能称出哪几种重量?每种重量各有几种可能方案?

考虑用母函数来解决问题:

我们如果x表示砝码。x的指数表示砝码的重量。这样:

1个1克的砝码能够用函数1 1*x^1表示, 1个2克的砝码能够用函数1 1*x^2表示。 1个3克的砝码能够用函数1 1*x^3表示, 1个4克的砝码能够用函数1 1*x^4表示。

上面这四个式子懂吗?

我们拿1 x^2来说,前面已经说过,x表示砝码,x的指数表示砝码的重量!初始状态时,这里就是一个质量为2的砝码。

那么前面的1表示什么?依照上面的理解,1事实上应该写为:1*x^0,即1代表重量为2的砝码没有取。

所以这里1 1*x^2 = 1*x^0 1*x^2,即表示2克的砝码有两种状态。不取或取,不取则为1*x^0。取则为1*x^2

接着讨论上面的1 x^2。这里x前面的系数有什么意义?

这里的系数表示状态数(方案数)

1 x^2。也就是1*x^0 1*x^2,也就是上面说的不取2克砝码。此时有1种状态;或者取2克砝码,此时也有1种状态。(分析!)

所以,前面说的那句话的意义大家能够理解了吧?

几种砝码的组合能够称重的情况,能够用以上几个函数的乘积表示:

(1 x)(1 x^2)(1 x^3)(1 x^4)

=(1 x x^2 x^4)(1 x^3 ^4 x^7)

=1 x x^2 2*x^3 2*x^4 2*x^5 2*x^6 2*x^7 x^8 x^9 x^10

从上面的函数知道:可称出从1克到10克,系数便是方案数。(!

。!经典!!。)

比如右端有2^x^5 项,即称出5克的方案有2种:5=3 2=4 1。相同。6=1 2 3=4 2;10=1 2 3 4。

故称出6克的方案数有2种。称出10克的方案数有1种 。


接着上面,接下来是另外一种情况:

另外一种:

求用1分、2分、3分的邮票贴出不同数值的方案数:

大家把这样的情况和第一种比較有何差别?第一种每种是一个,而这里每种是无限的。

以展开后的x^4为例。其系数为4。即4拆分成1、2、3之和的拆分方案数为4;

即 :4=1 1 1 1=1 1 2=1 3=2 2

这里再引出两个概念”整数拆分“和”拆分数“:

所谓整数拆分即把整数分解成若干整数的和(相当于把n个无差别的球放到n个无标志的盒子。盒子同意空。也同意放多于一个球)。

整数拆分成若干整数的和,办法不一,不同拆分法的总数叫做拆分数。

代码模板:

代码语言:javascript复制
#include <iostream>  
using namespace std;  
const int _max = 10001;   
// c1是保存各项质量砝码能够组合的数目  
// c2是中间量。保存没一次的情况  
int c1[_max], c2[_max];     
int main()  
{      
    int nNum;     
    int i, j, k;  
   
    while(cin >> nNum)  
    {  
        for(i=0; i<=nNum;   i)   // ---- ①  
        {  
            c1[i] = 1;  
            c2[i] = 0;  
        }  
        for(i=2; i<=nNum;   i)   // ----- ②  
        {  
   
            for(j=0; j<=nNum;   j)   // ----- ③  
                for(k=0; k j<=nNum; k =i)  // ---- ④  
                {  
                    c2[j k]  = c1[j];  
                }  
            for(j=0; j<=nNum;   j)     // ---- ⑤  
            {  
                c1[j] = c2[j];  
                c2[j] = 0;  
            }  
        }  
        cout << c1[nNum] << endl;  
    }  
    return 0;  
} 

我们来解释下上面标志的各个地方:(***********!!

!重点!

!!***********)

① 、首先对c1初始化,由第一个表达式(1 x x^2 ..x^n)初始化,把质量从0到n的全部砝码都初始化为1. ② 、 i从2到n遍历。这里i就是指第i个表达式。上面给出的另外一种母函数关系式里。每个括号括起来的就是一个表达式。 ③、j 从0到n遍历,这里j就是(前面i個表达式累乘的表达式)里第j个变量,(这里感谢一下seagg朋友给我指出的错误。大家能够看下留言处的讨论)。如(1 x)(1 x^2)(1 x^3),j先指示的是1和x的系数,i=2运行完之后变为 (1 x x^2 x^3)(1 x^3)。这时候j应该指示的是合并后的第一个括号的四个变量的系数。 ④ 、 k表示的是第j个指数,所以k每次增i(由于第i个表达式的增量是i)。 ⑤ 、把c2的值赋给c1,而把c2初始化为0。由于c2每次是从一个表达式中開始的。

发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/115889.html原文链接:https://javaforall.cn

0 人点赞