武汉理工大学算法分析与设计 课内实验

---

前言

《算法设计与分析》课程实验的目的是为了使学生在课堂学习的同时，通过一系列的实验，使学生加深理解和更好地掌握《算法设计与分析》课程教学大纲要求的内容。

主要仪器设备及耗材

1.安装了Windows 10操作系统的PC机1台
2.PC机系统上安装了Microsoft Visual Studio 2019开发环境

---

以下是本篇文章正文内容，下面案例可供参考

一、分治法的应用

问题描述

i. 中位数问题
 问题描述
设X[ 0 : n - 1]和Y[ 0 : n – 1 ]为两个数组，每个数组中含有n个已排好序的数。找出X和Y的2n个数的中位数。
 编程任务
利用分治策略试设计一个O (log n)时间的算法求出这2n个数的中位数。
 数据输入
由文件input.txt提供输入数据。文件的第1行中有1个正整数n（n<=200），表示每个数组有n个数。接下来的两行分别是X，Y数组的元素。
 结果输出
程序运行结束时，将计算出的中位数输出到文件output.txt中。

输入文件示例 输出文件示例
input.txt output.txt
3 14
5 15 18
3 14 21

 实现提示
比较两个序列的中位数大小，如果两个数相等，则该数为整个2n个数据的中位数，否则通过比较，分别减少两个序列的查找范围，确定查找的起止位置，继续查找。

ii. Gray码问题
 问题描述
Gray码是一个长度为2n的序列。序列中无相同的元素，每个元素都是长度为n位的串，相邻元素恰好只有一位不同。用分治策略设计一个算法对任意的n构造相应的Gray码。
 编程任务
利用分治策略试设计一个算法对任意的n构造相应的Gray码。
 数据输入
由文件input.txt提供输入数据n。
 结果输出
程序运行结束时，将得到的所有编码输出到文件output.txt中。

输入文件示例 输出文件示例
input.txt
3
output.txt
0 0 0
0 0 1
0 1 1
0 1 0
1 1 0
1 1 1
1 0 1
1 0 0

 实现提示
把原问题分解为两个子问题，分别对两个子问题的每个数组后一位加0和1。

设计思想

分治策略
（1） 了解用分治法求解的问题：当要求解一个输入规模为n，且n的取值相当大的问题时，如果问题可以分成k个不同子集合，得到k个不同的可独立求解的子问题，其中1<k≤n，而且子问题与原问题性质相同，原问题的解可由这些子问题的解合并得出。那末，对于这类问题分治法是十分有效的。
（2） 掌握分治法的一般控制流程。

DanC(p,q)
  global n，A[1:n]; integer m,p,q;  // 1<=p<=q<=n
  if Small(p,q) then return G(p,q);
  else m=Divide(p,q); // p<=m<q
      return Combine(DanC(p,m),DanC(m+1,q));
  endif
end DanC

（3） 实现典型的分治算法的编程与上机实验，验证算法的时间复杂性函数。

伪代码描述(中位数问题)

输入：由文件input.txt提供输入数据。文件的第1行中有1个正整数n（n<=200），表示每个数组有n个数。接下来的两行分别是X，Y数组的元素。
输出：将计算出的中位数输出到文件output.txt中。
1.如果数组x,y都只有一个元素，则中位数等于这两个数的平均数
2.如果数组x,y有两个元素，因为数组元素按升序排序，则比较较小数组的最后一个数和较大数组的最后一个数来判断中位数等于哪两个数的平均数
3.当数组元素数较大时，则通过减小查询范围，查找的起止位置，继续查找。如果是偶数，则除以2，如果是奇数，则除以2加1。如果得到的两个位置的数相等，则该数为整个2n个数据的中位数，如果小于，则从x数组后一位开始，继续调用该函数；大于同理。

时间性能分析（中位数问题）

由于每次都是将数组分成两部分寻找中位数，从n,n/2,n/4,n/8……一直到1，所以时间复杂度为O(log2n).

源码（中位数问题）

#include<iostream>
#include<fstream>
using namespace std;
int zhongweishu(int *a, int *b,int n) {
   

	if (n==1) return (*a+*b)/2;//x,y数组只有一个元素
	else if (n == 2) {
   
		if (*a >= *b) {
   
			 if(*(b+1)>*(a+1))return (*a + *(a + 1)) / 2;
			 else return (*a + *(b + 1)) / 2;
		}
		else {
   
			if (*(a + 1) > * (b + 1))return