Pure Mathematics
Vol. 12  No. 12 ( 2022 ), Article ID: 59882 , 8 pages
10.12677/PM.2022.1212240

Faà di Bruno公式在恒等式及计数上的应用

沈萍

重庆师范大学数学与科学学院,重庆

收稿日期:2022年11月27日;录用日期:2022年12月23日;发布日期:2022年12月30日

摘要

Faà di Bruno’ formular的相关研究一直处于停滞状态,之后,由于Faà di Bruno在其它方面的应用,Lukacs开始将Faà di Bruno用于数学统计学;罗曼用umbral calculus的定理对Faà di Bruno方程进行了再论证,Constantine利用Faà di Bruno扩展了关于集合划分的恒等式,并给出了它的概率说明,Chu利用Faà di Bruno获得了一系列的行列式,Chou,Hsu和Shiue利用Faà di Bruno构造了一种具有互逆性的函数,由此推导了一组恒等方程,然而并没有有关用Fáa di Bruno’s公式求组合恒等式以及进行对称群及集合分割上的组合数的研究存在。本文使用Fáa di Bruno’s公式求组合恒等式以及进行对称群及集合分割上的组合数的研究,通过Faà di Bruno得到了各种著名组合数的恒等式,含括Catalan数,第一类及第二类Stirling数,q-二项式系数···等,及Faà di Bruno在计数上的应用;在第二节中,我们利用Faà di Bruno得到了多种组合数的恒等式;在第三节中,我们由第一类及第二类Stirling数、错排数、Bell数的指数生成函数用Faà di Bruno导出的恒等式获得这些数的组合意义,得到一些限制置换中圈结构和限制集合分割的子集大小的结果。

关键词

恒等式,Fáa di Bruno’s公式,高阶微分

Application of Faà di Bruno’s Formula in Identity and Counting

Ping Shen

School of Mathematics and Science, Chongqing Normal University, Chongqing

Received: Nov. 27th, 2022; accepted: Dec. 23rd, 2022; published: Dec. 30th, 2022

ABSTRACT

The research on Faà di Bruno has been in a stagnant state. Later, because of the application of Faà di Bruno in other aspects, Lukacs began to use Faà di Bruno in mathematical statistics; Roman demonstrated the Faà di Bruno equation again by using the theorem of umbral calculus. Constantine extended the identity of set partition by using Faà di Bruno and gave its probability explanation. Chu obtained a series of determinants by using Faà di Bruno, Chou; Hsu and Shiue used Faà di Bruno to construct a kind of reciprocal function, and derived a set of identity equations from it. However, there is no research on finding combinatorial identities by using Fáa di Bruno’s formula and on the number of combinations on symmetric groups and set partitions. Faà di Bruno's formula was used to find the combinatorial identity, and the combinatorial numbers on symmetric groups and set partitions were studied. In this paper, the identities of various famous combinatorial numbers, including Catalan numbers, the first and second Stirling numbers, q binomial coefficients, etc., and the application of Faà di Bruno in counting were obtained through Faà di Bruno; In the second section, we use Faà di Bruno to obtain the identities of multiple combinatorial numbers; in the third section, we obtain the combined meanings of Stirling numbers, staggered numbers and Bell numbers from the exponential generating functions of the first and second types of Stirling numbers with the identities derived from Faà di Bruno, and the subset size of set partitions can be obtained.

Keywords:Identical Equation, Faà di Bruno’s Formula, Higher Order Differential

Copyright © 2022 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

现代数学可以分为两大类:一类是研究连续对象的,如分析、方程等;另一类就是研究离散对象的组合数学。微积分以及现代数学的发展,为现代工业的发展打下了坚实的基石 [1]。而组合数学的发展,为本世纪电脑的变革打下了坚实的基石。从十七世纪到十八世纪,组合数学和随机运算的发展,到了十九世纪,它已经彻底地颠覆了传统的概念,产生了一种全新的、充满活力的新的数学派别 [2]。因此,组合数学的发展,使传统的数理学科成为了以解析与代数为主导的学科。尤其是随着计算机越来越多地渗透到各行各业,除了数字运算之外,也出现了许多快速发展的非数字运算方法,这些方法大多采用了合成方法 [3]。

