HDU 1394 Minimum Inversion Number



http://acm.hdu.edu.cn/showproblem.php?pid=1394

用一种新的方法求逆序数——分治法（之前使用的是线段树）

原理说明

划分问题：把元素序列分为个数尽量相等的两半；

递归求解：统计i和j均在左面或均在右面的逆序对个数；

合并问题：统计i在左边，j在右边的逆序对个数。我们发现对于右边的每个j，统计左边比它大的个数f(j)，则所有f(j)之和就是答案。在归并排序中合并是由小到大进行的，当右边的a[j]复制到临时空间时，左边还没来得及复制到临时空间的那些数就是左边所有比a[j]大的数，此时只要计数左边元素的个数即可，由于左边的区间是[p,mid)，因此共有m-p个元素。时间复杂度为O(nlog(n))的。

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>

using namespace std;

const int N = 5000;
int a[N << 1],t[N << 1],b[N << 1], n, cnt;

void merge_sort(int x, int y){
    if(y - x > 1){
        int mid = x + (y - x)/2;
        int p = x, q = mid, i = x;
        merge_sort(x, mid);
        merge_sort(mid, y);
        while(p < mid || q < y){///左右有一个非空就继续合并
            if(q >= y || (p < mid && a[p] <= a[q]))///从左半数组复制到临时空间（右半数组为空，那么左半数组必不为空）
                t[i++] = a[p++];
            else
                {t[i++] = a[q++]; cnt += mid - p;}
        }
        for(int i=x; i<y; i++)  a[i] = t[i];
    }
}

int main(){
//    freopen("in.txt", "r", stdin);
    while(scanf("%d",&n) == 1){
        for(int i=1; i<=n; i++) {scanf("%d",&a[i]);b[i] = a[i];}
        cnt = 0;
        merge_sort(1,n+1);
        int ans = cnt;
        for(int i=1; i<=n; i++){
            cnt = cnt + n - 1 - 2*b[i];
            ans = min(ans, cnt);
        }
        printf("%d\n",ans);
    }
    return 0;
}