3286 火柴排队 2013年NOIP全国联赛提高组

http://codevs.cn/problem/3286/
题目描述 Description
涵涵有两盒火柴，每盒装有 n 根火柴，每根火柴都有一个高度。现在将每盒中的火柴各自排成一列，同一列火柴的高度互不相同，两列火柴之间的距离定义为：
，其中 ai表示第一列火柴中第 i 个火柴的高度，bi表示第二列火柴中第 i 个火柴的高度。
每列火柴中相邻两根火柴的位置都可以交换，请你通过交换使得两列火柴之间的距离最小。请问得到这个最小的距离，最少需要交换多少次？如果这个数字太大，请输出这个最小交换次数对 99,999,997 取模的结果。

输入描述 Input Description
共三行，第一行包含一个整数 n，表示每盒中火柴的数目。
第二行有 n 个整数，每两个整数之间用一个空格隔开，表示第一列火柴的高度。
第三行有 n 个整数，每两个整数之间用一个空格隔开，表示第二列火柴的高度。

输出描述 Output Description
输出共一行，包含一个整数，表示最少交换次数对 99,999,997 取模的结果。

样例输入 Sample Input
[Sample 1]
4
2 3 1 4
3 2 1 4
[Sample 2]
4
1 3 4 2
1 7 2 4

样例输出 Sample Output
[Sample 1]
1
[Sample 2]
2

数据范围及提示 Data Size & Hint
【样例1说明】
最小距离是 0，最少需要交换 1 次，比如：交换第 1 列的前 2 根火柴或者交换第 2 列的前 2 根火柴。
【样例2说明】
最小距离是 10，最少需要交换 2 次，比如：交换第 1 列的中间 2 根火柴的位置，再交换第 2 列中后 2 根火柴的位置。
【数据范围】
对于 10%的数据， 1 ≤ n ≤ 10；
对于 30%的数据，1 ≤ n ≤ 100；
对于 60%的数据，1 ≤ n ≤ 1,000；
对于 100%的数据，1 ≤ n ≤ 100,000，0 ≤火柴高度≤ 2^31 - 1。

这就要求逆序对的对数了。为什么呢？
请注意只能够移动相邻的两根火柴。
下面是证明：（来自百度NOIP贴吧）
1、一个序列有序，当且仅当其不存在逆序对
2、如果一个序列不是有序的，就必然存在一对相邻的两个数，这两个数构成一对逆序对（否则整个序列就是有序的了）
3、交换2中的两个数，可以使得整个序列的逆序对数恰好减少1
综上所述，命题得证

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
using namespace std;
struct sir
{
    int a,num;
};
sir a[100005],b[100005];
int Rank[100005],C[100005],c[100005];
long long cnt=0;
const int mod=99999997;
int cmp1(sir x,sir y)
{
    return x.a<y.a;
}
int cmp2(sir x,sir y)
{
    return x.num<y.num;
}
void Mergesort(int l,int r)//归并排序求逆序对
{
    int mid,i,j,tmp;
    if(r>l+1)
    {
        mid=(l+r)/2;//分治
        Mergesort(l,mid);
        Mergesort(mid,r);
        tmp=l;
        for(i=l,j=mid;i<mid&&j<r;)
        {
            if(C[i]>C[j])//逆序对，我们按升序排列
            {
                c[tmp++]=C[j++];
                cnt=(cnt+mid-i)%mod;
            }
            else{
                c[tmp++]=C[i++];
            }
        }
        //给一个区间排序不会影响这个区间中的数和另一个区间的数的逆序对个数
        if(j<r)
            for(;j<r;j++)c[tmp++]=C[j];
        else
            for(;i<mid;i++)c[tmp++]=C[i];
        for(i=l;i<r;i++)C[i]=c[i];
    }
}
int main()
{
    int n;
    cin>>n;
    for(int i=0;i<n;i++) scanf("%d",&a[i].a),a[i].num=i;
    for(int i=0;i<n;i++) scanf("%d",&b[i].a),b[i].num=i;
    sort(a,a+n,cmp1);
    sort(b,b+n,cmp1);
    for(int i=0;i<n;i++)a[i].a=b[i].a=i;
    sort(a,a+n,cmp2);
    sort(b,b+n,cmp2);
    for(int i=0;i<n;i++)Rank[a[i].a]=i;
    for(int i=0;i<n;i++)C[i]=Rank[b[i].a];
    Mergesort(0,n);
    cout<<cnt<<endl;
    return 0;
}