排序算法一:快速排序

快速排序的第一种实现(单指针移动,挖空填数)

快速排序由于排序效率在同为O(N*logN)的几种排序方法中效率较高,因此经常被采用,再加上快速排序思想----分治法也确实实用,因此很多软件公司的笔试面试,包括像腾讯,微软等知名IT公司都喜欢考这个,还有大大小的程序方面的考试如软考,考研中也常常出现快速排序的身影。

总的说来,要直接默写出快速排序还是有一定难度的,因为本人就自己的理解对快速排序作了下白话解释,希望对大家理解有帮助,达到快速排序,快速搞定

 

快速排序是C.R.A.Hoare于1962年提出的一种划分交换排序。它采用了一种分治的策略,通常称其为分治法(Divide-and-ConquerMethod)。

该方法的基本思想是:

1.先从数列中取出一个数作为基准数。

2.分区过程,将比这个数大的数全放到它的右边,小于或等于它的数全放到它的左边。

3.再对左右区间重复第二步,直到各区间只有一个数。

 

虽然快速排序称为分治法,但分治法这三个字显然无法很好的概括快速排序的全部步骤。因此我的对快速排序作了进一步的说明:挖坑填数+分治法

先来看实例吧,定义下面再给出(最好能用自己的话来总结定义,这样对实现代码会有帮助)。

 

以一个数组作为示例,取区间第一个数为基准数。

0

1

2

3

4

5

6

7

8

9

72

6

57

88

60

42

83

73

48

85

 

初始时,i = 0;  j = 9;   X = a[i] = 72

由于已经将a[0]中的数保存到X中,可以理解成在数组a[0]上挖了个坑,可以将其它数据填充到这来。

从j开始向前找一个比X小或等于X的数。当j=8,符合条件,将a[8]挖出再填到上一个坑a[0]中。a[0]=a[8]; i++;  这样一个坑a[0]就被搞定了,但又形成了一个新坑a[8],这怎么办了?简单,再找数字来填a[8]这个坑。这次从i开始向后找一个大于X的数,当i=3,符合条件,将a[3]挖出再填到上一个坑中a[8]=a[3]; j--;

  

数组变为:

0

1

2

3

4

5

6

7

8

9

48

6

57

88

60

42

83

73

88

85

 

i = 3;   j = 7;   X=72

再重复上面的步骤,先从后向前找,再从前向后找

从j开始向前找,当j=5,符合条件,将a[5]挖出填到上一个坑中,a[3] = a[5]; i++;

从i开始向后找,当i=5时,由于i==j退出。

此时,i = j = 5,而a[5]刚好又是上次挖的坑,因此将X填入a[5]。

  

数组变为:

0

1

2

3

4

5

6

7

8

9

48

6

57

42

60

72

83

73

88

85

 

可以看出a[5]前面的数字都小于它,a[5]后面的数字都大于它。因此再对a[0…4]和a[6…9]这二个子区间重复上述步骤就可以了。

  

  

对挖坑填数进行总结

1.i =L; j = R; 将基准数挖出形成第一个坑a[i]。

2.j--由后向前找比它小的数,找到后挖出此数填前一个坑a[i]中。

3.i++由前向后找比它大的数,找到后也挖出此数填到前一个坑a[j]中。

4.再重复执行2,3二步,直到i==j,将基准数填入a[i]中。

照着这个总结很容易实现挖坑填数的代码:

int AdjustArray(int s[], int l, int r) //返回调整后基准数的位置   
{   
    int i = l, j = r;   
    int x = s[l]; //s[l]即s[i]就是第一个坑   
    while (i < j)   
    {   
        // 从右向左找小于x的数来填s[i]   
        while(i < j && s[j] >= x)    
            j--;     
        if(i < j)    
        {   
            s[i] = s[j]; //将s[j]填到s[i]中,s[j]就形成了一个新的坑   
            i++;   
        }   
    
        // 从左向右找大于或等于x的数来填s[j]   
        while(i < j && s[i] < x)   
            i++;     
        if(i < j)    
        {   
            s[j] = s[i]; //将s[i]填到s[j]中,s[i]就形成了一个新的坑   
            j--;   
        }   
    }   
    //退出时,i等于j。将x填到这个坑中。   
    s[i] = x;   
    
    return i;   
}   

