快速排序的若干源码剖析(包括JAVA、Qt、算法导论)

 快排：快速排序，难度不是太大，但是是我经常把玩的小代码之一，引用快排代码创始人的一句话“要不是简练又明显没有错误，要不是复杂没有明显的错误”

 

起因：

主要是前一阵看刘汝佳的《算法艺术与信息学竞赛》P73页，碰巧我又亲手写了一遍那书上的测试了一下，发现竟然是错误的，然后我又测试了一下我以前写的快排的代码，发现也是错误的ORZ。然后我就翻了一下能找到的快排的源码。包括Qt、JAVA、算法导论，顺手用了一下系统的快排

简单说下

系统的快排 很快

Qt            绝对是另类的代码

JAVA        很符合一般人能想到的优化

算法导论   简练又明显没有错误

 

刘汝佳的《算法艺术与信息学竞赛》P73页

 

/* *选自刘汝佳的《算法艺术与信息学竞赛》P73页，此为书上的算法，经测试是 错误的快排 */ #include<stdio.h> long a[5]; long p(long l,long r){//可以发现这里面有点错误 long x=a[l],i=l-1,j=r+1,t; while(1){ do{//这种语句会出现很棘手的情况，使不该划分的被跳过去了 --j; } while(a[j]>x); do{ ++i; } while(a[i]>x); if(i<j){ t=a[i];a[i]=a[j];a[j]=t; } else{ return j; } } } void q(long l,long r){ long n; if(l<r){ n=p(l,r); q(l,n); q(n+1,r); } } int main(){ a[0]=5 ; a[1]=3 ; a[2]=2 ; a[3]=4 ; a[4]=1 ; q(0,4); int i; for(i=0;i<5;++i){ printf("%i ",a[i]); } return 0; }

 

直接用系统快排的速度惊人，并且也是编程的好习惯

/* *使用系统的快排,测试时会出现惊人的快 */ #include<cstdio> #include<cstdlib> using namespace std; int compare(const void*, const void*); int main() { int numsLength=5; int *nums=new int[numsLength]; nums[0]=5; nums[1]=7; nums[2]=2; nums[3]=4; nums[4]=1; for (int i=0;i<numsLength;++i) { printf("%i ",nums[i]); } printf("/n"); qsort(nums,numsLength,sizeof(nums[0]),compare); for (int i=0;i<numsLength;++i) { printf("%i ",nums[i]); } return 0; } int compare(const void* a, const void* b){ return *((const int *)a)-*((const int *)b); }

 

Qt的，不得不说，GOTO语句渐渐退出了人们的视线，但是Qt的快排从代码上来看，性能应该不错

 

/* *Qt4.6.1的排序，相当与众不同，当然，是正确的快排 */ void qSortHelper(RandomAccessIterator start, RandomAccessIterator end, const T &t, LessThan lessThan) { top: int span = int(end - start); if (span < 2) return; --end; RandomAccessIterator low = start, high = end - 1; RandomAccessIterator pivot = start + span / 2; if (lessThan(*end, *start)) qSwap(*end, *start); if (span == 2) return; if (lessThan(*pivot, *start)) qSwap(*pivot, *start); if (lessThan(*end, *pivot)) qSwap(*end, *pivot); if (span == 3) return; qSwap(*pivot, *end); while (low < high) { while (low < high && lessThan(*low, *end)) ++low; while (high > low && lessThan(*end, *high)) --high; if (low < high) { qSwap(*low, *high); ++low; --high; } else { break; } } if (lessThan(*low, *end)) ++low; qSwap(*end, *low); qSortHelper(start, low, t, lessThan); start = low + 1; ++end; goto top; }

 

JAVA的快排，优化很多，很JAVA

 

