[DLX]HDOJ4069 Squiggly Sudoku

题意：有9*9的格子

每个格子 由五部分组成：上(16)、右(32)、下(64)、左(128)、和该格的数值(0~9)  

                                若上下左右有分割格子的线 就加上相应的数， 该格的数值若为0，则是未知  1~9 则是已知

然后根据分割线 做数独(每行、每列、每宫都是1~9)

输出无解、多解或者一个解就输出那个解

这种数独与普通3*3的数独 的唯一区别就是 宫 的划分的方式不一样

16、32、64、128这几个数很特殊 分别是$2^4$、$2^5$、$2^6$、$2^7$  

也就是相应的二位数，若第4位为1，则上面有分割线

　　　　　　　　　　若第5位为1，则右边有分割线

　　　　　　　　　　若第6位为1，则下面有分割线

　　　　　　　　　　若第7位为1，则左边有分割线

那么只要随便dfs一下确定在哪个宫里就好啦

然后把模板里  本来宫的位置是(i/3)*3+(j/3)   换成   dfs搜出来的宫的位置就好啦~

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
//const int N=2e5+5;

const int N=9; //3*3数独
const int MaxN=N*N*N+10;    // 一格能填9个数  9*9格
const int MaxM=N*N*4+10;     // 9*9*4=(9+9+9)*9+9*9    (9+9+9)是9行 9列 9格 *9是9个数 9*9是81个格子
const int maxnode=MaxN*4+MaxM+10;
int g[MaxN];
int anss;
struct DLX
{
    int n, m, size;
    int U[maxnode], D[maxnode], R[maxnode], L[maxnode], Row[maxnode], Col[maxnode];
    int H[MaxN], S[MaxM];  // S: 各列节点数
    int ansd, ans[MaxN];
    void init(int _n, int _m)
    {
        n=_n;
        m=_m;
        for(int i=0; i<=m; i++)
        {
            S[i]=0;         //每一列元素个数
            U[i]=D[i]=i;//上下指针
            L[i]=i-1;      //←
            R[i]=i+1;     //→
        }
        R[m]=0;          //循环 最后一个指向第一个
        L[0]=m;           //第一个往前指向最后一个
        size=m;           // 节点总数
        for(int i=1; i<=n; i++)
            H[i]=-1;     //头指针
    }
    void Link(int r, int c)
    {
        S[Col[++size]=c]++;
        Row[size]=r;
        D[size]=D[c];
        U[D[c]]=size;
        U[size]=c;
        D[c]=size;
        if(H[r]<0)
            H[r]=L[size]=R[size]=size;
        else
        {
            R[size]=R[H[r]];
            L[R[H[r]]]=size;
            L[size]=H[r];
            R[H[r]]=size;
        }
    }
    void remove(int c)
    {
        L[R[c]]=L[c];
        R[L[c]]=R[c];
        for(int i=D[c]; i!=c; 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 resume(int c)
    {
        for(int i=U[c]; i!=c; i=U[i])
            for(int j=L[i]; j!=i; j=L[j])
                S[Col[U[D[j]]=D[U[j]]=j]]++;
        L[R[c]]=R[L[c]]=c;
    }
    void Dance(int d)
    {
        if(anss>1)
            return ;
        if(R[0]==0)
        {
            for(int i=0;i<d;i++)
                g[(ans[i]-1)/N]=(ans[i]-1)%N+1;
            anss++;
        }
        int c=R[0];
        for(int i=R[0]; i!=0; i=R[i])
            if(S[i]<S[c])
                c=i;
        remove(c);
        for(int i=D[c]; i!=c; i=D[i])
        {
            ans[d]=Row[i];
            for(int j=R[i]; j!=i; j=R[j])
                remove(Col[j]);
            Dance(d+1);
            for(int j=L[i]; j!=i; j=L[j])
                resume(Col[j]);
        }
        resume(c);
    }
} dlx;
int s[10][10];
int vis[10][10];
void dfs(int x,int y,int col)
{
    if(x<0||x>=9||y<0||y>=9||vis[x][y]!=-1)
        return ;
    if((s[x][y]&(1<<7))==0)
    {
        vis[x][y]=col;
        dfs(x,y-1, col);
    }
    if((s[x][y]&(1<<6))==0)
    {
        vis[x][y]=col;
        dfs(x+1,y, col);
    }
    if((s[x][y]&(1<<5))==0)
    {
        vis[x][y]=col;
        dfs(x,y+1, col);
    }
    if((s[x][y]&(1<<4))==0)
    {
        vis[x][y]=col;
        dfs(x-1,y, col);
    }
}

void palce(int &r, int &c1, int &c2, int &c3, int &c4, int i, int j, int k)
{
    r=(i*N+j)*N+k;  // 第几行
    c1=i*N+j+1;     //  第几个格子
    c2=N*N+i*N+k;   // 第i行上的k
    c3=N*N*2+j*N+k; // 第j列上的k
    c4=N*N*3+(vis[i][j])*N+k; // 某宫中的k;
}

int main()
{
    int t, ca=1;
    scanf("%d", &t);
    while(t--)
    {
        for(int i=0;i<N;i++)
            for(int j=0;j<N;j++)
                scanf("%d", &s[i][j]);
        int num=0;
        memset(vis, -1, sizeof(vis));
        for(int i=0;i<N;i++)
            for(int j=0;j<N;j++)
            {
                if(vis[i][j]==-1)
                    dfs(i, j, num++);
                if((s[i][j]-(1<<7))>=0)
                    s[i][j]-=1<<7;
                if((s[i][j]-(1<<6))>=0)
                    s[i][j]-=1<<6;
                if((s[i][j]-(1<<5))>=0)
                    s[i][j]-=1<<5;
                if((s[i][j]-(1<<4))>=0)
                    s[i][j]-=1<<4;
            }
        dlx.init(N*N*N, 4*N*N);
        for(int i=0; i<N; i++)
            for(int j=0; j<N; j++)
                for(int k=1; k<=9; k++)
                    if(s[i][j]==0 || s[i][j]==k)
                    {
                        int r, c1, c2, c3, c4;
                        palce(r, c1, c2, c3, c4, i, j, k);
                        dlx.Link(r, c1);
                        dlx.Link(r, c2);
                        dlx.Link(r, c3);
                        dlx.Link(r, c4);
                    }
        anss=0;
        dlx.Dance(0);
        printf("Case %d:
", ca++);
        if(anss==0)
            puts("No solution");
        else if(anss>1)
            puts("Multiple Solutions");
        else
        {
            for(int i=0;i<N;i++)
            {
                for(int j=0;j<N;j++)
                    printf("%d", g[i*N+j]);
                puts("");
            }
        }
    }
    return 0;
}