再写分治法的代码:

void quick_sort1(int s[], int l, int r)   
{   
    if (l < r)   
    {   
        int i = AdjustArray(s, l, r);//先成挖坑填数法调整s[]   
        quick_sort1(s, l, i - 1); // 递归调用    
        quick_sort1(s, i + 1, r);   
    }   
}   

这样的代码显然不够简洁,对其组合整理下:

//快速排序   
void quick_sort(int s[], int l, int r)   
{   
    if (l < r)   
    {   
        //Swap(s[l], s[(l + r) / 2]); //将中间的这个数和第一个数交换 参见注1   
        int i = l, j = r, x = s[l];   
        while (i < j)   
        {   
            while(i < j && s[j] >= x) // 从右向左找第一个小于x的数   
                j--;     
            if(i < j)    
                s[i++] = s[j];   
                
            while(i < j && s[i] < x) // 从左向右找第一个大于等于x的数   
                i++;     
            if(i < j)    
                s[j--] = s[i];   
        }   
        s[i] = x;   
        quick_sort(s, l, i - 1); // 递归调用    
        quick_sort(s, i + 1, r);   
    }   
}   

 

快速排序还有很多改进版本,如随机选择基准数,区间内数据较少时直接用另的方法排序以减小递归深度。有兴趣的筒子可以再深入的研究下。

 

注1,有的书上是以中间的数作为基准数的,要实现这个方便非常方便,直接将中间的数和第一个数进行交换就可以了。

快速派速的另一种实现(两个指针移动,交换的方法)

高快省的排序算法

有没有既不浪费空间又可以快一点的排序算法呢?那就是“快速排序”啦!光听这个名字是不是就觉得很高端呢。

假设我们现在对“6  1  2 7  9  3  4  5 10  8”这个10个数进行排序。首先在这个序列中随便找一个数作为基准数(不要被这个名词吓到了,就是一个用来参照的数,待会你就知道它用来做啥的了)。为了方便,就让第一个数6作为基准数吧。接下来,需要将这个序列中所有比基准数大的数放在6的右边,比基准数小的数放在6的左边,类似下面这种排列:

3  1  2 5  4  6  9 7  10  8

在初始状态下,数字6在序列的第1位。我们的目标是将6挪到序列中间的某个位置,假设这个位置是k。现在就需要寻找这个k,并且以第k位为分界点,左边的数都小于等于6,右边的数都大于等于6。想一想,你有办法可以做到这点吗?

排序算法显神威

方法其实很简单:分别从初始序列“6  1  2 7  9  3  4  5 10  8”两端开始“探测”。先从右往左找一个小于6的数,再从左往右找一个大于6的数,然后交换他们。这里可以用两个变量i和j,分别指向序列最左边和最右边。我们为这两个变量起个好听的名字“哨兵i”和“哨兵j”。刚开始的时候让哨兵i指向序列的最左边(即i=1),指向数字6。让哨兵j指向序列的最右边(即=10),指向数字。

094811yilrz1tkzkvlrriz.png

首先哨兵j开始出动。因为此处设置的基准数是最左边的数,所以需要让哨兵j先出动,这一点非常重要(请自己想一想为什么)。哨兵j一步一步地向左挪动(即j--),直到找到一个小于6的数停下来。接下来哨兵i再一步一步向右挪动(即i++),直到找到一个数大于6的数停下来。最后哨兵j停在了数字5面前,哨兵i停在了数字7面前。

095430axy0qkhxxkktkktk.png

095437kdandfxhbtokk2qh.png

现在交换哨兵i和哨兵j所指向的元素的值。交换之后的序列如下:

6  1  2  5  9 3  4  7  10  8

095448k1kevwlz41373e7k.png

095458ejza15wscjv7iw5c.png

到此,第一次交换结束。接下来开始哨兵j继续向左挪动(再友情提醒,每次必须是哨兵j先出发)。他发现了4(比基准数6要小,满足要求)之后停了下来。哨兵i也继续向右挪动的,他发现了9(比基准数6要大,满足要求)之后停了下来。此时再次进行交换,交换之后的序列如下:

