• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

计数排序(C语言版本)

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

让我们来谈谈数的排序思维

       计数排序假定待排序的全部元素都是介于0到K之间的整数。计数排序使用一个额外的数组countArray。当中第i个元素是待排序数组array中值等于i的元素的个数。然后依据数组countArray来将array中的元素排到正确的位置。

 算法的步骤例如以下:

  1. 找出待排序的数组中最大和最小的元素
  2. 统计数组中每一个值为i的元素出现的次数,存入数组countArray的第i
  3. 对全部的计数累加(从countArray中的第一个元素開始,每一项和前一项相加)
  4. 反向填充目标数组:将每一个元素i放在新数组的第countArray[i]项。每放一个元素就将countArray[i]减去1

稳定性和复杂度:

       计数排序是稳定的排序算法;平均时间复杂度、最优时间复杂度和最差时间复杂度都为O(n+k),空间复杂度为O(n+k)。当中,n为待排元素个数,k为待排元素的范围(0~k)。

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

#define RANDMAX 1000000
#define RANDMIN 900000

void getRandArray(int array[], int size);
void countSort(int array[], int size);
void isSorted(int array[], int size);

int main(int argc, char const *argv[])
{
    int size = 0;
    scanf("%d", &size);
    assert(size > 0);

    int *array = (int *)calloc(size, sizeof(int));
    getRandArray(array, size);

    clock_t begin;
    clock_t end;

    begin = clock();
    countSort(array, size);
    end = clock();
    //打印排序所花费的时间。在linux下单位为ms
    printf("%ld\n", (end - begin) / 1000);

    isSorted(array, size);
    free(array);

    return 0;
}


//利用伪随机数填充数组array。伪随机数的范围在RANDMIN~RANDMAX-1之间
void getRandArray(int array[], int size)
{
    assert(array != NULL && size > 0);

    srand((unsigned) time(NULL));
    int i = 0;
    for (i = 0; i < size; ++i) {
        array[i] = rand() % (RANDMAX - RANDMIN) + RANDMIN ;
    }
}

//从小到大进行排序
void countSort(int array[], int size)
{
    assert(array != NULL && size > 0);

    //计数数组,用于统计数组array中各个不同数出现的次数
    //由于数组array中的数属于0...RANDMAX-1之间,所以countArray的大小要够容纳RANDMAX个int型的值
    int *countArray = (int *) calloc(RANDMAX, sizeof(int));
    //用于存放已经有序的数列
    int *sortedArray = (int *) calloc(size, sizeof(int));

    //统计数组array中各个不同数出现的次数,循环结束后countArray[i]表示数值i在array中出现的次数
    int index = 0;
    for (index = 0; index < size; ++index) {
        countArray[array[index]]++;
    }

    //统计数值比index小的数的个数。循环结束之后countArray[i]表示array中小于等于数值i的元素个数
    for (index = 1; index < RANDMAX; ++index) {
        countArray[index] += countArray[index - 1];
    }

    for (index = size - 1; index >= 0; --index) {
        //由于数组的起始下标为0。所以要减一
        sortedArray[countArray[array[index]] - 1] = array[index];
        //这里减一是为了保证当有多个值为array[index]的元素时,后出现的值为array[index]的元素
        //放在后面,也就是为了保证排序的稳定性
        --countArray[array[index]];
    }

    memcpy(array, sortedArray, size * sizeof(int));
    free(sortedArray);
    free(countArray);
}

//推断数组array是否已经是有序的
void isSorted(int array[], int size)
{
    assert(array != NULL && size > 0);

    int unsorted = 0;
    int i = 0;
    for (i = 1; i < size; ++i) {
        if (array[i] < array[i - 1]) {
            unsorted = 1;
            break;
        }
    }

    if (unsorted) {
        printf("the array is unsorted!\n");
    } else {
        printf("the array is sorted!\n");
    }
}

       计数排序是非比較的排序算法,据说其排序速度要快于不论什么的比較排序算法(我还未验证,可是在排序100000000个10000以内的数时花费为606毫秒,而C语言的qsort函数则为10802毫秒,由此可见一斑),因为计数排序须要一个计数数组以及一个存放有序数列的数组。故计数排序对内存的要求比較高。

