codevs 2924 数独挑战
2924 数独挑战
“芬兰数学家因卡拉,花费3个月时间设计出了世界上迄今难度最大的数独游戏,而且它只有一个答案。因卡拉说只有思考能力最快、头脑最聪明的人才能破解这个游戏。”这是英国《每日邮报》2012年6月30日的一篇报道。这个号称“世界最难数独”的“超级游戏”,却被扬州一位69岁的农民花三天时间解了出来。
看到这个新闻后,我激动不已,证明我们OI的实力的机会来了,我们虽然不是思考能力最快、头脑最聪明的人,但是我们可以保证在1s之内解题。
好了废话不多说了……
数独是一种填数字游戏,英文名叫Sudoku,起源于瑞士,上世纪70年代由美国一家数学逻辑游戏杂志首先发表,名为Number Place,后在日本流行,1984年将Sudoku命名为数独,即“独立的数字”的省略,解释为每个方格都填上一个个位数。2004年,曾任中国香港高等法院法官的高乐德(Wayne Gould)把这款游戏带到英国,成为英国流行的数学智力拼图游戏。
玩家需要根据9×9盘面上的已知数字,推理出所有剩余位置(数据表示为数字0)的数字,并满足每一行、每一列、每一个粗线宫内的数字均含1-9,不重复。
现在给你一个数独,请你解答出来。每个数独保证有解且只有一个。
9行9列。
每个数字用空格隔开。0代表要填的数
行末没有空格,末尾没有回车。
输出答案。
排成9行9列。
行末没有空格,结尾可以有回车。
2 0 0 0 1 0 8 9 0
0 0 7 0 0 0 0 0 0
0 0 0 9 0 0 0 0 7
0 6 0 0 0 1 3 0 0
0 9 0 7 3 4 0 8 0
0 0 3 6 0 0 0 5 0
6 0 0 0 0 2 0 0 0
0 0 0 0 0 0 1 0 0
0 5 9 0 8 0 0 0 3
2 4 5 3 1 7 8 9 6
9 1 7 2 6 8 5 3 4
3 8 6 9 4 5 2 1 7
4 6 2 8 5 1 3 7 9
5 9 1 7 3 4 6 8 2
8 7 3 6 2 9 4 5 1
6 3 8 1 7 2 9 4 5
7 2 4 5 9 3 1 6 8
1 5 9 4 8 6 7 2 3
纯搜索,用empty[]数组记录下每个没有填数的位置,设n[i][l]=true数组表示第i行数l已经被填过,m[j][l]表示第j列数l已经被填过,k数组如下图所示
从第一个没有被填数的位置开始搜,填入一个数,更新m,n,k数组,注意回溯时也要更新m,n,k数组。
注意如果第i行j列没有填数,搜索时由于前面填的数的限制,什么数也不能填,那么此时一定要回退,直至这个格填上数再进行下一个,否则死循环。
#include<iostream>
#include<cstdlib>
#include<cstdio>
using namespace std;
int a[][],empty[],cnt;
bool n[][],m[][],k[][];
bool find(int i,int j,int l)//判断位置i,j填p是否合适
{
if(i<=&&j<=&&k[][l]) return false;
if(i<=&&j>=&&j<=&&k[][l]) return false;
if(i<=&&j>=&&k[][l]) return false;
if(i>=&&i<=&&j<=&&k[][l]) return false;
if(i>=&&i<=&&j>=&&j<=&&k[][l]) return false;
if(i>=&&i<=&&j>=&&k[][l]) return false;
if(i>=&&j<=&&k[][l]) return false;
if(i>=&&j>=&&j<=&&k[][l]) return false;
if(i>=&&j>=&&k[][l]) return false;
if(n[i][l]) return false;
if(m[j][l]) return false;
return true;
}
void mark(int i,int j,int l)//标记i,j位置填上了l
{
if(i<=&&j<=) k[][l]=true;
else if(i<=&&j>=&&j<=) k[][l]=true;
else if(i<=&&j>=)k[][l]=true;
else if(i>=&&i<=&&j<=) k[][l]=true;
else if(i>=&&i<=&&j>=&&j<=) k[][l]=true;
else if(i>=&&i<=&&j>=) k[][l]=true;
else if(i>=&&j<=) k[][l]=true;
else if(i>=&&j>=&&j<=) k[][l]=true;
else if(i>=&&j>=) k[][l]=true;
n[i][l]=true;
m[j][l]=true;
}
void back(int i,int j,int l)//位置i,j撤回l
{
if(i<=&&j<=) k[][l]=false;
else if(i<=&&j>=&&j<=) k[][l]=false;
else if(i<=&&j>=)k[][l]=false;
else if(i>=&&i<=&&j<=) k[][l]=false;
else if(i>=&&i<=&&j>=&&j<=) k[][l]=false;
else if(i>=&&i<=&&j>=) k[][l]=false;
else if(i>=&&j<=) k[][l]=false;
else if(i>=&&j>=&&j<=) k[][l]=false;
else if(i>=&&j>=) k[][l]=false;
n[i][l]=false;
m[j][l]=false;
}
void dfs(int now,int s)
{
if(s==)
{
for(int i=;i<=;i++)
{
for(int j=;j<=;j++) cout<<a[i][j]<<' ';
cout<<endl;
}
exit();
}
for(int i=now+;i<=cnt;i++)
{
bool ok=false;
if(!empty[i]) continue;//利用empty[i]是否为0来判断empty存的这个位置是否填上了数,0代表已经填上了
int x=empty[i]/,y=empty[i]%;
int h=empty[i];
for(int p=;p<=;p++)
if(find(x,y,p))
{
empty[i]=;
mark(x,y,p);//标记
a[x][y]=p;
dfs(i,s+);
back(x,y,p);//撤回标记
empty[i]=h;
}
if(empty[i]) return;//程序能运行的关键,没有这句就死循环了
}
}
int main()
{
for(int i=;i<=;i++)
for(int j=;j<=;j++)
{
cin>>a[i][j];
if(!a[i][j]) {empty[++cnt]=i*+j;continue;}
n[i][a[i][j]]=true;
m[j][a[i][j]]=true;
if(i<=&&j<=) k[][a[i][j]]=true;
else if(i<=&&j>=&&j<=) k[][a[i][j]]=true;
else if(i<=&&j>=) k[][a[i][j]]=true;
else if(i>=&&i<=&&j<=) k[][a[i][j]]=true;
else if(i>=&&i<=&&j>=&&j<=) k[][a[i][j]]=true;
else if(i>=&&i<=&&j>=) k[][a[i][j]]=true;
else if(i>=&&j<=) k[][a[i][j]]=true;
else if(i>=&&j>=&&j<=) k[][a[i][j]]=true;
else k[][a[i][j]]=true;
}
dfs(,-cnt);
}
for(int i=now+1;i<=cnt;i++)
{
bool ok=false;//**********
if(!empty[i]) continue;
int x=empty[i]/10,y=empty[i]%10;
int h=empty[i];
for(int p=1;p<=9;p++)
if(find(x,y,p))
{
ok=true;//********
empty[i]=0;
mark(x,y,p);
a[x][y]=p;
dfs(i,s+1);
back(x,y,p);
empty[i]=h;
}
if(!ok) return; //************
}
利用ok变量来判断empty[i]位置是否填上了数,但忽略了ok在回溯中没有更改,当这个位置被标记之后标记不会撤回。
更改方法:回溯时添加ok=false
错误原因:深搜回溯问题
错误1:代码及错误分析,警示
错误2:样例输出中,数与数之间有空格,这要是noip的话,100分就没了( ⊙ o ⊙ )!,注意注意!!!
上一个代码写的很丑,主要原因是k数组的调用,看了题解,发现k数组的第一维可以用如下式子计算:
(i-1)/3*3+(j-1)/3+1 代码长度优化版本:
#include<cstdio>
#include<cstring>
#include<cstdlib>
int a[][];
bool h[][],l[][],g[][];
void print()
{
for(int i=;i<=;i++)
{
for(int j=;j<=;j++)
printf("%d ",a[i][j]);
printf("\n");
}
exit();
}
void go(int x,int y)
{
if(a[x][y]!=)
{
if(x==&&y==) print();
if(y==) go(x+,); else go(x,y+);
}
if(a[x][y]==)
{
for(int i=;i<=;i++)
if(h[x][i]&&l[y][i]&&g[(x-)/*+(y-)/+][i])
{
a[x][y]=i;
h[x][i]=false;
l[y][i]=false;
g[(x-)/*+(y-)/+][i]=false;
if(x==&&y==) print();
if(y==) go(x+,); else go(x,y+);
a[x][y]=;
h[x][i]=true;
l[y][i]=true;
g[(x-)/*+(y-)/+][i]=true;
}
}
}
int main()
{
memset(h,true,sizeof(h));
memset(l,true,sizeof(l));
memset(g,true,sizeof(g));
for(int i=;i<=;i++)
for(int j=;j<=;j++)
{
scanf("%d",&a[i][j]);
if(a[i][j]>)
{
h[i][a[i][j]]=false;
l[j][a[i][j]]=false;
g[(i-)/*+(j-)/+][a[i][j]]=false;
}
}
go(,);
}
codevs 2924 数独挑战的更多相关文章
- codevs 2924 数独挑战 x(三种做法+超详细注释~)
2924 数独挑战 时间限制: 1 s 空间限制: 1000 KB 题目等级 : 钻石 Diamond 题目描述 Description “芬兰数学家因卡拉,花费3个月时间设计出了世界上迄今 ...
- codevs 2924 数独
数独挑战 Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://codevs.cn/problem/2924/ Description “芬兰数学家因 ...
- 数独挑战(codevs 2924)
2924 数独挑战 时间限制: 1 s 空间限制: 1000 KB 题目等级 : 钻石 Diamond 题解 查看运行结果 题目描述 Description “芬兰数学家因卡拉,花费3 ...
- codevs2924 数独挑战
2924 数独挑战 时间限制: 1 s 空间限制: 1000 KB 题目等级 : 钻石 Diamond 题解 题目描述 Description "芬兰数学家因卡拉,花费3个月时间设计出了世界 ...
- 洛谷P1784 数独
To 洛谷.1784 数独(类似(或者说相同)题:CODEVS.4966 简单数独(4*4数独) CODEVS.2924 数独挑战) 题目描述 数独是根据9×9盘面上的已知数字,推理出所有剩余空格的数 ...
- 洛谷 P1784 数独[DFS/回溯]
To 洛谷.1784 数独类似题:CODEVS.4966 简单数独(4*4数独) CODEVS.2924 数独挑战) 题目描述 数独是根据9×9盘面上的已知数字,推理出所有剩余空格的数字,并满足每一行 ...
- Dancing Links & Algorithm X
1 参考链接 http://www.cnblogs.com/steady/archive/2011/03/15/1984791.html#undefined http://en.wikipedia.o ...
- CODEVS 1174 靶形数独
题目描述 Description 小城和小华都是热爱数学的好学生,最近,他们不约而同地 迷上了数独游戏,好胜的他们想用数独来一比高低.但普通 的数独对他们来说都过于简单了,于是他们向Z 博士请教,Z ...
- 靶形数独(codevs 1174)
1174 靶形数独 2009年NOIP全国联赛提高组 时间限制: 4 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题解 查看运行结果 题目描述 Descri ...
随机推荐
- Xdebug文档(六) 分析PHP脚本
分析PHP脚本Xdebug内置分析器能让你找到脚本中的瓶颈并用额外的工具诸如KcacheGrind或WinCacheGrind工具可视化. 介绍 Xdebug分析器是分析PHP代码和判断瓶颈或确定代码 ...
- 基于WCF MSMQ 的企业应用解决方案
最近研究了一下基于MSMQ的WCF应用,从书上.网上查了很多资料,但始终没能彻底理解WCF-MSMQ的工作原理,也没能得到一个合理的应用解决方案.索性还是自己做个实验,探索一下吧.经过反复试验,颇有收 ...
- less hack 兼容
less hack 兼容 css做兼容是在所难免的,那么用less写css代码时怎样hack呢?倘若用css的方法直接在后面写上类似“\9”编译是要报错的.下面是我尝试的两个小方法仅供参考: 1. ...
- IE8兼容模式设置
设置---兼容性视图设置--添加此网站--在IE8中调试(防止调整IE内核后浏览器崩溃,360可通过设置极速模式-兼容模式 点击地址栏绿色图标)
- CSS3 Transitions, Transforms和Animation的使用
一.前言 CSS3动画相关的几个属性是:transition, transform, animation:分别理解为过渡,变换,动画.虽意义相近,但具体的功能和在CSS3中承担的工作有一定的差异. t ...
- SharePoint服务器端对象模型 之 使用CAML进展数据查询
SharePoint服务器端对象模型 之 使用CAML进行数据查询 一.概述 在SharePoint的开发应用中,查询是非常常用的一种手段,根据某些筛选.排序条件,获得某个列表或者某一些列表中相应的列 ...
- WeakReference在Handler中的应用
public class AutoActivity extends Activity { Handler handler = new Handler(){ public void handleMess ...
- IOS开发基础知识--碎片24
1:兼容字体大小6plue跟它以下的区别 #define FONT_COMPATIBLE_SCREEN_OFFSET(_fontSize_) [UIFont systemFontOfSize:(_fo ...
- GCD
Grand Central Dispatch(GCD)包含语言特点.运行库和系统增强功能,它提供系统的.综合的改进,进而去支持并发代码在iOS和OSX多核硬件上的执行. 子系统.CoreFoundat ...
- 软件工程--界面UI 的原型设计
经过今天的讨论, 我们组初步给出了设计方案. 安装了一晚上的 mockupBuilder 终于绘制了几份界面的原型图, 这里需要吐槽下 mockupBuilder, 这个软件很好用, 不过网页版和单机 ...