《算法导论》第7章 快速排序 个人笔记

第7章 快速排序

7.1 快速排序

void swap(int *a, int *b)  
{  
    int tmp = *a;  
    *a = *b;  
    *b = tmp;  
}  

int partition(int a[], int low, int high)  
{  
    int privotKey = a[low];
    while(low < high){
        while(low < high  && a[high] >= privotKey) --high;
        swap(&a[low], &a[high]);  
        while(low < high  && a[low] <= privotKey ) ++low;  
        swap(&a[low], &a[high]);  
    }
    return low;  
}

void quickSort(int a[], int low, int high){  
    if(low < high){  
        int privotLoc = partition(a, low, high);
        quickSort(a, low, privotLoc -1);
        quickSort(a, privotLoc + 1, high);
    }  
}  

7.2 快速排序的性能

快速排序的运行时间依赖于划分是否平衡。
最好情况：partition得到的两个子问题的规模都不大于n/2，递归式为

$$T(n)=2T(n/2)+\Theta(n)$$
根据主定理，上述递归式的解为 $T(n)=\Theta(n\lg n)$。

7.3 快排的随机版本

int randomized_partition(int a[], int low, int high){
    int i = random(low, high);
    swap(&a[low], &a[i]);
    partition(a, low, high);
}

void randomized_quickSort(int a[], int low, int high){
    if(low < high){  
        int privotLoc = randomized_partition(a, low, high);
        randomized_quickSort(a, low, privotLoc -1);
        randomized_quickSort(a, privotLoc + 1, high);
    } 
}

7.4 快排分析

最坏情况：

$$T(n)=max_{\{0\leq q\leq n-1\}}(T(n)+T(n-q-1))+\Theta(n)$$
猜测 $T(n)\leq cn^2$成立，代入上式，得：
$$T(n)=max_{\{0\leq q\leq n-1\}}(cq^2+c(n-q-1)^2)+\Theta(n)\\ =c*max_{\{0\leq q\leq n-1\}}(q^2+(n-q-1)^2)+\Theta(n)\\ \leq c*(n-1)^2/2+\Theta(n)\leq cn^2$$

期望运行时间：

定义：$z_i$是数组A中地i小的元素，$Z_{ij}=\{z_i,z_{i+1},...,z_j\}$，$X_{ij}=I\{z_i~compare~with~z_j \}$，$X=\sum_{i=1}^{n-1}\sum_{j=i+1}^{n}X_{ij}$

$$E[X]=\sum_{i=1}^{n-1}\sum_{j=i+1}^{n}E[X_{ij}]=\sum_{i=1}^{n-1}\sum_{j=i+1}^{n}P(z_i~compare~with~z_j)$$
$z_i$与 $z_j$会进行比较，当且仅当 $Z_{ij}$中将被选为主元的第一个元素是 $z_i$或 $z_j$。于是有 $P(z_i~compare~with~z_j)=\frac{2}{j-i+1}$，所以有
$$E[X]=\sum_{i=1}^{n-1}\sum_{j=i+1}^{n}\frac{2}{j-i+1}=\sum_{i=1}^{n-1}\sum_{k=1}^{n-i}\frac{2}{k+1} \\<\sum_{i=1}^{n-1}\sum_{k=1}^{n}\frac{2}{k}=\sum_{i=1}^{n-1}O(\lg n)=O(n\lg n)$$