快排序的非遞歸實現（原創）

本文簡要介紹快排序的非遞歸實現方式並給出了C語言版本的源代碼，非遞歸實現由於沒有函數調用的消耗，相對於遞歸方式有一定的優勢。
在快排序的遞歸實現中，函數quicksort(int nlow, int nhigh)遞歸的時候兩個參和high在改變。

非遞歸實現中可以用兩個棧int low[MAX_SIZE] , int high[MAX_SIZE]來模擬係統棧對nlow和nhigh的

存儲情況。非遞歸快排的具體實現比較簡單，需要說明的是，MAX_SIZE這裏等於32就可以了，因為

快速排序的空間複雜度是O(log2N)，而整數最大時2^32，因此兩個模擬棧的大小是32就可以了。

非遞歸快排序的源碼如下：

/*
 * Author: puresky
 * Date: 2010/01/08
 * Purpose: Nonrecursive quicksort
 */

#include <stdio.h>
#include <time.h>
#include <stdlib.h>

//Nonrecursive quicksort
void quicksort2(int a[], int n)
{
#define STACK_MAX 32 //the space complexity of quicksort is O(logN), so 32 is enough!
#define SWAP(a, b) {int temp = a; a = b; b = temp;}

        int low[STACK_MAX]; //low position stack
        int high[STACK_MAX]; // high position stack
        int top = -1; //stack top

        top++;
        low[top] = 0;
        high[top] = n - 1;

        while(top >= 0) // stack is not empty
        {
                int L = low[top];
                int R = high[top];
                top--;

                //printf("TOP:%d,L:%d,R:%d\n", top, L, R);

                if( L < R)
                {
                        //partion
                        SWAP(a[L], a[L + rand() % ( R - L + 1)]);
                       
                        int i;
                        int m = L;
                        for(i = L + 1; i <= R; ++i)
                        {       
                                if(a[i] < a[L])
                                {
                                        ++m;
                                        SWAP(a[m], a[i]);
                                }
                        }
                        SWAP(a[m], a[L]);
                       
                       //left part entering stack
                        if(L < m - 1)
                        {
                                top++;
                                low[top] = L;
                                high[top] = m - 1;
                        }
 
                        //rigth part entering stack
                        if(m + 1 < R)
                        {
                                top++;
                                low[top] = m + 1;
                                high[top] = R;
                        }
                }
        }
}

void print(int a[], int n)
{
        int i;
        for(i = 0; i < n; ++i)
                printf("%d ", a[i]);
        printf("\n");
}

int* init_array(int n)
{
        int i;
        int *a = (int*)malloc(sizeof(int) * n);
        for(i = 0; i < n; ++i)
                a[i] = rand() % n;
        return a;
}

void free_array(int *a)
{
        free(a);
}

int main(int argc, char** argv)
{
        srand(time(NULL));
        const int n = 100;
        int *a = init_array(n);
        quicksort2(a, n);
        print(a, n);
        free_array(a);
        system("pause");
        return 0;
}