当前位置:主页   - 电脑 - 程序设计 - C/C++
(算法导论习题解problem2.4)寻找一个序列中逆序对的数量
来源:网络   作者:   更新时间:2010-03-03
收藏此页】    【字号    】    【打印】    【关闭

  一个序列的逆序对是这样的两个元素, 对于序列A而言, i>j且A[i]<A[j], 于是A[i]和A[j]就形成一个逆序对.

  研究一个序列中逆序对的数量是有实际意义的, 对于插入排序而言, 它排序的时间与待排序序列的逆序对数量成正比.

  下面给出求出一个序列中逆序对数量的算法,类似于归并排序中使用的分治算法:一个序列的逆序对数量为它的左半部分逆序对的数量,加上右半部分逆序对的数量, 最后在合并的时候左半部分元素大于右半部分元素的数量.这几乎和归并算法的过程一模一样,只是在归并的时候加入了计数的操作, 完整的算法为:

#include <stdio.h>

void display(int array[], int size)
{
    int i;
    for (i = 0; i < size; ++i)
    {
        printf("%d ", array[i]);
    }

    printf("n");
}

int merge_inversion(int array[], int low, int middle, int high)
{
    int count = 0;
    int llen = middle - low + 1;
    int hlen = high - middle;
    int *L = (int *)malloc((llen + 1) * sizeof(int));
    int *H = (int *)malloc((hlen + 1) * sizeof(int));
    int i, j, k;

    for(i = 0; i < llen; i++)
        L[i] = array[low + i];
    L[i] = 99999;
    for(i = 0; i < hlen; i++)
        H[i] = array[middle + i + 1];
    H[i] = 99999;

    for(i = 0, j = 0, k = low; k < high + 1; k++)  
    {
        if(L[i] > H[j])
        {
            array[k] = H[j++];
            count += llen - i;
        }
        else
        {
            array[k] = L[i++];
        }
    }

    free(L);
    free(H);

    return count;
}

int count_inversion(int array[], int low, int high)
{
    int count = 0, middle;
    if(low < high)
    {
        middle = low + (high - low) / 2;
        count += count_inversion(array, low, middle);
        count += count_inversion(array, middle + 1, high);
        count += merge_inversion(array, low, middle, high);
    }

    return count;
}

int main()
{
    int array[]={2,8,3,6,1};
    printf("count of inversions is %dn",count_inversion(array, 0, 4));

    display(array, 5);

    return 0;
}

其它资源
来源声明

版权与免责声明
1、本站所发布的文章仅供技术交流参考,本站不主张将其做为决策的依据,浏览者可自愿选择采信与否,本站不对因采信这些信息所产生的任何问题负责。
2、本站部分文章来源于网络,其版权为原权利人所有。由于来源之故,有的文章未能获得作者姓名,署“未知”或“佚名”。对于这些文章,有知悉作者姓名的请告知本站,以便及时署名。如果作者要求删除,我们将予以删除。除此之外本站不再承担其它责任。
3、本站部分文章来源于本站原创,本站拥有所有权利。
4、如对本站发布的信息有异议,请联系我们,经本站确认后,将在三个工作日内做出修改或删除处理。
请参阅权责声明