[jobdu]数组中的逆序对

http://ac.jobdu.com/problem.php?pid=1348 数组中的逆序对也是个常见的题目,算法导论中也有一些描述,参考:http://www.cnblogs.com/wuyuegb2312/p/3156286.html。解法自然有朴素解法O(n^2)不表,然后有基于归并排序的,可以O(n*logn)。

1.在归并排序中,同样是对一个数组分为两段处理,在处理这两段时,并不会影响右段元素与左段元素的逆序关系,只有在归并时才会改变。
2.归并时的改变方式和插入排序是类似的:右段中取出元素放在左段其余所有元素前面时,相当于左段整体后移,后移的元素数就是这个逆序数。
3.由于归并排序使用的是分治法,将每次归并的逆序数累加,最后结果就是总的逆序数。并且,归并排序的时间复杂度是O(nlogn),优于插入排序。
根据以上的探讨,归并排序稍作修改,就获得了时间复杂度为O(nlogn)的寻找逆序对总数的算法了。

本来也到此结束了,后来又发现一个基于树状数组的解法(http://boj.haotui.com/thread-2792-1-1.html),不过还是太高级,就不深入研究了。树状数组略看了一下,和数字的和有关,比如可以解决多数求和的问题。“采用累加的方法还有一个局限,那就是,当修改掉数组中的元素后,仍然要你求数组中某段元素的和,就显得麻烦了。所以我们就要用到树状数组,他的时间复杂度为O(lgn)。”(http://www.cnblogs.com/zhangshu/archive/2011/08/16/2141396.html

具体代码:

1.C++中的数组新建和删除是 int* arr = new int[n]; delete arr[];
2.出现HEAP CORRUPTION DETECTED往往是操作new申请的内存溢出。比如我一开始调用了mergeCount(0, n)就错误的访问了a[n]。 http://blog.sina.com.cn/s/blog_511703010100lz33.html

3.借鉴了别人的写法,在一个循环里就把双指针移动写完了,主要是把到头后的判断加上,一开始使用||。
4.因为mergeSort的空间复杂度是o(n),所以也可以一开始就申请好数组,就不用在递归里不断new和delete了。
5.mergeCount的时候,当右边的小于左边的时,要注意计算对逆序数目。

#include <cstdio>

using namespace std;

long long tot;

int a[100005];

void mergeCount(int l, int r) // [l, m], [m+1, r]

{

    if (l >= r) return;

    int m = (l + r) / 2;

    mergeCount(l, m);

    mergeCount(m+1, r);

    int* tmp = new int[r-l+1];

    int i = l;

    int j = m+1;

    int k  = 0;

    while (i <= m || j <= r)

    {

        if (i > m)

        {

            tmp[k++] = a[j++];

        }

        else if (j > r)

        {

            tmp[k++] = a[i++];

        }

        else if (a[i] <= a[j])

        {

            tmp[k++] = a[i++];

        }

        else

        {

            tmp[k++] = a[j++];

            tot += (m-i+1);

        }

    }

    for (int j = 0; j <= (r - l); j++)

    {

        a[l+j] = tmp[j];

    }

    delete[] tmp;

}

   

int main()

{

    int n;  

    while(scanf("%d", &n) == 1)  

    {  

        for(int i=0; i<n; ++i)  

            scanf("%d", &a[i]);              

        tot = 0;  

        mergeCount(0, n-1);  

        printf("%lld\n", tot);  

    }  

}

  

你可能感兴趣的:(job)