2013年NOIP提高组 火柴排队

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

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

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

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

根据距离的定义，最小的距离应该是数列a的第 i 大数 与数列b的第 i 大数分别做差，离散化后就变成了 ai 与 bi 相等，即求数列 a 基于 数列 b 的逆序对个数。

#include<iostream>#include<cstdio>#include<algorithm>#include<cstring>using namespace std;const int size = 100010;int num[size];int templ[size];int temp[size];int ans = 0;int n;int a[size],b[size];int cmp[size];const int mod = 99999997;int EF(int l,int r,int k,int *tmp){    while(r - l > 1)    {        int mid = (l + r) >> 1;        if(tmp[k] > templ[mid])            l = mid;        else if (tmp[k] < templ[mid])            r = mid;        else            return mid;    }    return r;}void lsha(){    for(int i = 1 ; i <= n ; i ++)        templ[i] = a[i];    sort( templ+1 , templ+n+1 );    for(int i = 1 ; i <= n ; i ++)        a[i] = EF(0,n+1,i,a);}void lshb(){    for(int i = 1 ; i <= n ; i ++)        templ[i] = b[i];    sort( templ+1 , templ+n+1 );    for(int i = 1 ; i <= n ; i ++)        b[i] = EF(0,n+1,i,b);}void merge(int l,int r){    if(l == r)        return ;    int mid = (l + r) >> 1;    merge(l,mid) , merge(mid+1,r);    int ll = l , rr = mid + 1 , p = ll;    while(ll <= mid || rr <= r)    {        if(rr > r || (ll <= mid && cmp[a[ll]] <= cmp[a[rr]]))            temp[p++] = a[ll++];        else        {            temp[p++] = a[rr++];            ans += mid - ll + 1;            ans %= mod;        }    }    for(int i = l ; i <= r ; i ++)        a[i] = temp[i];}int main(){    int n;    scanf("%d",&n);    for(int i = 1 ; i <= n ; i ++)        scanf("%d",&a[i]);    for(int i = 1 ; i <= n ; i ++)        scanf("%d",&b[i]);    lsha();    memset(templ,0,sizeof(templ));    lshb();    for(int i = 1 ; i <= n ; i ++)        cmp[b[i]] = i;    merge(1,n);    printf("%d",ans%mod);    return 0;}

传送门 :
codevs 3286
tyvj 2515