生成函数(母函数)

参考文章:
生成函数(母函数)——目前最全的讲解
《小学生都能看懂的生成函数从入门到升天教程》《生成函数全家桶》
Acwing 进阶课程–生成函数

引入

任意给定一个无限长的序列a0,a1....an....a_{0},a_{1}....a_{n}....a0,a1....an....
定义函数g(x)=a0x0+a1x1+a2x2+...+anxn+...g(x)=a_{0}x^{0}+a_{1}x^{1}+a_{2}x^{2}+...+a_{n}x^{n}+...g(x)=a0x0+a1x1+a2x2+...+anxn+...
这是一个无穷级数,一般设x的取值范围为:-1<x<1。显然这个无穷级数收敛,我们可以在x取向正无穷时得到该无穷级数的值
我们称这个函数g(x)为序列的生成函数
我们不难看出
这么讲太抽象了,看一些例题就知道了

例题:

例题1 :砝码问题:

有三个砝码,分别是1g,2g,3g,求用这三个砝码一共能凑出多少个不同重量,并求每个重量的组成方案

用生成函数来解决这个问题:
我们将这个问题拆分成几个独立的子问题,然后利用乘法原理计算。在本题中,利用乘法原理来考虑每个砝码的选择方式。
在生成函数中,乘法原理可以映射到多项式乘法
对于每个砝码有选和不选两个方式
1个1g砝码可以看作是f1=x0+x1f1=x^0+x^1f1=x0+x1,x0x^0x0表示不取(重量为0),x1x^1x1表示只取一个
1个2g砝码可以看作是f2=x0+x2f2=x^0+x^2f2=x0+x2,x0x^0x0表示不取,x2x^2x2表示只取两个
1个3g砝码可以看作是f3=x0+x3f3=x^0+x^3f3=x0+x3,x0x^0x0表示不取,x3x^3x3表示只取三个
那么生成函数就是每个方案数的乘积 :
g(x)=f1∗f2∗f3=(1+x1)∗(1+x2)∗(1+x3)=1+x+x2+2x3+x4+x5+x6g(x)=f1*f2*f3=(1+x^1)*(1+x^2)*(1+x^3)=1+x+x^2+2x^3+x^4+x^5+x^6g(x)=f1f2f3=(1+x1)(1+x2)(1+x3)=1+x+x2+2x3+x4+x5+x6
最终算式中的1+x+x2+2x3+x4+x5+x61+x+x^2+2x^3+x^4+x^5+x^61+x+x2+2x3+x4+x5+x6就说明了0g的方案是1,3g的方案是2,sss克的方案数就是axsax^saxs的系数a,指数表示重量,系数表示方案数

例题2:砝码升级版

有三个砝码,分别是1g(共两个),2g(无数个),3g(共一个),求凑出4g重量的组成方案数
生成函数:
f1=(1+x+x2)f1=(1+x+x^2)f1=(1+x+x2)
f2=(1+x2+x4+x6+......)f2=(1+x^2+x^4+x^6+......)f2=(1+x2+x4+x6+......)
f3=(1+x3)f3=(1+x^3)f3=(1+x3)
f=f1∗f2∗f3=(1+x+x2)∗(1+x2+x4+x6+......)∗(1+x3)f=f1*f2*f3=(1+x+x^2)*(1+x^2+x^4+x^6+......)*(1+x^3)f=f1f2f3=(1+x+x2)(1+x2+x4+x6+......)(1+x3)
凑出4g重量的组成方案数有:
1,x4,11 ,x^4,11,x4,1
x,1,x3x,1,x^3x,1,x3
x2,x2,1x^2,x^2,1x2,x2,1
一共有三个方案,x4x^4x4的系数是3

例题3:苹果选择

现在有n种苹果,每种苹果无限个,问选k个苹果的方案数

方法一(隔板法):

