1. 程式人生 > >[DLX+bfs] hdu 4069 Squiggly Sudoku

[DLX+bfs] hdu 4069 Squiggly Sudoku

tor ios mes data true struct sed amp 9*9

題意:

給你9*9的矩陣。對於每一個數字。能減16代表上面有墻,能減32代表以下有墻。

。。

最後剩下的數字是0代表這個位置數要求,不是0代表這個數已知了。

然後通過墻會被數字分成9塊。

然後做數獨,這裏的數獨不是分成9個3*3的小塊而是通過墻分成的。

思路:

首先通過數字作出墻。

然後bfs求連通塊。dfs也能夠。目的是分塊。

然後就是dlx數獨模板題了。

這裏要註意的是假設找到答案2次就說明有多組解了。就應該停止返回了。不然會TLE。

代碼:

#include"stdio.h"
#include"algorithm"
#include"string.h"
#include"iostream"
#include"cmath"
#include"queue"
#include"map"
#include"vector"
#include"string"
using namespace std;
#define RN 9*9*9+5
#define CN 4*9*9+5
#define N 9*9*9*4+5
int wall[12][12][12][12];
int mp[12][12],used[12][12];
int dis[4][2]= {{0,1},{0,-1},{-1,0},{1,0}};
int kx,ff;
template<class T>inline void getd(T &x)
{
    int ch = getchar();
    bool minus = false;
    while(!isdigit(ch) && ch != ‘-‘)ch = getchar();
    if(ch == ‘-‘)minus = true, ch = getchar();
    x = ch - ‘0‘;
    while(isdigit(ch = getchar()))x = x * 10 - ‘0‘ + ch;
    if(minus)x = -x;
}

struct node
{
    int x,y;
};
void bfs(int x,int y)
{
    node cur,next;
    cur.x=x;
    cur.y=y;
    queue<node>q;
    q.push(cur);
    used[cur.x][cur.y]=kx;
    while(!q.empty())
    {
        cur=q.front();
        q.pop();
        for(int i=0; i<4; i++)
        {
            next.x=cur.x+dis[i][0];
            next.y=cur.y+dis[i][1];
            if(used[next.x][next.y]!=0 || wall[cur.x][cur.y][next.x][next.y]==1 ) continue;
            used[next.x][next.y]=kx;
            q.push(next);
        }
    }
    kx++;
    return ;
}

struct DLX
{
    int n,m,C;
    int U[N],D[N],L[N],R[N],Row[N],Col[N];
    int H[RN],S[CN],cnt,ans[RN];
    void init(int _n,int _m)
    {
        n=_n;
        m=_m;
        for(int i=0; i<=m; i++)
        {
            U[i]=D[i]=i;
            L[i]=(i==0?

m:i-1); R[i]=(i==m?0:i+1); S[i]=0; } C=m; for(int i=1; i<=n; i++) H[i]=-1; } void link(int x,int y) { C++; Row[C]=x; Col[C]=y; S[y]++; U[C]=U[y]; D[C]=y; D[U[y]]=C; U[y]=C; if(H[x]==-1) H[x]=L[C]=R[C]=C; else { L[C]=L[H[x]]; R[C]=H[x]; R[L[H[x]]]=C; L[H[x]]=C; } } void del(int x) { R[L[x]]=R[x]; L[R[x]]=L[x]; for(int i=D[x]; i!=x; i=D[i]) { for(int j=R[i]; j!=i; j=R[j]) { U[D[j]]=U[j]; D[U[j]]=D[j]; S[Col[j]]--; } } } void rec(int x) { for(int i=U[x]; i!=x; i=U[i]) { for(int j=L[i]; j!=i; j=L[j]) { U[D[j]]=j; D[U[j]]=j; S[Col[j]]++; } } R[L[x]]=x; L[R[x]]=x; } void dance(int x) { if(R[0]==0) { ff++; if(ff>=2) return ; cnt=x; for(int i=0; i<cnt; i++) { int tep=ans[i]-1; int a=tep/81,b=(tep%81)/9; mp[a+1][b+1]=tep%9+1; } return ; } int now=R[0]; for(int i=R[0]; i!=0; i=R[i]) { if(S[i]<S[now]) now=i; } del(now); for(int i=D[now]; i!=now; i=D[i]) { ans[x]=Row[i]; for(int j=R[i]; j!=i; j=R[j]) del(Col[j]); dance(x+1); if(ff>=2) return ; for(int j=L[i]; j!=i; j=L[j]) rec(Col[j]); } rec(now); return ; } } dlx; void getplace(int i,int j,int k,int &x,int &a,int &b,int &c) { x=(i-1)*81+(j-1)*9+k; a=81+(i-1)*9+k; b=81*2+(j-1)*9+k; c=81*3+(used[i][j]-1)*9+k; } int main() { int t,cas=1; cin>>t; while(t--) { memset(wall,0,sizeof(wall)); for(int i=1; i<=9; i++) { for(int j=1; j<=9; j++) { int x; getd(x); if(x-128>=0) { x-=128; wall[i][j][i][j-1]=1; } if(x-64>=0) { x-=64; wall[i][j][i+1][j]=1; } if(x-32>=0) { x-=32; wall[i][j][i][j+1]=1; } if(x-16>=0) { x-=16; wall[i][j][i-1][j]=1; } mp[i][j]=x; } } kx=1; memset(used,0,sizeof(used)); for(int i=1; i<=9; i++) for(int j=1; j<=9; j++) if(used[i][j]==0) bfs(i,j); dlx.init(9*9*9,4*9*9); for(int i=1; i<=9; i++) { for(int j=1; j<=9; j++) { int tep=(i-1)*9+j; int x,a,b,c; if(mp[i][j]==0) { for(int k=1; k<=9; k++) { getplace(i,j,k,x,a,b,c); dlx.link(x,tep); dlx.link(x,a); dlx.link(x,b); dlx.link(x,c); } } else { getplace(i,j,mp[i][j],x,a,b,c); dlx.link(x,tep); dlx.link(x,a); dlx.link(x,b); dlx.link(x,c); } } } ff=0; dlx.dance(0); printf("Case %d:\n",cas++); if(ff==0) puts("No solution"); else if(ff==1) { for(int i=1; i<=9; i++) for(int j=1; j<=9; j++) printf(j==9?"%d\n":"%d",mp[i][j]); } else puts("Multiple Solutions"); } return 0; }




[DLX+bfs] hdu 4069 Squiggly Sudoku