1. 引言
现代数学可以分为两大类:一类是研究连续对象的,如分析、方程等;另一类就是研究离散对象的组合数学。微积分以及现代数学的发展,为现代工业的发展打下了坚实的基石 [1]。而组合数学的发展,为本世纪电脑的变革打下了坚实的基石。从十七世纪到十八世纪,组合数学和随机运算的发展,到了十九世纪,它已经彻底地颠覆了传统的概念,产生了一种全新的、充满活力的新的数学派别 [2]。因此,组合数学的发展,使传统的数理学科成为了以解析与代数为主导的学科。尤其是随着计算机越来越多地渗透到各行各业,除了数字运算之外,也出现了许多快速发展的非数字运算方法,这些方法大多采用了合成方法 [3]。
在计算机科学、编码与密码学、物理等领域,都有着举足轻重的作用。化学、生物等学科都具有很大的实用价值,即计算科学 [4]。计算机的运算方法只是一个离散的物体,而不是一个数字,它的科学就在于它的结合,它给人一种感觉,就像它有思想一样,因此,它的计算方法就成为了计算机的中心 [5]。在商业经营、运输计划、作战指挥、财务分析等方面,美国有一间公司以组合数学为名;他们运用了综合数学的理论,以增加公司经营的收益,这个公司做的很好,另外,实验设计也是一个很实用的课题,其主要的数学理论是结合的。美国已经有一些专业公司,在此基础上,对其进行了研究 [6] [7]。近年来,一名德国知名的复合数学家运用复合数论对药品结进行了深入的探讨,为医药企业节约了成本。简而言之,组合数学的基础理论与实践,涉及计算机科学,生物学,化学;心理学,金融和基因工程学的最新的运用,其范围十分广阔。
本文首先研究使用Faà di Bruno公式求组合恒等式,然后进行对称群及集合分割上的组合数的研究,最后对论文进行总结,发觉Faà di Bruno公式的应用前景。由于本文研究偏理论研究,因此将不使用仿真例子进行研究,后续可在此基础上进行相关深入的研究。
2. 用Faà di Bruno求组合恒等式
对任意给定幂级数
,以
表示
的系数,即
。
假设f有反函数,考虑复合函数
和
,Chou,Hsu与Shiue将Faà di Bruno改写,得到下列两个式:
(1)
(2)
因为
和
在(1)和(2)中所处的地位恰好相反,称这两条恒等式为一组互逆关系式。特别地,考虑
,
,令
,若
,则由(1)和(2)有以下互逆关系:
(3)
(4)
定理1 对任意正整数n,有互逆关系
及
。
定理2. 对任意正整数n,r,有互逆关系
及
或
。
证明:令
,
,显然
,
故而
。
另一方面因为
观察
的系数可得
根据Faà di Bruno (4)可以得到
等号两边同乘n即得第一条恒等式;由Faà di Bruno的(3)可得第二条恒等式
注1. 除了上述定理外,连续整数的幂次和与
还有以下关系
定理3. 对任意整数n,r,有互逆关系。
由Faà di Bruno的(4)可得第二个恒等式。
3. 对称群及集合分割上的组合数
这一节,将由已知的第一类及第二类Stirling数、Bell数、错排数的指数生成函数出发,利用Faà di Bruno推导出这些数的组合意义,从而可用Faà di Bruno求出错排数及Bell数加上更多限制的生成函数。首先,先考虑第一类Stirling数
及错排数Dn。
对称群Sn中的每一个置换皆可唯一分解成互斥圈的乘积。若一个置换分解成互斥圈后,长度为i的圈有ki个
,则称这个置换是属于
型的置换。
下面是已知的事实:
给定非负整数
,则Sn中
型置换的个数为
(5)
第一类Stirling数的指数生成函数为
。
错排数Dn的指数生成函数为
。
定理4. 对任意正整数n,
。
证明:取
、
,显然
。由于
的
项系数为
。
因为
,对任意正整数i
因为
,只有
时
,故而,由Faà di Bruno(4)可知
也即
。
注4:定理4的恒等式可直接由错排数本身的定义得到,由(5)可知定理4中连加符号内的被加项分别为
型置换及
型置换的个数。
先观察定理中的置换,这些置换所对应的
为
内的元素,从而
必须满足
另一方面,对所有
,
代表长度为i的互斥圈的个数,所以
代表互斥圈的总数为r,而
则代表这个置换是
中的元素,因此
可以解释为
中互斥圈总数为r的置换的个数,于是得到了组合学中广为人知的
的组合意义:
中互斥圈总数为r的置换的个数。
而定理4中的置换对应的
需满足
其中因为
,这些置换分解成相斥圈后,都不具有长度为1的圈,换句话说就是没有固定点,因此
为Sn的所有置换中,不具有固定点的置换的个数,用比较贴近生活的言语表达即是:
有n个人参加宴会,每个人各戴一顶帽子,宴会进场时,n人皆将帽子交给服务生,宴会结束后n人拿回帽子,没有人拿到自己帽子的情况的个数。
这就是错排最初被提出时的问题。
取
、
,显然
。由于
的
项系数为
。
因为
由定理4证明的过程中,发现利用Faà di Bruno得到的等式可以更进一步的探讨错排的问题。
回顾当时的证明,取
的指数生成函数
、
,所以
,因为
中的
消去了
的t系数,
迫使
代表置换中长度为1的圈的个数为0,因为,可以得到此类置换的个数。
现在观察
,首先注意到因为
。
所以
为对称群元素个数
的指数生成函数,而由先前讨论可知,
乘上
后使得具有长度1的圈被排除而得到
的指数生成函数,因此猜测经由乘上e的幂次也能排除具有其他长度圈的置换。
定理5. 对任意正整数n,
证明:取
、
,显然
。对所有正整数i
,故而
。
利用Faà di Bruno的(3)得
。
等号两边同时乘以
得到
。
注8. 定理也可以直接由Bell数的定义把定理从r = 1到r = n的情况加起来得到,但这里由其生成函数出发,利用Faà di Bruno也能导出相同的结论。
定理中连加符号内的被加项为
。
可知其恰为将n元集合分割成ki个i元子集(1 ≤ I ≤ n)的方法数,定理中每一个被加项对应到的
都是
里的元素,即必须满足
所以可以将每一个被加项都视为一个分割,
表示这个分割的子集合总数,
则代表被分割的集合里有n个元素,因此可以解释成将n元相异集合分割成r个子集的方法数,于是得到一般熟知的
的组合意义—将n元相异集合分割成r个子集合的方法数。又因为
,立即可以得到Bn的组合意义—分割n元相异集合的方法数。
因为定理的启发,利用类似于定理中排除有某些特定长度圈的置换的方法,也可以用来限制分割集合时产生的子集大小,结果如下:
给定一集合
。一个n元集合若分割成子集后不具有元素个数为
的子集时,称此种分割为“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等人应用上面的一系列结果于统计之中。
在参变量多项式插值的研究中,需要估计一个复合函数
差商的界,王兴华和Floater分别用不同的方法独立的解决了该问题,并且,把该公式应用到行列式恒等式上,从而得到了一系列漂亮的组合恒等式。