上图是归并排序的过程。
#include<bits/stdc++.h>
using namespace std;
const int N=1e5+5;
int a[N];
int b[N];
int n;
void f(int l,int r){
if(l>=r)return;
int mid=(l+r)>>1;
f(l,mid);
f(mid+1,r);
int i=l,j=mid+1,k=l;
while(i<=mid&&j<=r){
if(a[i]<=a[j])b[k++]=a[i++];
else b[k++]=a[j++];
}
while(i<=mid)b[k++]=a[i++];
while(j<=r)b[k++]=a[j++];
for(int p=l;p<=r;p++)a[p]=b[p];
}
signed main(){
scanf("%d",&n);
for(int i=1;i<=n;i++)scanf("%d",&a[i]);
f(1,n);
for(int i=1;i<=n;i++)printf("%d ",a[i]);
}
上面是归并排序的代码。在排序的函数中,变量 i i i遍历的是靠左的那半个区间,而变量 j j j遍历的是靠右的区间,所以 i < j i<j i<j总是成立,当这个时候,如果 a [ i ] > a [ j ] a[i]>a[j] a[i]>a[j],那么就出现了一对逆序对,然而,因为这是归并排序,所以这两端区间都是有序的,所以当 a [ i ] > a [ j ] a[i]>a[j] a[i]>a[j]成立时, a [ i + 1 ] > a [ j ] a[i+1]>a[j] a[i+1]>a[j]也成立( i i i + + + 1 1 1 < m i d <mid <mid),所以, a [ j ] a[j] a[j]对整个答案做出的贡献就是 m i d − i + 1 mid-i+1 mid−i+1个逆序对。以此类推。
所以我们只需要在 a [ i ] > a [ j ] a[i]>a[j] a[i]>a[j]的时候用一个变量统计数量就行了
if(a[i]>a[j]){
b[k++]=a[j++];
ans+=mid-i+1;
}
完整代码:
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=5e5+5;
int n;
int a[N];
int b[N];
int ans;
void mergesort(int l,int r){
int mid=l+r>>1;
if(l>=r)return;
mergesort(l,mid);
mergesort(mid+1,r);
int i=l,j=mid+1,k=l;
while(i<=mid&&j<=r){
if(a[i]>a[j]){
b[k++]=a[j++];
ans+=mid-i+1;
}
else{
b[k++]=a[i++];
}
}
while(i<=mid)b[k++]=a[i++];
while(j<=r)b[k++]=a[j++];
for(int i=l;i<=r;i++)a[i]=b[i];
}
signed main(){
ios::sync_with_stdio(0);
cin>>n;
for(int i=1;i<=n;i++)cin>>a[i];
mergesort(1,n);
cout<<ans;
}