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 256
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 303
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 1538
<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 369
题意： 第三道数独题目了，还是深搜，这次的空白数量固定5个，开始一直WA，发现了题目中重要的描述，可能给出的puzzle本身就是错的。好吧。 然后就和前两道是一样的了。 代码： #include #include #include using namespace std; char square[10][10]; bool row[10][10], col[10][10], b
讯图GD460电脑DSP调音软件下载
最新发布
07-23
讯图GD460电脑DSP调音软件下载
港股历史逐笔十档订单簿分钟日数据下载
07-23
csv格式港股历史高频数据集（GB2312编码） ─────────────────────────────────── 【数据种类】逐笔成交｜订单簿快照｜分钟K线｜多周期聚合数据 【时间跨度】2005年至今完整历史数据 【数据精度】 - Level2逐笔数据：包含十档订单簿及逐笔成交明细（≥1TB） - 多频段K线：1/5/15/30/60分钟+日/周/月线（≥100GB） 【技术特性】 ● 存储方式：网盘加密压缩包 ● 数据规范：经时间戳对齐及异常值清洗 ● 格式兼容：原生csv适配Python(pandas.read_csv)及量化平台 【应用场景】 ◆ 高频策略验证：ARIMA/LSTM时序预测、做市商策略模拟、冲击成本测算 ◆ 风险建模：市场操纵模式识别、流动性风险压力测试 ◆ AI训练：深度强化学习(DRL)策略生成、OrderBook特征工程 ◆ 学术研究：市场微观结构分析、波动率聚类效应验证 【质量保障】 √ 订单簿重构误差＜0.1% √ Tick数据时延校准至交易所时钟 √ 提供完整的symbol-mapping表及除权处理说明 本数据集通过PCAOB审计标准校验，适用于《金融研究》等核心期刊论文的实证分析，满足资管机构合规性回测要求。 【版权说明】 数据源：银禾金融数据库，解释权归该数据库所有。
虚拟同步机离网模型与30KW PCS储能逆变技术研究：前端BUCK-BOOST电路与多功能的并网控制
07-23
内容概要：本文详细介绍了30KW PCS储能逆变器的虚拟同步机离网模型及其核心技术。首先阐述了离网与并网功能的融合，强调了逆变器在不同状态下的稳定性和可靠性。接着深入探讨了前端BUCK-BOOST电路的作用，解释了其对电压升降调节的关键意义。随后重点讲解了后端三电平逆变器的工作原理，特别是下垂控制加虚拟同步机控制策略的应用，实现了多台PCS的并联运行。此外，还讨论了C语言编程和S-function在DSP28XX系列芯片上的实现，提供了详细的代码和理论计算支持。最后，总结了该逆变器的速度与效率优势，展示了其在电力系统中的广泛应用前景。 适合人群：从事电力电子领域的工程师、研究人员及高校相关专业师生。 使用场景及目标：①理解和掌握30KW PCS储能逆变器的离网并网功能；②学习前端BUCK-BOOST电路和后端三电平逆变器的设计原理；③熟悉虚拟同步机控制技术和C语言编程在DSP芯片上的应用。 其他说明：本文基于大厂内部资料，提供了详尽的理论计算和实际问题解决方案，有助于提高产品研发和应用能力。
永磁同步电机直接转矩控制仿真的MATLAB 2010a版
07-23
内容概要：本文详细介绍了如何使用MATLAB 2010a进行永磁同步电机(PMSM)的直接转矩控制(DTC)仿真。首先，文章提供了PMSM的关键参数配置，强调了电感参数单位的选择。接着，逐步讲解了磁链观测器、滞环比较器以及电压矢量选择表的设计方法，并给出了具体的MATLAB代码实现。此外，还分享了一些实用的调参技巧，如调整磁链观测器、滞环宽度等，以优化仿真效果。最后展示了仿真结果图，包括电磁转矩跟踪、转速曲线和磁链轨迹，帮助读者更好地理解和验证DTC的工作原理。 适合人群：具有一定MATLAB基础并希望深入了解电机控制领域的工程师和技术爱好者。 使用场景及目标：适用于需要研究或应用永磁同步电机直接转矩控制技术的研究人员、学生或工程师。主要目标是掌握DTC的基本原理及其在MATLAB环境下的具体实现方法，同时学会如何根据实际情况调整参数以获得更好的性能。 阅读建议：由于文中涉及较多的技术细节和代码片段，在阅读过程中建议读者跟随作者提供的步骤亲自尝试编写代码，并结合理论知识深入理解每个模块的功能。对于遇到的问题可以通过查阅相关文献资料或者在线社区寻求解决方案。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值