数独游戏

你一定听说过“数独”游戏。

如图，玩家需要根据9×9盘面上的已知数字，推理出所有剩余空格的数字，并满足每一行、每一列、每一个同色九宫内的数字均含1-9，不重复。

数独的答案都是唯一的，所以，多个解也称为无解。

本图的数字据说是芬兰数学家花了3个月的时间设计出来的较难的题目。但对会使用计算机编程的你来说，恐怕易如反掌了。

本题的要求就是输入数独题目，程序输出数独的唯一解。我们保证所有已知数据的格式都是合法的，并且题目有唯一的解。

格式要求，输入9行，每行9个数字，0代表未知，其它数字为已知。
输出9行，每行9个数字表示数独的解。

例如：
输入（即图中题目）：
005300000
800000020
070010500
400005300
010070006
003200080
060500009
004000030
000009700

程序应该输出：
145327698
839654127
672918543
496185372
218473956
753296481
367542819
984761235
521839764

再例如，输入：
800000000
003600000
070090200
050007000
000045700
000100030
001000068
008500010
090000400

程序应该输出：
812753649
943682175
675491283
154237896
369845721
287169534
521974368
438526917
796318452

---

参考程序

#include<cstdio>
#include<cstring> 
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
#define set(a,b) memset(a,b,sizeof(a))
typedef long long ll;
#define N 9
int a[N][N],vis[N+1];

bool yanzheng(int x,int y){
    int i,j,yan[N+1]={0};
    for(i=0;i<N;i++){
        if(a[x][i]!=0&&yan[a[x][i]]){
            return false;
        }else{
            yan[a[x][i]]++;
        }
    }
    set(yan,0);
    for(i=0;i<N;i++){
        if(a[i][y]&&yan[a[i][y]]){
            return false;
        }else{
            yan[a[i][y]]++;
        }
    }
    set(yan,0);
    for(i=(x/3)*3;i<(x/3+1)*3 ;i++){
        for(j=(y/3)*3;j<(y/3+1)*3 ;j++){
            if(a[i][j]&&yan[a[i][j]]){
                return false;
            }else{
                yan[a[i][j]]++;
            }
        }
    }
    return true;
}
void dfs(int step){
    int i,j;
    if(step==81){
        for(i=0;i<N;i++){
            for(j=0;j<N;j++){
                cout<<a[i][j];
            }
            cout<<endl;
        }
        return ;
    }
    int x=step/N,y=step%N;
    if(a[x][y]){
        dfs(step+1);
    }
    else {
        for(int i=1;i<=N;i++){
            if(vis[i]<9){
                a[x][y]=i;
                if(yanzheng(x,y)){
                    vis[i]++;
                    dfs(step+1);
                    vis[i]--;
                    a[x][y]=0;
                }else{
                    a[x][y]=0;
                }
            }
        }
    }

}

int main()
{
    int i,j;
    char c[N+1];
    set(a,0);
    set(vis,0);
    for(i=0;i<N;i++){
        cin>>c;
        for(j=0;j<N;j++){
            a[i][j]=c[j]-'0';
            vis[c[j]-'0']++;
        }
    }
    dfs(0);
    return 0;
}