NOIP 2013 T2 火柴排队 ---->求逆序对

[NOIP2013T2]火柴排队

背景

noip2013day1

描述

涵涵有两盒火柴，每盒装有 n 根火柴，每根火柴都有一个高度。 现在将每盒中的火柴各
自 排成一列， 同一列火柴的高度互不相同， 两列火柴之间的距离定义为：∑(ai-bi)2，

其中 ai 表示第一列火柴中第 i 个火柴的高度， bi 表示第二列火柴中第 i 个火柴的高度。

每列火柴中相邻两根火柴的位置都可以交换，请你通过交换使得两列火柴之间的距离最
小。请问得到这个最小的距离，最少需要交换多少次？ 如果这个数字太大，请输出这个最
小交换次数对 99,999,997 取模的结果。

输入格式

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

输出格式

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

测试样例1

输入

4
2 3 1 4
3 2 1 4
输出

1
备注

对于 10%的数据， 1 ≤ n ≤ 10；
对于 30%的数据， 1 ≤ n ≤ 100；
对于 60%的数据， 1 ≤ n ≤ 1,000；
对于 100%的数据， 1 ≤ n ≤ 100,000， 0 ≤火柴高度≤ 231 − 1。

这里给出的是树状数组实现的求逆序对。。。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define mod 99999997 
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,ans;
int t[100005],rank[100005];
struct data{int v,num;}a[100005],b[100005];
inline bool cmp1(data a,data b)
{return a.v<b.v;}
inline int lowbit(int x){return x&(-x);}
inline int ask(int x)
{
    int tmp=0;
    for(int i=x;i;i-=lowbit(i))
        tmp+=t[i];
    return tmp;
}
void update(int x)
{
    for(int i=x;i<=n;i+=lowbit(i))
        t[i]++;
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++)
        {a[i].v=read();a[i].num=i;}
    for(int i=1;i<=n;i++)
        {b[i].v=read();b[i].num=i;}
    sort(a+1,a+n+1,cmp1);
    sort(b+1,b+n+1,cmp1);
    for(int i=1;i<=n;i++)
        rank[a[i].num]=b[i].num;
    for(int i=1;i<=n;i++)
    {
        update(rank[i]);
        ans=(ans+i-ask(rank[i]))%mod;
    }
    printf("%d",ans);
    return 0;
}