hdu 3111 DLX解数独

思路：裸的DLX解数独。关键是建图，感觉还不如写个dfs直接，DLX写这个的代码很烦。

#include<set>

#include<map>

#include<cmath>

#include<queue>

#include<cstdio>

#include<vector>

#include<string>

#include<cstdlib>

#include<cstring>

#include<iostream>

#include<algorithm>

#define pb push_back

#define mp make_pair

#define Maxn 400010

#define Maxm 200010

#define LL __int64

#define Abs(x) ((x)>0?(x):(-x))

#define lson(x) (x<<1)

#define rson(x) (x<<1|1)

#define inf 100000

#define lowbit(x) (x&(-x))

#define clr(x,y) memset(x,y,sizeof(x))

#define Mod 1000000007

using namespace std;

int L[Maxn],R[Maxn],D[Maxn],U[Maxn],S[Maxn],C[Maxn],X[Maxn],Q[Maxn],H[1010],id;

int g[10][10];

void init(int m)

{

    int i;

    for(i=0;i<=m;i++){

        D[i]=U[i]=i;

        L[i+1]=i;

        R[i]=i+1;

        S[i]=0;

    }

    R[m]=0;

    id=m+1;

    clr(H,-1);

}

void ins(int r,int c)

{

    D[id]=D[c];

    U[id]=c;

    U[D[c]]=id;

    D[c]=id;

    S[c]++;

    if(H[r]<0)

        H[r]=L[id]=R[id]=id;

    else{

        L[id]=H[r];

        R[id]=R[H[r]];

        L[R[H[r]]]=id;

        R[H[r]]=id;

    }

    C[id]=c;

    X[id++]=r;

}

void Remove(int c)

{

    int i,j;

    R[L[c]]=R[c];

    L[R[c]]=L[c];

    for(i=D[c];i!=c;i=D[i]){

        for(j=R[i];j!=i;j=R[j]){

            D[U[j]]=D[j];

            U[D[j]]=U[j];

            S[C[j]]--;

        }

    }

}

void Resume(int c)

{

    int i,j;

    R[L[c]]=c;

    L[R[c]]=c;

    for(i=D[c];i!=c;i=D[i]){

        for(j=R[i];j!=i;j=R[j]){

            U[D[j]]=j;

            D[U[j]]=j;

            S[C[j]]++;

        }

    }

}

bool dfs(int k)

{

    int i,j,c,temp;

    if(R[0]==0){

        for(i=0;i<k;i++){

            int r,k;

            temp=X[Q[i]];

            k=temp%9;

            if(!k) k=9;

            c=(temp%81)/9;

            if(temp%81==0) c=9;

            if(temp%81%9) c++;

            r=temp/81;

            if(temp%81) r++;

            g[r][c]=k;

        }

        return true;

    }

    temp=inf;

    for(i=R[0];i;i=R[i]){

        if(S[i]<temp){

            temp=S[i];

            c=i;

        }

    }

    Remove(c);

    for(i=D[c];i!=c;i=D[i]){

        Q[k]=i;

        for(j=R[i];j!=i;j=R[j])

            Remove(C[j]);

        if(dfs(k+1))

            return true;

        for(j=L[i];j!=i;j=L[j])

            Resume(C[j]);

    }



    Resume(c);

    return false;

}

void build()

{

    int i,j,k,b,r,c;

    init(324);

    for(i=1;i<=9;i++){

        for(j=1;j<=9;j++){

            if(g[i][j]){

                r=(i-1)*81+(j-1)*9+g[i][j];

                c=(i-1)*9+g[i][j];

                ins(r,c);

                c=81+(j-1)*9+g[i][j];

                ins(r,c);

                c=162+(i-1)*9+j;

                ins(r,c);

                c=81*3+(((i-1)/3)*3+(j+2)/3-1)*9+g[i][j];

                ins(r,c);

                continue;

            }

            for(k=1;k<=9;k++){

                r=(i-1)*81+(j-1)*9+k;

                c=(i-1)*9+k;

                ins(r,c);

                c=81+(j-1)*9+k;

                ins(r,c);

                c=162+(i-1)*9+j;

                ins(r,c);

                c=81*3+(((i-1)/3)*3+(j+2)/3-1)*9+k;

                ins(r,c);

            }

        }

    }

}

int main()

{

    int t,i,j,f=0;

    char str[20];

    scanf("%d",&t);

    while(t--){

        clr(g,0);

        for(i=1;i<=9;i++){

            scanf("%s",str);

            for(j=0;j<9;j++){

                if(str[j]!='?')

                    g[i][j+1]=str[j]-'0';

            }

        }

        if(t)

        scanf("%s",str);

        build();

        if(f)

            printf("---\n");

        f=1;

        if(!dfs(0)){

            printf("impossible\n");

            continue;

        }

        for(i=1;i<=9;i++){

            for(j=1;j<=9;j++){

                printf("%d",g[i][j]);

            }

            printf("\n");

        }

    }

    return 0;

}