------------------------------------------------------------------更新------------------------------------------------------------------

2014年10月20日

以上的实现有两个问题

       1、比較耗费内存:假设元素值的范围在9百万~1千万之间。上面的实现中countArray的大小须要1千万,而实际出现的不同值的元素就仅仅有1千万-9百万=1百万,也就是说用上面的实现会导致countArray中有90%的空间没有被利用,这个是非常浪费内存空间的!

尽管这个样例比較极端可是在详细应用面前还是须要注意!

       2、不能出现负数:假设元素值中有负数。就不能用上面的实现方法了,由于counArray的下标是待排数组的元素值,所以上面的实现不能排序包括负数的数列。

解决方法:

       1、第一个问题好解决:仅仅要首先一次遍历待排数组。找出最大值和最小值(眼下我所知的比較快的方法是时间复杂度为O(1.5n)的方法)。然后分配大小为(max-min+1)*sizeof(int)的空间给countArray即可了。至于统计的方法就得改变一点点了:

//从小到大进行排序,节省空间的版本号
void countSort2(int array[], int size)
{
    assert(array != NULL && size > 0);

    int min;
    int max;
    getMinAndMax(array, size, &min, &max);
    // printf("max - min + 1 = %d\n", max - min + 1);
    int countArraySize = max - min + 1;

    //计数数组,用于统计数组array中各个不同数出现的次数
    int *countArray = (int *) calloc(countArraySize, sizeof(int));
    //用于存放已经有序的数列
    int *sortedArray = (int *) calloc(size, sizeof(int));

    //统计数组array中各个不同数出现的次数。循环结束后countArray[i]表示数值i+min在array中出现的次数
    int index = 0;
    for (index = 0; index < size; ++index) {
        countArray[array[index] - min]++;
    }

    //统计数值比index小的数的个数。循环结束之后countArray[i]表示array中小于等于数值i+min的元素个数
    for (index = 1; index < countArraySize; ++index) {
        countArray[index] += countArray[index - 1];
    }

    for (index = size - 1; index >= 0; --index) {
        //由于数组的起始下标为0。所以要减一
        sortedArray[countArray[array[index] - min] - 1] = array[index];
        //这里减一是为了保证当有多个值为array[index]的元素时。后出现的值为array[index]的元素
        //放在后面,也就是为了保证排序的稳定性
        --countArray[array[index] - min];
    }

    memcpy(array, sortedArray, size * sizeof(int));
    free(sortedArray);
    free(countArray);
}

//待优化
void getMinAndMax(int array[], int size, int *min, int *max)
{
    assert(array != NULL && size > 0);
    if (size == 1) {
        *min = *max = array[0];
        return ;
    }

    *min = array[0];
    int index;
    for (index = 1; index < size; ++index) {
        if (array[index] < *min) {
            *min = array[index];
        }
    }
    *max = array[0];
    for (index = 1; index < size; ++index) {
        if (array[index] > *max) {
            *max = array[index];
        }
    }
}

       这样的方法能够改进在元素值范围比較小的数列的排序时间(在大部分情况时间较短)和空间复杂度。因为仅仅是为了验证想法的可行性,就没有在查找最值的函数上下功夫优化。等以后有时间了。再进一步优化。

       2、第二个问题:刚刚在洗漱的时候忽然想到了一个可行的方法:先一次遍历待排数组,找出负数中的最大值和最小值。正数中的最大值和最小值。创建两个计数数组,negativeCountArray用来统计待排数组中各个不同负数的出现个数,positiveCountArray用来统计待排数组中各个正数出现的个数,先统计总的负数的个数。然后再统计各个不同正数的个数,然后在往sortedArray中放元素的时候对正数和负数差别对待。这样就能让计数排序对包括负数的数列进行排序了。

代码例如以下:

