用MPI实现并行排序算法——奇偶交换排序(C++)

排序算法是算法领域的经典问题，如何在一个分布式系统上实现并行算法呢？

串行排序算法

冒泡排序是最经典的排序算法之一，算法按对比较元素的大小，a[0]与a[1]比较、a[1]与a[2]比较，以此类推，如果顺序不对就调换位置；每一次遍历数组，就可以将序列中待排序的最大值移动到最右边；n次遍历后就可以使数组有序

该算法虽然简单，但仔细想想，其固有的串行特征让其很难并行化，接下来我们考虑冒泡排序的变种之一——奇偶交换排序，该算法更加适合并行化

奇偶交换排序，其奇偶是指的奇数阶段和偶数阶段两个不同的阶段，在奇数阶段按下列顺序比较大小，如果顺序不对就调整位置：

在偶数阶段按下列顺序比较大小：

下面举一个例子，假设有如下序列：5 9 4 3

阶段0：比较和交换 (5,9)和(4,3)，获得序列 5 9 3 4

阶段1：比较和交换 (9,3)，获得序列5 3 9 4

阶段2：比较和交换 (5,3)和(9,4)，获得序列3 5 4 9

阶段3：比较和交换 (5,4)，获得序列 3 4 5 9

上述例子用了4个阶段将4个元素的序列有序化，实际操作中可能需要的次数更少

定理：对于n个元素的序列，作为奇偶交换排序的输入，那么至多经过n个阶段后，该序列一定能排好序

下面用伪代码来描述一下串行奇偶交换排序的算法流程：

for(int phase = 0; phase < n; phase++){
    if(phase % 2 == 0){        // 偶阶段，对0、1和2、3...等进行交换
        for(int i=0; i<n; i+=2){
            if(a[i-1] > a[i]){
                temp = a[i];
                a[i] = a[i-1];
                a[i-1] = temp;
            }
        }
    }
    else{                      // 奇阶段，对1、2和3、4...等进行交换 
        for(int i=1; i<n-1; i+=2){
            if(a[i] > a[i+1]){
                temp = a[i];
                a[i] = a[i+1];
                a[i+1] = temp;
            }
        }
    }
}

并行排序算法

可以很清楚的知道，奇偶交换排序更合适串行化，因为每个阶段的比较和交换操作是可以同时进行的

对于n个序列分配到p个进程的排序任务而言，假设n可以整除p；那么每个进程都分配到了多个数据，此时不再是简单的执行逆序对的交换操作了，(假设我们需要增序排序），而是将两个进程中的所有元素较小的一半保留至编号较小的进程，较大的一半数据保留至编号较大的进程

如何实现上述操作？第一个想法是排序，将需要交换数据的两个进程的数据合并起来排序，较小的前一半数据保留在编号较小的进程，较大的一半保留至编号较大的进程；如果原始的数据有序，那就可以使合并数据的效率更高，因此我们多了一个任务，在奇偶排序算法开始之前，将各个进程分到的数据使用串行的排序算法，如快排，借助c语言的sort函数可以快速实现

我们总结一下并行排序算法的流程：

1. 将各个进程中的数据用串行算法排序

2. 分奇偶阶段交换数据&