记录编号 81887 评测结果 AAAAAAAAAA
题目名称 [NOIP 2013]火柴排队 最终得分 100
用户昵称 Gravatarcstdio 是否通过 通过
代码语言 C++ 运行时间 0.115 s
提交时间 2013-11-19 14:02:29 内存使用 2.70 MiB
显示代码纯文本
#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<iomanip>
#include<vector>
using namespace std;
typedef long long ll;
const int SIZEN=100001;
const int MOD=99999997;
class MATCH{
public:
	int height;//高度
	int pos;//在原序列中的位置
};
bool operator < (MATCH x,MATCH y){
	return x.height<y.height;
}
int n;
MATCH A[SIZEN],B[SIZEN];
int agoal[SIZEN]={0};//其实就是A[i],B[i]是第几个
int lis[SIZEN]={0};
int ans=0;
bool flag[SIZEN]={0};
vector<int> atp;
ll inversion(int a[],int x,int y){
	if(x==y) return 0;
	int mid=(x+y)>>1;
	ll sum=inversion(a,x,mid)+inversion(a,mid+1,y);
	atp.clear();
	int i;
	int p1=x,p2=mid+1;
	for(i=x;i<=y;i++){//统计"比它靠后且比它小的"
		if(p1>mid) atp.push_back(a[p2++]);
		else if(p2>y){
			atp.push_back(a[p1++]);
			sum+=y-mid;//在这个数之前入队的后一半个数
		}
		else{
			if(a[p1]<=a[p2]){
				atp.push_back(a[p1++]);
				sum+=p2-mid-1;
			}
			else atp.push_back(a[p2++]);
		}
	}
	for(i=x;i<=y;i++) a[i]=atp[i-x];
	return sum%MOD;
}
void getlis(void){
	sort(A+1,A+1+n);
	sort(B+1,B+1+n);
	int i,temp;
	for(i=1;i<=n;i++) agoal[A[i].pos]=i;
	//将B的第agoal[i]大的数改成i
	for(i=1;i<=n;i++){
		temp=agoal[i];
		lis[B[temp].pos]=i;
	}
}
void read(void){
	scanf("%d",&n);
	int i;
	for(i=1;i<=n;i++){
		scanf("%d",&A[i].height);
		A[i].pos=i;
	}
	for(i=1;i<=n;i++){
		scanf("%d",&B[i].height);
		B[i].pos=i;
	}
}
int main(){
	freopen("MatchNOIP2013.in","r",stdin);
	freopen("MatchNOIP2013.out","w",stdout);
	read();
	getlis();
	printf("%lld\n",inversion(lis,1,n));
	return 0;
}