soj 1162. Sudoku[Special judge] 1317. Sudoku

最新推荐文章于 2020-01-23 10:29:56 发布
最新推荐文章于 2020-01-23 10:29:56 发布
阅读量1.4k 收藏 1
点赞数
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.youkuaiyun.com/jasison/article/details/8457637
本文探讨了数独求解问题的两种优化算法实现：一种是传统的暴力深搜，另一种是在剪枝的基础上改进的深搜算法。通过预处理计算每个待求解点的可能情况数并进行排序，显著减少了搜索深度，从而提高了求解效率。实验证明，剪枝后的算法在处理数独难题时，性能有了质的飞跃，从超时降低到几乎即时得出答案。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

两道数独求解的问题。

思路：

很显然，两道题目的思路都是，读入数据记录所有待求解的点，然后进行深搜即可。

但是不难看出的一点是，如果是单纯的暴力深搜，搜索的深度可能达到81，这样做可能会超时，于是要进行剪枝。

剪枝的方法很简单，预先算出进行下一步扩展的点的个数，然后进行排序，扩展点个数少的先搜索。

由于本题是求解数独，在最初的时候对各个待求解的点根据可能的情况数排序之后确定搜索的先后顺序即可。

方法：

记录每一行每一列每一个数字是否出现，记录所有带求解的点，深搜，剪枝。


1162 剪枝前代码：

#include <cstring>
#include <cstdio>

char square[10][10];
bool row[10][10], col[10][10], blo[10][10];
int x[85], y[85], todo, flag;
int block(int i, int j) {return (i-1)/3*3 + (j-1)/3 + 1;}
void dfs(int n)
{
    if (flag || n == todo) { flag = 1; return; }
    int xx = x[n], yy = y[n], bl = block(xx,yy);
    for (int temp = 1; temp <= 9; ++ temp)
    {
        if (!row[xx][temp] && !col[yy][temp] && !blo[bl][temp])
        {
            square[xx][yy] = temp+'0';
            row[xx][temp] = 1; col[yy][temp] = 1; blo[bl][temp] = 1;
            dfs(n+1); if (flag) return;
            row[xx][temp] = 0; col[yy][temp] = 0; blo[bl][temp] = 0;
        }
    }
}

int main()
{
    int t;
    scanf("%d", &t);
    while (t --)
    {
        flag = 0; todo = 0;
        memset(row, false , sizeof(row));
        memset(col, false , sizeof(col));
        memset(blo, false , sizeof(blo));
        for (int i = 1; i <= 9; ++ i)
        {
            scanf("%s", square[i]+1);
            for (int j = 1; j <= 9; ++ j)
            {
                if (square[i][j] == '0')
                {
                    x[todo] = i;
                    y[todo] = j;
                    todo ++;
                }
                else
                {
                    row[i][square[i][j]-'0'] = 1;
                    col[j][square[i][j]-'0'] = 1;
                    blo[block(i,j)][square[i][j]-'0'] = 1;
                }
            }
        }
        dfs(0);
        for (int i = 1; i <= 9; ++ i)
            printf("%s\n", square[i]+1);
    }
}
结果：Accepted 0.35s


剪枝后代码：

#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;

char square[10][10];
bool row[10][10], col[10][10], blo[10][10];
int todo, flag;
int block(int i, int j) {return (i-1)/3*3 + (j-1)/3 + 1;}
struct TODO
{
	int x, y, b, pNum;
	TODO() {}
	TODO(int _x, int _y): x(_x), y(_y), b(block(_x,_y)), pNum() {}
	void calc()
	{
		for (int i = 1; i <= 9; ++ i)
			if (!row[x][i] && !col[y][i] && !blo[b][i])
				pNum ++;
	}
};
bool cmp(const TODO &a, const TODO &b) {return a.pNum < b.pNum;}
TODO p[85];
void dfs(int n)
{
	if (flag || n == todo) { flag = 1; return; }
	int xx = p[n].x, yy = p[n].y, bl = p[n].b;
	for (int temp = 1; temp <= 9; ++ temp)
	{
		if (!row[xx][temp] && !col[yy][temp] && !blo[bl][temp])
		{
			square[xx][yy] = temp+'0';
			row[xx][temp] = 1; col[yy][temp] = 1; blo[bl][temp] = 1;
			dfs(n+1); if (flag) return;
			row[xx][temp] = 0; col[yy][temp] = 0; blo[bl][temp] = 0;
		}
	}
}