6  1  2 5  4  3  9  7 10  8

第二次交换结束,“探测”继续。哨兵j继续向左挪动,他发现了3(比基准数6要小,满足要求)之后又停了下来。哨兵i继续向右移动,糟啦!此时哨兵i和哨兵j相遇了,哨兵i和哨兵j都走到3面前。说明此时“探测”结束。我们将基准数6和3进行交换。交换之后的序列如下:

3  1 2  5  4  6  9 7  10  8

095506uz7e1uuukcblhkxv.png

095514cag5fumuqqg5jnsw.png

095530e0jf6p0y6aaaw2ir.png

到此第一轮“探测”真正结束。此时以基准数6为分界点,6左边的数都小于等于6,6右边的数都大于等于6。回顾一下刚才的过程,其实哨兵j的使命就是要找小于基准数的数,而哨兵i的使命就是要找大于基准数的数,直到i和j碰头为止。

OK,解释完毕。现在基准数6已经归位,它正好处在序列的第6位。此时我们已经将原来的序列,以6为分界点拆分成了两个序列,左边的序列是“3  1 2  5  4”,右边的序列是“9  7  10  8”。接下来还需要分别处理这两个序列。因为6左边和右边的序列目前都还是很混乱的。不过不要紧,我们已经掌握了方法,接下来只要模拟刚才的方法分别处理6左边和右边的序列即可。现在先来处理6左边的序列现吧。

左边的序列是“3  1  2 5  4”。请将这个序列以3为基准数进行调整,使得3左边的数都小于等于3,3右边的数都大于等于3。好了开始动笔吧

如果你模拟的没有错,调整完毕之后的序列的顺序应该是:

2  1  3  5  4

OK,现在3已经归位。接下来需要处理3左边的序列“2 1”和右边的序列“5 4”。对序列“2 1”以2为基准数进行调整,处理完毕之后的序列为“1 2”,到此2已经归位。序列“1”只有一个数,也不需要进行任何处理。至此我们对序列“2 1”已全部处理完毕,得到序列是“1 2”。序列“5 4”的处理也仿照此方法,最后得到的序列如下:

1  2  3 4  5  6 9  7  10  8

对于序列“9  7  10  8”也模拟刚才的过程,直到不可拆分出新的子序列为止。最终将会得到这样的序列,如下

1  2  3 4  5  6  7  8 9  10

到此,排序完全结束。细心的同学可能已经发现,快速排序的每一轮处理其实就是将这一轮的基准数归位,直到所有的数都归位为止,排序就结束了。下面上个霸气的图来描述下整个算法的处理过程。

232129ogop8gk0r8y7l70k.png

这是为什么呢?
快速排序之所比较快,因为相比冒泡排序,每次交换是跳跃式的。每次排序的时候设置一个基准点,将小于等于基准点的数全部放到基准点的左边,将大于等于基准点的数全部放到基准点的右边。这样在每次交换的时候就不会像冒泡排序一样每次只能在相邻的数之间进行交换,交换的距离就大的多了。因此总的比较和交换次数就少了,速度自然就提高了。当然在最坏的情况下,仍可能是相邻的两个数进行了交换。因此快速排序的最差时间复杂度和冒泡排序是一样的都是O(N2),它的平均时间复杂度为O(NlogN)。其实快速排序是基于一种叫做“二分”的思想。我们后面还会遇到“二分”思想,到时候再聊。先上代码,如下

  1. #include <stdio.h>
  2. int a[101],n;//定义全局变量,这两个变量需要在子函数中使用
  3. void quicksort(int left,int right)
  4. {
  5.     int i,j,t,temp;
  6.     if(left>right)
  7.        return;
  8.                                
  9.     temp=a[left]; //temp中存的就是基准数
  10.     i=left;
  11.     j=right;
  12.     while(i!=j)
  13.     {
  14.                    //顺序很重要,要先从右边开始找
  15.                    while(a[j]>=temp && i<j)
  16.                             j--;
  17.                    //再找右边的
  18.                    while(a[i]<=temp && i<j)
  19.                             i++;
  20.                    //交换两个数在数组中的位置
  21.                    if(i<j)
  22.                    {
  23.                             t=a[i];
  24.                             a[i]=a[j];
  25.                             a[j]=t;
  26.                    }
  27.     }
  28.     //最终将基准数归位
  29.     a[left]=a[i];
  30.     a[i]=temp;
  31.                             
  32.     quicksort(left,i-1);//继续处理左边的,这里是一个递归的过程
  33.     quicksort(i+1,right);//继续处理右边的 ,这里是一个递归的过程
  34. }
  35. int main()
  36. {
  37.     int i,j,t;
  38.     //读入数据
  39.     scanf("%d",&n);
  40.     for(i=1;i<=n;i++)
  41.                    scanf("%d",&a[i]);
  42.     quicksort(1,n); //快速排序调用
  43.                             
  44.     //输出排序后的结果
  45.     for(i=1;i<=n;i++)
  46.         printf("%d ",a[i]);
  47.     getchar();getchar();
  48.     return 0;
  49. }
