算法实验-1-5_自然合并排序算法

先举个例子看一下自然合并排序和合并排序的区别：

可以看出自然合并排序在最初时就有已经排好序的子数组，再利用合并排序的原理对于已经排好序的子数组进行两两合并。
一个解释：

合并排序 :
是将两个（或两个以上）有序表合并成一个新的有序表
自然合并排序：
就是在合并排序时，建立子数组段，通过数组已排序的（自然排序,不重新排序）子数组段。

问题：

自然合并排序算法是对合并排序算法的一种改进。设a[0:n-1]是无序数组，用一次对数组a的扫描可以找出其中自然排好序的子数组，然后将相邻的排好序的子数组段两两合并，继续合并相邻排好序的子数组段，直至将整个数组排好序。

思路：
自然合并排序是合并排序的改进，应用的原理是：对于一个给定的数组a，通常会有长度大于1的有序子序列，我们称之为自然有序子序列.自然归并排序是找出给定数组中的自然子序列，将其开始的下标位置保存起来，然后对两个子序列之间进行插入排序形成一个大的有序子序列，直到整个数组排序完成。
大概思路我们首先定义一个全局变量n，控制数组的长度。在给定数组长度的前提，考虑数组子序列数最遭的情况就是整个数组对我们要求的序列刚好成一个反序列，即子序列的个数就是数组元素的个数。所以定义两个长度为n的数组 a、b，一个保存数组元素a[i]，另一个保存自然子序列的起始下标值b[i]。
下面的代码一共是1个主函数4个子函数，主函数首先调用GetIndex()函数，获得自然分组后每个子串的起始坐标，以及自然分组的个数，再在GetIndex()中调用MergeSort()函数将子序列划分，再在MergeSort()中调用MergePass()函数，将相邻子数组进行合并，再在MergePass()函数中调用Merge()函数合并对应索引的数组。

---

上代码：

#include<iostream>
using namespace std;
#include <stdlib.h>

// 这个函数记录自然分组后每个子串的起始坐标
void GetIndex(int a[], int t[], int n, int &m){
 
    
    /*a为输入数组，t为自然分组后每个分组的起始值组成的数组，
	n为数组长度 ，m自然分组个数*/ 
    int j = 0;
    t[j++] = 0; 
    for(int i = 0; i < n-1; i++) {
        if(a[i+1] < a[i]) {
            t[j++] = i+1;
        }
    }
    m = j;      // m为自然分组的个数
}

//(3)
// 合并c[l:m]和c[m+1:r]到d[l:r]
void Merge(int *c,int *d,int l,int m,int r)

{
	/*
	
	*/
	int i=l,j=m+1,k=l;//左右子数组的当前处理索引，合并后的索引
	while(i<=m && j<=r) {
		if(c[i]<=c[j])  d[k++]=c[i++];
		else d[k++]=c[j++];
	}
	if(i>m)//左子数组全部并入数组d，同时右子数组还有值未并入数组d
		for(int q=j; q<=r; q++)
			d[k++]=c[q];
	else
		for(int q=i; q<=m; q++)
			d[k++]=c[q];
}

 //(2)
 // 将x中长度为s的相邻子数组合并，结果存入y中
void MergePass(int *x,int *y,int *t,int s,int n,int cnt)
{
/*
x为输入的数组，y为长度为n的数组 
t是分组后子数组的索引的数组，cnt为自然分组个数 
*/ 


// n: 数组元素数量

	int i=0;  //子数组起始位置
	while(i+2*s<=cnt) {   //s为子数组长度 ，cnt为自然分组个数 
		int r=((i+2*s)<cnt) ?  t[i+2*s]:n; 
		//
		Merge(x,y,t[i],t[i+s]-1,r-1);
		i+=2*s;
	}
	if(i+s<cnt) { //剩下的元素少于2s
		//剩下的元素多于s
		Merge(x,y,t[i],t[i+s]-1,n-1);
	} 
	else {
		//剩下的元素不多于s
		for(int j=t[i]; j<n; j++)
			y[j]=x[j];
	}
 
}


 //(1):相邻两个数组合并 
void MergeSort(int *a,int *t,int n,int cnt)
{
	/*a为原数组 
	t为自然分组后每个分组的起始值组成的数组，
	n为数组长度 ，cnt自然分组个数
	*/
	int b[n];
	int s=1;//子数组长度
	while(s<cnt) {
		
		//将长度为s的相邻子数组两两合并（合并成的新数组内部有序）
		//并将这样的结果保存到数组b
		MergePass(a,b,t,s,n,cnt);
		
		//子数组的长度增加一倍
		s+=s;
		//同样使相邻的两个（内部有序）子数组合并成一个长度为"2s"的内部有序新数组保存到数组a 
		MergePass(b,a,t,s,n,cnt);
		s+=s;
	}
}

int main()
{
	int n;
	cout<<"输入数组元素个数：";
	cin>>n;
	int a[n];
	int t[n]; 
	cout<<"输入"<<n<<"个元素：";
	for(int i=0;i<n;i++){
		cin>>a[i];
	} 
	//Print(a,n);
	int cnt;        // 记录自然分组个数
    GetIndex(a, t, n, cnt);
	MergeSort(a,t,n,cnt);
	for(int i=0;i<n;i++) cout<<a[i]<<" "; 
	return 0;
}