【学习笔记】概率生成函数

合集 - 学习笔记(42)

1.【学习笔记】Kruskal 重构树2022-11-25 2.【学习笔记】网络流2022-12-15 3.【学习笔记】高级数据结构2022-12-21 4.【学习笔记】线性基2022-12-27 5.【学习笔记】Link Cut Tree2022-12-28 6.【学习笔记】字符串后缀算法2022-12-30 7.【学习笔记】字符串回文算法2023-01-12 8.【学习笔记】组合数学2023-01-27 9.【学习笔记】多项式 1：基础操作2023-01-30 10.【学习笔记】多项式 2：集合幂级数2023-02-03 11.【学习笔记】多项式 3：多项式运算2023-02-11 12.【学习笔记】Prufer 序列2023-02-15 13.【学习笔记】多项式 4：生成函数2023-02-15 14.【学习笔记】DP 套 DP2023-02-19 15.【学习笔记】图的连通性2023-02-22 16.【学习笔记】差分约束2023-03-20 17.【学习笔记】长链剖分2023-04-22 18.【学习笔记】2-SAT2023-04-23 19.【学习笔记】根号算法2023-06-06 20.【学习笔记】Primal-Dual 原始对偶算法2023-06-15 21.【学习笔记】Bostan-Mori 算法2023-07-02 22.【学习笔记】狄利克雷卷积与高级筛法2023-07-03 23.【学习笔记】DP 优化 1：基础优化2023-07-03 24.【学习笔记】DP 优化 2：动态 DP2023-07-04 25.【学习笔记】李超线段树2023-07-09 26.【学习笔记】优化建图2023-07-11 27.【学习笔记】Segment Tree Beats2023-07-11 28.【学习笔记】插头 DP2023-07-13 29.【学习笔记】任意模数多项式乘法2023-07-16 30.【学习笔记】SG 函数与 SG 定理2023-08-05 31.【学习笔记】类欧几里得算法2023-08-06 32.【学习笔记】狄利克雷前/后缀和/差分2023-08-11 33.【学习笔记】DSU on Tree2023-08-22 34.【学习笔记】DP 优化 3：闵可夫斯基和优化 DP2023-09-01 35.【学习笔记】笛卡尔树2023-09-05 36.【学习笔记】Miller-Rabin 算法2023-10-24 37.【学习笔记】DP 优化 4：决策单调性2023-11-15 38.【学习笔记】DP 优化 5：wqs 二分优化 DP2023-12-04 39.【学习笔记】边分治2023-12-28 40.【学习笔记】KMP 相关算法2024-01-10

41.【学习笔记】概率生成函数2024-01-11

42.【学习笔记】离散对数和剩余2024-01-25

概述#

用生成函数刻画一些困难的概率期望问题，使用一些朴素的数学技巧来解出答案。

设 $F(x)$ 为概率生成函数，定义为：

F (x) = \sum_{i \geq 0} P (X = i) x^{i}

$F(x)=\sum_{i\ge 0}P(X=i)x^i$

容易发现 $F(1)=1$ 。

将 $F(x)$ 求导得到：

F^{'} (x) = \sum_{i \geq 0} i P (X = i) x^{i - 1}

$F'(x)=\sum_{i\ge 0}iP(X=i)x^{i-1}$

容易发现 $E(X)=F'(1)=1$ 。

同时根据高中数学知识可以得到：

\begin{aligned} V a r (X) & = \sum_{i \geq 0} (i - E (X))^{2} P (X = i) \\ = \sum_{i \geq 0} i^{2} P (X = i) - 2 i E (X) P (X = i) + E (X)^{2} P (X = i) \\ = \sum_{i \geq 0} i^{2} P (X = i) - E (X)^{2} \end{aligned}

$\begin{aligned} Var(X)&=\sum_{i\ge 0}(i-E(X))^2P(X=i)\\ &=\sum_{i\ge 0}i^2P(X=i)-2iE(X)P(X=i)+E(X)^2P(X=i)\\ &=\sum_{i\ge 0}i^2P(X=i)-E(X)^2 \end{aligned}$

而 $F(x)$ 的二阶导：

F^{″} (x) = \sum_{i \geq 0} i (i - 1) P (X = i) x^{i - 2}

$F''(x)=\sum_{i\ge 0}i(i-1)P(X=i)x^{i-2}$

于是 $Var(X)=F'(1)+F''(1)-(F'(1))^2$ 。

例题#

Luogu-P4548 CTSC 2006 歌唱王国 #

设 $f_i$ 表示在长度为 $i$ 时完成的概率， $g_i$ 表示在长度为 $i$ 时尚未完成的概率，记 $F(x),G(x)$ 为 $f,g$ 的生成函数，结合常数项 $f_0=0,g_0=1$ 可以得到：

F (x) + G (x) = x G (x) + 1

$F(x)+G(x)=xG(x)+1$

考虑将一个未完成的串后面直接接上目标串，此时可能出现没有完全接完就已经完成了的情况（原来就有一部分前缀），已经接上的部分既是前缀也是后缀，即 $\mathrm{border}$ ，剩下的部分就浪费了，所以可以写成：

\frac{x^{L} G (x)}{m^{| s |}} = \sum_{i = 1}^{| s |} \frac{x^{| s | - i} F (x)}{m^{| s | - i}} [s [1, i] \in B o r d e r (s)]

$\dfrac{x^LG(x)}{m^{|s|}}=\sum_{i=1}^{|s|} \dfrac{x^{|s|-i}F(x)}{m^{|s|-i}}[s[1,i]\in\mathrm{Border}(s)]$

因为要求 $F'(1)$ ，所以对第一个式子求导，得到：

F^{'} (x) + G^{'} (x) = G (x) + G^{'} (x)

$F'(x)+G'(x)=G(x)+G'(x)$

所以 $F'(1)=G(1)$ ，再把 $x=1$ 代入第二个式子。

G (1) = \sum_{i = 1}^{| s |} m^{i} [s [1, i] \in B o r d e r (s)]

$G(1)=\sum_{i=1}^{|s|}m^i[s[1,i]\in\mathrm{Border}(s)]$

所以答案就是 $m$ 的所有 $\mathrm{border}$ 次幂之和。

可见 PGF 解决问题的流程大致是设计一些数列的概率生成函数或普通生成函数，联立方程，求导，并将 $x=1$ 代入解决问题，上面部分已经阐明 $x=1$ 时的值有许多特殊之处。

P.S：伟大 APJ 指出， $F(x)$ 满足 $F(1)=1$ ，但是 $G(x)$ 不满足，所以只有 $F(x)$ 是概率生成函数。

参考资料#

《浅谈生成函数在掷骰子问题上的应用》 - 杨懋龙

作者：SoyTony

出处：https://www.cnblogs.com/SoyTony/p/17958685/Learning_Notes_about_Probability_Generating_Function

版权：本作品采用「署名-非商业性使用-相同方式共享 4.0 国际」许可协议进行许可。

posted @ 2024-01-11 15:29 SoyTony 阅读(258) 评论(0) 编辑收藏举报

刷新页面返回顶部

登录后才能查看或发表评论，立即登录或者逛逛博客园首页

相关博文：

· 【学习笔记】多项式 4：生成函数

· 【题解】Luogu-P5643 PKUWC 2018 随机游走

· 概率生成函数学习

· 概率生成函数（PGF）简记

· 2023.9.30 闲话

阅读排行：
· 单线程的Redis速度为什么快？
· 展开说说关于C#中ORM框架的用法！
· Pantheons：用 TypeScript 打造主流大模型对话的一站式集成库
· SQL Server 2025 AI相关能力初探
· 为什么退出登录或修改密码无法使 token 失效

公告

昵称： SoyTony
园龄： 3年2个月
粉丝： 144
关注： 92

+加关注

一言

私はあなたを爱して、私が地球の上のいかなる一人に比べてすべてあなたを爱するためです。

——『名探偵コナン』

2025年3月

日

一

二

三

四

五

六

Loading

SoyTony

【学习笔记】概率生成函数

概述#

例题#

Luogu-P4548 CTSC 2006 歌唱王国 #

参考资料#

公告

一言

搜索

合集 (8)

随笔分类 (863)

阅读排行榜

评论排行榜

推荐排行榜

Loading

SoyTony

【学习笔记】概率生成函数

概述#

例题#

Luogu-P4548 CTSC 2006 歌唱王国#

参考资料#

公告

一言

搜索

合集 (8)

随笔分类 (863)

阅读排行榜

评论排行榜

推荐排行榜

Luogu-P4548 CTSC 2006 歌唱王国 #