可以输入以下数据进行验证

1061279345108

运行结果是

12345678910

涨姿势环节

快速排序由 C. A. R. Hoare(东尼霍尔,Charles Antony Richard Hoare)在1960年提出,之后又有许多人做了进一步的优化。如果你对快速排序感兴趣可以去看看东尼霍尔1962年在Computer Journal发表的论文“Quicksort”以及《算法导论》的第七章。快速排序算法仅仅是东尼霍尔在计算机领域才能的第一次显露,后来他受到了老板的赏识和重用,公司希望他为新机器设计一个新的高级语言。你要知道当时还没有PASCAL或者C语言这些高级的东东。后来东尼霍尔参加了由Edsger Wybe Dijkstra(1972年图灵奖得主,这个大神我们后面还会遇到的到时候再细聊)举办的“ALGOL 60”培训班,他觉得自己与其没有把握去设计一个新的语言,还不如对现有的“ALGOL 60”进行改进,使之能在公司的新机器上使用。于是他便设计了“ALGOL 60”的一个子集版本。这个版本在执行效率和可靠性上都在当时“ALGOL 60”的各种版本中首屈一指,因此东尼霍尔受到了国际学术界的重视。后来他在“ALGOL X”的设计中还发明了大家熟知的“case”语句,后来也被各种高级语言广泛采用,比如PASCAL、C、Java语言等等。当然,东尼霍尔在计算机领域的贡献还有很多很多,他在1980年获得了图灵奖。

完整测试代码:

// 7-2.快速排序.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
#include <iostream>
#include <algorithm>
using namespace std;

void quick_sort1(int s[],int l,int r)
{
	if (l<r)
	{
		int i=l,j=r,x=s[l];
		while (i<j)
		{
			while (i<j&&s[j]>=x)
				j--;
			if(i<j)
				s[i++]=s[j];
			while(i<j&&s[i]<x)
				i++;
			if(i<j)
				s[j--]=s[i];
		}
		s[i]=x;
		quick_sort1(s,l,i-1);
		quick_sort1(s,i+1,r);
	}
}

void quick_sort2(int s[],int l,int r)
{
	if (l<r)
	{
		int i=l,j=r,x=s[l];
		while (i<j)
		{
			while(i<j&&s[j]>x)
				j--;
			while(i<j&&s[i]<=x)
				i++;
			if (i<j)
				swap(s[i],s[j]);
		}
		s[l]=s[i];
		s[i]=x;
		quick_sort2(s,l,i-1);
		quick_sort2(s,i+1,r);
	}
}
int _tmain(int argc, _TCHAR* argv[])
{
	int a[]={9,12,17,30,50,20,60,65,4,49};
	int b[]={9,12,17,30,50,20,60,65,4,49};
	quick_sort1(a,0,9);
	for (int i=0;i<10;i++)
	{
		cout<<a[i]<<" ";
	}
	cout<<endl;
	quick_sort2(b,0,9);
	for (int i=0;i<10;i++)
	{
		cout<<b[i]<<" ";
	}
	cout<<endl;
	return 0;
}
posted @ 2015-03-23 14:07  程序员修练之路  阅读(872)  评论(0编辑  收藏  举报