在计算机科学、编码与密码学、物理等领域,都有着举足轻重的作用。化学、生物等学科都具有很大的实用价值,即计算科学 [4]。计算机的运算方法只是一个离散的物体,而不是一个数字,它的科学就在于它的结合,它给人一种感觉,就像它有思想一样,因此,它的计算方法就成为了计算机的中心 [5]。在商业经营、运输计划、作战指挥、财务分析等方面,美国有一间公司以组合数学为名;他们运用了综合数学的理论,以增加公司经营的收益,这个公司做的很好,另外,实验设计也是一个很实用的课题,其主要的数学理论是结合的。美国已经有一些专业公司,在此基础上,对其进行了研究 [6] [7]。近年来,一名德国知名的复合数学家运用复合数论对药品结进行了深入的探讨,为医药企业节约了成本。简而言之,组合数学的基础理论与实践,涉及计算机科学,生物学,化学;心理学,金融和基因工程学的最新的运用,其范围十分广阔。

本文首先研究使用Faà di Bruno公式求组合恒等式,然后进行对称群及集合分割上的组合数的研究,最后对论文进行总结,发觉Faà di Bruno公式的应用前景。由于本文研究偏理论研究,因此将不使用仿真例子进行研究,后续可在此基础上进行相关深入的研究。

2. 用Faà di Bruno求组合恒等式

对任意给定幂级数 ϕ ( t ) = n = 0 a n t n ,以 [ ϕ n ] 表示 t n 的系数,即 ϕ ( t ) = n 0 [ ϕ n ] t n

假设f有反函数,考虑复合函数 ( f ϕ ) ( t ) ( f 1 ( f ϕ ) ) ( t ) ,Chou,Hsu与Shiue将Faà di Bruno改写,得到下列两个式:

[ f ϕ n ] = σ ( n ) ( D k f ) x = ϕ ( 0 ) [ ϕ 1 ] k 1 [ ϕ 2 ] k 2 [ ϕ n ] k n k 1 ! k 2 ! k n ! (1)

[ ϕ n ] = σ ( n ) ( D k f 1 ) x = f ϕ ( 0 ) [ f ϕ 1 ] k 1 [ f ϕ 2 ] k 2 [ f ϕ n ] k n k 1 ! k 2 ! k n ! (2)

因为 [ f ϕ n ] [ ϕ n ] 在(1)和(2)中所处的地位恰好相反,称这两条恒等式为一组互逆关系式。特别地,考虑 f ( x ) = e x ϕ ( t ) = n = 0 a n t n ,令 b n = [ f ϕ n ] ,若 ϕ ( 0 ) = 0 ,则由(1)和(2)有以下互逆关系:

b n = σ ( n ) a 1 k 1 a 2 k 2 a n k n k 1 ! k 2 ! k n ! (3)

a n = σ ( n ) ( 1 ) k 1 ( k 1 ) ! b 1 k 1 b 2 k 2 b n k n k 1 ! k 2 ! k n ! (4)

定理1 对任意正整数n,有互逆关系

C n = 1 n + 1 ( 2 n n ) = 1 n + 1 σ ( n ) 2 2 n k 1 k 1 k 1 ! 2 k 2 k 2 ! n k n k n !

4 n = 2 n σ ( n ) ( 1 ) k 1 ( k 1 ) ! i = 1 n ( 2 i i ) k i k i !

定理2. 对任意正整数n,r,有互逆关系

1 n + 2 n + + r n = n σ ( n ) ( 1 ) k 1 ( k 1 ) ! i = 1 n S ( r + i , r ) k i k i !

S ( n + r , r ) = σ ( n ) i = 1 n S ( 1 i + 2 i + + r i ) k i i k i k i ! S ( n , r ) = σ ( n r ) i = 1 n r ( 1 i + 2 i + + r i ) k i i k i k i !

证明:令 f ( x ) = e x ϕ ( t ) = ln 1 ( 1 t ) ( 1 2 t ) ( 1 r t ) ,显然 ϕ ( 0 ) = 0

( f ϕ ( t ) ) = 1 ( 1 t ) ( 1 2 t ) ( 1 r t ) = i r S ( i , r ) t i r

故而 b i = [ f ϕ i ] = S ( r + i , r )

另一方面因为

