算法导论学习笔记（1）——快排中hoarePartition的实现（问题已解决）

打算在学期结束之前（1月19日）看完，目前94页，总共六百页，不完全做习题及思考题（不细致考虑各种时间）。

之所以发表这个是公开一下，算是对自己的监督吧。以下进入正题

实现各种算法时都OK，在快速排序这里第一次遇到问题，算法出错，hoarePartition是快速排序中用到的一个比较古老的分块算法。

伪代码

HOARE-PARTITION(A,p,r)

{

　　x=A[p]

　　i=p-1

　　j=r+1

　　while(1)

　　{

　　　　do repeat j--

　　　　　　until A[j]<=x

　　　　     repeat i++

　　　　　　until A[i]>=x

　　　　if i<j

　　　　　　then exchange A[i] and A[j]

　　　　　　else return j

　　}

}

思维比较清晰的算法，当同时找到坐标x两端的异常值时，左右互换

C代码

 1 #include <stdio.h>
 2 
 3 void RandomDoSort(int A[],int begin,int end);
 4 int Hoare_Partition(int A[],int begin,int end);
 5 
 6 void ChangeTwoInt(int* a,int* b);
 7 
 8 void main()
 9 {
10     int A[12]={13,19,9,5,12,8,7,4,21,2,6,11};
11     RandomDoSort(A,0,11);
12     for (int i=0;i<12;i++)
13     {
14         printf("%d\n",A[i]);
15     }
16 }
17 
18 
19 void RandomDoSort(int A[],int begin,int end)
20 {
21     int q=0;
22     if(begin<end)
23     {
24         q=Hoare_Partition(A,begin,end);
25         DoSort(A,begin,q-1);
26         DoSort(A,q+1,end);
27     }
28 }
29 
30 int Hoare_Partition(int A[],int begin,int end)
31 {
32     int xtemp=A[begin];
33     int i=begin-1;
34     int j=end+1;
35     while (1)
36     {
37         j--;
38         if (A[j]<=xtemp)
39         {
40             while(A[i]<xtemp)
41             {
42                 i++;
43             }
44         }
45         if (i<j)
46         {
47             ChangeTwoInt(&A[i],&A[j]);
48         }
49         else
50         {
51             return j;
52         }
53     }
54 }
55 
56 void ChangeTwoInt(int* a,int* b)
57 {
58     int temp=*a;
59     *a=*b;
60     *b=temp;
61 }

最后结果出错，有一个数没排序上，睡觉去，明早起来调试解决这问题

 其实问题并不出在Partition这个函数里，而是在Dosort函数里。上面这个Dosort函数简单地修改了之前另一个Partition搭配的Dosort，两者思想有细微区别。

先来看之前的Partition。

int Partition(int A[],int begin,int end)
{
    int x=A[end];
    int i=begin-1;
    int temp=0;
    for(int j=begin;j<=end-1;j++)
    {
        if(A[j]<x)
        {
            i++;
            temp=A[j];
            A[j]=A[i];
            A[i]=temp;
        }    
    }
    temp=A[i+1];　　　　　　　　　　//将坐标轴元素放在排序好的位置上
    A[i+1]=A[end];
    A[end]=temp;
    
    return i+1;
}

再来看HoarePartition

int Hoare_Partition(int A[],int begin,int end)
{
    int xtemp=A[begin];
    int i=begin-1;
    int j=end+1;
    while (1)
    {
        j--;
        if (A[j]<=xtemp)
        {
            while(A[i]<xtemp)
            {
                i++;
            }
        }
        if (i<j)
        {
            ChangeTwoInt(&A[i],&A[j]);
        }
        else
        {
            return j;
        }
    }
}

两者的区别是：前面那种在Partition过程中就将坐标轴元素放在排序好的位置上，自然在DOsort的迭代过程中不用再考虑该元素，但Hoare则没有这样做，在Hoare匹配的Dosort迭代时自然要将该元素也纳入迭代的考虑范围。

来看两个Dosort

之前那一种

void DoSort(int A[],int begin,int end)
{
    int q=0;
    if(begin<end)
    {
        q=Partition(A,begin,end);
        DoSort(A,begin,q-1);　　　　//A[q]不需要再处理
        DoSort(A,q+1,end);
    }
}

Hoare对应的Dosort

void HoareDoSort(int A[],int begin,int end)
{
    int q=0;
    if(begin<end)
    {
        q=Hoare_Partition(A,begin,end);　　　　　　
        DoSort(A,begin,q);　　　　　　　　　　//需要对A[q]再处理
        DoSort(A,q+1,end);
    }
}

 

以上。