数独求解_matlab

subject.png

空的位置用0补,用空格分隔,将题目记录到txt文档中。

data.txt

0 0 0 0 0 0 0 0 0
0 2 3 0 0 0 7 8 0
1 0 0 4 0 6 0 0 9
4 0 0 0 5 0 0 0 1
9 0 0 0 0 0 0 0 6
0 6 0 0 0 0 0 9 0
0 0 5 0 0 0 8 0 0
0 0 0 3 0 1 0 0 0
0 0 0 0 9 0 0 0 0

检查某个数字能否填入某个空格,只检查那一横行。

check_row.m

function [result] = check_row(matrix, i, j, num)
    result = true;
    for col = 1 : 9
        if matrix(i, col) == num
            result = false;
            break;
        end
    end
end

检查某个数字能否填入某个空格,只检查那一列。

check_col.m

function [result] = check_col(matrix, i, j, num)
    result = true;
    for row = 1 : 9
        if matrix(row, j) == num
            result = false;
            break;
        end
    end
end

检查某个数字能否填入某个空格,只检查所在的3*3小方阵。

check_block.m

function [result] = check_block(matrix, i, j, num)
    result = true;
    
    % 计算此i,j所属小九宫格的坐标
    block_row = floor( (i-1) / 3 ) + 1;
    block_col = floor( (j-1) / 3 ) + 1;
    
    block_row_from = (block_row * 3) - 2;
    block_row_to = block_row * 3;
    block_col_from = (block_col * 3) - 2;
    block_col_to = block_col * 3;
    
    
    for row = block_row_from : block_row_to
        for col = block_col_from : block_col_to
            if matrix(row, col) == num
                result = false;
                break;
            end
        end
    end
    
end

检查某个数字能否填入某个空格,联合检查:所在行、所在列、所在的3*3小方阵。

check_all.m

function [result] = check_all(matrix, i, j, num)
    result = check_row(matrix, i, j, num) && check_col(matrix, i, j, num) &&  check_block(matrix, i, j, num);
end

绘制格子和数据的 function :

function [] = draw(matrix, isSubplot, rows, cols, rank)

if isSubplot == true
    subplot(rows, cols, rank);
    
matrix0 = dlmread('data.txt', '');
matrix0 = matrix0';
    
plot(0,0);
hold on
for i = 0 : 9
    if rem(i,3) == 0
        plot([0;9],[i,i],'linewidth',3,'color','black');
        plot([i;i],[0,9],'linewidth',3,'color','black');
    else
        plot([0;9],[i,i],'linewidth',1,'color','black');
        plot([i;i],[0,9],'linewidth',1,'color','black');
    end
end

matrix = matrix';
for i = 1 : 9
    for j = 1 : 9
        if matrix(i,j) ~= 0
            if matrix(i,j) == matrix0(i,j)
                text(i-0.5, 9.5-j, strcat('\color{red}',num2str(matrix(i,j))));
            else
                text(i-0.5, 9.5-j, num2str(matrix(i,j)));
            end
        end
    end
end

end

使用穷举解决方案。

solve.m

% 整个大九宫格的index编号(从左到右,从上到下):1 ~ 81
function solve(matrix, index)
    if index > 81       % 已经找到解
        disp(matrix);
        draw(matrix, true, 1, 2, 2);
    else
        % 计算此index对应i、j
        i = floor((index - 1) / 9) + 1;
        j = mod((index-1), 9) + 1;
        
        if matrix(i, j) ~= 0    % 如果是个已知数的话,直接跳到下一格
            solve(matrix, index+1);
        else                    % 未知数
            for num = 1 : 9                                % 试值
                if check_all(matrix, i, j, num) == true    % 可以
                    matrix(i, j) = num;                    % 赋值
                    solve(matrix, index+1);                % 跳到下一格
                end
            end
        end
    end
end

主程序:将文本读入到二维数组中,调用解决程序,从位置1开始。

sudoku.m

matrix = dlmread('data.txt', ' ');

disp(matrix);
draw(matrix, true, 1, 2, 1);

solve(matrix, 1);

运行结果:

>> sudoku
     0     0     0     0     0     0     0     0     0
     0     2     3     0     0     0     7     8     0
     1     0     0     4     0     6     0     0     9
     4     0     0     0     5     0     0     0     1
     9     0     0     0     0     0     0     0     6
     0     6     0     0     0     0     0     9     0
     0     0     5     0     0     0     8     0     0
     0     0     0     3     0     1     0     0     0
     0     0     0     0     9     0     0     0     0

     5     4     9     7     3     8     1     6     2
     6     2     3     9     1     5     7     8     4
     1     8     7     4     2     6     3     5     9
     4     3     8     6     5     9     2     7     1
     9     5     1     8     7     2     4     3     6
     7     6     2     1     4     3     5     9     8
     3     9     5     2     6     4     8     1     7
     2     7     6     3     8     1     9     4     5
     8     1     4     5     9     7     6     2     3
Sudoku.png
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 160,999评论 4 368
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 68,102评论 1 302
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 110,709评论 0 250
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 44,439评论 0 217
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 52,846评论 3 294
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,881评论 1 224
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 32,062评论 2 317
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,783评论 0 205
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,517评论 1 248
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,762评论 2 253
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 32,241评论 1 264
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,568评论 3 260
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 33,236评论 3 241
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,145评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,941评论 0 201
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,965评论 2 283
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,802评论 2 275

推荐阅读更多精彩内容

  • 个人学习批处理的初衷来源于实际工作;在某个迭代版本有个BS(安卓手游模拟器)大需求,从而在测试过程中就重复涉及到...
    Luckykailiu阅读 4,538评论 0 11
  • .bat脚本基本命令语法 目录 批处理的常见命令(未列举的命令还比较多,请查阅帮助信息) 1、REM 和 :: 2...
    庆庆庆庆庆阅读 7,945评论 1 19
  • 一、Python简介和环境搭建以及pip的安装 4课时实验课主要内容 【Python简介】: Python 是一个...
    _小老虎_阅读 5,620评论 0 10
  • 常言道“父爱如山、母爱如水”,若说母爱是能让人心头蕴热的温度,那么父爱就是一种让人高山仰止的高度。父爱是一...
    言之命至者阅读 1,237评论 0 5
  • 引用简书链接 @dynamic与@synthesize的区别 @property有两个对应的词,一个是@synth...
    YM_1阅读 5,333评论 0 4