归并排序，树状数组,求逆序数 (openjudge 7662)

题目链接

描述
在Internet上的搜索引擎经常需要对信息进行比较，比如可以通过某个人对一些事物的排名来估计他（或她）对各种不同信息的兴趣，从而实现个性化的服务。

对于不同的排名结果可以用逆序来评价它们之间的差异。考虑1,2,…,n的排列i1，i2，…，in，如果其中存在j,k，满足 j < k 且 ij > ik， 那么就称(ij,ik)是这个排列的一个逆序。

一个排列含有逆序的个数称为这个排列的逆序数。例如排列 263451 含有8个逆序(2,1),(6,3),(6,4),(6,5),(6,1),(3,1),(4,1),(5,1)，因此该排列的逆序数就是8。显然，由1,2,…,n 构成的所有n!个排列中，最小的逆序数是0，对应的排列就是1,2,…,n；最大的逆序数是n(n-1)/2，对应的排列就是n,(n-1),…,2,1。逆序数越大的排列与原始排列的差异度就越大。

现给定1,2,…,n的一个排列，求它的逆序数。

输入
第一行是一个整数n，表示该排列有n个数（n <= 100000)。
第二行是n个不同的正整数，之间以空格隔开，表示该排列。
输出
输出该排列的逆序数。
样例输入
6
2 6 3 4 5 1
样例输出
8
提示

1. 利用二分归并排序算法（分治）；
2. 注意结果可能超过int的范围，需要用long long存储。

这里我们用归并排序，来求逆序数，我们归并排序就是对每一半进行归并，当对于两个有序的序列如7 5 3 1 和8 6 4 2，我们依然按照归并排序从大到小排序的方法，然后记录当p2小于p1时，则p2到r2的所有数都小于p1位置的数，我们加上个(r2- p2+1)就行，最后的累加和就是逆序数

代码如下:

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;
typedef long long ll;
const int MAX = 100010;
ll Count;
void Merge(ll a[],int l,int m,int r,ll b[]){
    int pos,p1,p2;
    pos = 0,p1 = l,p2 = m+1;
    while(p1 <= m && p2 <= r){
        if(a[p2] >= a[p1]){
            b[pos++] = a[p2++];
        }
        else{
            Count += (r - p2 + 1);
            b[pos++] = a[p1++];
        }
    }
    while(p1 <= m) b[pos++] = a[p1++];
    while(p2 <= r) b[pos++] = a[p2++];
    pos = 0;
    for(int i=l;i<=r;i++){
        a[i] = b[pos++];
    }
}
void MergeSort(ll a[],int l,int r,ll b[]){
    if(l < r){
        int m = (l+r)/2;
        MergeSort(a,l,m,b);
        MergeSort(a,m+1,r,b);
        Merge(a,l,m,r,b);
    }
}

ll a[MAX],b[MAX];
int main(void){
    int n;
    Count = 0;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    MergeSort(a,1,n,b);

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

    return 0;
}

求逆序数还有一种利用树状数组的求法，不过要先离散化一下，代码也贴在了下面

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<queue>

using namespace std;
typedef long long ll;
#define lowbit(x) (x&-x)
const int MAX = 100010;
int n;
struct node{
    ll value;
    ll pos;
}q[MAX];
bool cmp(node a,node b){
    return a.value < b.value;
}
ll dat[MAX];
ll sum(int i){
    ll res = 0;
    while(i > 0){
        res += dat[i];
        i -= lowbit(i);
    }
    return res;
}
void add(int i){
    while(i <= n){
        dat[i] += 1;
        i += lowbit(i);
    }
}
int aa[MAX];
int main(void){

    memset(dat,0,sizeof(dat));
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%lld",&q[i].value);
        q[i].pos = i;
    }
    stable_sort(q+1,q+1+n,cmp);//这里不能用sort，因为sort是一个不稳定的排序，在元素相同时也有可能发生交换
    for(int i=1;i<=n;i++){
        aa[q[i].pos] = i;
    }
    ll ans = 0;
    for(int i=1;i<=n;i++){
        ans += (i - sum(aa[i]) - 1);
        add(aa[i]);
    }
    printf("%lld\n",ans);
    return 0;
}