并行算法的设计与分析15.pptVIP

并行算法的设计与分析 第 15 章 图论并行算法 15.2 图的传递闭包 15.2.1 传递闭包 1. 定义 给定有向图 G=(V,E), A =( a ij ) n × n 为 G 的邻接矩阵 , A 的传递闭包 A + =( b ij ) n × n , b ij 为 1 时当且仅当顶点 i 到顶点 j 之间有一条路径： 2. 串行传递闭包算法时间复杂度 : O( n 3 ) 示例 其他 有路径 到 或 j i j i b ij ? ? ? ? ? 0 1 0 1 0 0 A = 0 0 1 0 0 0 0 1 1 0 0 0 1 4 3 2 1 1 1 1 A + = 1 1 1 1 1 1 1 1 1 1 1 1 15.2 图的传递闭包 3. 基于布尔矩阵乘积的图传递闭包算法原理 令 B = A +I, I 为单位阵 , B =( b (1) ij ) n × n , b (1) ij =1 ? i = j 或 i ? j 有有向边 ? i ? j 有长为 0 或 1 的有向路径 . 对于 B 2 =( A +I) 2 =( b (2) ij ) n × n , b (2) ij = ∨ k =1~ n b (1) ik b (1) kj , ∨表示“逻辑或” 则有 , b (2) ij =1 ? 顶点 i ? 顶点 j 有长度 ≤ 2 的有向路径 . 类似地 , B k =( b ( k ) ij ) n × n , b ( k ) ij =1 ? i ? j 有长度 ≤ k 的有向路径 . 因此 , A + = B r , r ≥ n -1 所以 , B ? B 2 ? B 4 ? … ? B 2 = A + , 共有 次相乘 ? ? ) 1 log( ? n ? ? ) 1 log( ? n 15.2 图的传递闭包 4. 示例 1 4 3 2 (a) 1 1 0 0 0 1 1 0 0 0 1 1 1 0 0 1 (b) B=A+I= 1 1 1 0 0 1 1 1 1 0 1 1 1 1 0 1 (c) B 2 = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 (d) B 4 = =A + 15.2 图的传递闭包 5. SIMD-CC 上的传递闭包算法 ? 求图传递闭包等价于求图的连通性 ( 连通矩阵 ), 用 (A+I) 自乘 log( n -1) 次获得 ? 将 n 3 个处理器排成 n × n × n 的三维阵列 , P r 的坐标为 ( i , j , k ), r = i * n 2 + j * n + k ， P r 有三个寄存器 A ( i , j , k ), B ( i , j , k ) 和 C ( i , j , k ), 0≤ i , j , k ≤ n -1 初始时 : A (0, j , k ) ? a jk ， 0≤ j , k ≤ n -1; 结束时 : C (0, j , k ) 为 A + 的 ( j , k ) 元素 ? 算法 15.1 SIMD-CC 上的传递闭包（连通矩阵）算法 Begin // 输入 : A n × n , 输出 : C n × n (1) for j=0 to n-1 do in parallel A(0, j, j)=1 endfor // 形成 A