#NOIP2013#火柴排序(逆序对)

2284: 【NOIP2013】火柴排队

时间限制: 1 Sec  内存限制: 128 MB

题目描述

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

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

输入

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

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

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

输出

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

样例输入

样例一

Copy (如果复制到控制台无换行，可以先粘贴到文本编辑器，再复制)

4

2 3 1 4

3 2 1 4

样例二

Copy (如果复制到控制台无换行，可以先粘贴到文本编辑器，再复制)

4

1 3 4 2

1 7 2 4

样例输出

样例一

1

样例二

2

提示

【输入输出说明】

一：最小距离是 0，最少需要交换 1 次，比如：交换第 1 列的前 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³¹ − 1。

【解】

      两串火柴的高度为两个无序序列，因距离为两火柴高度之差的平方，所以只有当 1 串的第i高火柴和 2 串的第i高火柴放在同一编号位时，才符合题意要求。

      对两串火柴的高度进行从小到大排序。（排序前记录下原始编号）

      对两串火柴从小到大边号进行组合（可以将其中一串的编号作为数组下标，将另一串的编号对应存入数组）。

      在组合出的数据（数组）中查找逆序对，逆序对的数量即为题目的解。

      【证明】

       1.如果一个序列有序，必然不存在逆序对。

       2.若一个队列不是有序的，就必然存在一对相邻的两个数，它们会构成一对逆序对。

       3.交换 2 中的两个数，整个序列的逆序对数量会减少一。

     

       查找逆序对时的时间限制：

      1.普通查找：

for(i=0;i<n;i++){
	for(j=i+1;j<n;j++)
		if(c[j]<c[i])
			ans++;

       //时间复杂度为O(n^2) 可以通过70%的数据

      2.利用归并排序的思想：

while(i<=mid&&j<=r)
	if(c[i]<=c[j])d[k++]=c[i++];
	else {
		d[k++]=c[j++];
		ans+=mid+1-i;
		ans%=99999997;
    }
while(i<=mid)
	d[k++]=c[i++];
while(j<=r)
	d[k++]=c[j++];
for(i=l;i<=r;i++)
	c[i]=d[i];

      //时间复杂度为O(nlogn) AC

#include<cstdio>
#include<cstring>
#include<limits.h>
#include<cmath>
#include<algorithm>
using namespace std;
long long ans;
int c[1000005],d[1000005];
struct node{
int xx,yy;
}a[1000005],b[1000005];
int cmp(const void *a,const void *b){
return (*(node*)a).xx<(*(node*)b).xx;
}
void inversion(int *c,int l,int mid,int r){//逆序对
int i=l,j=mid+1,k=l;
while(i<=mid&&j<=r)
if(c[i]<=c[j])d[k++]=c[i++];
else {
d[k++]=c[j++];
ans+=mid+1-i;
ans%=99999997;
   }
while(i<=mid)
d[k++]=c[i++];
while(j<=r)
d[k++]=c[j++];
for(i=l;i<=r;i++)
c[i]=d[i];
}
void merge(int *c,int l,int r){//查找
if(l<r){
int mid=(l+r)/2;
merge(c,l,mid);
merge(c,mid+1,r);
inversion(c,l,mid,r);
}
}
int main(){
int i,n;
scanf("%d",&n);
for(i=0;i<n;i++){
scanf("%d",&a[i].xx);
a[i].yy=i;
}
for(i=0;i<n;i++){
scanf("%d",&b[i].xx);
b[i].yy=i;
}
qsort(a,n,sizeof(node),cmp);
qsort(b,n,sizeof(node),cmp);
for(i=0;i<n;i++)
c[b[i].yy]=a[i].yy;
merge(c,0,n-1);
printf("%Ild\n",ans);
}