软件工程第一次作业——数独的求解与生成
代码的GitHub地址:https://github.com/Liu-SD/sudoku
personal software process stages | 预估耗时 | 实际耗时 |
计划 | ||
估计这个任务需要多少时间 | 10 min | 10 min |
开发 | ||
需求分析(包括学习新技术) | 180 min | 190 min |
生成设计文档 | 0 min(没做设计文档) | 0 min |
设计复审(和同事审核设计文档) | 0 min(没有同事复审) | 0 min |
代码规范(为目前的开发制定合适的规范) | 30 min | 60 min |
具体设计 | 180 min | 180 min |
具体编码 | 240 min | 240 min |
代码复审 | 60 min | 180 min |
测试(自我测试,修改代码,提交修改) | 60 min | 50 min |
报告 | ||
测试报告 | 180 min | 180 min |
计算工作量 | 0 min(没有做这项工作) | 0 min |
事后总结,并提出过程改进计划 | 30 min | 30 min |
合计 | 1070 min | 1160 min |
解题思路描述
题目要求分为两个部分,一个是解数独,一个是自动生成数独。并且生成数独时对第一个数字做出了限制。所以可以认为自动生成数独是只有一个数字限制时的解数独。所以说这两者大部分的实现算法是相同的。解数独时是在找到第一个解的时候返回真,而生成数独是在找到第N个解时返回真。
因此问题转换为在R个限制条件(1<=R<=81)的情况下寻找到N(1<=N<=1000000)个解时返回真,否则返回假。解的输出可以把输出流作为参数传入,找到一个答案时将答案输出。这个问题是典型的回溯问题,在网上查资料了解到可以将问题转化为精确覆盖问题,使用DLX(dancing link)算法实现。DLX的重点在于十字链表的实现。
设计实现过程
实现算法的第一步需要使用十字链表表示稀疏矩阵。我把十字链表封装为cross_link类,这个类,类中包括build方法,建立空的矩阵。build函数在构造函数中调用。稀疏矩阵建立以后,head指向矩阵头,rows数组存放行的头指针,cols数组存放列的头指针。类中包括insert方法插入元素。delrow,delcol,recoverrow,recovercol四个方法删除或恢复整行或整列。这四个方法将是算法实现的关键。在删除行(列)时,修改每一行(列)中元素的邻居节点到另一侧,同时保留被删除行(列)的邻居指针,方便恢复。删除和恢复需要以栈的方式操作。即后删除的先恢复。
算法被封装到dlx类中,dlx类调用并维护cross_link类的实例。类中的find_one和find_many实现解数独和生成数独两个需求。find_one调用_find_one递归函数,find_many调用_find_many递归函数。
_find_one函数流程:
1. 如果head的右结点为空,即矩阵没有列,则将result数组结果输出,stack中的所有删除操作使用recover做逆操作,即恢复为最初始的矩阵。返回真。
2. 否则,找到元素最少的列min_col_index,对于这一列的每一个元素i:
将元素i所在行压如result栈中。
对于元素i所在行的每个元素j:
对于元素j所在列的每个元素p:
删除p所在的行,压如stack栈中。
删除j所在的列,压如stack栈中。
以stack和result栈的栈顶指针为参数递归调用自己,如果为真,返回真。
否则,把元素i从result栈弹出。
把这一次压如stack栈中的元素弹出并且recover函数恢复删除的行和列。
函数最后返回假。
_find_many函数和_find_one类似,在找到结果时将结果输出到文件,计数器加一,在计数器等于目标数目之前都一直返回假。
十字链表的搭建容易出错,而且增删操作对于之后的算法十分重要,因此单元测试的重点放在了十字链表类的测试。设计单元测试时,先测试十字链表的插入操作,每插入一个元素,就计算元素个数并且断言。之后测试链表的删除行(列)以及之后行(列)的复原。
程序性能改进
我在程序性能改进方面没有做很多的操作,仅仅将输出到文件的部分从逐个数字写入改为了一个数独盘的一次性写入。但这个改进大大提升了程序的性能。生成一百万个数独从之前的五分钟变为当前的十四秒。
于是,进一步的,我修改了从文件读入的方式,将之前的逐个数字读入改为逐行读入。
生成一百万个数独时,VS性能分析图:
执行次数最多的函数为:
由此可见递归调用_find_many的执行次数最多,对链表的四个操作十分频繁,对其进行优化十分重要。
上图是递归调用的函数关系图。
代码说明
十字链表类的关键代码为删除和恢复行或者列。
void cross_link::delrow(int r) {for (Cross i = rows[r]; i != NULL; i = i->right) { // 对于该行每一个元素cols[i->col]->count--; // 删除后该列的计数器减一if (i->up)i->up->down = i->down; // 在它上方存在元素的情况下,上方元素改为指向该元素下方元素if (i->down)i->down->up = i->up; // 下方元素同理} }void cross_link::recoverrow(int r) {for (Cross i = rows[r]; i != NULL; i = i->right) { // 对于该行每一个元素cols[i->col]->count++; // 恢复后该列的计数器加一if (i->up)i->up->down = i; // 在上方元素不空的情况下,将原本指向自己下方的上方元素指回自己if (i->down)i->down->up = i; // 下方元素同理} }
以上是行的删除和恢复一行元素的代码,删除恢复列的情形和一上代码类似。
dlx类的关键代码为函数的递归调用。以下是_find_many的代码:
1 bool dlx::_find_many(int stack_top, int result_pos, int N, int & n, std::ofstream &fout) 2 { 3 if (!head->right) { // 在链表不存在列的情况下, 4 int matrix[9][9]; 5 for (int i = 0; i < 81; i++) { // 将result数组里面的解翻译为9*9的矩阵 6 int j = result[i] - 1; 7 int val = j % 9 + 1; 8 int pos = j / 9; 9 int row = pos / 9; 10 int col = pos % 9; 11 matrix[row][col] = val; 12 } 13 char str[19 * 9 + 2] = { 0 }; 14 for (int i = 0; i < 9; i++) { // 用字符串记录矩阵并输出到文件 15 for (int j = 0; j < 9; j++) { 16 str[i * 19 + 2 * j] = matrix[i][j] + '0'; 17 str[i * 19 + 2 * j + 1] = ' '; 18 } 19 str[i * 19 + 18] = '\n'; 20 } 21 str[19 * 9] = '\n'; 22 str[19 * 9 + 1] = '\0'; 23 fout << str; 24 if (++n >= N) // 计数器加一,如果大于要求的数量,返回真,否则返回假 25 return true; 26 else 27 return false; 28 } 29 int min_col_count = 100; 30 int min_col_index = -1; 31 for (Cross p = head->right; p != NULL; p = p->right) { // 找到元素最少的列 32 if (min_col_count > p->count) { 33 min_col_count = p->count; 34 min_col_index = p->col; 35 } 36 } 37 for (Cross a = cols[min_col_index]->down; a != NULL; a = a->down) { // 对于该列的所有元素 38 result[result_pos++] = a->row; // 将该元素的行号压如result栈中 39 int new_stack_top = stack_top; 40 for (Cross b = rows[a->row]->right; b != NULL; b = b->right) { // 对于该元素所在行的所有元素 41 for (Cross c = cols[b->col]->down; c != NULL; c = c->down) { // 对于该元素所在列的所有元素 42 A->delrow(c->row); // 删除该元素所在行 43 stack[new_stack_top++] = c->row; // 并记录下删除操作(删除行时压入正的行号,删除列时压入负的列号) 44 } 45 A->delcol(b->col); // 删除该元素所在列 46 stack[new_stack_top++] = -b->col; // 记录删除操作 47 } 48 if (_find_many(new_stack_top, result_pos, N, n, fout)) // 调用下一级函数,如果找到的数独数量达到了需求,则返回真 49 return true; 50 for (int i = new_stack_top - 1; i >= stack_top; i--) { // 否则将压入stack栈中的删除操作弹出,并且做其逆操作 51 if (stack[i] > 0) 52 A->recoverrow(stack[i]); 53 else 54 A->recovercol(-stack[i]); 55 } 56 result_pos--; // 最后将部分解从result数组中弹出,进入下一轮循环 57 } 58 return false;在循环完所有情况还没有达到数量需求的情况下,返回假 59 }
以下是_find_one的代码:
1 bool dlx::_find_one(int stack_top, int result_pos) 2 { 3 if (!head->right) { // 找到一个解时记录答案并且直接返回真 4 for (int i = stack_top - 1; i >= 0; i--) { 5 if (stack[i] > 0) 6 A->recoverrow(stack[i]); 7 else 8 A->recovercol(-stack[i]); 9 } 10 return true; 11 } 12 int min_col_count = 100; 13 int min_col_index = -1; 14 for (Cross p = head->right; p != NULL; p = p->right) { 15 if (min_col_count > p->count) { 16 min_col_count = p->count; 17 min_col_index = p->col; 18 } 19 } 26 for (Cross a = cols[min_col_index]->down; a != NULL; a = a->down) { 27 result[result_pos++] = a->row; 28 int new_stack_top = stack_top; 29 for (Cross b = rows[a->row]->right; b != NULL; b = b->right) { 30 for (Cross c = cols[b->col]->down; c != NULL; c = c->down) { 31 A->delrow(c->row); 32 stack[new_stack_top++] = c->row; 33 } 34 A->delcol(b->col); 35 stack[new_stack_top++] = -b->col; 36 } 37 if (_find_one(new_stack_top, result_pos)) 38 return true; 39 for (int i = new_stack_top - 1; i >= stack_top; i--) { 40 if (stack[i] > 0) 41 A->recoverrow(stack[i]); 42 else 43 A->recovercol(-stack[i]); 44 } 45 result_pos--; 46 } 47 return false; 48 }
程序实现完成。_find_one和_find_many存在代码冗余的现象,下一步修改目标是将冗余部分抽取出来或是将两个函数合并为一个。
转载于:https://www.cnblogs.com/liusd/p/7581921.html
软件工程第一次作业——数独的求解与生成相关推荐
- 软件工程大作业——数独游戏
软件工程大作业--数独游戏1 一.PSP表格 二.问题分析 三.系统设计 四.具体实现 五.单元测试 六.程序性能及质量分析 七.GUI 八.总结 代码地址:https://github.com/fr ...
- 软件工程第一次作业(补充)
软件工程第一次作业补充 花20分钟写一个能自动生成小学四则运算题目的"软件",要求:除了整数以外,还要支持真分数的四则运算.将代码上传至coding.net, 并将地址发布至自己的 ...
- 软件工程第一次作业补充
软件工程第一次作业(2) 关于<构建之法>的5个问题 1)P28,2.1.3回归测试具体怎么操作? 2)P46讲到了软件工程师的成长,那么对于我们大学生来说,需要培养哪方面的品质? 3)P ...
- 广工软件工程第一次作业
这个作业属于哪个课程 广工软件工程学习 这个作业要求在哪里 软件工程第一次学习 这个作业的目标 评估当前的自己.展望未来中的所有问题和要求 其他参考文献 无 文章目录 git 自我评估 1. 个人简介 ...
- BUAA 软件工程 第一次作业
BUAA 软件工程2022 第一次作业 项目 内容 这个作业属于哪个课程 北航 2022 春季敏捷软件工程 这个作业的要求在哪里 作业说明链接 我在这个课程的目标是 了解并提高自己对软件工程的认识和实 ...
- 软件工程第一次作业:写一篇自己的博客
这个作业属于哪个课程 18级软件工程基础 这个作业要求在哪里 第一次个人作业:阅读与准备 我在这个课程的目标是 学会创建自己的博客以及Markdown的语法 其他参考文献 git优点缺点 其他参考文献 ...
- 2017年秋季学期软件工程第一次作业(曹洪茹)
作业一 在开始作业要求的正文之前,我先简单谈谈自此课开课以来,包括读了许多大牛写的博文之后的几点感悟和思考. 首先,作为一名有四年地方大学生活经验的军校研究生,我很激动也很庆幸在研究生阶段能遇到这么一 ...
- 软件工程——第一次作业
Part one:自我介绍 首先进行一些自我介绍,我叫贾雅杰,河北廊坊人,本科就读于南昌大学计算机科学与技术专业,现有编程能力稍弱,希望通过这学期的课程有所提高.第一次使用博客,希望大家多多提出宝贵的 ...
- 李昂 软件工程第一次作业
第1次个人作业 2.分析软件 第一部分 结缘计算机 1.计算机是你喜欢的领域吗?是你擅长的领域吗? 这两个问题,不得不说我的答案都是否.我很难喜欢计算机编程,更遑论软件硬件原理等等,感觉有些晦涩了,无 ...
- 高级软件工程第一次作业--准备
1) 回想一下你对计算机/软件工程专业的畅想 考研之所以选择计算机专业,是因为本科就是这个专业.不去跨专业,是因为觉得换个专业考,比起那些科班出身的人,考上的机率会更小,也有一部分原因是因为比起计 ...
最新文章
- java中final class的一点思考
- 2011寒假-操作系统学习笔记
- c++服务器开发学习--02--MySQL,Redis,ASIO,iocp,TrinityCore代码结构,c++对象模型
- 正则表达式不包含某个字符串_JMeter必知必会系列(18) JMeter正则表达式提取器疑难分析...
- How does setModel and getModel work in Fiori
- cannot find -lunwind-x86_64
- js横向滚动_seleniumJS处理滚动条
- c 调用上层类中函数_Matlab系列之函数嵌套
- quick: setup_mac.sh分析
- 珠海格力工厂一线员工待遇如何?
- pandas作图_pandas绘图
- 小米集团王嵋因错误表达致歉并请辞;亚马逊云服务出现中断,许多网站受到影响;deepin 深度系统更新发布|极客头条...
- 探秘综合布线产品质保问题
- string类型输入一行字符串,带空格
- 谷歌与IE浏览器兼容问题
- 学计算机专业独立显卡有必要吗,独立显卡驱动有什么用(显卡驱动有必要安装吗)...
- react 脚手架创建后暴漏配置文件 运行yarn eject 报错 (已解决)
- ChatGPT封杀潮,禁入学校,AI顶会特意改规则,LeCun:要不咱把小模型也禁了?...
- 详解Linux内核IO技术栈
- 1341:【例题】一笔画问题——欧拉(回)路