Sudoku

中英对译题目

Sudoku
数独
题目描述
In the game of Sudoku, you are given a large 9 × 9 grid divided into smaller 3 × 3 subgrids. For example,
在数独游戏中，你会得到一个9×9的大网格，分为较小的3×3子网格。例如

Given some of the numbers in the grid, your goal is to determine the remaining numbers such that the numbers 1 through 9 appear exactly once in

(1) each of nine 3 × 3 subgrids, (2) each of the nine rows, and (3) each of the nine columns.
给定网格中的一些数字，你的目标是确定剩余的数字，使数字1到9在

（1）九个3×3子网格中的每一个中（2）九行中的每行和（3）九列中的每列中都恰好出现一次。

输入输出样例

输入样例1：

.2738..1..1...6735.......293.5692.8...........6.1745.364.......9518...7..8..6534. ......52..8.4......3...9...5.1...6..2..7........3.....6...1..........7.4.......3. end

输出样例1：

527389416819426735436751829375692184194538267268174593643217958951843672782965341 416837529982465371735129468571298643293746185864351297647913852359682714128574936

【耗时限制】1000ms 【内存限制】256MB

---

常见剪枝技巧

1. 可行性剪枝： 在搜索过程中，及时对当前状态进行检查，如果发现分支已经无法生成有效的解就剪枝回溯。某些题目条件的范围限制是一个区间，此时可行性剪枝也被称为“上下界剪枝”。

2. 最优性剪枝： 在最优化问题的搜索中，如果当前的代价已经超过了之前搜到的最优解，则可以剪枝回溯。

3. 记忆化搜索： 在搜索过程中记录每个状态的搜索结果，在重复遍历一个状态时直接检索并返回，记忆化搜索时一种特殊的最优性剪枝。

4. 排除等效冗余： 在搜索过程中，如果我们能够判定从搜索树的当前节点上沿着某几条不同分支到达的子树是等效的，那么只需要对其中的一条分支执行搜索。

5. 优化搜索顺序： 在一些搜索问题中，搜索树的各个层次、分支之间的顺序不是一定的，不同的搜索顺序会产生不同的搜索树形态，其规模大小也可能相差甚远。

---

题目正式讲解

搜索：依次搜索每个空位置，判断应该填写的数字，如果某位置无可选数字则回溯。

如何判断某位置是否可以填某数字？

① 枚举：枚举该位置所在的行、列和宫中的所有数字，通过数组标记法判断是否冲突。

② 常数优化：3个9位二进制数a, b, c分别标记每行、每列、每宫是否已经有1-9中的数字。a) 用二进制从右往左的1-9位，分别代表1-9。b) 用二进制0，表示该数字被占用，1表示该数组未被占用。b) a & b & c的结果中，二进制位为1的表示可选的数。

                                    

性能优化：优化搜索顺序，优先搜索可用数字最少的位置，可以极大的降低搜索量。

---

为了方便大家拿代码

楼主把带注释的代码截了个图

一下是没有注释的代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <sstream>
using namespace std;
const int N=15;
string s;
int h[N],l[N],g[N],idle[1<<N],num[1<<N],cnt;
inline int GID(int x,int y){return x/3*3+y/3;}
inline int lowbit(int x){return x&-x;}
void turn(int i,int k){
    int x=i/9,y=i%9,z=GID(x,y);
    h[x]^=1<<(k-1);
    l[y]^=1<<(k-1);
    g[z]^=1<<(k-1);
}
bool dfs(int cur){
    if(cur==cnt) return true;
    int k=10,t;
    for(int i=0;i<81;i++){
        if(s[i]=='.'){
            int x=i/9,y=i%9,z=GID(x,y);
            int w=h[x]&l[y]&g[z];
            if(idle[w]<k){k=idle[w];t=i;}
        }
    }
    int x=t/9,y=t%9,z=GID(x,y);
    int w=h[x]&l[y]&g[z];
    while(w){
        int now=num[lowbit(w)];
        s[t]=now+'0';turn(t,now);
        if(dfs(cur+1)) return true;
        turn(t,now);s[t]='.';
        w-=lowbit(w);
    }
    return 0;
}
void solve(){
    for(int i=0;i<9;i++) h[i]=l[i]=g[i]=(1<<9)-1;
    cnt=0;
    for(int i=0;i<81;i++){
        if(s[i]!='.') turn(i,s[i]-'0');
        else ++cnt;
    }
    if(dfs(0))
        for(int i=0;i<81;i++) cout<<s[i];
    cout<<endl;
}
inline int count(int x){
    int cnt=0;
    while(x){
        cnt++;
        x&=x-1;
    }
    return cnt;
}
int main(){
    int maxs=(1<<9)-1;
    for(int i=0;i<=maxs;i++) idle[i]=count(i);
    for(int i=0;i<9;i++) num[1<<i]=i+1;
    while(cin>>s&&s[0]!='e') solve();
    return 0;
}