组合数学之母函数问题

母函数问题是组合数学中非常经典的问题,大概是本科二年级的课程,很有意思的一门课,当然也是很精深的一门课。

定义

对于序列a0,a1,a2,…构造函数G(x):

则称函数G(x)是序列a0,a1,a2,…的母函数。

很明显,根据二项展开式,很容易知道(1+x)^n是序列C(n,0),C(n,1),…,C(n,n)的母函数。如果已知序列a0,a1,a2,…则对应的母函数G(x)便可根据定义给出。反之,如若已经求得序列的母函数G(x),则该序列也随之确定。 将序列a0,a1,a2,…记为{an}

Example

  • 砝码问题

    经典的母函数问题就是砝码问题,即若有1克、2克、3克、4克的砝码各一枚,能称出哪几种重量?各有几种可能方案?

    利用母函数,可以轻松求解该问题。

    考虑构造母函数,用x的指数表示称出的质量,则:

    1个1克的砝码可以用函数1+x表示,

    1个2克的砝码可以用函数1+x^2表示,

    1个3克的砝码可以用函数1+x^3表示,

    1个4克的砝码可以用函数1+x^4表示,

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

    (1+x)( 1+x^2)( 1+x^3)( 1+x^4),将该式子展开,x的指数表示可以称出的重量,其系数表示该种重量的方案数。

  • 邮票问题

    砝码问题指定了每种砝码的数量,还有一个经典的问题就是邮票问题,即求用1分、2分、3分的邮票贴出不同数值的方案数。

    与砝码问题所不同的是,邮票允许重复,所以构造的母函数也就不同,

    这一类问题,常常涉及到“整数拆分”的概念。所谓的“整数拆分”就是把一个整数分解成若干个整数的和,其实就相当于把n个无区别的球放到n个无标志的盒子,盒子允许空着,也允许放多个球。整数拆分成若干个整数的和,不同拆分法的总数叫做拆分数。

    以上述的邮票问题为例,将上式展开后,以x^4为例,其系数为4,即将4拆分成1、2、3的拆分数是4,即4=1+1+1+1=1+1+2=1+3=2+2.

编码实现

对于母函数这一类问题,其实求解的思路很清晰,最关键的问题就是如何模拟多项式的展开。

以邮票问题为例,编码实现。输入为一个整型数值n,输出用1分、2分、3分…n分的邮票组合出数值n的方案数。

/*
    Name:
    Copyright:
    Author: Lijiansong
    Date: 06/08/15 17:38
    Description:
    组合数学,母函数问题,求用1分、2分、3分、4分等的邮票贴出数值n的方案数。
*/
#include<iostream>
using namespace std;
const int num=10000;
// c1表示各项质量砝码可以组合的数目
// c2是中间量,表示每一次的情况
int c1[num+1],c2[num+1];
int main(){
    int n,i,j,k;
    while(cin>>n){//n表示待组合的数值
        for(i=0;i<=n;++i){
    //1+x+x^2+x^3+...+x^n初始化,将质量从0-n的所有砝码初始化为1
            c1[i]=1;
            c2[i]=0;
        }
        for(int i=2;i<=n;++i){//i从2遍历到n,表示第i个表达式
            for(j=0;j<=n;++j)
//j从0遍历到n,j表示在第i个表达式里的第j个变量
//k表示的是第j个指数,k每次递增i(因为第i个表达式的增量是i)
                    c2[j+k]+=c1[j];
                 }
                 for(j=0;j<=n;++j){
//将c2的值赋给c1,而把c2初始化为0,因为c2每次是从一个表达式中开始的
                    c1[j]=c2[j];
                    c2[j]=0;
                 }
        }
        cout<<c1[n]<<endl;
    }
    return 0;
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

时间: 2024-08-11 17:07:50

组合数学之母函数问题的相关文章

组合数学 随记

组合数学 1.母函数 母函数(生成函数): 生成函数有普通型生成函数和指数型生成函数两种(本题是普通型). 形式上,普通型母函数用于解决多重集的组合问题, 指数型母函数用于解决多重集的排列问题. 母函数还可以解决递归数列的通项问题(例如使用母函数解决斐波那契数列,Catalan数的通项公式). 1.普通母函数:  /*hdu2082*/ 构造母函数G(x), G(x) = a0 + a1*x + a2* + a3* +....+ an*,  则称G(x)是数列a0,a1…an的母函数. 通常普通

组合数学 - 母函数 --- 模板 + 详解

Square Coins Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 8341    Accepted Submission(s): 5674 Problem Description People in Silverland use square coins. Not only they have square shapes but

组合数学 - 母函数的变形 --- hdu 1171:Big Event in HDU

Big Event in HDU Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 24002    Accepted Submission(s): 8458 Problem Description Nowadays, we all know that Computer College is the biggest department

【组合数学】 03 - 母函数和递推关系

1. 母函数 1.1 母函数的定义 计数问题的结果一般可以表示为自然数集上的函数\(f(n)\),比如组合数\(\binom{n}{k}\)可以看成是关于\(k\)的函数.孤立的通项公式很难发现数值之间的内在联系,从而丢失了结果的整体性.本节介绍一下欧拉提出的母函数思想,它是计数问题的一个基本工具. 你一定知道组合数\(\binom{n}{k}\)其实是二项式\((1+x)^n\)的系数,换个角度想就是多项式\((1+x)^n\)中包含了整个数列\(\{\binom{n}{k}\}\).利用单个

组合数学 - 母函数 + 模板总结

// Memory Time // 1347K 0MS // by : Snarl_jsb // 2014-09-19-18.23 #include<algorithm> #include<cstdio> #include<cstring> #include<cstdlib> #include<iostream> #include<vector> #include<queue> #include<stack>

组合数学 - 母函数的运用 --- hdu 1709 :The Balance

The Balance Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 5706    Accepted Submission(s): 2311 Problem Description Now you are asked to measure a dose of medicine with a balance and a number o

组合数学 - 母函数的运用 --- 模板题

Holding Bin-Laden Captive! Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 15064    Accepted Submission(s): 6750 Problem Description We all know that Bin-Laden is a notorious terrorist, and he h

组合数学 - 母函数的运用 + 模板 --- hdu : 2082

找单词 Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 4093    Accepted Submission(s): 2933 Problem Description 假设有x1个字母A, x2个字母B,..... x26个字母Z,同时假设字母A的价值为1,字母B的价值为2,..... 字母Z的价值为26.那么,对于给定的字母,可以找到

组合数学 - 母函数 + 模板题 : 整数拆分问题

Ignatius and the Princess III Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 13129    Accepted Submission(s): 9285 Problem Description "Well, it seems the first problem is too easy. I will let