记录编号 |
81887 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[NOIP 2013]火柴排队 |
最终得分 |
100 |
用户昵称 |
cstdio |
是否通过 |
通过 |
代码语言 |
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;
}