寻找两个等长的有序序列的中位数

王道课后习题-3

题目：一个长度为L(L>=1)的升序序列S，处在第[L/2]个位置的数称为S的中位数。例如，若序列S1=(11,13,15,17,19),则S1的中位数是15，两个序列的中位数就是含它们所有元素的升序序列的中位数。例如，若S2=(2,4,6,8,20)，则S1和S2的中位数是11。现在有两个等长升序序列A和B，试设计一个在时间和空间两方面都尽可能高效的算法，找出两个序列A和B的中位数。要求：

1）给出算法的基本设计思想。

2）根据设计思想，采用C或C++或Java语言描述算法，关键之处给出注释。

3）说明你所设计算法的时间复杂度和空间复杂度。

解答：

1）算法的基本设计思想如下：

        采用分治法，分别求两个升序序列A、B的中位数，设为a和b，求序列A、B的中位数过程如下：

①若a=b，则a或b即为所求中位数，算法结束。

②若a<b，则舍弃序列A中较小的一半，同时舍弃B中较大的一半，要求两次舍弃的长度相等。

③若a>b，则舍弃序列A中较大的一半，同时舍弃B中较小的一半，要求两次舍弃的长度相等。

在保留的两个升序序列中，重复过程①、②、③，直到两个序列中均只含一个元素时为止，较小者即为所求的中位数。

2）使用C语言描述算法如下：

int M_Search(int A[],int B[],int n){
    int s1,d1,m1,s2,d2,m2;
    s1=0;d1=n-1;
    s2=1;d2=n-1;   
    while(s1 != d1 || s2 != d2 ){    //即L>1时，当A，B中各剩一个数时退出循环
        m1=(s1+d1)/2;
        m2=(s2+d2)/2;
        if(A[m1]==B[m2]){
            return A[m1];            //满足设计思路里的条件若a=b时
        }
        if(A[m1]<B[m2]){             //若a<b时
            if((s1+d1)%2==0){        //若元素个数为奇数时
                s1=m1;
                d2=m2;
            }else{                   //若元素个数为偶数时
                s1=m1+1;             //从下中位开始
                d2=m2;
            }
        }elae{                       //若a>b时
            if((s1+d1)%2==0){        //若元素个数为奇数时
                d1=m1;
                s2=m2;
            }else{                   //若元素个数为偶数时
                d1=m1;
                s2=m2+1;             //从下中位开始
            }
        }
    }
    return A[s1]<B[s2]?A[s1]:B[s2];  //比较二者大小输出较小的数
}
            

        这段代码主要就是分情况讨论，先按算法设计里的分成三种a和b比大小的情况，再在不相等的情况中根据元素个数奇偶情况来讨论，主要需要注意的是判断A、B的舍弃部分，以及起始索引在元素个数为偶数时舍弃较小部分那个序列的起始索引应从下中位开始，也就是我们向下取整求出的中位+1。

3）算法的时间复杂度为O()，空间复杂度为O(1)。