【NOIP 2013 DAY.1】火柴排队【codevs 3286】

分析：贪心策略。第一行第一小对第二行第一小、第一行第二小对第二行第二小。。。类推。

即：排序，求排序的次数。

（归并排序求逆序对）【记录交换的次数即是答案】*推荐使用归并。本题最优解法。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>

using namespace std;

const int mod=99999997;
int a[100005],b[100005];
int l[100005],r[100005];
int s[100005],t[100005];
int rr[100005];
int sum;
int n;

void m_sort(int x,int y)
{
if(y-x>1)
{
int m=x+(y-x>>1);
int p=x,q=m,i=x;
m_sort(x,m);
m_sort(m,y);
while(p<m || q<y)
{
if(q>=y || (p<m && s[p]<=s[q])) t[i++]=s[p++];
else
{
t[i++]=s[q++];
sum+=m-p;
}
}
for(i=x;i<y;i++) s[i]=t[i];
}
sum%=mod;
}

int main()
{
scanf("%d",&n);
for(int i=0;i<n;i++)
{
scanf("%d",a+i);
l[i]=a[i];
}
for(int i=0;i<n;i++)
{
scanf("%d",b+i);
r[i]=b[i];
}
sort(l,l+n);
sort(r,r+n);
for(int i=0;i<n;i++)
{
int locB=lower_bound(r,r+n,b[i])-r;
rr[locB]=i;
}
for(int i=0;i<n;i++)
{
int locA=lower_bound(l,l+n,a[i])-l;
s[rr[locA]]=i;
}
m_sort(0,n);
printf("%d",sum);
return 0;
}

（线段树/树状数组求逆序对）

#include<iostream>

#include<cstdlib>

#include<cstdio>

#include<algorithm>

using namespace std;

struct stone

{

int num,c;

} l1[100010],l2[100010];

int n,c[100010];

int tree[100010];

int lowbit(int x)

{

return (-x)&x;

}

bool cmp(stone x,stone y)

{

return x.c<y.c;

}

void insert(int x)

{

while(x<=n)

{

tree[x]++;

x+=lowbit(x);

}

}

int get(int x)

{

int sum=0;

while(x>0)

{

sum+=tree[x];

x-=lowbit(x);

}

return sum;

}

int main()

{

cin>>n;

for(int i=1;i<=n;i++)

{

scanf("%d",&l1[i].c);

l1[i].num=i;

}

for(int i=1;i<=n;i++)

{

scanf("%d",&l2[i].c);

l2[i].num=i;

}

sort(l1+1,l1+n+1,cmp);

sort(l2+1,l2+n+1,cmp);

long long cnt=0;

for(int i=1;i<=n;i++)

{

c[l1[i].num]=l2[i].num;

}

//for(int i=1;i<=n;i++) cout<<c[i]<<" ";

for(int i=1;i<=n;i++)

{

insert(c[i]);

cnt+=i-get(c[i]);

cnt=cnt%99999997;

}

cout<<cnt;

return 0;

}

std：：

#include<cstdio>
#include<algorithm>
#define MAXN 100000
using namespace std;
struct node
{
int x,cnt;
}A[MAXN+10],B[MAXN+10];
int aa[MAXN+10];
int T[MAXN+10];
int ans=0;
bool cmp(node a,node b)
{
return a.x<b.x;
}
void marge(int l,int mid,int r)
{
int i=l,j=mid+1,k;
for(k=l;k<=r;k++)
if((i<=mid)&&((j>r)||(aa[i]<=aa[j])))
T[k]=aa[i++];
else{
T[k]=aa[j++];
ans+=mid-i+1;
ans%=99999997;
}
for(i=l;i<=r;i++)
aa[i]=T[i];
}
void margef(int l,int r)
{
if(l<r)
{
int mid=(l+r)/2;
margef(l,mid);
margef(mid+1,r);
marge(l,mid,r);
}
}
int main()
{
freopen("match.in","r",stdin);
freopen("match.out","w",stdout);
int n,i;
scanf("%d",&n);
for(i=1;i<=n;i++)
{
scanf("%d",&A[i].x);
A[i].cnt=i;
}
for(i=1;i<=n;i++)
{
scanf("%d",&B[i].x);
B[i].cnt=i;
}
sort(A+1,A+1+n,cmp);
sort(B+1,B+1+n,cmp);
for(i=1;i<=n;i++)
aa[A[i].cnt] = B[i].cnt;
margef(1,n);
printf("%d\n",ans%99999997);
return 0;
}

my ac code：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#define maxn 100010
int n;
int aa[maxn],T[maxn];
typedef struct{
int val,id;
}node;
node a[maxn],b[maxn];
int ans;
bool cmp(node A,node B)
{
return A.val<B.val;
}
void marge(int l,int r)
{
if(l==r)return;
int mid=(l+r)/2;
marge(l,mid);
marge(mid+1,r);
int i=l,j=mid+1,k;
for(k=l;k<=r;k++)
{
if((i<=mid)&&((j>r)||aa[i]<=aa[j]))T[k]=aa[i++];
else
{
T[k]=aa[j++];
ans+=mid-i+1;
ans%=99999997;
}
}
for(i=l;i<=r;i++)aa[i]=T[i];
}
int main()
{
freopen("match.in","r",stdin);
freopen("match.out","w",stdout);
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i].val);
a[i].id=i;
}
for(int i=1;i<=n;i++)
{
scanf("%d",&b[i].val);
b[i].id=i;
}
std::sort(a+1,a+1+n,cmp);
std::sort(b+1,b+1+n,cmp);
for(int i=1;i<=n;i++)aa[a[i].id]=b[i].id;
marge(1,n);
printf("%d\n",ans%99999997);
return 0;
}