我们先不用生成函数,用组合数来做:
我们设第i种苹果选了xi个,那么有∑i=1nxi=k,xi>=0\sum_{i=1}^{n}x_{i}=k,x_{i}>=0i=1nxi=k,xi>=0
我们用隔板法来做,需要xi大于0x_{i}大于0xi0
我们令ai=xi+1a_{i}=x_{i}+1ai=xi+1
有:∑i=1nai=k+n,ai>0\sum_{i=1}^{n}a_{i}=k+n,a_{i}>0i=1nai=k+n,ai>0
这个我们可以看作是有k+n个小球,现在有n个桶,必须保证每个桶至少有一个球,方案数为:Ck+n−1n−1C_{k+n-1}^{n-1}Ck+n1n1
∑i=1nai=k+n,ai>0\sum_{i=1}^{n}a_{i}=k+n,a_{i}>0i=1nai=k+n,ai>0方程组解的个数就是方案数

隔板法讲解:
这是隔板法的模板,做法就是看作有n-1个隔板,k+n个球有k+n-1个间隙,将这个n-1个隔板插入到k+n-1个间隙中,这样就会分割出n个空间,每个空间至少有一个小球。方案数就是:k+n-1个间隙中选出n-1个位置放隔板,Ck+n−1n−1C_{k+n-1}^{n-1}Ck+n1n1

方法二(生成函数):

f1=(1+x1+x2+x3+....)=1−xn1−x=11−x,n取向无穷大f1=(1+x^1+x^2+x^3+....)=\frac{1-x^{n}}{1-x}=\frac{1}{1-x},n取向无穷大f1=(1+x1+x2+x3+....)=1x1xn=1x1,n
f2=11−xf2=\frac{1}{1-x}f2=1x1

f=f1∗f2∗...∗fn=1(1−x)n=a0+a1∗x+...+akk+..=∑k∞Ck+n−1n−1xkf=f1*f2*...*fn=\frac{1}{(1-x)^n}=a_{0}+a_{1}*x+...+a_{k}^k+..=\sum_{k}^{∞}C_{k+n-1}^{n-1}x^{k}f=f1f2...fn=(1x)n1=a0+a1x+...+akk+..=kCk+n1n1xk
我们刚才用隔板法求出每种方案数:
所以:ak=Ck+n−1n−1a_{k}=C_{k+n-1}^{n-1}ak=Ck+n1n1

总结:

生成函数大多用来解决有限或者无线物品的组合方案
生成函数⇔原问题分解为若干个子问题⇔求每个子问题的生成函数⇔所有子问题生成函数乘起来就是总方案数生成函数\Leftrightarrow原问题分解为若干个子问题\Leftrightarrow求每个子问题的生成函数\Leftrightarrow所有子问题生成函数乘起来就是总方案数

生成函数常用公式:

高等数学基本上都学过

普通型生成函数:

单下标序列普通型生成函数:

A(x)=∑n=0∞anxnA(x)=\sum_{n=0}^{∞}a_{n}x^nA(x)=n=0anxn

双下标序列普通型生成函数:

A(x1,x2)=∑n=0∞a(n,m)x1nx2mA(x_{1},x_{2})=\sum_{n=0}^{∞}a_{(n,m)}x_1^nx_2^mA(x1,x2)=n=0a(n,m)x1nx2m

普通型生成函数常用公式:

11−x=∑i=0∞xi\frac{1}{1-x}=\sum_{i=0}^{∞}x^i1x1=i=0xi

指数型生成函数:

二项式定理:

(x+y)n=∑i=1nCniaibn−i(x+y)^n=\sum_{i=1}^{n}C_{n}^ia^ib^{n-i}(x+y)n=i=1nCniaibni
(x+1)n=∑i=0nCnixi(x+1)^{n}=\sum_{i=0}^{n}C_{n}^ix^i(x+1)n=i=0nCnixi
1(1−x)n=∑i=0∞Cn+i−1ixi\frac{1}{(1-x)^n}=\sum_{i=0}^{∞}C_{n+i-1}^{i}x^{i}(1x)n1=i=0Cn+i1ixi

例题:

HDU1028
HDU1085
洛谷P2000
BZOJ3028

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/316138.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

