算法与并行编程.doc

算法和并行编程 （Algorithms） 1. 在下面三个问题中任选一个，给出相应的并行求解算法，具体说明使用的数据结构和分解方法，并分析算法的时间复杂度。 a)图的连通分支（connected components of graphs） b)图的极大独立集（Maximal independent set） c)单源/全源最短路（single (all pair) shortest path）。 2. 在下面三个问题中任选一个，用MPI，Pthread或OpenMP中的一种编写出程序。 a)矩阵相乘的并行程序 b)π计算 c)求数组的最小值或最大值。 算法和并行编程 （Algorithms）答案 1. 在下面三个问题中任选一个，给出相应的并行求解算法，具体说明使用的数据结构和分解方法，并分析算法的时间复杂度。 a)图的连通分支（connected components of graphs） 答：对图进行深度优先遍历可以找出连通分支： 并行求解算法：将图G的邻接矩阵划分为p个部分，每个部分分配给一个进程。 第一步，每个进程Pi计算子图Gi的深度优先生成森林。 第二步，生成森林按对合并，直到只剩下一个生成森林。 合并方法：假设将A合并到B，对A中的每条边（u，v），确定B中两个顶点是否处在同一棵树上，如果不是，则将两棵树合并。 算法的时间复杂度：Θ(n2/p) ------------------------------------------------------------------------------------------------------------------- 答：Perform DFS on the graph to get Connected Components: The n * n adjacency matrix is partitioned into p blocks. Partition the graph across processors and run independent connected component algorithms on each processor. second, spanning forests are merged pairwise until only one spanning forest remains. Each processor can compute its local spanning forest in time Θ(n2/p). Merging is done by embedding a logical tree into the topology. There are log p merging stages, and each takes time Θ(n). Thus, the cost due to merging is Θ(n log p). During each merging stage, spanning forests are sent between nearest neighbors. Recall that Θ(n) edges of the spanning forest are transmitted. The parallel run time of the connected-component algorithm is For a cost-optimal formulation p = O(n / log n). b)图的极大独立集（Maximal independent set） 答：图的极大独立集算法的思想： Luby算法：计算一个图的顶点V的最大独立集I： 集合I初始设置为空，候选顶点集合C=V 一个唯一的随机数分配给C中每一个顶点，如果某个顶点的随机数小于所有与它相邻顶点的随机数，则将它加入到I中；更新集合C，除去加入到I中的顶点以及与它们相邻的所有顶点 重复②至C为空 并行求解算法： 使用三个数组，一个存储进入最大独立集I中的节点，一个存储候选集C，一个存储随机数R。将候选集C在p个进程中划分。每个进程为从C中分配给它的顶点产生一个随机数。当所有顶点都产生随机数后，就可以来确定哪些顶点包含在I中。 算法的时间复杂度： O(log|V|) ------------------------------------------------------------------------------------------------------------------- 答：Maximal independent set算法的思想：We us