二分的应用——归并排序

S · B · Y · X

于 2023-08-27 14:30:30 发布

阅读量388

点赞数

分类专栏：算法文章标签：排序算法算法

本文链接：https://blog.csdn.net/XieYilang/article/details/132522879

版权

算法专栏收录该内容

7 篇文章 1 订阅

订阅专栏

前几天我们讲了二分，那么今天我们就来讲讲它的应用——归并排序。

1.归并排序

归并排序的思想

归并排序其实就是把一个数列二分成很多份，在合并起来。具体怎么拆分和合并，我们等一下会讲。

因为归并排序是用了二分的思想，所以它的时间复杂度为 $\Theta (log_2 n)$ 。

归并排序是稳定的。

拆分

拆分这一步其实很简单，就是两两的分开，直接看代码就好了。

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.例题

洛谷 P1908 逆序对

题目描述

猫猫 TOM 和小老鼠 JERRY 最近又较量上了，但是毕竟都是成年人，他们已经不喜欢再玩那种你追我赶的游戏，现在他们喜欢玩统计。

最近，TOM 老猫查阅到一个人类称之为“逆序对”的东西，这东西是这样定义的：对于给定的一段正整数序列，逆序对就是序列中 ai>aj 且 i<j 的有序对。知道这概念后，他们就比赛谁先算出给定的一段正整数序列中逆序对的数目。注意序列中可能有重复数字。

Update:数据已加强。

输入格式

第一行，一个数 n，表示序列中有 n个数。

第二行 n 个数，表示给定的序列。序列中每个数字不超过 109109。

输出格式

输出序列中逆序对的数目。

输入输出样例

输入 #1

6
5 4 2 6 3 1

输出 #1

说明/提示

对于 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);
}