//从小到大进行排序。节省空间,可排序负数的版本号
void countSort3(int array[], int size)
{
    assert(array != NULL && size > 0);

    int negativeMin;
    int negativeMax;
    int positiveMin;
    int positiveMax;
    _getMinAndMax(array, size, &negativeMin, &negativeMax, &positiveMin, &positiveMax);
    int positiveCountArraySize = positiveMax - positiveMin + 1;
    int negativeCountArraySize = negativeMax - negativeMin + 1;

    //计数数组,用于统计数组array中各个不同数出现的次数
    int *positiveCountArray = (int *) calloc(positiveCountArraySize, sizeof(int));
    int *negativeCountArray = (int *) calloc(negativeCountArraySize, sizeof(int));
    //用于存放已经有序的数列
    int *sortedArray = (int *) calloc(size, sizeof(int));

    //统计数组array中各个不同的数出现的次数,循环结束后positiveCountArray[i]表示数值i+positiveMin>=0在array中出现的次数
    //negativeCountArray[i]表示数值i+negativeMin<0在array中出现的次数
    int index = 0;
    for (index = 0; index < size; ++index) {
        if (array[index] < 0) {
            negativeCountArray[array[index] - negativeMin]++;
        } else {
            positiveCountArray[array[index] - positiveMin]++;
        }
    }

    //统计数值比index+negativeMin小的负数的个数。循环结束之后negativeCountArray
    for (index = 1; index < negativeCountArraySize; ++index) {
        negativeCountArray[index] += negativeCountArray[index - 1];
    }
    int negativeCount = negativeCountArray[index - 1];//负数的总数

    //统计数值比index小的非负数数的个数,循环结束之后positiveCountArray[i]表示array中小于等于数值i的值为非负数的元素个数
    for (index = 1; index < positiveCountArraySize; ++index) {
        positiveCountArray[index] += positiveCountArray[index - 1];
    }

    for (index = size - 1; index >= 0; --index) {
        if (array[index] < 0) {
            sortedArray[negativeCountArray[array[index] - negativeMin] - 1] = array[index];
            --negativeCountArray[array[index] - negativeMin];
        } else {
            //因为非负数总是比负数来的大,所以在把正数放到sortedArray时要在下标上加上负数的总个数
            sortedArray[positiveCountArray[array[index] - positiveMin] - 1 + negativeCount] = array[index];
            --positiveCountArray[array[index] - positiveMin];
        }
    }

    memcpy(array, sortedArray, size * sizeof(int));
    free(sortedArray);
    free(positiveCountArray);
    free(negativeCountArray);
}

//待优化
void _getMinAndMax(int array[], int size, int *negativeMin, int *negativeMax, int *positiveMin, int *positiveMax)
{
    assert(array != NULL && size > 0);
    if (size == 1) {
        *negativeMin = *negativeMax = *positiveMin = *positiveMax = array[0];
        return ;
    }

    int firstNegative = 0;
    int firstPositive = 0;
    int index;
    for (index = 0; index < size; ++index) {
        if (array[index] < 0) {
            firstNegative = array[index];
        }
    }
    for (index = 0; index < size; ++index) {
        if (array[index] > 0) {
            firstPositive = array[index];
        }
    }

    *negativeMin = firstNegative;
    *positiveMin = firstPositive;
    for (index = 0; index < size; ++index) {
        if (array[index] < 0 && array[index] < *negativeMin) {
            *negativeMin = array[index];
        }
        if (array[index] > 0 && array[index] < *positiveMin) {
            *positiveMin = array[index];
        }
    }
    *negativeMax = firstNegative;
    *positiveMax = firstPositive;
    for (index = 0; index < size; ++index) {
        if (array[index] < 0 && array[index] > *negativeMax) {
            *negativeMax = array[index];
        }
        if (array[index] > 0 && array[index] > *positiveMax) {
            *positiveMax = array[index];
        }
    }
}

       因为仅仅是为了验证想法的可行性,就没有在查找最值的函数上下功夫优化,等以后有时间了。再进一步优化。

其它八种排序算法的博客:

通用9内部排序的类型(C语言)


版权声明:本文博客原创文章,博客,未经同意,不得转载。


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
C#Random()函数详解发布时间:2022-07-13
下一篇:
c++jsoncpp解析json发布时间:2022-07-13
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap