首页 > 其他分享 >生成函数的一个科技——指数族

生成函数的一个科技——指数族

时间:2023-01-27 18:55:39浏览次数:31  
标签:frac 函数 text sum 生成 科技 mathcal Hand Card

浅谈指数族 \(\text{Exponential Family}\)

本文参考了 《Generating Functionology》Chapter 3 的内容

0.引入

某天,有人(不是我)问了如下问题(不用考虑任何置换),

  • 如何求第一类、第二类斯特林数的生成函数?
  • 如何求大小为 \(n\) 连通标号图的生成函数?
  • 如何求大小为 \(n\) 的二分图数量的生成函数?
  • ...

学习了 \(\text{Exponential Family}\) ,我们可以轻松解决上述问题。

1.一些定义

  • 一个 \(\text{Card} \space C(S,p)\) 是一个包含有限集 \(S \subset \N_+\) 和一个描述 \(S\) 的图片(包括图)的二元组。

    它的权值为 \(|S|\) ,若 \(S = [n]\) ,则该 \(\text{Card}\) 被称为标准的

  • 一个 $\text{Hand} \space $ 为一些 \(\text{Card}\) 的集合,满足对于给定的 \(n\) ,有 \(\sum |S_i| = n\) 且 \(\cup S_i = [n]\) ,即它们可以恰好凑成 \([n]\) 。

    它的权值为 \(\sum |S_i|\) 。

  • 一个 \(\text{Deck}\) 是若干权值相同的**标准 \(\text{Card}\) **集合,记 \(d_i= |D_i|\) ,

    \(D(x)\) 为 \(\{d_i\}\) 的 \(EGF\) 。

  • 一个指数族(乱翻译的)\(\mathcal{F}\) 则包含了 \(D_1,D_2,...\) 。

\(h(n,k)\) 表示权值为 \(n\) 且恰好有\(k\) 个数 \(\text{Card}\) 的 \(\text{Hand}\) 个数,

令 $$ \large \mathcal{H}(x,y)=\sum_{n,k \geq 0}h(n,k)\frac{xn}{n!}yk $$ 。

如果省略 \(y\) 或 \(k\) ,例如 \(h(n)\) 、\(\mathcal{H}(x)\) 表示所有合法的 \(y\) 之和。

2. 一些性质

The exponential formula

对于一个 \(\mathcal{F}\) ,有

