NOIP2013火柴排队

题目描述

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

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

输入

输入文件为 match.in。

共三行，第一行包含一个整数 n，表示每盒中火柴的数目。

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

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

输出

输出文件为 match.out。

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

样例输入

【输入输出样例 1】
4
2 3 1 4
3 2 1 4
【输入输出样例 2】4
1 3 4 2
1 7 2 4

样例输出

【输入输出样例 1】
1
【输入输出样例 2】
2

提示

【输入输出样例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 ≤火柴高度≤ 231 − 1。

首先看到题目中的，思考一下什么时候sum取到最小值，随便取几个数据代入就可以得出：sum取最小值当且仅当数列{an},{bn}已经按相同顺序排序时。

其实这一个也是有数学依据的。

根据排序不等式：

（以下内容摘自百度百科）排序不等式表述如下，设有两组数a1,a2,……an和b1,b2,……bn，当满足a1≤a2≤……≤an,b1≤b2≤……≤bn则有a1bn+a2bn-1+……+anb1≤a1bt1+a2bt2+……+anbtn≤a1b1+a2b2+anbn式中t1，t2，……，tn是1，2，……，n的任意一个排列，当且仅当a1=a2=……=an或b1=b2=……=bn时成立。一般为了便于记忆，常记为：反序和≤乱序和≤同序和.

数学描述很复杂，总结起来就一句话，两数列对应项乘积之和，当且仅当相同顺序排序时最小，当且仅当相反顺序排序时最大。

 反正可以想象出来两个数列对应项在数列中位置相同的时候，sum最小。

那么原问题就转化为：把数列b的数字大小关系移动成 数列a的大小关系所需要的最小移动次数。

将数列a，b离散化，找出b中每一个数字对应在a中的位置，然后问题就转化为将b数列冒泡排序所需要的最小交换次数。

也就是求逆序对的个数。

为什么这样？

根据冒泡排序的原理，当一个数列为顺序数列的时候，逆序对个数为0。每当后一个小于前一个数字时，逆序对个数增加一个。也就是每次交换前后两个数字，逆序对个数减少一个，那么最少的交换次数应该就是逆序对的个数（将逆序对个数减为0的最少交换次数）。

下面是代码：

#include <iostream>
#include <algorithm>

using namespace std;
const int N=100009;
int n,a[N],b[N],cnt;
void lisan(int *a);
void Ini();
void merg(int l,int mid,int r);
void st(int l,int r);
int main()
{
    Ini();
    st(1,n);
    cout<<cnt<<endl;
    return 0;
}
void Ini()
{
    int m,t[N];
    cin>>n;
    for(int i=1;i<=n;i++)cin>>a[i];
    lisan(a);
    for(int i=1;i<=n;i++)cin>>b[i];
    lisan(b);
    for(int i=1;i<=n;i++)t[a[i]]=i;
    for(int i=1;i<=n;i++)b[i]=t[b[i]];
}
void lisan(int *g)
{
    int t[N],m;
    for(int i=1;i<=n;i++)t[i]=g[i];
    sort(t+1,t+1+n);
    m=unique(t+1,t+1+n)-t-1;
    for(int i=1;i<=n;i++)
        g[i]=lower_bound(t+1,t+1+m,g[i])-t;
}
void st(int l,int r){
    if(l==r)return;
    int mid=(l+r)>>1;
    st(l,mid);
    st(mid+1,r);
    merg(l,mid,r);
}
void merg(int l,int mid,int r){
    int i=l,j=mid+1;
    for(int k=l;k<=r;k++){
        if(j>r||i<=mid&&b[i]<=b[j])a[k]=b[i++];
        else a[k]=b[j++],cnt=(cnt+mid-i+1)%99999997;
    }
    for(int k=l;k<=r;k++){
        b[k]=a[k];
    }
}