CICC科普栏目|编程中的基本数据结构与算法思想
编程的关键在于选择数据结构和算法,数据结构用于描述问题,算法用于描述解决问题的方法和步骤。
描述问题的数据除了各数据元素本身,还要考虑各元素的逻辑关系,主要是一对一的线性关系,一对多的树型关系和多对多的图形关系。另外,内存中对各数据元素的存储只有顺序存储和链式存储两种方式,所以数据结构还要考虑数据的存储结构,并考虑逻辑结构与数据结构如何有效地结合到一起。
用算法描述问题,当问题比较复杂时,通常的思路是分而治之,并辅以适当的数据结构。
1 分治法Divide and Conquer
分治法通常描述为以下三步:
Divide the problem into more subproblems(分解问题为众多的子问题);
Conuqe(solve) the subproblems(解决各子问题);
Combine(merge) the solution of subproblems(if need)(合并各子问题的解(如果需要)).
如用分治法来计算2^10?
2^10=2^5*x^5=2^2*x^3*x^5=32*32=1024
相对于顺序查找,二分查找有更高的效率,前提是二分查找需要事先排好序:
int binarySearchLoop(int arr[], int len, int findData)
if(arr==NULL || len =0)
return -1;
int start = 0;
int end = len-1;
while(start=end)
int mid = start+(end-start)/2;
if(arr[mid] == findData)
return mid;
else if(findData arr[mid])
end = mid-1;
else
start = mid+1;
return -1;
2 枚举法也是一种暴力缩小问题规模的算法
简单的枚举算法也是可以优化的,即尽可能缩小搜索的空间,如判断质数:
质数(prime number)又称素数,有无限个。质数定义为在大于1的自然数中,除了1和它本身以外不再有其他因数。
判断质数的函数:
int isPrime(int n)
if(n= 1)// 小于等于1的整数不可能是素数
return 0;
if(n == 2); // 2 是素数
return 1;
if(n%2 == 0); // 能被2整除的其他整数都不是素数
return 0;
int limit = (int)sqrt((double)n)+1;
for(int i = 3; i = limit; i=i+2)
if(n % i == 0)
return 0;
return 1;
isPrime没有必要枚举所有的因子。
I 只要发现任何一个大于1小于n的因子,就能停下来报告n不是素数。
II 如果n能被2整除,直接报告n不是素数。如果n不能被2整除,那么它也不可能被4或6或其他偶数整除。因此,isPrime只需要检查2和奇数(由3开始,步长为2)。但注意有个特例,2能被2整除,但2是素数。
III 如果n不是素数,则必有一个因子小于√ n 。因此不需要检查到n为止。只需检查到√ n(n=√ n*√ n) 。
因为如果n能被2~n-1之间任一整数整除,其二个因子必定有一个小于或等于√n,另一个大于或等于√n。例如24可以表示为:2*12、3*8、4*6,前面的因子小于√24,后面的因子大于√24,检验出了小因子,即可判断n是否为素数,就像逻辑运算的短路求值。
3 程序的模块化
分治法在程序思想中的应用就是实现程序的模块化,包括面向过程的函数化和面向对象的对象化。
许多原因都促使我们将应用程序分解成函数,下面仅列举其中三个:
函数一般小而具体。用一系列函数来写程序,胜于一气呵成写完整个程序。这称为“分而治之”,使你的精力一次集中在一个函数上。
包含许多小函数的应用程序比单一的长程序更容易阅读和调试。
函数可以重用。函数写好后可在程序的其他任何地方调用。这减少了编码量,提高了开发效率。
4 函数调用与栈
首先讨论一个从a点出发去f点,然后回到a点的问题(中间的b、c、d、e都有多个分岔口):
a→b2→c1→d3→e2→f,每个分岔口都有一个信封,告诉你应该走哪一个分支,为了能够正确地回到起点a,正确的做法是拿到一个信封后,即将这个信封叠在上一次拿到的信封的上面,回去时,依次从上面拿取信封,按提示即可正确返回。
其做法就是依次放入,依次取出,信封之间是顺序关系,只在一端操作,也就是不管是放入还是取出都不在中间操作。这样一种思路在计算机上用数据来描述就是后进先出的栈,函数的调用、返回,递归、回溯算法都需要使用栈这种数据结构(由程序员或递归时由编译器来实现)。
在C++中,函数不能嵌套定义,但可以嵌套调用,在函数调用时,编译器需要确保在逐级调用后能够回归到最初的调用点,编译器会隐式实现一个堆栈,用来保存每一级函数调用时的函数返回地址和局部变量,依次入栈和出栈。
C++也支持递归函数的递归调用,同样是由编译器隐式地实现了一个堆栈。
5 深度搜索与广度搜索
如果将上述的问题稍微扩展一点,要从源点到目标点,中间的节点可能有多个分叉,这样的问题可以用一个树或图来描述。
而探路的方法可以分为两种,一种是深度优先搜索(下一点、下一点……回溯……),一种是广度优先搜索(下一点的全部分叉、下一点的全部分叉……):
5.1 深度优先搜索用栈(stack)来实现,整个过程可以想象成一个倒立的树形:
1)把根节点压入栈中。
2)每次从栈中弹出一个元素,搜索所有在它下一级的元素,把这些元素压入栈中。并把这个元素记为它下一级元素的前驱。
3)找到所要找的元素时结束程序。
4)如果遍历整个树还没有找到,结束程序。
5.2 广度优先搜索使用队列(queue)来实现,整个过程也可以看做一个倒立的树形:
1)把根节点放到队列的末尾。
2)每次从队列的头部取出一个元素,查看这个元素所有的下一级元素,把它们放到队列的末尾。并把这个元素记为它下一级元素的前驱。(取出的元素也可以保存到一个队列)
3)找到所要找的元素时结束程序。
4)如果遍历整个树还没有找到,结束程序。
广度优先搜索相对于深度优先搜索,因为是逐层探索的,可以确保以较少的点到达目标点,缺点是存储量较大。
6 递归算法
递归就是某个函数直接或间接的调用自身。
语法形式上: 在一个函数的运行过程中, 调用这个函数自己:
直接调用: 在fun中直接执行fun;
间接调用: 在fun1中执行fun2; 在fun2中又执行fun1 ;
问题的求解过程是划分成许多相同性质的子问题的求解,而小问题的求解过程可以很容易的求出。这些子问题的解就构成里原问题的解。
待求解问题的解可以描述为输入变量x的函数f(x)。
通过寻找函数g( ),使得f(x) = g(f(x-1))。
且已知f(0)的值, 就可以通过f(0)和g( )求出f(x)的值。
扩展到多个输入变量x, y, z等, x-1也可以推广到 x - x1 , 只要递归朝着 “出口” 的方向即可。
递归算法分解出的子问题与原问题之间是纵向的, 同类的关系(枚举分解出的子问题之间是横向的, 同类的关系)。
递归的三个要点:
递归式:如何将原问题划分成子问题;
递归出口:递归终止的条件, 即最小子问题的求解,可以允许多个出口;
界函数:问题规模变化的函数, 它保证递归的规模向出口条件靠拢。
如一个求阶乘的递归程序,给定n, 求阶乘n!
阶乘的栈:
二分搜索的递归实现:
int binarySearchRecursion(int arr[], int findData, int start, int end)
if(arr==NULL || startend)
return -1;
int mid = start+(end-start)/2;
if(arr[mid] == findData)
return mid;
else if(findData arr[mid])
binarySearchRecursion(arr, findData, start, mid-1);
else
binarySearchRecursion(arr, findData, mid+1, end);
7 归并排序
归并排序(merge sort)是建立在归并操作上的一种有效的排序算法。该算法是分治法(Divide and Conquer)的一个非常典型的应用。将已有序的子序列合并,得到完全有序的序列;即先使每个子序列有序,再使子序列段间有序。若将两个有序表合并成一个有序表,称为2-路归并(2-way or binary merges sort)。
归并排序在1945年由冯·诺伊曼首次提出。
2-路归并的基本思路就是将数组分成二组A,B,如果这二组组内的数据都是有序的,那么就可以很方便的将这二组数据进行排序。如何让这二组组内数据有序?
可以将A,B组各自再分成二组。依次类推, 当分出来的小组只有一个数据时,可以认为这个小组组内已经达到了有序,然后再合并相邻的二个小组就可以了。这样 通过先递归的分解数列,再合并数列就完成了归并排序。
归并排序的效率是比较高的,设数列长为N,将数列分开成小数列一共要logN步,每步都是一个合并有序数列的过程,时间复杂度可以记为O(N),故一共为O(N*logN)。因为归并排序每次都是在相邻的数据中进行操作,所以归并排序在O(N*logN)的几种排序方法(快速排序,归并排序,希尔排序,堆排序)也是效率比较高的。
归并排序的实现分为 递归实现与 非递归(迭代)实现。递归实现的归并排序是算法设计中分治策略的典型应用,我们将一个大问题分割成小问题分别解决,然后用所有小问题的答案来解决整个大问题。非递归(迭代)实现的归并排序首先进行是两两归并,然后四四归并,然后是八八归并,一直下去直到归并了整个数组。
7.1 归并排序 分解
可以看到这种结构很像一棵完全二叉树, 分阶段可以理解为就是递归拆分子序列的过程,递归深度为log2n。
7.2 归并排序 合并相邻有序子序列
再来看看 并阶段,我们需要将两个已经有序的子序列合并成一个有序序列,比如上图中的最后一次合并,要将[4,5,7,8]和[1,2,3,6]两个已经有序的子序列,合并为最终序列[1,2,3,4,5,6,7,8],来看下实现步骤。
申请空间,使其大小为两个已经排序序列之和,该空间用来存放合并后的序列;
设定两个指针,最初位置分别为两个已经排序序列的起始位置;
比较两个指针所指向的元素,选择相对小的元素放入到合并空间,并移动指针到下一位置;temp[index++] = A[i] = A[j] ? A[i++] : A[j++];
重复步骤3直到某一指针到达序列尾;
将另一序列剩下的所有元素直接复制到合并序列尾;
7.3 归并排序 动图演示
7.4 归并排序 代码
8 回溯法和分书问题
回溯算法实际上是一个类似枚举的搜索尝试过程,主要是在搜索尝试过程中寻找问题的解,当发现已不满足求解条件时,就“回溯“返回,尝试别的路径。可以参考一下走迷宫的过程,一开始会随机选择一条道路前进,一直到走不通之后就会回头直到找到另外一条没有试过的道路前进。实际上,走迷宫的算法就是回溯法的经典问题。
回溯法实际上 也是一种试错的思路,通过不断尝试解的组合来达到求解可行解和最优解的目的。虽然都有穷搜的概念蕴含其中,但是回溯法和穷举查找法是不同的。对于一个问题的所有实例,穷举法注定都是非常缓慢的,但应用回溯法至少可以期望对于一些规模不是很小的实例,计算机在可接受的时间内对问题求解。
许多复杂的规模的问题都可以使用回溯法,有”通用解题方法”的美称。分书问题和八皇后都是典型的回溯法问题。
分书问题能够较有代表性地表现数据描述、递归、回溯的算法思路。
有编号为0,1,2,3,4的5本书,准备分给5个人A,B,C,D,E,写一个程序,输出所有皆大欢喜的分书方案。
每个人的阅读兴趣用一个二维数组like描述:
Like[i][j] = true i喜欢书j
Like[i][j] = false i不喜欢书j
设计一个函数trynext(int i)给第i个人分书。
用一个一维数组take表示某本书分给了某人。take[j]=i+1;//把第j本书分配给第i个人
依次尝试把书j分给人i。
如果第i个人不喜欢第j本书,则尝试下一本书,如果喜欢,并且第j本书尚未分配,则把书j分配给i。
如果i是最后一个人,则方案数加1,输出该方案。否则调用trynext(i+1)为第i+1个人分书。
如果对第i个人枚举了他喜欢的所有的书,都没有找到可行的方案,那就回到前一个状态i-1,让i-1把分到的书退回去,重新找喜欢的书,再递归调用函数,寻找可行的方案。
#include iostream
#include conio.h
using namespace std;
int like[5][5]={
{0,0,1,1,0},
{1,1,0,0,1},
{0,1,1,0,1},
{0,0,0,1,0},
{0,1,0,0,1}};
int take[5]={0,0,0,0,0};//记录每一本书的分配情况
int n;//n表示分书方案数
void trynext(int i);
int main
n=0;
trynext(0);
getch;
return 0;
//对第 i 个人进行分配
void trynext(int i)
int j,k;
for(j=0;j5;j++)
if(like[i][j]take[j]==0)
take[j]=i+1;//把第j本书分配给第i个人
if(i==4)//第5个人分配结束,也即所有的书已经分配完毕,可以将方案进行输出
n++;
cout"第"n"种分配方案"endl;
for(k=0;k5;k++)
cout"第"k"本书分配给"(char)(take[k]+'A'-1)endl;
coutendl;
else
trynext(i+1);//递归,对下一个人进行分配
take[j]=0;//回溯,寻找下一种方案
当like矩阵的值为
附归并排序的代码:
#include stdio.h
#include stdlib.h
#include limits.h
// 分类 -------------- 内部比较排序
// 数据结构 ---------- 数组
// 最差时间复杂度 ---- O(nlogn)
// 最优时间复杂度 ---- O(nlogn)
// 平均时间复杂度 ---- O(nlogn)
// 所需辅助空间 ------ O(n)
// 稳定性 ------------ 稳定
// 合并两个已排好序的数组A[left...mid]和A[mid+1...right]
void Merge(int A[], int left, int mid, int right)
int len = right - left + 1;
int *temp = new int[len]; // 辅助空间O(n)
int index = 0;
int i = left; // 前一数组的起始元素
int j = mid + 1; // 后一数组的起始元素
while (i = mid j = right)
temp[index++] = A[i] = A[j] ? A[i++] : A[j++]; // 带等号保证归并排序的稳定性
while (i = mid)
temp[index++] = A[i++];
while (j = right)
temp[index++] = A[j++];
for (int k = 0; k len; k++)
A[left++] = temp[k];
// 递归实现的归并排序(自顶向下)void MergeSortRecursion(int A[], int left, int right){if (left == right) // 当待排序的序列长度为1时,递归开始回溯,进行merge操作return;int mid = (left + right) / 2;MergeSortRecursion(A, left, mid); //左半部分排好序MergeSortRecursion(A, mid + 1, right); //右半部分排好序Merge(A, left, mid, right); //合并左右部分}// 非递归(迭代)实现的归并排序(自底向上)void MergeSortIteration(int A[], int len){int left, mid, right;// 子数组索引,前一个为A[left...mid],后一个子数组为A[mid+1...right]for (int i = 1; i len; i *= 2) // 子数组的大小i初始为1,每轮翻倍{left = 0;while (left + i len) // 后一个子数组存在(需要归并){mid = left + i - 1;right = mid + i len ? mid + i : len - 1;// 后一个子数组大小可能不够Merge(A, left, mid, right);left = right + 1; // 前一个子数组索引向后移动}}}int main{int A1[] = { 6, 5, 3, 1, 8, 7, 2, 4 }; // 从小到大归并排序int A2[] = { 6, 5, 3, 1, 8, 7, 2, 4 };int n1 = sizeof(A1) / sizeof(int);int n2 = sizeof(A2) / sizeof(int);MergeSortRecursion(A1, 0, n1 - 1); // 递归实现MergeSortIteration(A2, n2); // 非递归实现printf("递归实现的归并排序结果:");for (int i = 0; i n1; i++){printf("%d ", A1[i]);}printf("");printf("非递归实现的归并排序结果:");for (i = 0; i n2; i++){printf("%d ", A2[i]);}printf("");system("pause");return 0;}
- END -
来源:机器人网
数学中国站长多年收集的超400G数学建模资料免费分享,已经做了详细的分类,使用非常方便!
链接:
提取码:5V7D
复制内容打开移动云盘PC客户端,操作更方便哦
关注公众号了解更多
会员申请 请在公众号内回复“个人会员”或“单位会员
欢迎关注中国指挥与控制学会媒体矩阵
CICC官方网站
CICC官方微信公众号
《指挥与控制学报》官网
国际无人系统大会官网
中国指挥控制大会官网
全国兵棋推演大赛
全国空中智能博弈大赛
搜狐号
一点号