ϕ ( t ) = ln 1 ( 1 t ) ( 1 2 t ) ( 1 r t ) = i = 1 r ( ln ( 1 i t ) ) = i = 1 r j 1 i j j t j = j 1 i = 1 r i j j t j

观察 t n 的系数可得

a i = [ ϕ n ] = i = 1 r i n n = 1 n + 2 n + + r n n

根据Faà di Bruno (4)可以得到

1 n + 2 n + + r n n = σ ( n ) ( 1 ) k 1 ( k 1 ) ! S ( r + 1 , r ) k 1 S ( r + 2 , r ) k 2 S ( r + n , r ) k n k 1 ! k 2 ! k n !

等号两边同乘n即得第一条恒等式;由Faà di Bruno的(3)可得第二条恒等式

S ( n , r ) = σ ( n r ) i = 1 n r ( 1 i + 2 i + + r i ) k i i k i k i !

注1. 除了上述定理外,连续整数的幂次和与 S ( n , k ) 还有以下关系

1 n + 2 n + + r n = k = 1 n k ! S ( n , k ) ( r + 1 k + 1 )

定理3. 对任意整数n,r,有互逆关系。

由Faà di Bruno的(4)可得第二个恒等式。

3. 对称群及集合分割上的组合数

这一节,将由已知的第一类及第二类Stirling数、Bell数、错排数的指数生成函数出发,利用Faà di Bruno推导出这些数的组合意义,从而可用Faà di Bruno求出错排数及Bell数加上更多限制的生成函数。首先,先考虑第一类Stirling数 S ( n , k ) 及错排数Dn

对称群Sn中的每一个置换皆可唯一分解成互斥圈的乘积。若一个置换分解成互斥圈后,长度为i的圈有ki ( i = 1 , 2 , , n ) ,则称这个置换是属于 1 k 1 2 k 2 n k n 型的置换。

下面是已知的事实:

给定非负整数 k 1 , k 2 , , k n ,则Sn 1 k 1 2 k 2 n k n 型置换的个数为

n ! 1 k 1 k 1 ! 2 k 2 k 2 ! n k n k n ! (5)

第一类Stirling数的指数生成函数为 n r s ( n , r ) t n n ! = 1 r ! ( ln ( 1 + t ) ) r

错排数Dn的指数生成函数为 n 0 D n t n n ! = 1 e t ( 1 t )

定理4. 对任意正整数n, D n = k 2 + k 3 + + k n = k 2 k 2 + 3 k 3 + + n k n = n n ! 2 k 2 k 2 ! n k n k n !

证明:取 f ( x ) = ln ( x ) ϕ ( t ) = 1 e t ( 1 t ) ,显然 f ϕ ( 0 ) = 0 。由于 ϕ ( t ) t n 项系数为 [ ϕ n ] = D n n !

因为 f ϕ ( t ) = ln ( 1 t ) t = t 2 2 + t 3 3 + ,对任意正整数i

