数独(sudoku)的生成与破解

数独(sudoku)的生成与破解
最近在网上比较流行的智力游戏。笔者本人也玩过，可以下个模拟游戏试试，简单的还可以，太难就无从下手了。虽然偶脑子不好使，但偶是计算机科班出身，怕你不成，老规矩，编程破解。

首先，偶是第一次做数独的程序，可能程序不强，以后有时间再改进改进。望高手评析。

还是把数独游戏的规则说一说吧，或许你是刚刚听说这个名字的朋友。数独(sudoku)，起源于瑞士，于1970 年代由美国的一家数学逻辑游戏杂志首先发表，当时名为Number Place。及后在日本大力推广下得以发扬光大，于1984 年取名“数独”，即“独立的数字”的省略，在一个9x9的方格中，有81个小方格组成，然后又分9个大块，每块由3x3的方格组成，就是中国的九宫图，大九宫里面再套9个小九宫，共九九八十一个小格子，游戏开始前会有一些格子上写好了数，你需要在剩下的格子里填数，真到把所有格子填满，并且要求，任何一行或一列或者一个小九宫中没有相同的数字，当然你只能用1-9之间的9个数字。如下图就是一个数独。希望我解释清楚了，如果你还不清楚，用google搜索一下相关资料，点这里http://www.google.com/search?q=sudoku%20%E6%95%B0%E7%8B%AC&hl=zh-CN&lr=&nxpt=10.1471893728569764719035

好啦，言归正传。简单来说，我的破解法非常简单，就是超级无敌强硬大搜索，呵呵，别拿你想数独的那套来让计算机想，它比你可笨得多了，它只会照程序做事，做得快而已，快就是它的本事，其它的一概不会。

核心算法：深度优先搜索（其它形式的搜索也可以）

数据结构：如果用递归的形式写深搜，定义在函数dfs里的所有变量都可以看成是这里的数据结构，因为它们自动地被系统压入栈内，所以，省了，你唯一要做的就是一个二维数组，存放当前数独的状态。

当然有了这些，偶还不敢动手做，如果你做过马遍历的程序，大概会有点怕，那才8x8，这里是9x9，不来点‘启发式’谁敢动手写程序，有可能一个数独来几千几万个解，一个解要搜80层上下（估计），懂得树这个数据结构的人就会明白，80层是什么概念，1-9有9个数字，就是9叉树，至少是9^80量级的代价，什么？计算机反正算得快？也不行，再快的计算机遇到指数复杂度的程序也得变回傻子！谢天谢地，棋局尺寸是固定不变的，我们需要做的就是，剪枝。

偶的启发式思想来源于偶想数独的思路，数独之所以难是因为可行情况太多，把这种不确定性降低就会使它变得简单。某个格子上可以填的数字的个数就称为它的不确定度吧，首先，如果一开始空格比较少，那空格上的不确定度就小，数独也就相对容易，同时，随着填上去数字增多，剩余空格上的不确定度也会降低，如果某个格子上的不确定度降到1，那这个格子可以先确定下来，如果降到了0，哦，非常遗憾，在前面的填数中一定是填错了，剪枝也发生在这里，你不得不回退。

详细地说，如果把每个空格做为分枝，那要优先选择哪个分枝进行搜索呢？对每个空格计算它的权值，也就是它的不确定度，然后从中选出最小的一个进行搜索，同时放弃其它空格在这一层上的搜索！也就是说对剩余的空格交给子层处理。这样在某一个结点上的处理就包括两步：1，选择最佳空格，2，遍历这个空格的所有可行值，填入空格，递归。

OK，看程序：

#ifndef SUDOKU_RICK_0701_
#define SUDOKU_RICK_0701_
class CSudoku
{
 int map[9][9];
 int solves;
 int check(int,int,int*);
 void dfs();
public:
 CSudoku(int n=40);// 随机生成数独,n越大越难
 CSudoku(int *data);// 人工指定数独
 virtual ~CSudoku();
 void display();// 输出数独
 void resolve();// 解数独
};
#endif

#include "sudoku.h"
#include "stdio.h"
#include "stdlib.h"
#include "time.h"

CSudoku::CSudoku(int n)
{
 int i,j,k,m,mark[10],temp,blanks=0;
 srand(time(0));
 for(i=0;i<9;++i)
 {
  for(j=0;j<9;++j)
   map[i][j]=j+1;
  for(j=0;j<9;++j)
  {
   int a=rand()%9;
   int b=rand()%9;
   temp=m