数独游戏——数模处理(MATLAB整数规划)

一、问题描述

  数独是源自18世纪瑞士的一种数学游戏，盘面是个九宫，每一宫又分为九个小格。在这八十一格中根据已知数字和解题条件，利用逻辑和推理，在其他的空格上填入1-9的数字，使1-9中的每个数字在每一行、每一列和每一宫中都只出现一次。

二、模型建立

  将 $9\times 9$ 的数独矩阵转化为 $9\times 9\times 9$ 的三维矩阵 ，三维矩阵中每一个元素的值为0或1。例如，若原数独矩阵的第一行第二列的数字为7，则三维矩阵第7层的对应位置为1，其它层该位置皆为0，总共个 $9\times 9\times 9$ 变量，用 $X(i,j,k)$ 表示第 $k$ 层矩阵的第 $i$ 行第 $j$ 列的元素。

  约束条件：

1. 9层矩阵的同一位置只可以出现一个1：
   $$\sum _{k=1}^{9}X(i,j,k)=1,i=1,2,\cdots ,9,j=1,2,\cdots ,9$$

2. 任何一层的任何一行只能有一个1：
   $$\sum _{j=1}^{9}X(i,j,k)=1,i=1,2,\cdots ,9,k=1,2,\cdots ,9$$

3. 任何一层的任何一列只能有一个1：
   $$\sum _{i=1}^{9}X(i,j,k)=1,j$$