[ f ϕ n ] = { 0 , i = 1 1 i , i 1

因为 [ f ϕ 1 ] = 0 ,只有 k 1 = 0 [ f ϕ 1 ] k 1 [ f ϕ 2 ] k 2 [ f ϕ n ] k n k 1 ! k 2 ! k n ! 0 ,故而,由Faà di Bruno(4)可知

D n n ! = σ ( n ) [ f ϕ 1 ] k 1 [ f ϕ 2 ] k 2 [ f ϕ n ] k n k 1 ! k 2 ! k n ! = k 2 + k 3 + + k n = k 2 k 2 + 3 k 3 + + n k n = n 1 2 k 2 k 2 ! n k n k n !

也即 D n = k 2 + k 3 + + k n = k 2 k 2 + 3 k 3 + + n k n = n n ! 2 k 2 k 2 ! n k n k n !

注4:定理4的恒等式可直接由错排数本身的定义得到,由(5)可知定理4中连加符号内的被加项分别为 1 k 1 2 k 2 n k n 型置换及 1 0 2 k 2 n k n 型置换的个数。

先观察定理中的置换,这些置换所对应的 ( k 1 , k 2 , , k n ) σ ( n , r ) 内的元素,从而 ( k 1 , k 2 , , k n ) 必须满足

{ k 1 + k 2 + + k n = r k 1 + 2 k 2 + + n k n = n

另一方面,对所有 1 i n k i 代表长度为i的互斥圈的个数,所以 k 1 + k 2 + + k n = r 代表互斥圈的总数为r,而 k 1 + 2 k 2 + + n k n = n 则代表这个置换是 S n 中的元素,因此 σ ( n , r ) n ! 1 k 1 k 1 ! 2 k 2 k 2 ! n k n k n ! 可以解释为 S n 中互斥圈总数为r的置换的个数,于是得到了组合学中广为人知的 | S ( n , r ) | 的组合意义: S n 中互斥圈总数为r的置换的个数。

而定理4中的置换对应的 ( k 1 , k 2 , , k n ) 需满足

{ k 1 = 0 k 1 + k 2 + + k n = k ( 1 k n ) k 1 + 2 k 2 + + n k n = n

其中因为 k 1 = 0 ,这些置换分解成相斥圈后,都不具有长度为1的圈,换句话说就是没有固定点,因此 k 2 + k 3 + + k n = k 2 k 2 + 3 k 3 + + n k n = n n ! 2 k 2 k 2 ! n k n k n ! 为Sn的所有置换中,不具有固定点的置换的个数,用比较贴近生活的言语表达即是:

有n个人参加宴会,每个人各戴一顶帽子,宴会进场时,n人皆将帽子交给服务生,宴会结束后n人拿回帽子,没有人拿到自己帽子的情况的个数。

这就是错排最初被提出时的问题。

f ( x ) = ln ( x ) ϕ ( t ) = 1 e t ( 1 t ) ,显然 f ϕ ( 0 ) = 0 。由于 ϕ ( t ) t n 项系数为 [ ϕ n ] = D n n !

因为 f ϕ ( t ) = ln ( 1 t ) t = t 2 2 + t 3 3 +

由定理4证明的过程中,发现利用Faà di Bruno得到的等式可以更进一步的探讨错排的问题。

回顾当时的证明,取 D n 的指数生成函数 ϕ ( t ) = 1 e t ( 1 t ) f ( x ) = ln ( x ) ,所以 f ϕ ( t ) = ln ( 1 t ) t = t 2 2 + t 3 3 + ,因为 ln 1 e t ( 1 t ) 中的 ln 1 e t 消去了 ln 1 1 t 的t系数, [ f ϕ 1 ] = 0 迫使 k 1 = 0 代表置换中长度为1的圈的个数为0,因为,可以得到此类置换的个数。

现在观察 ϕ ( t ) ,首先注意到因为 1 1 t = 1 + t + t 2 + t 3 + = 1 + t + 2 ! t 2 2 ! + 3 ! t 3 3 ! +

所以 1 1 t 为对称群元素个数 | S n | 的指数生成函数,而由先前讨论可知, 1 1 t 乘上 e t 后使得具有长度1的圈被排除而得到 D n 的指数生成函数,因此猜测经由乘上e的幂次也能排除具有其他长度圈的置换。

定理5. 对任意正整数n,

B n = σ ( n ) n ! ( 1 ! ) k 1 k 1 ! ( 2 ! ) k 2 k 2 ! ( n ! ) k n k n !

证明:取 f ( x ) = e x ϕ ( t ) = e t 1 ,显然 ϕ ( 0 ) = 0 。对所有正整数i

[ ϕ i ] = 1 i ! ,故而 [ f ϕ 1 ] = e [ e t 1 n ] = B n n !

利用Faà di Bruno的(3)得 B n n ! = σ ( n ) 1 ( 1 ! ) k 1 k 1 ! ( 2 ! ) k 2 k 2 ! ( n ! ) k n k n !

等号两边同时乘以 n ! 得到 B n = σ ( n ) n ! ( 1 ! ) k 1 k 1 ! ( 2 ! ) k 2 k 2 ! ( n ! ) k n k n !

注8. 定理也可以直接由Bell数的定义把定理从r = 1到r = n的情况加起来得到,但这里由其生成函数出发,利用Faà di Bruno也能导出相同的结论。

定理中连加符号内的被加项为 n ! ( 1 ! ) k 1 k 1 ! ( 2 ! ) k 2 k 2 ! ( n ! ) k n k n !

可知其恰为将n元集合分割成ki个i元子集(1 ≤ I ≤ n)的方法数,定理中每一个被加项对应到的 ( k 1 , k 2 , , k n ) 都是 σ ( r ) 里的元素,即必须满足

{ k 1 + k 2 + + k n = r k 1 + 2 k 2 + + n k n = n

所以可以将每一个被加项都视为一个分割, k 1 + k 2 + + k n = r 表示这个分割的子集合总数, k 1 + 2 k 2 + + n k n = n 则代表被分割的集合里有n个元素,因此可以解释成将n元相异集合分割成r个子集的方法数,于是得到一般熟知的 S ( n , r ) 的组合意义—将n元相异集合分割成r个子集合的方法数。又因为 σ ( n ) = k = 0 n σ ( n , k ) ,立即可以得到Bn的组合意义—分割n元相异集合的方法数。

因为定理的启发,利用类似于定理中排除有某些特定长度圈的置换的方法,也可以用来限制分割集合时产生的子集大小,结果如下:

给定一集合 P N 。一个n元集合若分割成子集后不具有元素个数为 j P 的子集时,称此种分割为“P免除分割”。

4. 总结

代数组合学是组合数学中最重要的分支之一,它有着丰富的内容和广泛的应用,著名的Faà di Bruno公式,Bell多项式以及对称群循环指标是代数组合学中重要的组成部分,三者息息相关,密不可分 [8]。Faà di Bruno公式起源可以追溯到十九世纪,其主要是解决复合函数的高阶求导问题,很多数学家曾对此作出过贡献:Scott,Meyer,Hoppe,Faà di Bruno等等。Faà di Bruno公式有很多种表达形式,最常见的是用Bell多项式来表示,然而,事实上,Bell多项式本来是集合分拆的数学工具,与Faà di Bruno公式并无直接联系,后来,Riordan发现Bell多项式非常适合对Faà di Bruno的表达。

同样运用Faà di Bruno公式,王兴华证明了对于复多项式,Halley族迭代的所有额外不动点都是斥性的。该结果比之前的发现要早六年解决其中的问题。Faà di Bruno公式也被应用到曲线插值的误差估计上,Floater和他的合作者研究了参数选择在插值逼近阶上的影响,证明了弦长参数值的三次及三次以上多项式能够饱和逼近曲线,井且给出了相应的算法.该算法也能够被用来估计曲线的长度 [9]。另外,王兴华通过利用Bell多项式和对称群的循环指标,建立了同时获得多项式所有根的并行和区间迭代.该结果成为文献的主要部分,另外,UlrichAbel和Mircea Ivan把Bell多项式运用到了区间趋向于零时差商的澈分中值定理的渐进展开上面。

一个复合函数的高阶导数是什么?至今最好的答案就是Faà di Bruno公式,从Gousat和ValleePoussin的著作中可以看到,曾经Faà di Bruno公式被认为是实分析中的一个结果。而Riordan和Comtet把它看作是代数组合学的一部分,Faà di Bruno公式广泛地出现在整数拆分、数学统计、矩阵理论、有限差分的微积分学、科学计算、对称函数等各个领域中。

Faà di Bruno公式有着悠久的历史。Faà di Bruno曾经把他的公式发表在期刊上,两者都是在1855年的12月份,文章较短,本质上基本相同。现在常见的Faà di Bruno公式有三种形式分别是集合分拆形式,Bell多项式形式和行列式形式。事实上,根据Lukacs关于函数高阶导数的显式表达问题,最早在1810年Lacriox的关于微积分的论文中已经提到 [10]。除了Faà di Bruno公式,曾经已有不少数学工作者给出了各种类型的表示,Johnson在他的文章中列出了四种有关复合函数高阶导数的表达式,即分别是T.A.,Scott,Meyer和Hoppe公式。自从Faà di Bruno公式第一次出现以后,对它的各种证明不断涌现,而且,Faà di Bruno公式被广泛地应用到数学的各个领域。同时也出现了与其相关的一系列恒等式,随着问题的深化Faà di Bruno公式也进一步发展为多变量的情形和q的形式。

对称群的循环指标多项式是表现力极强的解析工具,有很大一类超越函数可以用对称群的循环指标多项式予以表示。王兴华利用对称群的循环指标多项式建立了一族求多项式所有零点的并行迭代。后来,王兴华给出了数值分析中经典的Hermite插值的对称群循环指标表示,可以看到Hermite插值可以被对称群循环指标表达得十分简单明了。人们对对称群的循环指标多项式重视程度还是不够,没有真正认识到它的重要性,比如拿十分优秀的数学手册来说,尽管其第24章组合分析中已经介绍了对称群的循环指标多项式,但是很少把它作为一种工具使用对于王兴华等人的结果完全可以列入第25章插值,数值微分和数值积分中。

差商是数值分析中很重要的一部分内容,特别是对样条理论的发展一直扮演著主要角色众所皆知,在B样条的递归关系中,它是一个基本工具。商有很多种定义方法,通常把它当作是牛顿插值的系数或者直接通过递归来定义。差商有多种表示方式,常见的有行列式比表示、Peano核表示方式、围道积分形式和Genocchi-Hermite表示形式。由于差商在数值分析中的重要性,所以一直以来许多作者对其性质和应用的研究从未间断过早期,对差商性质的应用多见于数值微分上,如上世纪九十年代,Powers等人研究了两种微分中值定理的渐进性质。第一种是与差商相关的形式,第二种是推广的Flett中值定理。2004年,UlrichAbel等人根据差商推广了该定理。同样利用差商的性质,Horwitz研究了另一种类型的中值的极限行为近来,杨专门对差商的性质进行了讨论,并把它们应用到基于第二类Chebyshev多项式零点的数值积分公式上。Floater在中给出了差商展开的误差,在此基础上,C.deBoor给出了差商的差商展开进一步,王兴华与其合作者得到了数值差商的一些结果,用一组节点的差商来表示同一个函数的给定节点的差商,当给定节点比较密集因而在其上的差商难以直接计算时,把它用另一组由自己选取的较为疏朗的节点上的差商表示出来,以便能顺利的进行数值计算。有趣的是,Adell等人应用上面的一系列结果于统计之中。

在参变量多项式插值的研究中,需要估计一个复合函数 h = f g 差商的界,王兴华和Floater分别用不同的方法独立的解决了该问题,并且,把该公式应用到行列式恒等式上,从而得到了一系列漂亮的组合恒等式。

文章引用

沈 萍. Faà di Bruno公式在恒等式及计数上的应用
Application of Faà di Bruno’s Formula in Identity and Counting[J]. 理论数学, 2022, 12(12): 2231-2238. https://doi.org/10.12677/PM.2022.1212240

参考文献

  1. 1. 李海侠. 组合数学中的容斥原理及其应用实例[J]. 内江科技, 2020, 41(3): 3.

  2. 2. 李丹. 组合数学中链与反链的教学探索及实践[J]. 软件导刊, 2020, 19(12): 71-73.

  3. 3. 姚娟, 于喜志. 基于点云数据特征组合数学模型的海量数据统计方法研究[J]. 现代职业教育, 2020(23): 2.

  4. 4. 牟丽丽. 代数组合学中的若干问题[D]: [博士学位论文]. 大连: 大连理工大学, 2016.

  5. 5. Gray, W.S. and Duffaut Espinosa, L.A. (2015) A Faà di Bruno Hopf Algebra for Analytic Nonlinear Feedback Control Systems. 149-217. https://doi.org/10.4171/143-1/4

  6. 6. Kahn, J., Steger, A. and Sudakov, B. (2021) Combinatorics. Oberwolfach Reports, 17, 6-89. https://doi.org/10.4171/OWR/2020/1

  7. 7. Naprienko, S. (2021) Combinatorics of Iwahori Whittaker Func-tions.

  8. 8. Settepanella, S. and Yamagata, S. (2021) Combinatorics of Discriminantal Arrangements. 10.48550/arXiv.2101.00544

  9. 9. Pinsky, R.G. (2021) A Survey of the Bridge between Combinatorics and Probability. https://doi.org/10.54550/ECA2021V1S3S3

  10. 10. Galvez-Carrillo, I., Kock, J. and Tonks, A. (2014) Groupoids and Faa di Bruno Formulae for Green Functions in Bialgebras of Trees. Advances in Mathematics, 254, 79-117. https://doi.org/10.1016/j.aim.2013.12.015

期刊菜单