数独问题数学建模.docx
- 文档编号:3155978
- 上传时间:2022-11-18
- 格式:DOCX
- 页数:17
- 大小:149.32KB
数独问题数学建模.docx
《数独问题数学建模.docx》由会员分享,可在线阅读,更多相关《数独问题数学建模.docx(17页珍藏版)》请在冰豆网上搜索。
数独问题数学建模
数独问题
摘要
本文是对数独问题进行求解。
结合数独生成的特点,立足于题中数独建模和求解的要求,建立了数独难度分析
函数和整数规划模型。
对于问题一,首先研究数独难度的影响因素,通过综合分析数独的特点结构,得出
可以在常数时间内计算出来以衡量数独的难易程度。
通过计算可知
,根据数独难度的划分得到如下结论:
数独难度系数为4,达到了极难的程度。
对于问题二,我们通过对此数独的分析和讨论,利用穷举法,通过matlab软件编程求解,最终得出答案,如表1所示。
对于问题三,我们利用回溯法思想,建立求解模型,具体算法一般采用如下步骤:
1).在此数独初盘选择一个空单元格;
2).取这个单元格中一个可能的候选数;
3).将这个候选数填入单元格中,迭代完成数独;
4).若这个候选数推导得到一个无效数独终盘,返回此单元格取其他候选数;
对于问题四采用整数规划模型,采用三维0-1变量的方法,运用lingo软件编程求解。
最终得到答案,如表1所示。
表1数独问题的唯一解
8
1
2
7
5
3
6
4
9
9
4
3
6
8
2
1
7
5
6
7
5
4
9
1
2
9
3
1
5
4
2
3
7
8
9
6
3
6
9
8
4
5
7
2
1
2
8
7
1
6
9
5
3
4
5
2
1
9
7
4
3
6
8
4
3
8
5
2
6
9
1
7
7
9
6
3
1
8
4
5
2
关键词:
数独数独难度分析穷举法回溯法整体规划
1问题的重述
前段时间芬兰一位数学家号称设计出全球最难的“数独游戏”,并刊登在报纸上,让大家去挑战。
该数独如下图所示:
数独是根据9×9盘面上的已知数字,推理出所有剩余空格的数字,并满足每一行、每一列、每一个粗线宫内的数字均含1-9,且不重复。
每一道合格的数独谜题都有且仅有唯一答案,推理方法也以此为基础,任何无解或多解的题目都是不合格的。
根据以上描述,试完成以下问题:
1.分析此数独的难度;
2.用穷举算法求解数独;
3.设计此数独求解的较优的算法;
4.建立数独求解模型并给出此数独的答案。
2模型的基本假设
1该数独问题存在唯一解。
3符号说明
表示空单元格候选数
表示候选数数
的加权函数
表示数独空单元格中的候选数数目函数
表示该数独的空格处
表示该数独难度的函数
表示数k是否填入数独方中的(i,j)处
表示往空格处填入0后数独方中(i,j)处的数
表示经过求解后数独方中(i,j)处的数
4模型的建立与求解
4.1问题1
4.1.1数独难度的影响因素
通过对数独的分析与研究,数独难度与数独候选数、逻辑推理方法、搜索步数、空格数以及空格的分布情况都有密切的关系。
通过大量的计算观察发现,用到的逻辑与推理方法越复杂,那么在数独中出现的候选数越多;反之,在数独题中出现的候选数越多,解决数独题所用到的逻辑推理方法一般也越难。
解答一个数独所用到的搜索步数越多,数独中的候选数越多。
反之,一般情况下也成立。
另外数独中的空格数以及空格的分布情况与候选数也有同样类似的关系。
综合这几个影响数独难度的几种因素,分析候选数和空格数为主要影响因素,再根据其构造加权规范函数
,计算数值来衡量数独难度。
4.1.2
函数的建立
加权规范函数建立在候选数列表的基础上。
根据候选数列表,计算出每一个空单元格中的候选数数目,将候选数数目与其相对应的加权函数结合起来,计算加权规范函数
。
定义单元格
中的候选数数目函数
为
,这个函数仅适用
于数独P中的空单元格,而数独
中的空单元格可以表示为
(1)
有
个候选数的候选数数目函数
(2)
我们赋予它相应的加权函数
,从而得到加权函数
(3)
不能准确的反映数独难度,
受数独中空单元格数目影响很大,呈正向关系,如在数独中删除单元格,数独空单元格数增加,导致
增加,即空格数越多,
越大,然而这并不符合所有的数独。
为了排除这一影响,将加权函数
规范,得到加权规范函数:
(4)
根据以上的分析,对于某单元格
,其候选数数
越大,其对应的加权函数
也应越大。
我们采用指数函数
”计算数独
的
,其中
为某空单元格的候选数数目。
计算发现,
与数独难度是正相关的,即
越大,数独的难度越大。
4.1.3
函数的求解
根据题中给出的数独,按照数独游戏应该满足的条件,可以得到该数独的空格处的候选数列表,如下表2所示:
表2空格处候选数列表
8
1246
24569
2347
12357
1234
13569
4579
12459
124
3
6
12578
1248
1589
45789
14579
1456
7
456
348
9
1348
2
458
13456
123469
5
2469
2389
2368
7
1689
2489
12469
12369
12368
269
2389
4
5
7
289
1269
24679
2468
24679
1
268
2689
5689
3
24569
23457
234
1
23479
237
2349
359
6
8
23467
2346
8
5
2367
23469
39
1
2379
23567
9
2567
2378
123678
12368
4
257
2357
根据表2,把所得变量带入式(3)得:
(5)
由
代入式(4)式中得:
(6)
4.1.4数独难度的划分
根据计算所得
大小,我们将数独题难度分为四个区间,分别表示简单、中等、难、极难。
为方便表示,我们用1、2、3、4来表示难度系数。
1)若
,数独简单,有较多候选数的空单元格很少,此时数独题用一些简单的直观法就可以解决,用1表示。
2)若
,数独有一定难度,要解决此数独要用到候选数法中的一些简单方法,且与直观法结合起来推理,用2表示。
3)若
,数独比较难,内部逻辑结构复杂,将直观法与候选数法结合起来一般可以解决问题,用3表示。
4)若
,这个数独很难,内部的逻辑结构相当复杂,将直观法与候选数法结合起来不一定可以解决问题,甚至有时候需要对某些空单元格进行猜测,用4表示。
在这里,我们将(0,1)粗略分为四个区间,用来相对表示数独的相对难度。
根据数独难度的划分,由式(6)可得此数独难度系数为4,达到了极难的程度。
4.2问题2
4.2.1算法的介绍
本问中需要的是用穷举法对数独问题进行求解,首先介绍一下穷举法:
穷举法,或称为暴力破解法,是基于计算机特点而进行解题的思维方法。
一般是在一时找不出解决问题的更好途径(即从数学上找不到求解的公式或规则)时,可以根据问题中的部分条件(约束条件)将所有可能解的情况列举出来,然后通过逐个验证是否符合整个问题的求解要求,而得到问题的解。
这样解决问题的方法我们称之为穷举算法。
穷举算法特点是算法简单,但运行时所花费的时间量大。
因此,我们在用穷举方法解决问题时,应尽可能将明显的不符合条件的情况排除在外,以尽快取得问题的解。
4.2.2求解的思想
结合本问中需要用穷举法解决数独问题,最终算出上面所给出的数独问题的解。
针对此数独问题,在此先介绍一下我的算法思想:
1)建立一个堆栈来存放数据;
2)根据每行、每列和一个小九宫中不能出现相同的数字的规则来找出所有空格中的所有可能值;
3)从可能值中选取一个可能项最少的并提取一个出来,若还有可能值就将其放入堆栈中去,若提出的值不满足条件则从堆栈中再提取一个值来继续求解直到找到满足条件的解;
举个例子吧,对这一数独问题,可以很快找到第八行第七列的可能值为3和9,其它空格的可能值都超过了三个,现取出3出来进行尝试,那么放入堆栈中的是9和其它的可能值,还有a(数独值),然后一直按这种方法进行下去,要是遇到不满足则从堆栈中重新拿出一个值来,直到结果满足结束循环。
下面列了一个流程图,如下流程图所示。
首先进行对程序中的所用符号进行说明,先将数独中问题的初始值(空格为0)存入数组a,将所有空格中的可能值存入数组y。
4.2.3问题的求解
根据该流程图进行编程,并在matlab中实现,具体程序见附录1。
经过2分钟左右求解得到最终结果,如下表4所示:
否是
是否
否
是
利用穷举法编程求解的流程图
4.3问题3
根据问题2中处理方法,发现穷举算法的特点是算法简单,但运行时所花费的时间较长。
因此,我们在此基础上进行改进,尽可能将明显的不符合条件的情况排除在外,以尽快取得问题的解。
在求解数独的过程中,遍历此数独所有可能的搜索树,直至找到数独的解为止!
在这个过程中,我们采用回溯法进行求解。
回溯法是一种搜索算法,其基本思路是:
在一个问题中,根据题意给出的边界条件划定出所有可能解的范围(称为可能解),根据题意确定出约束条件。
利用程序顺次在所有可能解中,搜索时按照深度搜索的方式进行。
即在第一层选定一个满足约束条件的解,然后以该可能解为出发点,搜索第二层的一个可能解(试探)。
如果搜索到第二层的一个可能解,则继续搜索第三层的一个可能解。
依次类推,直到所有层的可能解都被找到,则得到了该问题的一个完整解。
如果第二层所有的可能解都不满足约束条件,则返回第一层,放弃原有的可能解,使用第一层的下一个可能解(回溯)。
以此类推,寻找第二层的一个可能解。
具体算法一般采用如下步骤:
1)在此数独初盘选择一个空单元格;
2)取这个单元格中一个可能的候选数;
3)将这个候选数填入单元格中,迭代完成数独;
4)若这个候选数推导得到一个无效数独终盘,返回此单元格取其他候选数;
由于回溯法是在不断地试探和回溯中运算,因此也可以称为试探法或者试探—回溯法。
从上面的描述中可知,回溯法得到的问题的解只是根据不同的初始条件获得的第一个完全满足所有约束条件的解,因此该解的获得和初始条件有关。
如果想要获得该问题的全部解,则需要遍历所有的可能的初始条件,也就是遍历所有的第一层的可能解。
回溯法相对于其他穷举的特点在于,不必把问题的每一层的所有的可能解都遍历一遍,只要当前的可能解不满足约束条件就抛弃该解,寻求下一个可能解,而不必求解其余的下层解。
当当前层的所有可能解都不满足约束条件,」则回溯到上一层,抛弃上一层的当前可能解。
从以上分析中结合数独问题的规则,得出数独问题的约束条件为:
l)每一格的数值范围仅限于l一9。
2)每一格内的数字在当前行不允许重复。
3)每一格内的数字在当前列不允许重复。
4)每一格内的数字在当前小九宫内不允许重复
此算法最坏运行时间是
,其中。
n为数独中数字未确定的单元格数,k是数字已经确定的单元格数。
但是将这个方法稍作修改就可以在线性时间解决一个数独题目,即首先选择数独中候选数个数最少的空单元格,再在其中选择候选数进行试验一检验。
虽然回溯法有超多项式的最坏时间,但是这种算法可以判断任何数独是否具有唯一解,进而解决这个数独,难度很大的数独求解经常用到此法。
4.4问题4
4.4.1问题的分析
根据对数独问题的分析,可知数独是根据9×9盘面上的已知数字,推理出所有剩余空格的数字,并满足每一行、每一列、每一个粗线宫内的数字均含1-9,并且不能重复。
为了叙述方便,对一个空的9阶数独方的行、列、区进行重新编码
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 问题 数学 建模
![提示](https://static.bdocx.com/images/bang_tan.gif)