1. 引言
图的完美匹配计数问题是一个具有很强应用背景的问题,主要涉及到量子化学和统计物理领域。在量子化学领域,完美匹配被称为Kekulé结构,而Kekulé结构是描述分子稳定性的一个重要指标;在统计物理领域,完美匹配则被称为Dimer构型,而Dimer构型是描述在格点模型中粒子分布的一种方式[1] [2]。因此,研究完美匹配的计数问题不仅有助于理解这些领域的复杂现象,还可能为解决实际问题提供理论依据和计算方法。同时,对图的完美匹配计数问题的研究不仅在实践中具有重要意义,而且在图论的其它领域也有广泛的应用,如匹配算法、网络流等领域[3] [4]。
然而,Valiant证明了即使是二部图的完美匹配计数问题也是NP-难的[5]。科学家Kasteleyn首次提出运用图的Pfaffian定向来研究图的完美匹配计数问题[6]。并指出如果图G有一个Pfaffian定向,那么其完美匹配数就能在多项式时间内算出。Kasteleyn的一个经典定理就是证明了所有的平面图都是Pfaffian的[7]。利用Pfaffian定向的方法,Kasteleyn得到了平面四角格子和环面格子图完美匹配数的精确表达式[6]。Little推广了Kasteleyn的结果并且证明如果一个图G不包含
的细分,则G有一个Pfaffian定向[8]。尽管有许多等价条件证明图Pfaffian定向的存在性[8] [9],但没有多项式时间算法来检查图G的一个给定定向是否为Pfaffian定向。Lu和Wu构造了能嵌入到莫比乌斯带,Klein瓶上的四方形网格的一个Pfaffian定向,并给出了这些图的完美匹配数的显式表达式[10]。运用Pfaffian方法还可以简化某些图类的完美匹配数的计算[11]-[13]。Yan和Zhang利用圈C4,路P4与树的乘积图的Pfaffian性得到了一些图的完美匹配数的显式表达式[14]。Propp介绍了许多关于图的完美匹配计数问题的应用背景及其研究进展[15]。更多相关的内容可以参考文献[16]-[18]。
除了利用Pfaffian定向来计算图的完美匹配数外,还可以利用图自身的特殊结构求其匹配数。比如在1985年,张福基等人研究了树状多六边形(简写为TPH)的完美匹配数,得出TPH图的完美匹配计数相当于树状苯型渺位稠合多环芳香体系的Kekulé结构计数[19]。在1997年,张和平等人研究四角系统完美匹配数,证明非基本四角系统的完美匹配数等于一些基本子四角系统完美匹配数的乘积[20]。在2005年,林泓等人研究了若干四角系统的完美匹配计数,利用组合递推法给出了几类四角系统的完美匹配数的显式表达式[21]。在2014年,唐保祥等人利用划分、求和,再嵌套递推的方法给出了几类特殊图完美匹配数的显式表达式[22]。马京成等人研究了3-正则Halin图的完美匹配计数问题,给出了3-正则Halin图的完美匹配数的一个递推公式,同时,他们还给出了两类特殊的Halin图的完美匹配数下界的计算公式[23]。2023年,叶银珠等人利用Dong等人给出的偶数条边树的线图完美匹配数的表达式,得到一类繁星树线图的完美匹配数的最大值和最小值[24] [25]。更多关于完美匹配数的计算问题,请参考文献[26]-[28]。
富勒烯是一种完全由碳组成的中空分子,形状有球形,椭球形,柱形和管状。而管状富勒烯图就是两端有封口的单壁碳纳米管或巴基管[29]。富勒烯由于其特殊的机械、物理、化学性能在工程材料、催化、吸附分离、储能器件电极材料等诸多领域得到了广泛应用[30]。在1998年,Došlić研究了富勒烯图的完美匹配数的下界,证明n个顶点的富勒烯图至少有
个完美匹配[31]。在2001年,张和平等人改进了Došlić等人的结果,得出n个顶点的富勒烯图的完美匹配数的下界为
[32]。在2009年,Kardoš等人利用四色定理得到在一般情况下的具有n个顶点且不包含非平凡5-边割的富勒烯图至少有
个完美匹配[33]。在2012年,鞠阳研究了富勒烯邻接矩阵积和式的界估计结果及其与富勒烯完美匹配数之间的定量关系,得出富勒烯完美匹配数的一个新的下界估计:
[34]。在2016年,王芳研究了富勒烯图完美匹配数和共振型个数的比较,得出对所有的富勒烯图,富勒烯图中共振型个数严格小于完美匹配个数[35]。后来在2020年,钱进研究了富勒烯的完美匹配数的下界及其稳定性,得出比之前更优的富勒烯的完美匹配数的下界估计[36]。Lovász和Plummer曾经提出一个关于图的完美匹配计数问题的一个猜想:任意2-边连通3-正则图都有指数多个完美匹配[37]。Esperet等人在2011年证明了n个顶点的2-边连通3-正则图至少有
个完美匹配[38],从而解决了Lovász和Plummer提出的猜想。
本文研究的管状富勒烯图是由六个同心的六边形层组成,两端都由一个六边形以及与这个六边形相邻的六个五边形面构成的顶盖封口,记作F。本文证明了以下结果,验证了Lovász和Plummer猜想的正确性。
定理1 对于管状富勒烯图F,F的完美匹配数为
(其中s为F的顶点数)。
本文的结构安排如下:第一部分,介绍了完美匹配计数问题的研究背景,研究进展,富勒烯的完美匹配数的研究进展,以及本文的主要结论。第二部分为准备工作,给出本文需要用的一些符号和概念。第三部分,利用图的结构和递推公式证明本文主要结论。
2. 准备工作
设
和
分别是图G的顶点集和边集,分别用
和
表示图G的顶点集和边集的个数。设
,
,如果
,通过从G中删除e而保持顶点和其余边不变来获得一个
条边的图,得到的图用
表示。类似地,如果
,通过从G中删除顶点v以及所有与v关联的边来得到一个
个顶点的图,得到的图用
表示。对
,设
是通过从G中删除在H中的元素而获得的子图。
图G的途径是指一个点边交替序列
,其中
是图的一条边(
),称W是从
到
的一条途径,或一条
途径。顶点
分别称为W的起点和终点,而
称为W的内部顶点,W中的边数k称为W的长度。若
,则W称为闭途径。若途径W中,
互不相同,则W称为迹。若
,则W称为闭迹。若在迹W中,
互不相同,则W称为一条路。
的路称为一个圈,记为
。长度为奇数的圈称为奇圈,长度为偶数的圈称为偶圈。
对于一个图G,如果能将图G边不相交(除端点处相交)的画在平面上,则称图G是可平面图,这个画出的无边相交的图称作图G的平面嵌入图。G的平面嵌入图将平面划分为若干连通的闭区域。这些闭区域称为G的面。有限的区域称为有限面或内部面,无限的区域称为无限面或外部面,包围面的边称为该面的边界。如果两个面有一条共同的边,则这两个面相邻。一个面与其边界上的顶点和边是关联的,如果G是2-连通可平面图,则G的每个面的边界都是圈,为方便起见,我们经常用面的边界来代替面。
Figure 1. Tubular fullerene graphs F
图1. 一类管状富勒烯图F
特别地,富勒烯图中的五边形面和六边形面分别称为5长圈和6长圈。
图G的匹配是一组两两不交的边构成的集合。如果M是一个匹配,则称M的每条边的两个端点在M下匹配,并且与M的一条边关联的每个顶点称被M覆盖。一个完美匹配是覆盖图的所有顶点的匹配。如果图G有两个完美匹配M1和M2,M1和M2如果有一条边不同,则称M1和M2是图G的两个不同的完美匹配。图G的完美匹配数就是图G的所有不同的完美匹配的数目,记为
。
3. 主要结果及其证明
管状富勒烯图F是由同心的六边形层,两端分别加一个帽子构成,每个帽子又由一个六边形面以及六个五边形邻面构成,除了两端的帽子,其余的每个同心六边形层都由六个相邻的六边形构成(如图1所示为该类管状富勒烯图的平面嵌入图)。我们对F做一些记号,记
为F的所有同心层,这里H0和Hn是F的最里层和最外层,都由一个六边形面构成。而H1是和H0相邻的由六个五边形面构成的一层,即
构成一端的帽子,
构成另一端的帽子,对于
,六边形层
与
与
均相邻,对
,我们设
。对
上的顶点顺时针进行标号
,则C1和Cn是6长圈,
是12长圈,且
,
,
Figure 2. The label of tubular fullerene graphs
图2. 标号的管状富勒烯图
如图2所示。连接
和
的边称为横跨边(比如
是C1和C2之间的一条横跨边),
和
之间的所有横跨边构成集合,记作
。
本文中我们都沿用如图2所示的对管状富勒烯图F的标号,由于F由
层同心层构成,故F有
个顶点,12个五边形面,
个六边形面。记该类管状富勒烯图为
。
定理1的证明思路:设
表示管状富勒烯图
的所有完美匹配构成的集合,
表示所有完美匹配的个数,C1和C2之间的横跨边集
,因为圈C1上有偶数个点,所以对于任意一个完美匹配
,
不可能为奇数,即
可能取值为
。下面我们分别求出在每一种情况下的
的完美匹配数,进而求和得到
的大小。
引理2 设
,则
。
证明:对任意的
,当
时,6长圈
上的点只能与C1上的点匹配,则圈C1上的完美匹配要么是
,要么是
,即圈C1上有两种完美匹配。记
是从图
中删除圈C1上的顶点后获得的子图,如图3所示。
观察图
,从12长圈
上的顶点中选一个2度点,比如
,则
要么与
匹配,要么与
匹配。当
与
匹配时,C2上的完美匹配为
;当
与
匹配时,C2上的完美匹配为
,则圈C2上也有两种完美匹配。再把从图
中删除C2上的点得到的子图,记作
。
对
进行类似对
的分析,依此类推,每个圈
上都有两种完美匹配,最后剩一个6长圈
,
也有两种完美匹配。则
(1)
Figure 3. Graphs
图3. 图
引理3 记
,则
。
证明:对任意的
,当
时,
或
或
(
;
,当
时,
;当
时,
)。对上述
,设
;
;
;则
(
且
),
。故
(2)
设
,
,
,
,
,
。
见图4,图5所示图
,
,
,记
为子图
的完美匹配的数目,对
,
,类似。
Claim 3.1
。
证明:因为
(
,当
时,
;当
时,
)。所以圈C1上的点
不可能被
覆盖(
),这与
是完美匹配矛盾。比如当
时,
,则圈C1上的点
只能与圈C2上的点
匹配,但是边
,所以圈C1上的点
不可能被
覆盖,则当
且
时,
。
Claim 3.2
。
证明:因为
(
,当
时,
;当
时,
)。当
时,即
,C1上剩余四个点
被
覆盖的方式只有一种,即
。因为
,
,所以此时
(如图4所示)。
当
时,即
,C1上剩余四个点
被
覆盖的方式只有一种,即
,并且从
中删除C1上的点及
后获得的子图与
同构。因而,此时
。类似地,当
时,
。
综上,当
且
时,
。
Claim 3.3
。
证明:观察图
(图4),边集
,
为图
所有完美匹配构成的集合。对边集E2中是完美匹配边的数目
进行分类,其中
。因为圈C2在图
上剩余偶数个点,所以对于任意一个完美匹配
,
可能取值为
。又因
有一个悬挂点
(即度为1的点),因而点
必须与点
匹配,即
。因而
。当
或6时,因
,所以E2中除
以外剩余的五条边
中至少有三条边属于
。当
(
)属于
时,C2上剩余点
不会被
覆盖,与
Figure 4. Graphs
图4. 图
是
的完美匹配矛盾。当
属于
时,C2上剩余点
不会被
覆盖。同样与
是
的完美匹配矛盾。由此矛盾可知
且
。
当
时,由于
,故
或
或
或
或
。
当
时,C2上剩余点
被
覆盖的方式唯一,即
。从
中删除C2上的点及点
后剩下的图同构于
。因而此时
的完美匹配数等于
的完美匹配数,即
。由对称性,当
时,
。
当
时,C2上剩余点
被
覆盖的方式唯一,即
。从
中删除C2上的点及点
后剩下的图同构于
。因而此时
的完美匹配数等于
的完美匹配数,即
。由对称性,当
时,
。
当
时,C2上剩余点
被
覆盖的方式唯一,即
。从
中删除C2上的点及点
后剩下的图同构于
。因而此时
的完美匹配数等于
的完美匹配数,即
。
综上,
(3)
Claim 3.4
。
证明:因为
(
,当
时,
;当
时,
)。当
时,即
时,C1上剩余四个点
被
覆盖的方式只有一种,即
。从
中删除C1上的点及
后剩下的子图为
。故此时
。由对称性,当
或
时,
。当i的取值分别为4,5或6时,与i的取值分别为1,2或3时的情况重复。
所以当
且
时,
有
种完美匹配。
Claim 3.5
。
证明:观察图
(图5),边集
,设
为图
所有完美匹配构成的集合。对边集E2中是完美匹配边的数目
进行分类,其中
。因为圈C2在图
上剩余偶数个点,所以对于任意一个完美匹配
,
不可能为奇数,
可能取值为
。因为在图
中,由点
诱导的子图是两个四长路:
和
,分别记作P1和P2。所以对任意的
,
。另一方面,若横跨边
中有两条相继的边
(
)属于
,则路P1上剩余一点
不会被
覆盖,矛盾。该矛盾说明
中不能有两条相继的边属于
,同理,
中也不能有两条相继的边属于
。这意味着
。
记
,
,则
,
。由前面的分析可知:
,
(其中
)。所以当
时,
且
。即
,但P1上剩余三个点
不可能同时被
覆盖,矛盾。该矛盾说明
。
当
时,
且
。由对称性,我们又分为以下三种情况:
情况1:
(
,当
时,
)。
因为
且
。故
的取值为
或
。当
时,C2上剩余点
被
覆盖的方式唯一,即
。从
中删除C2上的点及
后剩下的图同构于
。因而此时
的完美匹配数等于
的完美匹配数,即
。由对称性,当
时,
。
情况2:
(
,当
时,
)。
因为
且
。故
的取值为
或
或
或
。当
时,C2上剩余点
被
覆盖的方式唯一,即
。从
中删除C2上的点及点
后剩下的图同构于
。因而此时
的完美匹配数等于
的完美匹配数,即
。由对称性,当
或
或
时,
。
情况3:
(
,当
时,
)。
因为
且
。故
的取值为
或
或
。当
时,C2上剩余点
被
覆盖的方式唯一,即
。从
中删除C2上的点及点
后剩下的图同构于
。因而此时
的完美匹配数等于
的完美匹配数,即
。由对称性,当
或
时,
。
综上,
(4)
观察
(如图5)进行类似Claim 3.3,Claim 3.5的分析,则
(5)
由(3)~(5)式得
Figure 5. Graphs
(left), graphs
(right)
图5. 图
(左),图
(右)
(6)
(7)
再由(6)~(7)式得
(8)
由(3) – 2 × (5),(4) – (5)分别得
(9)
(10)
由(8),(9)得
(11)
其中(11)式的特征方程的根为
,
,
。易知
,
。得出递推式(11)式的通解为
(12)
由(8),(10)得
(13)
线性递推式(13)式的特征方程的根与递推式(11)式的特征方程的根相同,易知
,
,
。因此递推式(13)式的通解为
(14)
综上,当
时,图
的完美匹配数为
(15)
引理4 记
,则
。
证明:对任意的
,当
时,
或
或
(
,当
时,
;当
时,
)。对上述
,设
,
,
。则
(
且
),
。故
(16)
设
,
,
,
,
,
。
见图6,图7所示图
,
,
。记
为子图
的完美匹配的数目,对
,
类似。
Claim 4.1
。
证明:因为
(
,当
时,
;当
时,
)。所以此时
必须与
匹配。与
矛盾,因而
。比如当
时,
,则C1上剩余点
必须分别与
匹配,此时
,与
矛盾。
Claim 4.2
。
证明:因为
(
,当
时,
;当
时,
)。所以此时
必须与
匹配。与
矛盾,因而
。比如当
时,
,则C1上剩余点
必须分别与
匹
配,此时
,与
矛盾。
Claim 4.3
。
证明:因为
(
,当
时,
;当
时,
)。当
时,
,C1上剩余两个点
被
覆盖的方式只有一种,即
。因为
,
。从
中删除
及点
后剩下的图即为
。因而此时
,由对称性,当
时,
,所以
。
Figure 6. Graphs
图6. 图
Claim 4.4
。
证明:观察图
(如图6),边集
,设
为图
所有完美匹配构成的集合。对边集E2中是完美匹配边的数目
进行分类,其中
。因为圈C2在图
上剩余偶数个点,所以对于任意一个完美匹配
,
不可能为奇数,
可能取值为
。
对图
,点
必须与点
匹配,点
必须与点
匹配,点
必须与点
匹配,则边集
,所以
可能取值为4或6。又当
时,即
,则圈C2上点
不能被
覆盖,矛盾。从而
。
对任意
,当
时,由于
,故
或
或
。当
时,则圈C2上点
被
覆盖的方式只有一种,即
。从
中删除
中边关联的点以及点
后剩下的图同构于
,因而此时
。由对称性,当
时,
。当
时,则圈C2上点
被
覆盖的方式只有一种,即
。从
中删除
中边关联的点以及点
后剩下的图同构于
,因而此时
。
综上,
(17)
Claim 4.5
。
证明:观察图
(如图7),边集
,设
为图
所有完美匹配构成的集合。对边集E2中是完美匹配边的数目
进行分类,其中
。因为圈C2在图
上剩余偶数个点,所以对于任意一个完美匹配
,
不可能为奇数,
可能取值为
。且在图
中,圈C2被分成两个2-长路
和两个孤立点
(见图7左)。
所以对任意
,
且
中恰好有一条边属于
,同理,
Figure 7. Graphs
(left), graphs
(right)
图7. 图
(左),图
(右)
中恰好有一条边也属于
,为此,我们有
且
或
或
或
。当
时,圈C2上点
被
覆盖的方式只有一种,即
。从
中删除
中边关联的点以及点
后剩下的图同构于
,因而此时
。当
时,圈C2上点
被
覆盖的方式只有一种,即
。从
中删除
中边关联的点以及点
后剩下的图同构于
,因而此时
。由对称性,当
时,
。当
时,圈C2上点
被
覆盖的方式只有一种,即
。从
中删除
中边关联的点以及点
后剩下的图同构于
,因而此时
。
综上,
(18)
观察
(如图7)进行类似Claim 4.4,Claim 4.5的分析,可得
(19)
由(17)~(19)式得
(20)
(21)
再由(20)~(21)式得
(22)
由(17),(22)联立得
(23)
其中(23)式的特征方程的根为
,
,
。易知
,
。得出递推式(23)式的通解为
(24)
故对任意的
,当
时,图
的完美匹配数为
(25)
引理5 令
,则
。
证明:当
时,即C1和C2之间的横跨边
是完美匹配
中的边。从图
中删除横跨边
中的边关联的点得到图
(如图8左所示),则C2和C3之间的横跨边
也是完美匹配
中的边。再从图
中删除横跨边E2中的边关联的点得到图
,则C3和C4之间的横跨边E3也是完美匹配
中的边。依此类推,最后
和
之间的横跨边
也是完美匹配
中的边。则
时,图
的完美匹配如图8右所示。
即当
时,图
的完美匹配数为
(26)
Figure 8. Graphs
(left), a perfect matching of graphs
(right)
图8. 图
(左),图
的一种完美匹配(右)
定理1的证明:结合引理2~5,可得
(27)
设s是F的顶点数,则
,即
代入等式(27),可得
(28)
即图F的完美匹配数为
(其中s为F的顶点数)。
4. 结论
在本文中,我们得出一类管状富勒烯图F的完美匹配数的计算公式,验证了Lovász和Plummer的任意2-边连通3-正则图都有指数多个完美匹配的猜想的正确性。对于其他管状富勒烯图,我们可以尝试运用本文中的方法求出其完美匹配数的计算公式。在推导文中的递推公式时过程比较繁琐,后续可以寻求更好的方法解决此类问题。
基金项目
本文研究得到了国家自然科学基金项目(No.11801148, 11626089)以及河南理工大学博士基金项目(No.B2014-060)的资助。
NOTES
*通讯作者。