\[ \large \mathcal{H}(x,y)= \exp\{yD(x)\} $$. 为了证明上述定理,首先我们介绍如下引理: #### The Fundamental Lemma of Labeled Counting 令 $\mathcal{F^{\\'}} ,F^{\\'\\'}$ 为两个指数族,且它们的 $\text{Hand}$ 互不相同,$\mathcal{F} = \mathcal{F^{\\'}} \oplus F^{\\'\\'}$ 为它们的合并,则有 $\mathcal{H}(x,y)=\mathcal{H}^{\\'}(x,y)\mathcal{H}^{\\''}(x,y)$。 证明: 注意到 $ \large h(n,k)=\sum_{n_1,k_1 \geq 0} \binom{n}{n_1} h^{\\'}(n_1,k_1)h^{\\''}(n-n_1,k-k_1)$ 其意义为选定 $n$ 的一个 $n_1$ 大小子集,然后再从其中任选 $k_1$ 个 $\text{Card}$ 。 而右式 $ = [\frac{x^n}{n!}y^k] \mathcal{H}^{\\'}(x,y)\mathcal{H}^{\\''}(x,y) $ . $\square$ 接下来我们考虑仅有一个 $D_i$ 非零,的情况,则 $h(n,k)$ 只能在 $ h(ki,k) $ 处有值。 不妨令 $n = ik$,由于所有 $\text{Card}$ 的都是从 $D_i$ 当中选择,因此最后需要除以 $k!$ 。 也就是说: $$\large h(ki,k)=d_i^k\frac{n!}{i!^k}\frac{1}{k!}$$ , 从而有 \]

\mathcal{H}(x,y) &= \sum_{n,k \geq 0}h(n,k)\frac{xn}{n!}yk \
&= \sum_k h(n,k) \frac{x{ik}}{n!}yk \
&= \sum_k \frac{d_ikxyk}{k!(i!)k} = \exp{\frac{y\times d_ix^i}{i!}}

\[ 然后,我们再考虑将不同的 $D_i$ 合并在一起,不妨认为仅由 $D_i$ 产生的 $\text{Hand}$ 的生成函数记作 $\mathcal{H}_i(x,y)$ ,有引理知,我们有: \]

\mathcal{H}(x,y)=\prod_{i\geq1}H_i(x,y)= \exp{y\sum_{i \geq 1} \frac{d_ix^i}{i!}}=\exp

\[$\square$ ## 3. 应用 #### 求第一类斯特林数 容易发现大小为 $i$ 的圆排列 $d_i = (i-1)!$ ,从而 $D(x)=\sum_{i\geq1}\frac{x^n}{n}=\log\frac{1}{1-x}$, 因此 $\mathcal{H}(x,y)=e^{yD(x)}=\frac{1}{(1-x)^y}$。 #### 求第二类斯特林数 大小为 $i$ 的子集只有 $1$ 个,因此 $D(x) = e^x-1$ , 从而 $\mathcal{H}(x,y)=e^{y(e^x-1)}$。 > $ \mathcal{H} (x)=e^{e^x-1}$ 其实就是著名的贝尔数的 $EGF$。 #### 求连通图个数 如果一个 $\text{Hand}$ 表示一个图, 那么一个 $\text{Deck}$ 就表示一个连通图。 而显然,$\mathcal{H}(x)=\sum_{n\geq0}2^{\binom{n}{2}}\frac{x^n}{n!}$, 因此 $D(x)=\log \mathcal{H}(x)$, #### 二分图数量 首先考虑如果黑白两色可以区分, 则 $h_i=\sum_{k\geq0}\binom{n}{k}2^{k(n-k)}$ , 从而有 $D(x)=\log \mathcal{H}_0(x)$ , 但是不考虑颜色,我们会有 $\mathcal{H}(x)=\exp \frac{D(x)}{2}=\sqrt{\mathcal{H}(x)}$。 > 思考:如何计算每个点度数都为 2 的图的个数? ##\]

标签:frac,函数,text,sum,生成,科技,mathcal,Hand,Card
From: https://www.cnblogs.com/LitDarkness/p/17069167.html

相关文章

  • 数据分析自动化 生成报告可视化
    数据分析是指用适当的统计分析方法对收集来的大量数据进行分析,将它们加以汇总和理解并消化,以求最大化地开发数据的功能,发挥数据的作用。数据分析过程自动化随着信息化......
  • hdu:"红色病毒"问题(指数型母函数用e^x指数函数来计算)
    ProblemDescription医学界发现的新病毒因其蔓延速度和Internet上传播的”红色病毒”不相上下,被称为”红色病毒”,经研究发现,该病毒及其变种的DNA的一条单链中,胞嘧啶,......
  • 递归函数的全局变量使用技巧
    递归函数的全局变量使用技巧我希望提取以下数组中每个path的值放入一个数组letarr=[{path:'a',b:2,children:[{......
  • 为什么pollard rho算法要用x^2+c生成随机序列
    \(f(x)=x^2+c\)满足:\(\forallx\equivy\pmodn\),\(f(x)\equivf(y)\pmodn\).用\(f(x)\)生成一列数\(a_0=x,a_n=f^{(n)}(x)=f(a_{n-1})\).这意味着对要分解......
  • 生成函数推导组合恒等式
    上接https://www.cnblogs.com/juruo-zzt/p/15369446.html可能循环论证了!范德蒙德卷积\[\sum_i\binomni\binomm{m-i}=\binom{n+m}{n}\]\[[x^n](x+1)^{n}(x+1)^m=[x^n......
  • MySQL基础篇(运算符、排序分页、多表查询、函数)
    MySQL基础篇​​数据库概述​​​​数据库与数据库管理系统​​​​数据库与数据库管理系统的关系​​​​Mysql介绍​​​​RDBMS与非RDBMS​​​​关系型数据库(RDBMS)......
  • 读Java8函数式编程笔记02_流
    1. 外部迭代1.1. for循环是一个封装了迭代的语法糖1.1.1. 本质上来讲是一种串行化操作1.2. 很难抽象出不同操作2. 内部迭代2.1. 内部迭代中的相应接口:Stream......
  • 利用python函数调用ffmpeg批量进行转码
    本人学习python没几天,代码也没记住,写个函数到处查笔记,东拼西凑的。累……但是最终还是搞定了。欢迎高手指导,谢谢!单个文件转码#学会如何在python调用bat文件importos,......
  • 函数递归
    目录函数的递归调用介绍回溯与递推函数递归调用介绍函数不仅可以嵌套定义,还可以嵌套调用,即在调用一个函数的过程中,函数内部又调用另一个函数,而函数的递归调用指的是在......
  • 根据swagger生成无侵入性的文档
    所需工具:1、swagger生成的json文件2、工具网站:https://editor.swagger.io/操作步骤:1、将swagger.json文件的内容粘贴到工具网站左侧。2、从上方的工具栏依次选择“Gen......