/* *JDK1.6的排序，主要是快排，长度为7一下使用冒泡，7以上使用快排，分割点是用多点取中的方法,当然，是正确的快排 */ #include<cstdio> using namespace std; void swap(long x[], int a, int b) ; void sort1(long x[], int off, int len); void vecswap(long x[], int a, int b, int n); int med3(long x[], int a, int b, int c); int main() { int numsLength=5; long *nums=new long[numsLength]; nums[0]=5; nums[1]=7; nums[2]=2; nums[3]=4; nums[4]=1; sort1(nums, 0, numsLength); for (int i=0;i<numsLength;++i) { printf("%i ",nums[i]); } return 0; } void sort1(long x[], int off, int len) { if (len < 7) { for (int i=off; i<len+off; i++) for (int j=i; j>off && x[j-1]>x[j]; j--) swap(x, j, j-1); return; } // Choose a partition element, v int m = off + (len >> 1); // Small arrays, middle element if (len > 7) { int l = off; int n = off + len - 1; if (len > 40) { // Big arrays, pseudomedian of 9 int s = len/8; l = med3(x, l, l+s, l+2*s); m = med3(x, m-s, m, m+s); n = med3(x, n-2*s, n-s, n); } m = med3(x, l, m, n); // Mid-size, med of 3 } long v = x[m]; // Establish Invariant: v* (<v)* (>v)* v* int a = off, b = a, c = off + len - 1, d = c; while (true) { while (b <= c && x[b] <= v) { if (x[b] == v) swap(x, a++, b); b++; } while (c >= b && x[c] >= v) { if (x[c] == v) swap(x, c, d--); c--; } if (b > c) break; swap(x, b++, c--); } // Swap partition elements back to middle int s, n = off + len; s = a-off< b-a?a-off:b-a ; vecswap(x, off, b-s, s); s = d-c< n-d-1?d-c: n-d-1; vecswap(x, b, n-s, s); // Recursively sort non-partition-elements if ((s = b-a) > 1) sort1(x, off, s); if ((s = d-c) > 1) sort1(x, n-s, s); } /** * Swaps x[a] with x[b]. */ void swap(long x[], int a, int b) { long t = x[a]; x[a] = x[b]; x[b] = t; } /** * Swaps x[a .. (a+n-1)] with x[b .. (b+n-1)]. */ void vecswap(long x[], int a, int b, int n) { for (int i=0; i<n; i++, a++, b++) swap(x, a, b); } /** * Returns the index of the median of the three indexed longs. */ int med3(long x[], int a, int b, int c) { return (x[a] < x[b] ? (x[b] < x[c] ? b : x[a] < x[c] ? c : a) : (x[b] > x[c] ? b : x[a] > x[c] ? c : a)); }

 

压轴的《算法导论》，不会令你失望的，代码相当简练，在随机情况下性能也可以

 

 /* *选自《算法导论》P85页，此为书上的算法，经测试是 正确的快排，并且代码长度相当短，在其他算法相关的书上竟几乎没有再出现 */ #include<cstdio> using namespace std; int quicksort(int *nums,int s,int e); int main() { int numsLength=5; int *nums=new int[numsLength]; nums[0]=5; nums[1]=7; nums[2]=2; nums[3]=4; nums[4]=1; for (int i=0;i<numsLength;++i) { printf("%i ",nums[i]); } printf("/n"); quicksort(nums,0,numsLength-1); for (int i=0;i<numsLength;++i) { printf("%i ",nums[i]); } return 0; } int quicksort(int *nums,int s,int e) { if (s>=e) { return 0; } int x=nums[e];//此处可以找个随机点与最后一个点交换下值，实现随机化 int cutPoint=s-1; int temp; for (int i=s;i<e;++i) {//注意此处把最后一个点给空过去了 if (nums[i]<=x) { ++cutPoint; temp=nums[cutPoint];nums[cutPoint]=nums[i];nums[i]=temp; } } temp=nums[cutPoint+1];nums[cutPoint+1]=nums[e];nums[e]=temp;//最后才把分割值这个点给放到分割点了 quicksort(nums,s,cutPoint); quicksort(nums,cutPoint+2,e); return 0; }