前几天我们讲了二分,那么今天我们就来讲讲它的应用——归并排序。
1.归并排序
归并排序的思想
归并排序其实就是把一个数列二分成很多份,在合并起来。具体怎么拆分和合并,我们等一下会讲。
因为归并排序是用了二分的思想,所以它的时间复杂度为。
归并排序是稳定的。
拆分
拆分这一步其实很简单,就是两两的分开,直接看代码就好了。
void msort(int l,int r){
if(l==r) return;
int m=(l+r)/2;
msort(l,m);
msort(m+1,r);
}
合并
接下来就是合并,我们举个例子。
假设有两数列{2,3}和{1,4}。
首先我们定一个指针p和另外两个指针i和j,再来一个记录合并好之后的数列的数组k。
首先a[i]>a[j],于是我们就把a[j]放进k里。
k[p++]=a[i++];
现在i=1,j=2;
a[i]<a[j],所以把a[i]放进k里,然后i++;
k[p++]=a[i++];
a[i]还是<a[j],于是
k[p++]=a[i++];
结束了吗?还没有!
j还没到终点!
k[p++]=a[j++];
那么这个时候,可以看到,我们的k数组已经处于有序的状态,所以这就是合并了。
void msort(int l,int r){
if(l==r) return;
int m=(l+r)/2;
msort(l,m);
msort(m+1,r);
int i=l,j=m+1,p=l; //三个指针
while(i<=m&&j<=r){
if(a[i]<=a[j]) k[p++]=a[i++]; //如果a[i]<=a[j],把a[i]放进k里,i指针向后移一位。
else k[p++]=a[j++]; //否则把a[j]放进k里,j指针向后移一位。
}
while(i<=m) k[p++]=a[i++]; //如果i还没到终点,则把i后面的全部堆进k数组里。
while(j<=r) k[p++]=a[j++]; //如果j还没到终点,则把j后面的全部堆进k数组里。
for(int o=l;o<=r;o++) a[o]=k[o]; //把排序好的数组放回原数组里。
}
前面合并的时候,已经保证两个后面要合并的数组是有序的,所以直接合并就好了。
完整代码
#include<bits/stdc++.h>
using namespace std;
int a[100010],k[100010],n,ans;
void msort(int l,int r){
if(l==r) return;
int m=(l+r)/2;
msort(l,m);
msort(m+1,r);
int i=l,j=m+1,p=l;
while(i<=m&&j<=r){
if(a[i]<=a[j]) k[p++]=a[i++];
else k[p++]=a[j++];
}
while(i<=m) k[p++]=a[i++];
while(j<=r) k[p++]=a[j++];
for(int o=l;o<=r;o++) a[o]=k[o];
}
int main(){
scanf("%d",&n);
for(int i=0;i<n;i++) scanf("%d",&a[i]);
msort(0,n-1);
for(int i=0;i<n;i++) printf("%d ",a[i]);
}
2.例题
题目描述
猫猫 TOM 和小老鼠 JERRY 最近又较量上了,但是毕竟都是成年人,他们已经不喜欢再玩那种你追我赶的游戏,现在他们喜欢玩统计。
最近,TOM 老猫查阅到一个人类称之为“逆序对”的东西,这东西是这样定义的:对于给定的一段正整数序列,逆序对就是序列中 ai>aj 且 i<j 的有序对。知道这概念后,他们就比赛谁先算出给定的一段正整数序列中逆序对的数目。注意序列中可能有重复数字。
Update:数据已加强。
输入格式
第一行,一个数 n,表示序列中有 n个数。
第二行 n 个数,表示给定的序列。序列中每个数字不超过 109109。
输出格式
输出序列中逆序对的数目。
输入输出样例
输入 #1
6 5 4 2 6 3 1
输出 #1
11
说明/提示
对于 25% 的数据,n≤2500
对于 50% 的数据,n≤4×10^4。
对于所有数据,n≤5×10^5
请使用较快的输入输出
题解
这道题乍一看好像用冒泡排序可以做,每交换一个数就cnt++嘛。但是这个数据范围不允许我们这么做。
这时候,我们就要用到归并排序。
我们在合并的时候,遇到a[i]>a[j]时,那么就肯定是产生了逆序对,而且数量为m-i+1。
因为你的a[i]到a[m]是有序的,所以a[i]>a[j]就说明a[i+1],a[i+2],a[i+3],一直到a[m]都是要大于a[j]的,所以这一次产生的逆序对就为m-i+1个。
2 3
i
1 4
j
(2,3都与1构成了逆序对,因为a[i]到a[m]都大于等于a[i],而a[i]>a[j],所以a[i...m]>a[j],ans+=m-i+1)
通过这个方法,我们就可以在较短的时间内找出所有的逆序对。
(值得注意的是,由于逆序对数量太多,得用long long来存ans)
完整代码
#include<bits/stdc++.h>
using namespace std;
long long a[500010],k[500010],n,ans;
void msort(long long l,long long r){
if(l==r) return;
long long m=(l+r)/2;
msort(l,m);
msort(m+1,r);
long long i=l,j=m+1,p=l;
while(i<=m&&j<=r){
if(a[i]<=a[j]) k[p++]=a[i++];
else k[p++]=a[j++],ans+=m-i+1;
}
while(i<=m) k[p++]=a[i++];
while(j<=r) k[p++]=a[j++];
for(long long o=l;o<=r;o++) a[o]=k[o];
}
int main(){
scanf("%lld",&n);
for(long long i=0;i<n;i++) scanf("%lld",&a[i]);
msort(0,n-1);
printf("%lld\n",ans);
}
3.总结
归并排序其实就是对二分的应用,算是一个比较快的排序方法,而且还稳定。所以总结来说归并排序还是个不错的排序方法。就是代码有点麻烦。
关于二分,请看:【C++算法讲解】整数二分