目录
实验目的
理解分治法的算法思想,阅读实现书上已有的部分程序代码并完善程序,加深对分治法 的算法原理及实现过程的理解。
实验内容
一、用分治法实现一组无序序列的两路合并排序和快速排序。要求清楚合并排序及快速排 序的基本原理,编程实现分别用这两种方法将输入的一组无序序列排序为有序序列后输出。
二、采用基于“五元中值组取中值分割法”(median-of-median-of-five partitioning)的线性 时间选择算法,找出 N 个元素集合 S 中的第 k 个最小的元素,使其在线性时间内解决。
实验步骤
一、两路合并排序和快速排序
1、排序是数据处理中常用的重要手段,是指将一个元素序列调整为按指定关键字值的递增 (或递减)次序排列的有序序列。
用分治法求解排序问题的思路是,按某种方式将序列分成两个或多个子序列,分别进行排序,再将已排序的子序列合并成一个有序序列。 合并排序和快速排序是两种典型的符合分治策略的排序算法。
2、两路合并排序
#include <iostream>
using namespace std;
int n ;
int *arr;
void Merge(int left,int mid ,int right) {
int *temp = new int[right-left+1];
int i=left,j=mid+1,k=0;
while(i<=mid && j<=right){
if(arr[i]<arr[j]) temp[k++]=arr[i++];
else temp[k++] = arr[j++];
}
while(i<=mid) temp[k++]=arr[i++];
while(j<=right) temp[k++]=arr[j++];
for(i=0;i<k;i++){
arr[left++]=temp[i];
}
}
void MergeSort(int left,int right) {
if(left<right) {
int mid = (left+right)/2;
MergeSort(left,mid);
MergeSort(mid+1,right);
Merge(left,mid,right);
}
}
void Input() {
cout << "请输入待排序的个数:";
cin >> n;
arr = new int [n];
cout << "请输入待排序的数组:";
for(int i=0; i<n; i++) {
cin >> arr[i];
}
}
void Output() {
for(int i=0; i<n; i++) {
cout << arr[i] << ' ';
}
}
int main() {
Input();
MergeSort(0,n-1);
Output();
return 0;
}
3、快速排序
#include <iostream>
using namespace std;
#define INF 2147483647
int n;
int *arr;
void Swap(int i,int j){
int temp = arr[i];
arr[i]=arr[j];
arr[j]=temp;
}
int Partition(int left,int right) {
int i=left,j=right+1;
do {
do i++;
while (arr[i]<arr[left]);
do j--;
while (arr[j]>arr[left]);
if(i<j) Swap(i,j);
} while(i<j);
Swap(left,j);
return j;
}
void QuickSort(int left,int right) {
if(left<right) {
int j = Partition(left,right);
QuickSort(left,j-1);
QuickSort(j+1,right);
}
}
void Input() {
cout << "请输入待排序的个数:";
cin >> n;
arr = new int [n+1];
cout << "请输入待排序的数组:";
for(int i=0; i<n; i++) {
cin >> arr[i];
}
arr[n] = INF;
}
void Output() {
for(int i=0; i<n; i++) {
cout << arr[i] << ' ';
}
}
int main() {
Input();
QuickSort(0,n);
Output();
return 0;
}
二、线性时间选择算法——寻找第 k 个最小元
#include <iostream>
using namespace std;
int n,m;
int *arr;
void Input() {
cout << "请输入待排序的个数:";
cin >> n;
arr = new int [n];
cout << "请输入待搜索的数组:";
for(int i=0; i<n; i++) {
cin >> arr[i];
}
cout << "请输入要查找的第几个最小元:";
cin >> m;
}
void Swap(int i,int j) {
int temp = arr[i];
arr[i]=arr[j];
arr[j]=temp;
}
void InsertSort(int left,int right) {
for(int i=left+1; i<=right; i++) {
int key=arr[i];
int j=i-1;
while(j>=0 && arr[j]>key) {
arr[j+1]=arr[j];
j--;
}
arr[j+1]=key;
}
}
int Partition(int left,int right) {
int i=left,j=right+1;
do {
do i++;
while(arr[left]>arr[i]);
do j--;
while(arr[left]<arr[j]);
if(i<j) Swap(i,j);
} while(i<j);
Swap(left,j);
return j;
}
int Ceil(int x,int y) {
return x%y==0?x/y:x/y+1;
}
int Select(int k,int left,int right,int r) {
int n=right-left+1;
if(n<=r) {//若问题足够小,使用直接插入排序,取其中得第k小元素,其下标为;left+k-1
InsertSort(left,right);
return left+k-1;
}
for(int i=1; i<=n/r; i++) {
InsertSort(left+(i-1)*r,left+i*r-1); //二次取中规则求每组的中间值
Swap(left+i-1,left+(i-1)*r+Ceil(r,2)-1);//将每组的中间值集中存放在子表前部
}
int j=Select(Ceil(n/r,2),left,left+(n/r)-1,r);//求二次中间值,其下标为j
Swap(left,j); // 二次中间值为主元,并换至left处
j=Partition(left,right); //对表(子表)进行分划操作
if(k==j-left+1) return j; //返回第K小元素下标
else if(k<j-left+1) return Select(k,left,j-1,r);//在左子表求第k小元素
else return Select(k-(j-left+1),j+1,right,r);//在右子表求第k-(j-left+1)小元素
}
int main() {
Input();
int index = Select(m,0,n-1,5);
cout << arr[index];
return 0;
}