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里的函数进行了测试。