int main()
{
	int t;
	scanf("%d", &t);
	while (t --)
	{
		flag = 0; todo = 0;
		memset(row, false , sizeof(row));
		memset(col, false , sizeof(col));
		memset(blo, false , sizeof(blo));
		for (int i = 1; i <= 9; ++ i)
		{
			scanf("%s", square[i]+1);
			for (int j = 1; j <= 9; ++ j)
			{
				if (square[i][j] == '0')
					p[todo++] = TODO(i,j);
				else
				{
					row[i][square[i][j]-'0'] = 1;
					col[j][square[i][j]-'0'] = 1;
					blo[block(i,j)][square[i][j]-'0'] = 1;
				}
			}
		}
		for (int i = 0; i < todo; ++ i) p[i].calc();
		sort(p, p+todo, cmp);
		dfs(0);
		for (int i = 1; i <= 9; ++ i)
			printf("%s\n", square[i]+1);
	}
}
结果：Accepted 0.02s


1317剪枝前代码：

#include <cstring>
#include <cstdio>
#define N 15
char square[N][N];
bool row[N][N], col[N][N], blo[N][N];
int x[85], y[85], ans[N][N], todo, num;
inline int block(int i, int j) {return (i-1)/3*3 + (j-1)/3 + 1;}
void dfs(int n)
{
    if (n == todo)
    {
        num ++;
        if (num == 1)
            for (int i = 1; i <= 9; ++ i)
                for (int j = 1; j <= 9; ++ j)
                    ans[i][j] = square[i][j]-'0';
        return;
    }
    int xx = x[n], yy = y[n], bl = block(xx,yy);
    for (int temp = 1; temp <= 9; ++ temp)
    {
        if (!row[xx][temp] && !col[yy][temp] && !blo[bl][temp])
        {
            square[xx][yy] = temp+'0';
            row[xx][temp] = 1; col[yy][temp] = 1; blo[bl][temp] = 1;
            dfs(n+1);
            row[xx][temp] = 0; col[yy][temp] = 0; blo[bl][temp] = 0;
        }
    }
}

int main()
{
    int t;
    scanf("%d", &t);
    for (int Case = 1; Case <= t; ++ Case)
    {
        if (Case != 1) printf("\n");
        num = 0; todo = 0;
        memset(row, false , sizeof(row));
        memset(col, false , sizeof(col));
        memset(blo, false , sizeof(blo));
        for (int i = 1; i <= 9; ++ i)
        {
            scanf("%s", square[i]+1);
            for (int j = 1; j <= 9; ++ j)
            {
                if (square[i][j] == '_')
                {
                    x[todo] = i;
                    y[todo] = j;
                    todo += 1;
                }
                else
                {
                    row[i][square[i][j]-'0'] = 1;
                    col[j][square[i][j]-'0'] = 1;
                    blo[block(i,j)][square[i][j]-'0'] = 1;
                }
            }
        }
        dfs(0);
        if (num == 0) printf("Puzzle %d has no solution\n", Case);
        else if (num == 1)
        {
            printf("Puzzle %d solution is\n", Case);
            for (int i = 1; i <= 9; ++ i)
            {
                for (int j = 1; j <= 9; ++ j)
                    printf("%d", ans[i][j]);
                printf("\n");
            }
        }
        else printf("Puzzle %d has %d solutions\n", Case, num);
    }
}
结果：Time Limit Exceeded 


剪枝后代码：

#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;

#define N 12
char square[N][N], ans[N][N];
bool row[N][N], col[N][N], blo[N][N];
int todo, num;
int block(int i, int j) {return (i-1)/3*3 + (j-1)/3 + 1;}
struct TODO
{
	int x, y, b, pNum;
	TODO() {}
	TODO(int _x, int _y): x(_x), y(_y), b(block(_x,_y)), pNum(0) {}
	void calc()
	{
		for (int i = 1; i <= 9; ++ i)
			if (!row[x][i] && !col[y][i] && !blo[b][i])
				pNum ++;
	}
};
bool cmp(const TODO &a, const TODO &b) {return a.pNum < b.pNum;}
TODO p[85];
void dfs(int n)
{
	if (n == todo)
	{
		num ++;
		if (num == 1) memcpy(ans, square, sizeof(square));
		return;
	}
	int xx = p[n].x, yy = p[n].y, bl = p[n].b;
	for (int temp = 1; temp <= 9; ++ temp)
	{
		if (!row[xx][temp] && !col[yy][temp] && !blo[bl][temp])
		{
			square[xx][yy] = temp+'0';
			row[xx][temp] = 1; col[yy][temp] = 1; blo[bl][temp] = 1;
			dfs(n+1);
			row[xx][temp] = 0; col[yy][temp] = 0; blo[bl][temp] = 0;
		}
	}
}

