0 引言
多处理器系统已经成为解决科学与工程计算领域重大挑战性问题的必备工具,系统互连网络的结构性质对系统的性能起着不可忽视的作用。在构建实际多处理器系统时,除了要考虑系统的性能之外,还要充分考虑系统的制造成本。这使得系统互连网络的诸多设计要求之间包含着矛盾,选择网络时只能在不同要求中寻找一种均衡。超立方体网络
[1 ] 和星图网络
[2 ] 是较早被设计出的互连网络,已经被广泛应用于实际系统的构建。当今多处理器系统的应用范围越来越广,为了满足日益多元化的计算需求,又有许多互连网络被相继提出,如
k 元
n 方体网络
[3 ] 、(
n ,
k )-星图网络
[4 ] 和cactus-based网络
[5 ] 等。作为星图网络的推广,(
n ,
k )-星图网络不仅保留了星图网络的许多优良性质(如对称性、可嵌入性等),而且在规模的选取上更加灵活
[4 ] 。近来人们对(
n ,
k )-星图的网络性质进行了一系列的研究。Chiang等
[6 ] 确定了(
n ,
k )-星图的直径和平均距离;Cheng等
[7 ] 研究了(
n ,
k )-星图的强匹配排除问题;Chang等
[8 -9 ] 分别在Preparata-Metze-Chien(PMC)模型与比较诊断模型下分析了(
n ,
k )-星图的条件诊断度;Feng
[10 ] 在点故障模型下给出了使得(
n ,
k )-星图中不存在(
n -
m ,
k -
m )-星图子网络所需要破坏的最小点数的上下界;Kumar等
[11 ] 对(
n ,
k )-星图的定向直径进行了研究。
在实际中,一些特定计算任务(如计算量不大但可靠性要求较高的任务)的调度算法只需要系统指派某个子网络来执行,这样不仅能使系统资源得到有效利用,而且可以规避大规模网络有效性差的缺点。基于这一考虑,学者们广泛关注互连网络的无故障子网络存在概率估计问题。Das等
[12 ] 根据
n 维超立方体的结构性质给出了
n 维超立方体中存在无故障子超立方体的概率的计算公式。Chang等
[13 ] 给出了不同点可靠性下无故障(
n -1)维子超立方体的存在概率的计算方法;与Das等
[12 ] 的方法相比,Chang等
[13 ] 提出的方法在确保结果精度的同时具有更优的计算效率。近年来,采用Chang等
[13 ] 提出的方法,许多互连网络的无故障子网络存在概率估计问题得到了研究
[14 -19 ] 。
(
n ,
k )-星图的无故障子网络存在概率估计研究可以为基于(
n ,
k )-星图网络构建的实际系统的任务调度算法设计提供理论参考。Li等
[20 ] 分析了(
n ,
k )-星图中无故障(
n -1,
k -1)-星图子网络的存在概率;但(
n ,
k )-星图中任意规模子网络的存在概率估计尚未研究,这在一定程度上限制了(
n ,
k )-星图网络在实际系统中的应用。基于此,本文将在不同点可靠性下评估(
n ,
k )-星图中无故障(
n -
m ,
k -
m )-星图子网络(
m ≥ 1 )的存在概率。
1 基本概念和性质
本文中没有解释说明而直接运用的图论术语及记号参见文献[
21 ]。假定正整数
i 和
j 且满足
j ≤ i ,令
i 表示
1,2 , ⋯ , i ,
A i j 表示
i ! / ( i - j ) ! 。(
n ,
k )-星图(
n 和
k 为整数且满足
1 ≤ k ≤ n - 1 ),记作
S n , k ,其顶点集为
V S n , k = u 1 u 2 ⋯ u k u 1 , u 2 , ⋯ , u k ∈ n 且它 们两 两互 不相 同 。 两个顶点
u 1 u 2 ⋯ u k 和
v 1 v 2 ⋯ v k 相邻当且仅当以下两个条件之一成立:
(1) 存在一个整数 2 ≤ i ≤ k 使得 u 1 = v i , v 1 = u i 且对于所有的 2 ≤ s ≤ k s ≠ i 均有 u s = v s ;
(2) u 1 ≠ v 1 且对于所有的 2 ≤ s ≤ k 均有 u s = v s 。
由
S n , k 的定义可知,
S n , k 是一个具有
A n k = n ! / n - k ! 个顶点的
n - 1 -正则图。
S 4,2 和
S 4,3 如
图1 所示。
给定正整数 n , k , m 满足 1 ≤ k ≤ n - 1 , 1 ≤ m ≤ k - 1 。设 i 1 , i 2 , ⋯ , i m 为满足 2 ≤ i 1 < i 2 < ⋯ < i m ≤ k 的 m 个整数。设 a i 1 , a i 2 , ⋯ , a i m 为 n 中 m 个两两互不相同的数,记点集
M = b 1 b 2 ⋯ b i 1 - 1 a i 1 b i 1 + 1 ⋯ b i 2 - 1 a i 2 b i 2 + 1 ⋯ a i m b i m + 1 ⋯ b k : b 1 , b 2 , ⋯ , b i 1 - 1 , b i 1 + 1 , ⋯ , b i 2 - 1 , b i 2 + 1 , ⋯ , b i m + 1 , ⋯ , b k ∈ n \ a i 1 , a i 2 , ⋯ , a i m 且它 们两 两互 不相 同
(特别地,若 i m = k ,则 b i m + 1 ⋯ b k 为空字符串)。显然, S n , k 中由 M 导出的子图与 S n - m , k - m 同构。为了方便表述,用长度为 k 的字符串 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 来表示 S n , k 中由 M 导出的子图(本文称这种表示为 S n , k 的子网络简化表示),这里 X t = X X ⋯ X ︸ t 。将 S n , k 中与 S n - m , k - m 同构的子图称为 S n - m , k - m 子网络。
引理1[10 ] S n , k 中共有
k - 1 m n ! / n - m ! 个不同的
S n - m , k - m 子网络且任意一个
S n - m , k - m 子网络可用某个
X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 唯一表示,其中
1 ≤ m ≤ k - 1 。
2 无故障 S n - m , k - m 子网络的存在概率估计
2.1 理论分析
在实际中,许多多处理器系统在构建通信链路时会采用冗余设计的策略以确保通信的稳定(该情形下系统的链路故障往往可以忽略不计)
[17 ] 。基于这一背景,本文对概率故障条件下(假定不考虑
S n , k 中的边故障,
S n , k 中各个顶点具有相同的可靠性且发生故障相互独立)
S n , k 中无故障
S n - m , k - m 子网络的存在概率进行估计。
根据引理1可知, S n , k 中不同的 S n - m , k - m 子网络共有 k - 1 m n ! / n - m ! 个。将 S n , k 中第 i 个 S n - m , k - m 子网络不存在故障点记为事件 A i ,这里 1 ≤ i ≤ k - 1 m n ! / n - m ! 。记 S n , k 中顶点的可靠性概率为 p 。令 R n , k m ( p ) 表示 S n , k 中存在无故障 S n - m , k - m 子网络的概率(这里 1 ≤ m ≤ k - 1 ), P r · 表示事件 · 发生的概率,则:
R n , k m ( p ) = ∑ i = 1 k - 1 m n ! / ( n - m ) ! P r ( A i ) + ( - 1 ) ∑ 1 ≤ i < j ≤ k - 1 m n ! / ( n - m ) ! P r ( A i ⋂ A j ) + ( - 1 ) 2 ∑ 1 ≤ i < j < l ≤ k - 1 m n ! / ( n - m ) ! P r ( A i ⋂ A j ⋂ A l ) + ⋯ + ( - 1 ) k - 1 m n ! / ( n - m ) ! - 1 P r ( A 1 ⋂ A 2 ⋂ ⋯ ⋂ A k - 1 m n ! / ( n - m ) ! ) 。 (1)
定理1 设 n , k , m 为正整数且满足 1 ≤ k ≤ n - 1 和 1 ≤ m ≤ k - 1 ,则
R n , k m ( p ) ≤ m i n k - 1 m n ! / n - m ! p A n - m k - m , 1 。
证明 注意到 S n , k 中各个顶点具有相同的可靠性且发生故障相互独立。由 V S n - m , k - m = A n - m k - m 可知,对于 1 ≤ i ≤ k - 1 m n ! / n - m ! 有 P r ( A i ) = p A n - m k - m 成立。根据式(1),
R n , k m ( p ) ≤ ∑ i = 1 k - 1 m n ! / n - m ! P r ( A i ) = k - 1 m n ! / n - m ! p A n - m k - m 。
又因为 R n , k m ( p ) ≤ 1 ,故 R n , k m ( p ) ≤ m i n k - 1 m n ! / n - m ! p A n - m k - m , 1 。定理1证毕。
引理2 设 u 和 l 为正整数且满足 u < l , x 1 , x 2 , ⋯ , x u 和 y 1 , y 2 , ⋯ , y u 均为 u 个 l 中两两互不相同的整数序列。当 x 1 , x 2 , ⋯ , x u 给定时,满足 y 1 ≠ x 1 , y 2 ≠ x 2 , ⋯ , y u ≠ x u 的 y 1 , y 2 , ⋯ , y u 的选取方法有 A l u - ∑ i = 1 u - 1 i - 1 u i A l - i u - i 种。
证明 记事件 B 表示从 l 中选取两两互不相同的 y 1 , y 2 , ⋯ , y u ,事件 B e ( 1 ≤ e ≤ u )表示从 l 中选取两两互不相同的 y 1 , y 2 , ⋯ , y u 且满足 y e = x e 。令 f · 表示事件 · 相应的不同选取方法的种数。此时,满足 y 1 ≠ x 1 , y 2 ≠ x 2 , ⋯ , y u ≠ x u 的 y 1 , y 2 , ⋯ , y u 的不同选取方法有 f B - f B 1 ⋃ B 2 ⋃ ⋯ ⋃ B u 种。由 B 的定义可知, f B = A l u 。由 B e ( 1 ≤ e ≤ u )的定义和容斥原理可知,
f B 1 ⋃ B 2 ⋃ ⋯ ⋃ B u = ∑ 1 ≤ i ≤ u f B i + ( - 1 ) ∑ 1 ≤ i < j ≤ u f B i ⋂ B j + ( - 1 ) 2 ∑ 1 ≤ i < j < q ≤ u f B i ⋂ B j ⋂ B q + ⋯ + - 1 u - 1 f B 1 ⋂ B 2 ⋂ ⋯ ⋂ B u = - 1 1 - 1 u 1 A l - 1 u - 1 + - 1 2 - 1 u 2 A l - 2 u - 2 + ⋯ + - 1 u - 1 u u A l - u 0 = ∑ i = 1 u - 1 i - 1 u i A l - i u - i 。
综上所述,当 x 1 , x 2 , ⋯ , x u 给定时,满足 y 1 ≠ x 1 , y 2 ≠ x 2 , ⋯ , y u ≠ x u 的 y 1 , y 2 , ⋯ , y u 的选取方法有 A l u - ∑ i = 1 u - 1 i - 1 u i A l - i u - i 种。引理2证毕。
由式(1) 可知,通过厘清两个 S n - m , k - m 子网络的不同相交情形来计算式(1) 右端前两项之和可以得出 R n , k m ( p ) 的下界的理论公式。
定理2 设 n , k , m 为正整数且有 1 ≤ k ≤ n - 1 和 1 ≤ m ≤ k - 1 ,令
R L = k - 1 m n ! / n - m ! p A n - m k - m - 1 2 ∑ i = 0 m - 1 k - 1 i k - 1 - i m - i k - 1 - m m - i A n m A n - m m - i p 2 A n - m k - m - A n - ( 2 m - i ) k - ( 2 m - i ) + ∑ i = 0 m k - 1 i k - 1 - i m - i k - 1 - m m - i ( ( A n m ) 2 - A n m A n - m m - i ) p 2 A n - m k - m ,
则 R n , k m ( p ) ≥ m a x R L , 0 。
证明 注意到 R n , k m ( p ) ≥ 0 。根据式(1) 和定理1的证明,
R n , k m ( p ) ≥ m a x k - 1 m n ! / n - m ! p A n - m k - m - ∑ 1 ≤ i < j ≤ k - 1 m n ! / n - m ! P r ( A i ⋂ A j ) , 0 。
下面证明 k - 1 m n ! / n - m ! p A n - m k - m - ∑ 1 ≤ i < j ≤ k - 1 m n ! / n - m ! P r ( A i ⋂ A j ) = R L 。根据引理1, S n , k 中任意两个 S n - m , k - m 子网络可表示为 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 和 X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m ,其中 a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 均为 m 个 n 中两两互不相同的整数。根据 i 1 , i 2 , ⋯ , i m 和 j 1 , j 2 , ⋯ , j m 的不同取值,讨论以下情形。
情形 t 1 ≤ t ≤ m - 1 i 1 , i 2 , ⋯ , i m ⋂ j 1 , j 2 , ⋯ , j m = t 。
此时显然有 2 m - t ≤ k - 1 。不妨设 i g 1 , i g 2 , ⋯ , i g t ∈ j 1 , j 2 , ⋯ , j m 。注意到 i g 1 , i g 2 , ⋯ , i g t 为 2,3 , ⋯ , k 中 t 个两两互不相同的整数,有 k - 1 t 种不同的选法。当 i g 1 , i g 2 , ⋯ , i g t 选定时, i 1 , i 2 , ⋯ , i m \ i g 1 , i g 2 , ⋯ , i g t 中 m - t 个元素可在 2,3 , ⋯ , k \ i g 1 , i g 2 , ⋯ , i g t 中选取,有 k - 1 - t m - t 种不同的选法; j 1 , j 2 , ⋯ , j m \ i g 1 , i g 2 , ⋯ , i g t 中 m - t 个元素可在 2,3 , ⋯ , k \ i 1 , i 2 , ⋯ , i m 中选取,有 k - 1 - m m - t 种不同的选法。由于 i 1 , i 2 , ⋯ , i m 和 j 1 , j 2 , ⋯ , j m 的选取不受选取顺序的影响,故该情形下 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 有 1 2 k - 1 t k - 1 - t m - t k - 1 - m m - t 种不同的选法。当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时,考虑以下子情形。
子情形 t - 1 - l 1 ≤ l ≤ m - t a i g 1 = b i g 1 , a i g 2 = b i g 2 , ⋯ , a i g t = b i g t 且
a i 1 , a i 2 , ⋯ , a i m \ a i g 1 , a i g 2 , ⋯ , a i g t ⋂ b j 1 , b j 2 , ⋯ , b j m \ b i g 1 , b i g 2 , ⋯ , b i g t = l 。
根据 S n , k 的子网络简化表示的定义可知,
V X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m ⋂ V X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m = ∅ ,
即这一情形下两个 S n - m , k - m 子网络共有 2 A n - m k - m 个顶点。
a i 1 , a i 2 , ⋯ , a i m 为 n 中两两互不相同的整数,有 A n m 种不同的选法;由 a i g 1 = b i g 1 , a i g 2 = b i g 2 , ⋯ , a i g t = b i g t 可知,当 a i g 1 , a i g 2 , ⋯ , a i g t 选定时 b i g 1 , b i g 2 , ⋯ , b i g t 是确定的。注意到 a i 1 , a i 2 , ⋯ , a i m \ a i g 1 , a i g 2 , ⋯ , a i g t 与 b j 1 , b j 2 , ⋯ , b j m \ b i g 1 , b i g 2 , ⋯ , b i g t 有 l 个元素相同,不妨设 a x 1 , a x 2 , ⋯ , a x l ∈ b j 1 , b j 2 , ⋯ , b j m \ b i g 1 , b i g 2 , ⋯ , b i g t 。 a x 1 , a x 2 , ⋯ , a x l 可在 a i 1 , a i 2 , ⋯ , a i m \ a i g 1 , a i g 2 , ⋯ , a i g t 中选取,有 m - t l 种不同的选法;当 a x 1 , a x 2 , ⋯ , a x l 选定时满足 b j 1 , b j 2 , ⋯ , b j m \ b i g 1 , b i g 2 , ⋯ , b i g t 中有 l 个元素与 a x 1 , a x 2 , ⋯ , a x l 相同的选法有 A m - t l 种。 b j 1 , b j 2 , ⋯ , b j m \ b j g 1 , b j g 2 , ⋯ , b j g t , a x 1 a x 2 , ⋯ , a x l 中 m - t - l 个元素可在 n \ a i 1 , a i 2 , ⋯ , a i m 中选取,有 A n - m m - t - l 种不同的选法。因此,当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时, a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 有 A n m m - t l A m - t l A n - m m - t - l 种不同的选法。综上可知,在子情形 t - 1 - l 下 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 和 X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m 共有 1 2 k - 1 t k - 1 - t m - t k - 1 - m m - t A n m m - t l A m - t l A n - m m - t - l 种不同的选法。
子情形 t - 1 - m - t + 1 a i g 1 = b i g 1 , a i g 2 = b i g 2 , ⋯ , a i g t = b i g t 且
a i 1 , a i 2 , ⋯ , a i m \ a i g 1 , a i g 2 , ⋯ , a i g t ⋂ b j 1 , b j 2 , ⋯ , b j m \ b i g 1 , b i g 2 , ⋯ , b i g t = 0 。
令 w 1 , w 2 , ⋯ , w 2 m - t = i 1 , i 2 , ⋯ , i m ⋃ j 1 , j 2 , ⋯ , j m 且满足 w 1 < w 2 < ⋯ < w 2 m - t , c w x 满足 c w x = a w x , w x ∈ i 1 , i 2 , ⋯ , i m b w x , w x ∈ j 1 , j 2 , ⋯ , j m 。根据 S n , k 的子网络简化表示的定义可知,
V X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m ⋂ V X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m = V X w 1 - 1 c w 1 X w 2 - w 1 - 1 c w 2 ⋯ c w 2 m - t X k - w 2 m - t ,
即这一情形下两个 S n - m , k - m 子网络共有 2 A n - m k - m - A n - 2 m - t k - 2 m - t 个顶点。
a i 1 , a i 2 , ⋯ , a i m 为 n 中两两互不相同的整数,有 A n m 种不同的选法;由 a i g 1 = b i g 1 , a i g 2 = b i g 2 , ⋯ , a i g t = b i g t 可知,当 a i g 1 , a i g 2 , ⋯ , a i g t 选定时 b i g 1 , b i g 2 , ⋯ , b i g t 是确定的。 b j 1 , b j 2 , ⋯ , b j m \ b i g 1 , b i g 2 , ⋯ , b i g t 中 m - t 个元素可在 n \ a i 1 , a i 2 , ⋯ , a i m 中选取,有 A n - m m - t 种不同的选法。因此,当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时, a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 有 A n m A n - m m - t 种不同的选法。综上可知,在子情形 t - 1 - m - t + 1 下 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 和 X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m 共有 1 2 k - 1 t k - 1 - t m - t k - 1 - m m - t A n m A n - m m - t 种不同的选法。
子情形 t - 2 - s 1 ≤ s ≤ t - 1 a i g 1 = b i g 1 , a i g 2 = b i g 2 , ⋯ , a i g t = b i g t 中有且仅有 s 个等式成立。
根据 S n , k 的子网络简化表示的定义可知,
V X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m ⋂ V X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m = ∅ ,
即这一情形下两个 S n - m , k - m 子网络共有 2 A n - m k - m 个顶点。
a i 1 , a i 2 , … , a i m 为 n 中两两互不相同的整数,有 A n m 种不同的选法。由于 a i g 1 = b i g 1 , a i g 2 = b i g 2 , ⋯ , a i g t = b i g t 中有且仅有 s 个等式成立,可知成立的 s 个等式有 t s 种不同的选法;不妨设 a i g 1 = b i g 1 , a i g 2 = b i g 2 , ⋯ , a i g s = b i g s , a i g s + 1 ≠ b i g s + 1 , a i g s + 2 ≠ b i g s + 2 , ⋯ , a i g t ≠ b i g t ,则当 a i g 1 , a i g 2 , ⋯ , a i g s 选定时 b i g 1 , b i g 2 , ⋯ , b i g s 是确定的;由引理2可知,当 a i g s + 1 , a i g s + 2 , ⋯ , a i g t 选定时,符合 b i g s + 1 ≠ a i g s + 1 , b i g s + 2 ≠ a i g s + 2 , ⋯ , b i g t ≠ a i g t 的 b i g s + 1 , b i g s + 2 , ⋯ , b i g t 的选取方式有 A n - s t - s - ∑ i = 1 t - s - 1 i - 1 t - s i A n - s - i t - s - i 种; b j 1 , b j 2 , ⋯ , b j m \ b i g 1 , b i g 2 , ⋯ , b i g t 中 m - t 个元素可在 n \ b i g 1 , b i g 2 , ⋯ , b i g t 中选取,有 A n - t m - t 种不同的选取方式。因此,当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时, a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 有 A n m t s A n - s t - s - ∑ i = 1 t - s - 1 i - 1 A n - s - i t - s - i A n - t m - t 种不同的选法。综上可知,在子情形 t - 2 - s 下 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 和 X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m 共有 1 2 k - 1 t k - 1 - t m - t k - 1 - m m - t A n m t s A n - s t - s - ∑ i = 1 t - s - 1 i - 1 t - s i A n - s - i t - s - i A n - t m - t 种不同的选法。
子情形 t - 2 - t a i g 1 ≠ b i g 1 , a i g 2 ≠ b i g 2 , ⋯ , a i g t ≠ b i g t 。
根据 S n , k 的子网络简化表示的定义可知,
V X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m ⋂ V X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m = ∅ ,
即这一情形下两个 S n - m , k - m 子网络共有 2 A n - m k - m 个顶点。
a i 1 , a i 2 , ⋯ , a i m 为 n 中两两互不相同的整数,有 A n m 种不同的选法。由引理2可知,当 a i g 1 , a i g 2 , ⋯ , a i g t 选定时,满足 b i g 1 ≠ a i g 1 , b i g 2 ≠ a i g 2 , ⋯ , b i g t ≠ a i g t 的 b i g 1 , b i g 2 , ⋯ , b i g t 的选取方法有 A n t - ∑ i = 1 t - 1 i - 1 t i A n - i t - i 种; b j 1 , b j 2 , ⋯ , b j m \ b i g 1 , b i g 2 , ⋯ , b i g t 中 m - t 个元素可在 n \ b i g 1 , b i g 2 , ⋯ , b i g t 中选取,有 A n - t m - t 种选法。因此,当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时, a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 有 A n m A n t - ∑ i = 1 t - 1 i - 1 t i A n - i t - i A n - t m - t 种不同的选法。综上可知,在子情形 t - 2 - t 下 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 和 X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m 共有 1 2 k - 1 t k - 1 - t m - t k - 1 - m m - t A n m A n t - ∑ i = 1 t - 1 i - 1 t i A n - i t - i A n - t m - t 种不同的选法。
情形 m i 1 , i 2 , ⋯ , i m ⋂ j 1 , j 2 , ⋯ , j m = m 。
此时有 m ≤ k - 1 ,且 i 1 = j 1 , i 2 = j 2 , ⋯ , i m = j m 。显然, a i 1 = b j 1 , a i 2 = b j 2 , ⋯ , a i m = b j m 不能同时成立,这意味着 V X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m ⋂ V X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m = ∅ ,即这一情形下两个 S n - m , k - m 子网络共有 2 A n - m k - m 个顶点。
注意到 i 1 , i 2 , ⋯ , i m 为 2,3 , ⋯ , k 中 m 个两两互不相同的整数,有 k - 1 m 种不同的选法;当 i 1 , i 2 , ⋯ , i m 选定时 j 1 , j 2 , ⋯ , j m 是确定的。当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时,考虑以下子情形。
子情形 m - s 1 ≤ s ≤ m - 1 a i 1 = b j 1 , a i 2 = b j 2 , ⋯ , a i m = b j m 中有且仅有 s 个等式成立。
a i 1 , a i 2 , ⋯ , a i m 为 n 中两两互不相同的整数,有 A n m 种不同的选法。由于 a i 1 = b j 1 , a i 2 = b j 2 , ⋯ , a i m = b j m 中有且仅有 s 个等式成立,可知成立的 s 个等式有 m s 种不同的选法;不妨设 a i 1 = b j 1 , a i 2 = b j 2 , ⋯ , a i s ≠ b j s , a i s + 1 ≠ b j s + 1 , a i s + 2 ≠ b j s + 2 , ⋯ , a i m ≠ b j m ,则当 a i 1 , a i 2 , ⋯ , a i s 选定时 b j 1 , b j 2 , ⋯ , b j s 是确定的;由引理2可知,当 a i s + 1 , a i s + 2 , ⋯ , a i m 选定时,满足 b j s + 1 ≠ a i s + 1 , b j s + 2 ≠ a i s + 2 , ⋯ , b j m ≠ a i m 的 b j s + 1 , b j s + 2 , ⋯ , b j m 的选取方法有 A n - s m - s - ∑ i = 1 m - s - 1 i - 1 m - s i A n - s - i m - s - i 种。又因为 a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 的选取不受选取顺序的影响,故当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时, a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 有 1 2 A n m m s A n - s m - s - ∑ i = 1 m - s - 1 i - 1 m - s i A n - s - i m - s - i 种不同的选取方式。综上可知,在子情形 m - s 下 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 和 X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m 共有 1 2 k - 1 m A n m m s A n - s m - s - ∑ i = 1 m - s - 1 i - 1 m - s i A n - s - i m - s - i 种不同的选法。
子情形 m - m a i 1 ≠ b j 1 , a i 2 ≠ b j 2 , ⋯ , a i m ≠ b j m 。
a i 1 , a i 2 , ⋯ , a i m 为 n 中两两互不相同的整数,有 A n m 种不同的选法;由引理2可知,当 a i 1 , a i 2 , ⋯ , a i m 选定时,符合 b j 1 ≠ a i 1 , b j 2 ≠ a i 2 , ⋯ , b j m ≠ a i m 的 b j 1 , b j 2 , ⋯ , b j m 的选取方式有 A n m - ∑ i = 1 m - 1 i - 1 m i A n - i m - i 种。注意到 a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 的选取不受选取顺序的影响。因此,当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时, a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 有 1 2 A n m A n m - ∑ i = 1 m - 1 i - 1 m i A n - i m - i 种不同的选取方式。综上可知,在子情形 m - m 下 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 和 X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m 共有 1 2 k - 1 m A n m A n m - ∑ i = 1 m - 1 i - 1 m i A n - i m - i 种不同的选法。
情形 m + 1 i 1 , i 2 , ⋯ , i m ⋂ j 1 , j 2 , ⋯ , j m = 0 。
此时有 2 m ≤ k - 1 ,且 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 两两互不相同。 i 1 , i 2 , ⋯ , i m 为 2,3 , ⋯ , k 中 m 个两两互不相同的整数,有 k - 1 m 种不同的选法; j 1 , j 2 , ⋯ , j m 可以在 2,3 , ⋯ , k \ i 1 , i 2 , ⋯ , i m 中选取,有 k - 1 - m m 种不同的选取方式。由于 i 1 , i 2 , ⋯ , i m 和 j 1 , j 2 , ⋯ , j m 的选取不受选取顺序的影响,故该情形下 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 有 1 2 k - 1 m k - 1 - m m 种不同的选法。当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时,考虑以下子情形。
子情形 m + 1 - l a i 1 , a i 2 , ⋯ , a i m ⋂ b j 1 , b j 2 , ⋯ , b j m = l 1 ≤ l ≤ m 。
根据 S n , k 的子网络简化表示的定义可知,
V X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m ⋂ V X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m = ∅ ,
即这一情形下两个 S n - m , k - m 子网络共有 2 A n - m k - m 个顶点。
a i 1 , a i 2 , ⋯ , a i m 为 n 中两两互不相同的整数,有 A n m 种不同的选法;注意到 a i 1 , a i 2 , ⋯ , a i m 与 b j 1 , b j 2 , ⋯ , b j m 有 l 个元素相同,不妨设 a x 1 , a x 2 , ⋯ , a x l ∈ b j 1 , b j 2 , ⋯ , b j m 。 a x 1 , a x 2 , ⋯ , a x l 可在 a i 1 , a i 2 , ⋯ , a i m 中选取,有 m l 种不同的选法;当 a x 1 , a x 2 , ⋯ , a x l 选定时,满足 b j 1 , b j 2 , ⋯ , b j m 中有 l 个元素与 a x 1 , a x 2 , ⋯ , a x l 相同的选法有 A m l 种。 b j 1 , b j 2 , ⋯ , b j m \ a x 1 , a x 2 , ⋯ , a x l 中 m - l 个元素可在 n \ a i 1 , a i 2 , ⋯ , a i m 中选取,有 A n - m m - l 种不同的选法。因此,当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时, a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 有 A n m m l A m l A n - m m - l 种不同的选法。综上可知,在子情形 m + 1 - l 下 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 和 X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m 共有 1 2 k - 1 m k - 1 - m m A n m m l A m l A n - m m - l 种不同的选法。
子情形 m + 1 - m + 1 a i 1 , a i 2 , ⋯ , a i m ⋂ b j 1 , b j 2 , ⋯ , b j m = 0 。
令 w 1 , w 2 , ⋯ , w 2 m = i 1 , i 2 , ⋯ , i m ⋃ j 1 , j 2 , ⋯ , j m 且满足 w 1 < w 2 < ⋯ < w 2 m , c w x 满足 c w x = a w x , w x ∈ i 1 , i 2 , ⋯ , i m b w x , w x ∈ j 1 , j 2 , ⋯ , j m 。根据 S n , k 的子网络简化表示的定义可知, V X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m ⋂ V X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m = V X w 1 - 1 c w 1 X w 2 - w 1 - 1 c w 2 ⋯ c w 2 m X k - w 2 m ,即这一情形下两个 S n - m , k - m 子网络共有 2 A n - m k - m - A n - 2 m k - 2 m 个顶点。
a i 1 , a i 2 , ⋯ , a i m 为 n 中两两互不相同的整数,有 A n m 种不同的选法; b j 1 , b j 2 , ⋯ , b j m 为 n \ a i 1 , a i 2 , ⋯ , a i m 中两两互不相同的整数,有 A n - m m 种选法。因此,当 i 1 , i 2 , ⋯ , i m , j 1 , j 2 , ⋯ , j m 选定时, a i 1 , a i 2 , ⋯ , a i m 和 b j 1 , b j 2 , ⋯ , b j m 有 A n m A n - m m 种不同的选法。综上可知,在子情形 m + 1 - m + 1 下 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 和 X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m 共有 1 2 k - 1 m k - 1 - m m A n m A n - m m 种不同的选法。
由上述分析可知,
i 1 , i 2 , ⋯ , i m 和
j 1 , j 2 , ⋯ , j m 的所有情形都已考虑,各个情形下
X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 和
X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m 的相交情况见
表1 。
因此,
∑ 1 ≤ i < j ≤ k - 1 m n ! / n - m ! P r ( A i ⋂ A j ) = 1 2 ∑ t = 1 m - 1 k - 1 t k - 1 - t m - t k - 1 - m m - t A n m A n - m m - t p 2 A n - m k - m - A n - ( 2 m - t ) k - ( 2 m - t ) + 1 2 k - 1 m k - 1 - m m A n m A n - m m p 2 A n - m k - m - A n - 2 m k - 2 m + 1 2 ∑ t = 1 m - 1 k - 1 t k - 1 - t m - t k - 1 - m m - t A n m 2 - 1 2 ∑ t = 1 m - 1 k - 1 t k - 1 - t m - t k - 1 - m m - t A n m A n - m m - t p 2 A n - m k - m + 1 2 k - 1 m A n m A n m - 1 p 2 A n - m k - m + 1 2 k - 1 m k - 1 - m m A n m 2 - 1 2 k - 1 m k - 1 - m m A n m A n - m m p 2 A n - m k - m = 1 2 ∑ i = 0 m - 1 k - 1 i k - 1 - i m - i k - 1 - m m - i A n m A n - m m - i p 2 A n - m k - m - A n - ( 2 m - i ) k - ( 2 m - i ) + ∑ i = 0 m k - 1 i k - 1 - i m - i k - 1 - m m - i ( ( A n m ) 2 - A n m A n - m m - i ) p 2 A n - m k - m ,
即有 k - 1 m n ! / n - m ! p A n - m k - m - ∑ 1 ≤ i < j ≤ k - 1 m n ! / n - m ! P r ( A i ⋂ A j ) = R L 成立。定理2证毕。
2.2 有效性分析
要得出
R n , k m ( p ) 的精确值需要厘清不同数目
S n - m , k - m 子网络的所有相交情况。然而,
S n - m , k - m 子网络之间的相交情况会随着子网络数目的增大变得十分复杂。因此,本文仅给出了
R n , k m ( p ) 的上下界。下面对2.1小节得出的
R n , k m ( p ) 的上下界进行有效性分析。记
R n , k m ( p ) ¯ = m i n k - 1 m n ! / n - m ! p A n - m k - m = R U , 1 ,
R n , k m ( p ) ̲ = m a x R L , 0 。以不同规模的
S n , k S 7,4 , S 8,4 , S 7,5 , S 7,6 作为实验对象,对
S n , k 中存在无故障
S n - m , k - m (
m = 2,3 )子网络的概率的上界
R n , k m ( p ) ¯ 和下界
R n , k m ( p ) ̲ 进行数值计算,分别得出了满足
R L ≥ 0 、
R U ≤ 1 和
R n , k m ( p ) ¯ - R n , k m ( p ) ̲ ≤ δ (
δ = 0.01 )的
p 的取值范围,具体结果见
表2 。
从
表2 可以看出,当
m 不变时,满足
R L ≥ 0 、
R U ≤ 1 和
R n , k m ( p ) ¯ - R n , k m ( p ) ̲ ≤ δ 的
p 的取值范围均随着
S n , k 的规模的增大而变得越来越大;当
S n , k 的规模给定时,
m 的值越小,满足
R L ≥ 0 、
R U ≤ 1 和
R n , k m ( p ) ¯ - R n , k m ( p ) ̲ ≤ δ 的
p 的取值范围越大。这意味着,对于较大规模的
S n , k ,顶点可靠性
p 在很大范围内取值都可以满足
R L ≥ 0 、
R U ≤ 1 和
R n , k m ( p ) ¯ - R n , k m ( p ) ̲ ≤ δ ;进一步地,当上界
R n , k m ( p ) ¯ 和下界
R n , k m ( p ) ̲ 的差值较小(如
R n , k m ( p ) ¯ - R n , k m ( p ) ̲ ≤ δ )时,可以将
1 2 ( R n , k m ( p ) ¯ + R n , k m ( p ) ̲ ) 看作
R n , k m ( p ) 的较为精确的估计值。
3 无故障 S n - m , k - m 子网络搜寻算法
3.1 算法描述
给出在带有点故障的 S n , k 中搜寻无故障 S n - m , k - m 子网络的算法有助于基于 S n , k 构建的多处理器系统的任务分配算法的设计,具有一定的现实意义。
算法1 无故障 S n - m , k - m 子网络搜寻算法
输入:(n ,k )-星图 S n , k 和 S n - m , k - m 子网络的规模参数 n , k , m ( 1 ≤ k ≤ n - 1 , 1 ≤ m ≤ k - 1 )以及 S n , k 中一个故障点集 F 。
输出: S n , k 中无故障 S n - m , k - m 子网络的集合 S 。
1. 令 S ≔ X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m : 2 ≤ i 1 < i 2 < ⋯ < i m ≤ k , a i 1 , a i 2 , ⋯ , a i m ∈ n 且它 们两 两互 不相 同
2. while F ≠ ∅ do
3. 选择一个点 b 1 b 2 ⋯ b k ∈ F
4. S ≔ S \ X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m : 2 ≤ j 1 < j 2 < ⋯ < j m ≤ k
5. F ≔ F \ b 1 b 2 ⋯ b k
6. end while
7. 返回 S 。
算法1的正确性可由以下讨论得出:
由引理1, S n , k 中任一个 S n - m , k - m 子网络可用某个 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m 唯一表示,其中 1 ≤ m ≤ k - 1 , 2 ≤ i 1 < i 2 < ⋯ < i m ≤ k 。可知, S n , k 中所有 S n - m , k - m 子网络的集合可表示为 X i 1 - 1 a i 1 X i 2 - i 1 - 1 a i 2 ⋯ a i m X k - i m : 2 ≤ i 1 < i 2 < ⋯ < i m ≤ k , a i 1 , a i 2 , ⋯ , a i m ∈ n 且它 们两 两互 不相 同 。根据 S n , k 的子网络简化表示的定义,点 b 1 b 2 ⋯ b k ∈ V S n , k 发生故障会使得 X j 1 - 1 b j 1 X j 2 - j 1 - 1 b j 2 ⋯ b j m X k - j m : 2 ≤ j 1 < j 2 < ⋯ < j m ≤ k 中 k - 1 m 个不同的 S n - m , k - m 子网络遭到破坏。
3.2 无故障 S n - m , k - m 子网络搜寻算法实例
本节使用MATLAB实现算法1。分别以
S 5,3 和
S 5,4 为例,在随机生成不同的故障点集
F 之后,利用算法1搜寻
S n , k 中不同规模的无故障
S n - m , k - m 子网络(
m = 1,2 ),部分结果见
表3 和
表4 。由实验结果可以看出,若仅考虑点故障的发生且
S n , k 中所有故障点已知,则利用算法1可以搜寻出
S n , k 中所有的无故障
S n - m , k - m 子网络。
4 基于蒙特卡罗的 R n , k m ( p ) 的近似评估方法
4.1 算法描述
当基于 S n , k 构建的多处理器系统需指派 S n - m , k - m 子网络执行用户任务时,相应算法的设计依赖于精确的无故障 S n - m , k - m 子网络的存在概率估计结果。注意到,当 R n , k m ( p ) 的上下界的差值较大时,利用上界 R n , k m ( p ) ¯ 和下界 R n , k m ( p ) ̲ 不能较为准确地评估 R n , k m ( p ) 。下面给出一个基于蒙特卡罗的 R n , k m ( p ) 的近似计算方法。
算法2 R n , k m ( p ) 的基于蒙特卡罗的近似评估算法
输入:(n ,k )-星图 S n , k 和 S n - m , k - m 子网络的规模参数 n , k , m ( 1 ≤ k ≤ n - 1 , 1 ≤ m ≤ k - 1 )、顶点可靠性 p 和模拟次数 N 。
输出: S n , k 中无故障 S n - m , k - m 子网络的存在概率 R n , k m ( p ) 。
1. 把 S n , k 中 n ! / n - k ! 个顶点分别记为 T [ 0 ] , T [ 1 ] , ⋯ , T [ n ! / n - k ! - 1 ]
2. for i = 1 to N do
3. 生成[0,1]内服从均匀分布的 n ! / n - k ! 个随机数 γ 0 , γ 1 , ⋯ , γ n ! / n - k ! - 1
4. set F ≔ ∅
5. for j = 0 to n ! / n - k ! - 1 do
6. if γ j < 1 - p then
7. F ≔ F ⋃ T [ j ]
8. end if
9. end for
10. 以规模参数 n , k , m 及故障点集 F 为输入调用算法1,将得到的无故障 S n - m , k - m 子网络的集合记为 S
11. if S > 0 then
12. set x i ≔ 1
13. else
14. set x i ≔ 0
15. end if
16. end for
17. set R n , k m ( p ) ≔ ∑ i = 1 N x i / N
18. 返回 R n , k m ( p ) 。
算法2的正确性可由以下讨论得出:
将生成的[0,1]内服从均匀分布的 n ! / n - k ! 个随机数 γ 0 , γ 1 , ⋯ , γ n ! / n - k ! - 1 与 S n , k 的顶点 T [ 0 ] , T [ 1 ] , ⋯ , T [ n ! / n - k ! - 1 ] 一一对应,把符合 γ j < 1 - p 的对应顶点记成故障点,剩余顶点记成非故障点,则这一网络状态下 S n , k 的顶点可靠性为 p 。在该网络状态下借助算法1可以得到 S n , k 中无故障 S n - m , k - m 子网络的集合 S ;重复以上模拟过程N 次,若存在无故障 S n - m , k - m 子网络(即 S > 0 )的累计模拟次数为 ∑ i = 1 N x i ,则此时 ∑ i = 1 N x i / N 可以看作 R n , k m ( p ) 的估计值,且模拟次数N 越大该估计值越精确。
4.2 实验分析
以不同规模的
S n , k S 7,4 , S 8,4 , S 7,5 , S 7,6 作为实验对象,本节借助MATLAB实现算法2对
S n , k 中无故障
S n - m , k - m 子网络(
m = 2,3 )的存在概率进行估计,并将所得的近似评估结果与定理1和定理2给出的上下界进行对比。对于某一选定的
S n , k ,若
S n , k 的顶点可靠性
p 在[0,1]内给定时,执行算法2(设置模拟次数
N = 10 000 )对
S n , k 中无故障
S n - m , k - m 子网络的存在概率
R n , k m ( p ) 进行估算,并根据定理1和定理2对
R n , k m ( p ) 的上界
R n , k m ( p ) ¯ 和下界
R n , k m ( p ) ̲ 进行计算,具体实验结果如
图2 所示。
通过
图2 可以看出,对于任意的
S n , k ,随着顶点可靠性
p 逐渐减小,利用算法2得出的
R n , k m ( p ) 的近似评估结果与
R n , k m ( p ) 的上下界基本吻合,这表明当上界
R n , k m ( p ) ¯ 和下界
R n , k m ( p ) ̲ 的差满足精度要求时,利用上下界的平均值
1 2 ( R n , k m ( p ) ¯ + R n , k m ( p ) ̲ ) 和算法2都能得出较为精确的
R n , k m ( p ) ;当上下界之间的差值不满足精度要求时,利用上界
R n , k m ( p ) ¯ 和下界
R n , k m ( p ) ̲ 无法较为准确的估计
R n , k m ( p ) ,此时利用算法2得出的结果是较为精确的(模拟次数
N 越大结果越精确)。
5 结论
对于 1 ≤ k ≤ n - 1 和 1 ≤ m ≤ k - 1 ,本文在概率故障条件下估计了 S n , k 中无故障 S n - m , k - m 子网络的存在概率 R n , k m ( p ) ,得出了 R n , k m ( p ) 的上下界,给出了在仅考虑点故障的 S n , k 中搜寻无故障 S n - m , k - m 子网络的算法,并基于蒙特卡罗模拟提出了 R n , k m ( p ) 的近似评估方法。结果表明,当 R n , k m ( p ) 的上下界的差值满足评估精度时,上界和下界的平均值可以看作 R n , k m ( p ) 的较为精确的估计值,且该平均值与近似评估方法得出的结果是一致的;当 R n , k m ( p ) 的上下界的差值不满足评估精度时,利用基于蒙特卡罗的近似方法可以得出较为精确的评估结果。
值得注意的是,为了达到预期评估精度,基于蒙特卡罗的无故障 S n - m , k - m 子网络存在概率的近似评估方法往往需设置较大的模拟次数,这可能会导致评估效率大大降低,其原因如下:当 S n , k 的顶点可靠性为 p 时,算法2的每次模拟过程中生成的故障点集 F 满足 | F | ≈ n ! ( n - k ) ! ( 1 - p ) (这是由于 | V ( S n , k ) | = n ! ( n - k ) ! 且 1 - | F | | V ( S n , k ) | 可以看作 p 的估计值),这表明在算法2的步骤10中调用算法1时算法1的步骤2—步骤6循环体的执行次数约为 n ! ( n - k ) ! ( 1 - p ) ;若 n ! ( n - k ) ! ( 1 - p ) 和模拟次数 N 均较大时,利用算法2估计 R n , k m ( p ) 的效率会较低。高效的无故障子网络存在概率的近似评估方法值得进一步研究。