work-10
0. 问题描述
见老师博客
1.架构简介
经过软件工程的课程,我将学到的很多知识应用到了这次作业中首先,我从架构上来讲解下我的这次作业。
由于各个语言优势不相同,例如在C++课上我们讲到了C++的尴尬地位(但其实是一个折衷)。所以在这次作业中我也用到了类似思想,例如对于需要效率并不高的UI我采用了Python的PYQT界面的自动生成与方便的类型封装快速帮我解决了程序编写问题。但是对于效率依赖较高的的算法部分,我认为使用C会更快速解决问题,算法实现效率更好。
所以我采用的是界面由python编写,然后使用gcc编译出的DLL文件进行调用。以此发挥各个语言的优势。而在软件工程中我认为更为重要的一点就是封装,经过我了解,我才用了dll+.h文件的代码交流方式,对于代码的实现细节使用方并不需要了解,更为重要的是,这也是对代码创作者版权的一种保护。
2界面简介
界面简洁明了,左边的网格可供输入,也可以通过load一个矩阵文件来进行,通过点击GO!表示确认右边的矩阵,点击Reset将程序恢复到初始状态。三个选项钩分别表示不定形,垂直相连,水平相连的三个参数。然后4个播放功能键意义也很明了,分别表示播放,从头来,上一步,下一步。
在逻辑层面上封闭了非法操作,如当执行到最后一个状态,>>和>都会被锁定,保护程序的正确性。
3.算法简介
3.1 original.h
#ifndef ORIGINAL_H_INCLUDED
#define ORIGINAL_H_INCLUDED
#define SOLVE_CHOSEN -10001
/*
该文件包含了一系列求对象最大子对象的函数
*/
int maxsumline(int *p,int size);
/*
求一个序列p中包含的最大子序列,并将被选择的位置置为SOLVE_CHOSEN
*/
int maxsumcycle(int *p,int size);
/*
求一个循环序列p中包含的最大子序列,并将被选择的位置置为SOLVE_CHOSEN
循环序列代表着该序列首尾相连。
*/
int maxsumblock(int a[],int n,int m,int cycle,int expand);
/*
求一个矩阵的最大子矩阵,并将被选择的位置置为SOLVE_CHOSEN
n,m分别是它的长宽,cycle为1时代表该矩阵的左右是相连的,expand为1时代表着该矩阵的上下是相连的,并将
*/
#endif
问题分析
在一维情况下我们已经分析得到了基于长度n的O(n)时间复杂度的算法。那么我们可以先考虑在二维情况下是否可以得到基于长度n宽度m的O(m)时间复杂度的算法。如我在作业1里分析。设s[x][y]为以坐标(0,0)为左上角,(x,y)为右下角的点所形成的的矩形的加和。以(a,b)(x,y)构成的矩形的值为,(s[x][y] - s[a-1][y])-(s[x][b-1] - s[a-1][b-1]),不具备一维时的单调性,只能通过在此枚举一行。时间复杂度为O(m*n*n),无法达到最好的O(m*n)。
也就是说对于普通的问题,我们只需要枚举2行的组合即先枚举i再枚举小于等于i的j,加和j-i的区间,1维处理就可以了。
而对于-v的垂直相连参数,是很容易转化为普通问题的,普通问题中只考虑j-i的区间,而这里再考虑下i-n与0-j的区间就可以了,时间复杂度也为O(m*n*n)。
而对于-h的水平相连参数,我们可以从转化出的一维问题中考虑。对于1维情况下如果收尾相连应该如何处理。一开始我考虑的是复制一遍贴在右边,但其实实现起来限制过于复杂。如果选择了超越经线0的矩形其实就是踢掉了中间的一块矩形,于是只需要找到最小的矩形,然后用正行的权减去它,与普通解想比较即可。
3.2 atype.h
#ifndef ATYPE_H_INCLUDED
#define ATYPE_H_INCLUDED
#define SOLVE_CHOSEN -10001
#include <stdio.h>
#define N 32
#define M 32
#include <math.h>
#include <stdlib.h>
#include <time.h>
void SAA(int v,float T,float r,float Tmin);
/*
退火算法,v温度取的次数,T表示温度,r表示冷却速率,Tmin表示最终温度
*/
int deal(int o,int p,int c[],int vertical,int horizontal);
/*
使用的是退火算法求出一个相对较优解,o,p代表着矩阵c的长宽,vertical为1表示矩阵垂直相连,horizontal为1表示矩阵水平相连。返回得到块的值。
并将选中的点置为SOLVE_CHOSEN
执行过以次deal后才可以使用之后的单步功能
*/
int get_situation(int idx);
/*
返回第i次的退火解
*/
#endif
对于-a参数,这个历史遗留问题。经过并查集缩点之后(将正权点加合在一起作为一个点,并认为它的负权为0),建立一个Graph,然后通过随机挪动其中的一个点,即选还是不选作为状态转移可否的依据。然后再通过SAA的当前温度给出一个概率,表示新状态没有原来好但依然接受的概率。
3.3返回中间结果
在最最开始老师刚说的时候,我确实是想边运行边输出,但当老师说了可以先预处理一遍然后再进行帧操作,我决定使用这个方法,即在execute的时候将中间结果保存到一个文件里,然后进行运算。
4代码优化与注释
4.1注释见Github
4.2代码优化
4.2.1命名优化
首先有一个极为严重的问题往往被我们忽视,那就是build-in function,这是c语言内嵌函数。我在代码编写中写出了一个warning提示我声明y1这个变量在c里是个build-in function,但我一开始忽略了该问题,结果导致文件无法输出,至此我也不知道这两者之间的联系。不过我认为可以说明一些问题,那就是对于自己不了解的一些warning恰恰应该给出足够的重视,而对自己了然的一些warning或许可以采取忽视的态度,而不能因为连这个warning都不知道是什么意思就忽略掉了。
命名上由于C是面向过程的,很多地方写的东西只是该层函数所需要的,但我对于全局变量的把握依然不高,但是算法上的东西真的有必要去追求完美的拓展性吗,我觉得这个问题有待商榷,更多的我觉得只需要帮列在.h文件里的那些函数命名规范好就够了,因为其他人根本就不需要看里面的东西呀。
4.2.2 MagicNumber
把所有的魔数都转化成了一个常数类型。
4.2.3功能分化
把原本放在一起的两个dll分别以atype.dll与original.dll的两个文件进行完整剥离。另外在两个.c文件里我写了main函数用来测试。
4.2.4代码覆盖率
非常有意义的一项测试,我只是尝试了几个较小的数据,对dll里的函数进行了测试。