int main()
{
	int t;
	scanf("%d", &t);
	for (int Case = 1; Case <= t; ++ Case)
	{
		if (Case != 1) printf("\n");
		num = 0; todo = 0;
		memset(row, false , sizeof(row));
		memset(col, false , sizeof(col));
		memset(blo, false , sizeof(blo));
		for (int i = 1; i <= 9; ++ i)
		{
			scanf("%s", square[i]+1);
			for (int j = 1; j <= 9; ++ j)
			{
				if (square[i][j] == '_')
					p[todo++] = TODO(i,j);
				else
				{
					row[i][square[i][j]-'0'] = 1;
					col[j][square[i][j]-'0'] = 1;
					blo[block(i,j)][square[i][j]-'0'] = 1;
				}
			}
		}
		for (int i = 0; i < todo; ++ i) p[i].calc();
		sort(p, p+todo, cmp);
		dfs(0);
		if (num == 0) printf("Puzzle %d has no solution\n", Case);
		else if (num == 1)
		{
			printf("Puzzle %d solution is\n", Case);
			for (int i = 1; i <= 9; ++ i)
				printf("%s\n", ans[i]+1);
		}
		else printf("Puzzle %d has %d solutions\n", Case, num);
	}
}
结果：Accepted 0.3s

jasison
关注
SoJ0825.github.io
02-14
【标题】"SoJ0825.github.io" 指的是一个个人或组织在GitHub上托管的静态网站。在GitHub上，用户可以利用GitHub Pages服务创建并发布自己的个人网页或者项目页面，这个URL通常以用户名.github.io的形式呈现。在这个...
sicily 1317 Sudoku
anhong585258的博客
01-25 257
Description Sudoku is a placement puzzle. The goal is to enter a symbol in each cell of a grid, most frequently a9x9grid made up of3x3subgrids. Each row, column and subgrid must contain...
sicily 1317. Sudoku
Dijason的专栏
11-29 1273
题目链接sicily 1317 这题咋一看以为不简单，因为第一感觉是时间复杂度很大，可想想写起来应该不难，于是飞快的写了个dfs。 结果是证明了我的感觉，超时了。确实是写起来容易的代码来对付时间复杂度很大很大的是会超时的。 只能做些优化了。看了别人的做法，也想到了之前马的周游一样的优化方法。从能够填写的数字个数最少的位置开始 深搜。原因是产生的树上端都是小，然后搜到下面的位置的话，一旦不符合很快
1317
ailinyingai的博客
01-23 304
def run(n): def no_zero(num): return False if '0' in str(num) else True for m in range(n): if no_zero(m) and no_zero(n-m): return [m, n-m]
sicily1317_数独
yaozhiyi的专栏
06-24 1504
深度优先搜索即可题目：ProblemSudoku is a placement puzzle. The goal is to enter a symbol in each cell of a grid, most frequently a 9 x 9 grid made up of 3 x 3 subgrids. Each row, column and subgrid must contain only one instance of each symbol. Sudoku initially becam
深入探讨SoJ0825.github.io的HTML实现
2. SoJ0825.github.io：这个名字表明这个GitHub页面的仓库名为“SoJ0825”，这是用户名或组织名，后面跟着“.github.io”格式表明这是一个个人或组织的GitHub页面。用户通过这种命名的仓库可以在GitHub上创建一个...
soj.rar_soj 1488_soj 1822
09-24
在给定的压缩包"soj.rar"中，包含了三道来自SOJ（ ACM/ICPC Online Judge）的编程竞赛题目：1019、1488 和 1822。这些题目主要涉及动态规划这一算法设计与分析的重要方法。动态规划是一种通过将复杂问题分解为更小的...
soj.rar_OJ PHP_oj_oj系统_soj
09-23
【标题】"soj.rar_OJ PHP_oj_oj系统_soj" 提供的是一个在线判题系统（Online Judge，简称OJ）的PHP实现，名为"soj"。这个系统主要用于编程竞赛或教育场景，它能自动评估参赛者提交的代码是否符合预设测试用例，从而...
soj.rar_1876_sicily 1022_sicily 1934_sicily soj IP
09-21
标题中的"soj.rar_1876_sicily 1022_sicily 1934_sicily soj IP" 提到了几个关键元素，这些元素指向了一个编程竞赛或者作业集中的多个问题。"soj"通常代表"Software Judge"，这是一个用于在线编程竞赛或教育平台的...
Sicily1317-Sudoku-位运算暴搜
denghuang8508的博客
12-30 205
最终代码地址：https://github.com/laiy/Datastructure-Algorithm/blob/master/sicily/1317.c 这题博主刷了1天，不是为了做出来，AC之后在那死磕性能... 累积交了45份代码，纪念一下- - 以上展示了从1.25s优化到0.03s的艰苦历程... 来看题目吧，就是一个数独求解的题： 1317. ...
sicily 1317
alfredtofu的专栏
01-08 1539
<br />这题不难，跟马的周游问题一样<br />开始没优化，就直接穷举（看到10秒嘛。。），结果TL，后来优化一下，对可能出现的数字少的点先进行枚举。<br />思路就这样，可是我做了四个小时。。。。。。。。。最后我认真认真的看一次题目，发现如果有多个解，solution要加s。。。。这基本的复数形式我错了。。。o(︶︿︶)o 唉<br /> <br />#include <iostream> #include <cstring> using namespace std; char a[9][9]
SOJ 1016
Tsinting的博客
02-27 422
1016. Boundaries on A New Kind of Science Constraints Time Limit: 10 secs, Memory Limit: 32 MB Description Stephen Wolfram in a “New Kind of Science” describes a special kind of cellular automat
「笔记」深搜的剪枝技巧
Loceaner的博客
12-24 84
深搜的剪枝技巧 剪枝方法 1.优化搜索顺序 在不同的问题中，搜索树的各个层次、各个分支之间的顺序不是固定的，不同的搜索顺序会产生不同的搜索树形态，其规模大小也相差甚远。 2.排除等效冗余 在搜索过程中，若能判断从搜索树当前节点上沿某几条不同分支到达的子树是相同的，那么只需对其中一条分支执行搜索。 3.可行性剪枝 可行性剪枝也叫上下界剪枝，其是指在搜索过程中，及时对当前状态进行检查，若发现分支已无...
soj 1664. Su-Su-Sudoku
jasison
01-05 370
题意： 第三道数独题目了，还是深搜，这次的空白数量固定5个，开始一直WA，发现了题目中重要的描述，可能给出的puzzle本身就是错的。好吧。 然后就和前两道是一样的了。 代码： #include #include #include using namespace std; char square[10][10]; bool row[10][10], col[10][10], b
MATLAB GUI设计平台：多算法雷达一维恒虚警检测CFAR可视化界面，实现噪声波形与目标检测回波-检测门限波形图的可视化
07-25
内容概要：本文介绍了利用MATLAB GUI设计平台构建的一个多算法雷达一维恒虚警检测(CFAR)可视化界面。该界面允许用户选择不同的噪声类型（如高斯白噪声、瑞利噪声）、目标类型（如点目标、扩展目标），并选择多种CFAR算法（如OS-CFAR、CA-CFAR）。用户还可以手动输入关键参数，如噪声功率、目标回波功率和门限值。系统能够实时生成并显示噪声波形和目标检测的回波-检测门限波形图，帮助用户直观地评估不同算法在各种噪声条件下的表现。 适合人群：从事雷达信号处理研究的专业人士、高校师生及相关领域的科研工作者。 使用场景及目标：适用于需要对比和优化雷达系统中CFAR算法性能的研究项目，旨在提升雷达系统的检测精度和适应能力。 其他说明：通过这个可视化工具，研究人员可以更加便捷地探索不同参数配置下各CFAR算法的效果，为实际工程应用提供理论支持和技术指导。
49页PPTAI赋能制造业-250624165915.pptx
07-25
【49页PPT】AI赋能制造业_250624165915.pptx
Delphi 12.3控件之BCompare-zh-5.1.2.31185.rar
07-25
Delphi 12.3控件之BCompare-zh-5.1.2.31185.rar
MATLABSimulink中小电流系统单相接地故障仿真的建模与分析 电力系统仿真 v2.1
最新发布
07-25
利用MATLAB/Simulink构建小电流系统单相接地故障仿真模型的方法。主要探讨了三种不同的中性点接地方式（不接地、经消弧线圈接地、经小电阻接地）下，故障特征的变化情况。文中不仅提供了具体的模型搭建步骤，还展示了如何通过Simulink进行波形数据的采集与分析，如零序电流、零序电压及其五次谐波波形图的获取方法。此外，作者分享了一些实用技巧，比如正确配置Powergui模块以确保仿真精度。 适用人群：电力系统研究人员、高校师生以及从事电力工程相关工作的技术人员。 使用场景及目标：适用于需要深入理解不同类型中性点接地方式对电网稳定性影响的研究项目；也可作为教学案例帮助学生掌握电力系统故障诊断的基本原理和技术手段。 其他说明：文中涉及大量MATLAB/Simulink的具体操作细节，对于初学者来说可能具有一定难度，建议有一定编程基础后再行尝试。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值