数独问题
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
（9）式表示每一列只能有一个数k；
（10）式表示每一个区内只能有一个数k；
（11）式表示 只能取0或1；
（12）上所有式中
4.4.3模型的求解
条件：
（1）每一（i，j）处只能填一个数 k；
（2）第i 行只有一个 k；
（3）第j 列只有一个 k；
（4）第一个区只有一个k；
4.4.2.3该9阶数独方的数学模型
由上述分析可建立该9阶数独方的整数规划模型如下：
注：
（7）式表示数独方中每一（i，j）处只能填入一个数k；
（8）式表示每一行只能有一个数k；
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来表示难度系数。
在这里，我们将(0，1)粗略分为四个区间，用来相对表示数独的相对难度。根据数独难度的划分，由式(6)可得此数独难度系数为4，达到了极难的程度。
4.2问题2
4.2.1算法的介绍
本问中需要的是用穷举法对数独问题进行求解，首先介绍一下穷举法：穷举法，或称为暴力破解法，是基于计算机特点而进行解题的思维方法。一般是在一时找不出解决问题的更好途径（即从数学上找不到求解的公式或规则）时，可以根据问题中的部分条件（约束条件）将所有可能解的情况列举出来，然后通过逐个验证是否符合整个问题的求解要求，而得到问题的解。这样解决问题的方法我们称之为穷举算法。穷举算法特点是算法简单，但运行时所花费的时间量大。因此，我们在用穷举方法解决问题时，应尽可能将明显的不符合条件的情况排除在外，以尽快取得问题的解。
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
举个例子吧，对这一数独问题，可以很快找到第八行第七列的可能值为3和9，其它空格的可能值都超过了三个，现取出3出来进行尝试，那么放入堆栈中的是9和其它的可能值，还有a（数独值），然后一直按这种方法进行下去，要是遇到不满足则从堆栈中重新拿出一个值来，直到结果满足结束循环。下面列了一个流程图，如下流程图所示。首先进行对程序中的所用符号进行说明，先将数独中问题的初始值(空格为0)存入数组a，将所有空格中的可能值存入数组y。
4.1.2 函数的建立
加权规范函数建立在候选数列表的基础上。根据候选数列表，计算出每一个空单元格中的候选数数目，将候选数数目与其相对应的加权函数结合起来，计算加权规范函数 。
定义单元格 中的候选数数目函数 为 ，这个函数仅适用
于数独P中的空单元格，而数独 中的空单元格可以表示为
(1)
有 个候选数的候选数数目函数
回溯法相对于其他穷举的特点在于，不必把问题的每一层的所有的可能解都遍历一遍，只要当前的可能解不满足约束条件就抛弃该解，寻求下一个可能解，而不必求解其余的下层解。当当前层的所有可能解都不满足约束条件，」则回溯到上一层，抛弃上一层的当前可能解。从以上分析中结合数独问题的规则，得出数独问题的约束条件为:
l)每一格的数值范围仅限于l一9。
数独问题
摘要
本文是对数独问题进行求解。结合数独生成的特点，立足于题中数独建模和求解的要求，建立了数独难度分析 函数和整数规划模型。
对于问题一，首先研究数独难度的影响因素，通过综合分析数独的特点结构，得出 可以在常数时间内计算出来以衡量数独的难易程度。通过计算可知 ，根据数独难度的划分得到如下结论：数独难度系数为4，达到了极难的程度。
表示该数独的空格处
表示该数独难度的函数
表示数k是否填入数独方中的（i，j）处
表示往空格处填入0后数独方中（i，j）处的数
表示经过求解后数独方中（i，j）处的数
4模型的建立与求解
4.1问题1
4.1.1数独难度的影响因素
通过对数独的分析与研究，数独难度与数独候选数、逻辑推理方法、搜索步数、空格数以及空格的分布情况都有密切的关系。通过大量的计算观察发现，用到的逻辑与推理方法越复杂，那么在数独中出现的候选数越多;反之，在数独题中出现的候选数越多，解决数独题所用到的逻辑推理方法一般也越难。解答一个数独所用到的搜索步数越多，数独中的候选数越多。反之，一般情况下也成立。另外数独中的空格数以及空格的分布情况与候选数也有同样类似的关系。综合这几个影响数独难度的几种因素，分析候选数和空格数为主要影响因素，再根据其构造加权规范函数 ，计算数值来衡量数独难度。
1)若 ，数独简单，有较多候选数的空单元格很少，此时数独题用一些简单的直观法就可以解决，用1表示。
2)若 ，数独有一定难度，要解决此数独要用到候选数法中的一些简单方法，且与直观法结合起来推理，用2表示。
3)若 ，数独比较难，内部逻辑结构复杂，将直观法与候选数法结合起来一般可以解决问题，用3表示。
4)若 ，这个数独很难，内部的逻辑结构相当复杂，将直观法与候选数法结合起来不一定可以解决问题，甚至有时候需要对某些空单元格进行猜测，用4表示。
1）在此数独初盘选择一个空单元格；
2）取这个单元格中一个可能的候选数；
3）将这个候选数填入单元格中，迭代完成数独；
4）若这个候选数推导得到一个无效数独终盘，返回此单元格取其他候选数；
由于回溯法是在不断地试探和回溯中运算，因此也可以称为试探法或者试探—回溯法。从上面的描述中可知，回溯法得到的问题的解只是根据不同的初始条件获得的第一个完全满足所有约束条件的解，因此该解的获得和初始条件有关。如果想要获得该问题的全部解，则需要遍历所有的可能的初始条件，也就是遍历所有的第一层的可能解。
: 表示数k不填入（i，j）处
4.4.2.2 目标函数与约束
规划问题的目标函数一般都是最大或最小型的，又因根据题意可知该 9阶数独方的可行解唯一，且该可行解就是最优解，所以确定目标函数为单位数1。
通过对数独这个问题的分析可以得出，所有剩余空格的数字，并满足每一行、每一列、每一个粗线宫内的数字均含1-9，并且不能重复。因此可以得到如下约束
4.4.2模型的建立
由于数 要么填入该9阶数独方矩阵的（i，j）处，要么不填入，只有这两种状态，因而可用三维0和1 变量来表示。又该9 阶数独方的解唯一，可将这唯一解作为目标函数，行、列、区的约束条件可用三维的0和1 变量表出。
4.4.2.1 决策变量
令 : 表示数k填入 n 阶数独方矩阵的（i，j）处
(2)
我们赋予它相应的加权函数 ，从而得到加权函数
(3)
不能准确的反映数独难度， 受数独中空单元格数目影响很大，呈正向关系，如在数独中删除单元格，数独空单元格数增加，导致 增加，即空格数越多， 越大，然而这并不符合所有的数独。为了排除这一影响，将加权函数 规范，得到加权规范函数：
(4)
根据以上的分析，对于某单元格 ，其候选数数 越大，其对应的加权函数 也应越大。我们采用指数函数 ”计算数独 的 ，其中 为某空单元格的候选数数目。计算发现， 与数独难度是正相关的，即 越大，数独的难度越大。
2)每一格内的数字在当前行不允许重复。
3)每一格内的数字在当前列不允许重复。
4)每一格内的数字在当前小九宫内不允许重复
此算法最坏运行时间是 ，其中。n为数独中数字未确定的单元格数，k是数字已经确定的单元格数。但是将这个方法稍作修改就可以在线性时间解决一个数独题目，即首先选择数独中候选数个数最少的空单元格，再在其中选择候选数进行试验一检验。虽然回溯法有超多项式的最坏时间，但是这种算法可以判断任何数独是否具有唯一解，进而解决这个数独，难度很大的数独求解经常用到此法。
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，且不重复。 每一道合格的数独谜题都有且仅有唯一答案，推理方法也以此为基础，任何无解或多解的题目都是不合格的。
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
回溯法是一种搜索算法，其基本思路是:在一个问题中，根据题意给出的边界条件划定出所有可能解的范围(称为可能解)，根据题意确定出约束条件。利用程序顺次在所有可能解中，搜索时按照深度搜索的方式进行。即在第一层选定一个满足约束条件的解，然后以该可能解为出发点，搜索第二层的一个可能解(试探)。如果搜索到第二层的一个可能解，则继续搜索第三层的一个可能解。依次类推，直到所有层的可能解都被找到，则得到了该问题的一个完整解。如果第二层所有的可能解都不满足约束条件，则返回第一层，放弃原有的可能解，使用第一层的下一个可能解(回溯)。以此类推，寻找第二层的一个可能解。具体算法一般采用如下步骤：
4.2.3 问题的求解
根据该流程图进行编程，并在matlab中实现，具体程序见附录1。经过2分钟左右求解得到最终结果，如下表4所示：
否 是
是 否
否Leabharlann Baidu
是
利用穷举法编程求解的流程图
4.3问题3
根据问题2中处理方法，发现穷举算法的特点是算法简单，但运行时所花费的时间较长。因此，我们在此基础上进行改进，尽可能将明显的不符合条件的情况排除在外，以尽快取得问题的解。在求解数独的过程中，遍历此数独所有可能的搜索树，直至找到数独的解为止!在这个过程中，我们采用回溯法进行求解。
4.1.3 函数的求解
根据题中给出的数独，按照数独游戏应该满足的条件，可以得到该数独的空格处的候选数列表，如下表2所示：
表2空格处候选数列表
8
1246
24569
2347
12357
1234
13569
4579
1345679
12459
124
3
6
12578
1248
1589
45789
14579
1456
7
456
4.4问题4
4.4.1问题的分析
根据对数独问题的分析，可知数独是根据9×9盘面上的已知数字，推理出所有剩余空格的数字，并满足每一行、每一列、每一个粗线宫内的数字均含1-9，并且不能重复。
为了叙述方便，对一个空的 9阶数独方的行、列、区进行重新编码。数独方的行序为从上到下记作1～9， 列序为从左到右记作1～9；然后把每一个粗线宫作为一个区，区的行序为从上到下记作1～3， 列序为从左到右记作1～3；区行的行序为从上到下记作1～3， 区列的列序为从左到右记作1～3。具体的编码方式如下图1所示：
对于问题二，我们通过对此数独的分析和讨论，利用穷举法，通过matlab软件编程求解，最终得出答案，如表1所示。
对于问题三，我们利用回溯法思想，建立求解模型，具体算法一般采用如下步骤：
1).在此数独初盘选择一个空单元格；
2).取这个单元格中一个可能的候选数；
3).将这个候选数填入单元格中，迭代完成数独；
根据以上描述，试完成以下问题：
1. 分析此数独的难度；
2. 用穷举算法求解数独；
3. 设计此数独求解的较优的算法；
4.建立数独求解模型并给出此数独的答案。
2模型的基本假设
1该数独问题存在唯一解。
2该数独难度等级的划分具有一定的科学性。
3符号说明
表示空单元格候选数
表示候选数数 的加权函数
表示数独空单元格中的候选数数目函数
4.2.2求解的思想
结合本问中需要用穷举法解决数独问题，最终算出上面所给出的数独问题的解。针对此数独问题，在此先介绍一下我的算法思想：
1）建立一个堆栈来存放数据；
2）根据每行、每列和一个小九宫中不能出现相同的数字的规则来找出所有空格中的所有可能值；
3）从可能值中选取一个可能项最少的并提取一个出来，若还有可能值就将其放入堆栈中去，若提出的值不满足条件则从堆栈中再提取一个值来继续求解直到找到满足条件的解；
（10）式表示每一个区内只能有一个数k；
（11）式表示 只能取0或1；
（12）上所有式中
4.4.3模型的求解
条件：
（1）每一（i，j）处只能填一个数 k；
（2）第i 行只有一个 k；
（3）第j 列只有一个 k；
（4）第一个区只有一个k；
4.4.2.3该9阶数独方的数学模型
由上述分析可建立该9阶数独方的整数规划模型如下：
注：
（7）式表示数独方中每一（i，j）处只能填入一个数k；
（8）式表示每一行只能有一个数k；
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来表示难度系数。
在这里，我们将(0，1)粗略分为四个区间，用来相对表示数独的相对难度。根据数独难度的划分，由式(6)可得此数独难度系数为4，达到了极难的程度。
4.2问题2
4.2.1算法的介绍
本问中需要的是用穷举法对数独问题进行求解，首先介绍一下穷举法：穷举法，或称为暴力破解法，是基于计算机特点而进行解题的思维方法。一般是在一时找不出解决问题的更好途径（即从数学上找不到求解的公式或规则）时，可以根据问题中的部分条件（约束条件）将所有可能解的情况列举出来，然后通过逐个验证是否符合整个问题的求解要求，而得到问题的解。这样解决问题的方法我们称之为穷举算法。穷举算法特点是算法简单，但运行时所花费的时间量大。因此，我们在用穷举方法解决问题时，应尽可能将明显的不符合条件的情况排除在外，以尽快取得问题的解。
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
举个例子吧，对这一数独问题，可以很快找到第八行第七列的可能值为3和9，其它空格的可能值都超过了三个，现取出3出来进行尝试，那么放入堆栈中的是9和其它的可能值，还有a（数独值），然后一直按这种方法进行下去，要是遇到不满足则从堆栈中重新拿出一个值来，直到结果满足结束循环。下面列了一个流程图，如下流程图所示。首先进行对程序中的所用符号进行说明，先将数独中问题的初始值(空格为0)存入数组a，将所有空格中的可能值存入数组y。
4.1.2 函数的建立
加权规范函数建立在候选数列表的基础上。根据候选数列表，计算出每一个空单元格中的候选数数目，将候选数数目与其相对应的加权函数结合起来，计算加权规范函数 。
定义单元格 中的候选数数目函数 为 ，这个函数仅适用
于数独P中的空单元格，而数独 中的空单元格可以表示为
(1)
有 个候选数的候选数数目函数
回溯法相对于其他穷举的特点在于，不必把问题的每一层的所有的可能解都遍历一遍，只要当前的可能解不满足约束条件就抛弃该解，寻求下一个可能解，而不必求解其余的下层解。当当前层的所有可能解都不满足约束条件，」则回溯到上一层，抛弃上一层的当前可能解。从以上分析中结合数独问题的规则，得出数独问题的约束条件为:
l)每一格的数值范围仅限于l一9。
数独问题
摘要
本文是对数独问题进行求解。结合数独生成的特点，立足于题中数独建模和求解的要求，建立了数独难度分析 函数和整数规划模型。
对于问题一，首先研究数独难度的影响因素，通过综合分析数独的特点结构，得出 可以在常数时间内计算出来以衡量数独的难易程度。通过计算可知 ，根据数独难度的划分得到如下结论：数独难度系数为4，达到了极难的程度。
表示该数独的空格处
表示该数独难度的函数
表示数k是否填入数独方中的（i，j）处
表示往空格处填入0后数独方中（i，j）处的数
表示经过求解后数独方中（i，j）处的数
4模型的建立与求解
4.1问题1
4.1.1数独难度的影响因素
通过对数独的分析与研究，数独难度与数独候选数、逻辑推理方法、搜索步数、空格数以及空格的分布情况都有密切的关系。通过大量的计算观察发现，用到的逻辑与推理方法越复杂，那么在数独中出现的候选数越多;反之，在数独题中出现的候选数越多，解决数独题所用到的逻辑推理方法一般也越难。解答一个数独所用到的搜索步数越多，数独中的候选数越多。反之，一般情况下也成立。另外数独中的空格数以及空格的分布情况与候选数也有同样类似的关系。综合这几个影响数独难度的几种因素，分析候选数和空格数为主要影响因素，再根据其构造加权规范函数 ，计算数值来衡量数独难度。
1)若 ，数独简单，有较多候选数的空单元格很少，此时数独题用一些简单的直观法就可以解决，用1表示。
2)若 ，数独有一定难度，要解决此数独要用到候选数法中的一些简单方法，且与直观法结合起来推理，用2表示。
3)若 ，数独比较难，内部逻辑结构复杂，将直观法与候选数法结合起来一般可以解决问题，用3表示。
4)若 ，这个数独很难，内部的逻辑结构相当复杂，将直观法与候选数法结合起来不一定可以解决问题，甚至有时候需要对某些空单元格进行猜测，用4表示。
1）在此数独初盘选择一个空单元格；
2）取这个单元格中一个可能的候选数；
3）将这个候选数填入单元格中，迭代完成数独；
4）若这个候选数推导得到一个无效数独终盘，返回此单元格取其他候选数；
由于回溯法是在不断地试探和回溯中运算，因此也可以称为试探法或者试探—回溯法。从上面的描述中可知，回溯法得到的问题的解只是根据不同的初始条件获得的第一个完全满足所有约束条件的解，因此该解的获得和初始条件有关。如果想要获得该问题的全部解，则需要遍历所有的可能的初始条件，也就是遍历所有的第一层的可能解。
: 表示数k不填入（i，j）处
4.4.2.2 目标函数与约束
规划问题的目标函数一般都是最大或最小型的，又因根据题意可知该 9阶数独方的可行解唯一，且该可行解就是最优解，所以确定目标函数为单位数1。
通过对数独这个问题的分析可以得出，所有剩余空格的数字，并满足每一行、每一列、每一个粗线宫内的数字均含1-9，并且不能重复。因此可以得到如下约束
4.4.2模型的建立
由于数 要么填入该9阶数独方矩阵的（i，j）处，要么不填入，只有这两种状态，因而可用三维0和1 变量来表示。又该9 阶数独方的解唯一，可将这唯一解作为目标函数，行、列、区的约束条件可用三维的0和1 变量表出。
4.4.2.1 决策变量
令 : 表示数k填入 n 阶数独方矩阵的（i，j）处
(2)
我们赋予它相应的加权函数 ，从而得到加权函数
(3)
不能准确的反映数独难度， 受数独中空单元格数目影响很大，呈正向关系，如在数独中删除单元格，数独空单元格数增加，导致 增加，即空格数越多， 越大，然而这并不符合所有的数独。为了排除这一影响，将加权函数 规范，得到加权规范函数：
(4)
根据以上的分析，对于某单元格 ，其候选数数 越大，其对应的加权函数 也应越大。我们采用指数函数 ”计算数独 的 ，其中 为某空单元格的候选数数目。计算发现， 与数独难度是正相关的，即 越大，数独的难度越大。
2)每一格内的数字在当前行不允许重复。
3)每一格内的数字在当前列不允许重复。
4)每一格内的数字在当前小九宫内不允许重复
此算法最坏运行时间是 ，其中。n为数独中数字未确定的单元格数，k是数字已经确定的单元格数。但是将这个方法稍作修改就可以在线性时间解决一个数独题目，即首先选择数独中候选数个数最少的空单元格，再在其中选择候选数进行试验一检验。虽然回溯法有超多项式的最坏时间，但是这种算法可以判断任何数独是否具有唯一解，进而解决这个数独，难度很大的数独求解经常用到此法。
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，且不重复。 每一道合格的数独谜题都有且仅有唯一答案，推理方法也以此为基础，任何无解或多解的题目都是不合格的。
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
回溯法是一种搜索算法，其基本思路是:在一个问题中，根据题意给出的边界条件划定出所有可能解的范围(称为可能解)，根据题意确定出约束条件。利用程序顺次在所有可能解中，搜索时按照深度搜索的方式进行。即在第一层选定一个满足约束条件的解，然后以该可能解为出发点，搜索第二层的一个可能解(试探)。如果搜索到第二层的一个可能解，则继续搜索第三层的一个可能解。依次类推，直到所有层的可能解都被找到，则得到了该问题的一个完整解。如果第二层所有的可能解都不满足约束条件，则返回第一层，放弃原有的可能解，使用第一层的下一个可能解(回溯)。以此类推，寻找第二层的一个可能解。具体算法一般采用如下步骤：
4.2.3 问题的求解
根据该流程图进行编程，并在matlab中实现，具体程序见附录1。经过2分钟左右求解得到最终结果，如下表4所示：
否 是
是 否
否Leabharlann Baidu
是
利用穷举法编程求解的流程图
4.3问题3
根据问题2中处理方法，发现穷举算法的特点是算法简单，但运行时所花费的时间较长。因此，我们在此基础上进行改进，尽可能将明显的不符合条件的情况排除在外，以尽快取得问题的解。在求解数独的过程中，遍历此数独所有可能的搜索树，直至找到数独的解为止!在这个过程中，我们采用回溯法进行求解。
4.1.3 函数的求解
根据题中给出的数独，按照数独游戏应该满足的条件，可以得到该数独的空格处的候选数列表，如下表2所示：
表2空格处候选数列表
8
1246
24569
2347
12357
1234
13569
4579
1345679
12459
124
3
6
12578
1248
1589
45789
14579
1456
7
456
4.4问题4
4.4.1问题的分析
根据对数独问题的分析，可知数独是根据9×9盘面上的已知数字，推理出所有剩余空格的数字，并满足每一行、每一列、每一个粗线宫内的数字均含1-9，并且不能重复。
为了叙述方便，对一个空的 9阶数独方的行、列、区进行重新编码。数独方的行序为从上到下记作1～9， 列序为从左到右记作1～9；然后把每一个粗线宫作为一个区，区的行序为从上到下记作1～3， 列序为从左到右记作1～3；区行的行序为从上到下记作1～3， 区列的列序为从左到右记作1～3。具体的编码方式如下图1所示：
对于问题二，我们通过对此数独的分析和讨论，利用穷举法，通过matlab软件编程求解，最终得出答案，如表1所示。
对于问题三，我们利用回溯法思想，建立求解模型，具体算法一般采用如下步骤：
1).在此数独初盘选择一个空单元格；
2).取这个单元格中一个可能的候选数；
3).将这个候选数填入单元格中，迭代完成数独；
根据以上描述，试完成以下问题：
1. 分析此数独的难度；
2. 用穷举算法求解数独；
3. 设计此数独求解的较优的算法；
4.建立数独求解模型并给出此数独的答案。
2模型的基本假设
1该数独问题存在唯一解。
2该数独难度等级的划分具有一定的科学性。
3符号说明
表示空单元格候选数
表示候选数数 的加权函数
表示数独空单元格中的候选数数目函数
4.2.2求解的思想
结合本问中需要用穷举法解决数独问题，最终算出上面所给出的数独问题的解。针对此数独问题，在此先介绍一下我的算法思想：
1）建立一个堆栈来存放数据；
2）根据每行、每列和一个小九宫中不能出现相同的数字的规则来找出所有空格中的所有可能值；
3）从可能值中选取一个可能项最少的并提取一个出来，若还有可能值就将其放入堆栈中去，若提出的值不满足条件则从堆栈中再提取一个值来继续求解直到找到满足条件的解；