逆序对

求解逆序对的两种方向

1. 归并排序法

   归并排序是应用分治的思想实现排序的过程。

   归并排序是将数列a[l,h]分成两半a[l,mid]和a[mid+1,h]分别进行归并排序，然后再将这两半合并起来。在合并的过程中（设l<=i<=mid，mid+1<=j<=h），当a[i]<=a[j]时，并不产生逆序数；当a[i]>a[j]时，在前半部分中比a[i]大的数都比a[j]大，将a[j]放在a[i]前面的话，逆序数要加上mid+1-i。因此，可以在归并排序中的合并过程中计算逆序数.
   因为k = i + j - mid - 1 => mid+1-i = j - k

/**
    题意:给你长度为n个不同的数，问最少需要交换多少次使得升序 只能交换相邻的两个 其实就是冒泡排序的过程
    题解:这题可以采用归并排序来实现求逆序对。其实归并排序的副产品就是逆序对
*/
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
using namespace std;
const int maxn = 5e5+10;
typedef long long ll;
ll num[maxn],temp[maxn];
ll ans;
void mergenum(int low,int mid,int high)
{
    int i = low,j = mid + 1,k = low;
    while(i<=mid && j<=high) {
        if(num[i] <= num[j]) temp[k++] = num[i++];
        else ans += j - k,temp[k++] = num[j++];
    }
    while(i<=mid) temp[k++] = num[i++];
    while(j<=high) temp[k++] = num[j++];
    for(int i=low;i<=high;i++) num[i] = temp[i];
}
void mergesort(int low,int high)
{
    if(low < high){
        int mid = (low + high)>>1;
        mergesort(low,mid);
        mergesort(mid+1,high);
        mergenum(low,mid,high);
    }
}
int main()
{
    int n;
    while(~scanf("%d",&n),n)
    {
        for(int i=0;i<n;i++) scanf("%lld",&num[i]);
        ans = 0;
        mergesort(0,n-1);
        printf("%lld\n",ans);
    }
    return 0;
}

2.树状数组法

众所周知，线段树和树状树可以用来解决区间操作问题，就是因为这两个算法区间操作的时间复杂度很低O(logN)，才让这种方法具有可行性。

首先先来看一个序列 6 1 2 7 3 4 8 5，此序列的逆序数为5+3+1=9。冒泡法可以直接枚举出逆序数，但是时间复杂度太高O(n^2)。冒泡排序的原理是枚举每一个数组，然后找出这个数后面有多少个数是小于这个数的，小于它逆序数+1。仔细想一下，如果我们不用枚举这个数后面的所有数，而是直接得到小于这个数的个数，那么效率将会大大提高。

总共有N个数，如何判断第i+1个数到最后一个数之间有多少个数小于第i个数呢？不妨假设有一个区间 [1,N]，只需要判断区间[i+1,N]之间有多少个数小于第i个数。如果我们把总区间初始化为0，然后把第i个数之前出现过的数都在相应的区间把它的值定为1，那么问题就转换成了[i+1,N]值的总和。再仔细想一下，区间[1,i]的值+区间[i+1,N]的值=区间[1,N]的值(i已经标记为1)，所以区间[i+1,N]值的总和等于N-[1,i]的值！因为总共有N个数，不是比它小就是比它(大或等于)。

现在问题已经转化成了区间问题，枚举每个数，然后查询这个数前面的区间值的总和，i-[1,i]既为逆序数。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
typedef struct node
{
    int x,tag;
}point;
const int maxn = 5e5+10;
point num[maxn];
int c[maxn];
int position[maxn];
int lowbit(int x) { return x & -x;}
void add(int pos,int x) {
    for(int i=pos;i<maxn;i += lowbit(i)) c[i] += x;
}
ll sum(int pos){
    ll ans = 0;
    for(int i=pos;i;i -= lowbit(i)) ans += c[i];
    return ans;
}
bool cmp(const point &a,const point &b){
    return a.x < b.x;
}
void init()
{
    memset(c,0,sizeof(c));
}
int main()
{
    int n;
    while(~scanf("%d",&n),n)
    {
        init();
        for(int i=1;i<=n;i++) scanf("%d",&num[i].x),num[i].tag = i;
        sort(num+1,num+1+n,cmp);
        for(int i=1;i<=n;i++) position[num[i].tag] = i;          ///表示这个数放在第几个位置
        ll ans = 0;
        for(int i=1;i<=n;i++) {///按照顺序放入树状数组中
            add(position[i],1);
            ans += i - sum(position[i]);///sum(position[i])表示这个[1,position]有多少个,即<=position[i]的有多少个,于是>position[i]有i - sum(position[i])个(感觉类似于插入排序)
        }
        printf("%lld\n",ans);
    }
}