AtCoder Beginner Contest 234

A - Weird Function B - Longest Segment C - Happy New Year! D - Prefix K-th Max E - Arithmetic Number F - Reordering G - Divide a Sequence 写个函数 int f(int x){return x*x2*x3;} int main() { int t;scanf("%d", &t);cout<<f(f(f(t)t)f(f(t…

分享一个.NET平台开源免费跨平台的大数据分析框架.NET for Apache Spark

今天早上六点半左右微信群里就看到张队发的关于.NET Spark大数据的链接https://devblogs.microsoft.com/dotnet/introducing-net-for-apache-spark/ &#xff0c;正印证了“微软在不断通过.NET Core补齐各领域开发&#xff0c;真正实现一种语言的跨平台”这句话。那么我们今天就…

Codeforces Round #760 (Div. 3)

E. Singers’ Tour F. Reverse G. Trader Problem 推推式子就行了。 int a[N]; int main() {int t;scanf("%d", &t);while(t --){int n;LL sum 0;scanf("%d", &n);for(int i 1;i < n;i ) scanf("%d", ai), sum a[i];a[0] a[n]…

acwing3132. 食物(BZOJ3028)

acwing3132. 食物 题意&#xff1a; 你当然要帮他计算携带 N 件物品的方案数。 承德汉堡&#xff1a;偶数个。 可乐&#xff1a;0 个或 1 个。 鸡腿&#xff1a;0 个&#xff0c;1 个或 2 个。 蜜桃多&#xff1a;奇数个。 鸡块&#xff1a;4 的倍数个。 包子&#xff1a;0 个…

持续畅销20年的《C#高级编程》出第11版了!

TA是谁&#xff1f;Wrox精品红皮书&#xff0c;引领无数程序员进入程序开发殿堂&#xff0c;C#专家级指南&#xff0c;是经验丰富的程序员提高效率的更快捷方式&#xff0c;连续畅销20年&#xff0c;累计销量超30万册。TA出生名门&#xff1a;TA战绩辉煌&#xff1a;2019新的征…

cfF. Boring Queries

cfF. Boring Queries 题意&#xff1a; n个数组a[]&#xff0c;q个询问&#xff0c;每次询问区间[l,r]的lcm值 题目要求强制在线 1<n<1e5 1<a<2e5 1<q<1e5 题解&#xff1a; 添加链接描述 添加链接描述 添加链接描述 我们一般求lcm都是直接通过ab/gcd(a…

Educational Codeforces Round 119 (Rated for Div. 2)

D. Exact Change E. Replace the Numbers G. Subsequences Galore 因为1和2的数量最大值不是很多&#xff0c;多了的话可以用3代替&#xff0c;那么枚举1和2的数量然后二分3的数量 int a[110], n; bitset<10> bit; bool ch(int x) {for(int i 1;i < n;i ){int num …

.NET微服务体系结构中为什么使用Ocelot实现API网关

为什么要使用API网关而不是直接通信&#xff1f;在微服务架构中&#xff0c;客户端应用程序通常需要使用来自多个微服务的功能。如果直接执行该消费&#xff0c;则客户端需要处理多个微服务端点以进行呼叫。当应用程序发展并引入新的微服务或更新现有的微服务时会发生什么&…

P2000 拯救世界

P2000 拯救世界 题意&#xff1a; 为了拯救世界&#xff0c;小 a 和 uim 决定召唤出 kkksc03 大神和 lzn 大神。根据古籍记载&#xff0c;召唤出任何一位大神&#xff0c;都需要使用金木水火土五种五行神石来摆一个特定的大阵。而在古籍中&#xff0c;记载是这样的&#xff1…

Codeforces Round #762 (Div. 3)

E. MEX and Increments F. Let’s Play the Hat? G. Unusual Minesweeper H. Permutation and Queries 用个优先队列模拟。 map<int,int>ma; priority_queue<int> q;int main() {int t;scanf("%d", &t);while(t --){int n;scanf("%d", …

基于Jenkins Pipeline的ASP.NET Core持续集成实践

最近在公司实践持续集成&#xff0c;使用到了Jenkins的Pipeline来提高团队基于ASP.NET Core API服务的集成与部署&#xff0c;因此这里总结一下。一、关于持续集成与Jenkins Pipeline1.1 持续集成相关概念互联网软件的开发和发布&#xff0c;已经形成了一套标准流程&#xff0c…

踩不出足迹(牛客练习赛88 )

踩不出足迹(牛客练习赛88 ) 题意&#xff1a; 长度为n的数组a&#xff0c;每个数是一个k位二进制 定义一下操作&#xff1a; 令第一次得到的结果为 a1a_1a1​。你需要从第二个数开始&#xff0c;每次可以选择与上一次得到的结果异或或者同或起来。 问最大值是多少&#xff1f…

Codeforces Round #766 (Div. 2)

D. Not Adding E. Not Escaping F. Not Splitting 直接枚举就行了&#xff0c;原本还想的是素倍数&#xff0c;但是素倍数也不行。 bool dis[N];int main() {int n, x;scanf("%d", &n); int ans -n;while(n --)scanf("%d", &x), dis[x] 1;for(…

编程语言之父谈语言设计,龟叔大赞TypeScript

争论哪门编程语言孰优孰劣&#xff0c;长期以来都是程序员乐此不疲的“娱乐活动”。之所以说是娱乐活动&#xff0c;因为这些争论到最后往往只是各自在发泄情绪&#xff0c;再则就是&#xff0c;脱离使用场景去讨论所谓哪门语言更好并没意义。但如果让编程语言作者坐在一起讨论…

P2656 采蘑菇

P2656 采蘑菇 题意&#xff1a; 有n个点&#xff0c;m个单向边&#xff0c;每个边都有边权&#xff0c;如果经过这个边&#xff0c;可以获得其边权&#xff0c;而其边权会变成原来的p倍(0.1<p<0.8)&#xff0c;向下取整 从s点出发&#xff0c;问最多可以采到的蘑菇 题…

Codeforces Round #764 (Div. 3)

A. Plus One on the Subset B. Make AP C. Division by Two and Permutation D. Palindromes Coloring E. Masha-forgetful F. Interacdive Problem G. MinOr Tree 就是最小值逐步增加到最大值的过程。 int main() {int t;scanf("%d", &t);for(int _ 1;_ <…

你必须知道的 SmartSql

介绍SmartSql MyBatis Cache(Memory | Redis) R/W Splitting Dynamic Repository Diagnostics ......简洁、高效、高性能、扩展性、监控、渐进式开发&#xff01;她是如何工作的&#xff1f;SmartSql 借鉴了 MyBatis 的思想&#xff0c;使用 XML 来管理 SQL &#xff0c;并…

CF785D Anton and School - 2

CF785D Anton and School - 2 题意&#xff1a; 给定一个长度≤210^5由(和)组成的字符串&#xff0c;问有多少个子串&#xff08;可以不连续&#xff09;&#xff0c;前半部分是由(组成后半部分由)组成. 题解&#xff1a; 怎么括号匹配能出这么多题 如何才能不重不漏的选出…

OsharpNS轻量级.net core快速开发框架简明入门教程

OsharpNS官方资源项目地址&#xff1a;https://github.com/i66soft/osharp-ns20演示地址&#xff1a;https://www.osharp.org 直接使用QQ登录可以查看效果文档地址&#xff1a;https://docs.osharp.org 正在完善中....发布博客&#xff1a;https://www.cnblogs.com/guomingfeng…

cf1491C. Pekora and Trampoline

cf1491C. Pekora and Trampoline 题意&#xff1a; 有n个蹦床&#xff0c;每个蹦床有它的弹力值bib_{i}bi​,从i蹦床起跳可以落到ibiib_{i}ibi​的位置上(前提是这个位置有蹦床)&#xff0c;跳完后&#xff0c;蹦床的弹力值会减1&#xff0c; 题解&#xff1a; 对于一个蹦床…