hdu 4210 Su-domino-ku

题目连接：http://acm.hdu.edu.cn/showproblem.php?pid=4210

Dancing Links + dfs。

额，好像又有一段时间没更新博客了，最近比较忙，上周校赛终于圆满结束了，我这一颗心也终于能放下来了。

上周还有一些另外收获，先保密吧 ^_^

这周本来不准备去武汉参加比赛的，想在五一的时候好好出去玩玩，但由于种种原因吧， 老师还是决定让我们这队去。哎，压力又是不小啊。

这几天抓紧集训一下，希望能找回一些状态。

上周就看过这道题目，当时没仔细读，以为给出的骨牌没什么用，就是一个很裸的数独问题，今天写好了代码，发现结果和给出的样例不太一样。

又读了读题目，发现原来需要判断能否用骨牌拼接并且每种骨牌只能用一次。这就增加了题目难度，每找到一种合适的数独，就需要再判断是否满足

骨牌的要求，我是用dfs来判断是否满足要求的。

code:(好长时间没写Dancing Links，又查了好多次模板)

# include<stdio.h>
# include<string.h>
# include<stdlib.h>
# define RR 750
# define CC 350
# define V RR*CC
int U[V],D[V];
int L[V],R[V];
int C[V],ROW[V];
int H[RR],S[CC];
int size;
int map[10][10];
int hash1[RR],hash2[RR],hash[RR],OK[90];
int visit[12][12],vis[12][12];
int dir[2][2]={0,1,1,0};
void Link(int r,int c)
{
    S[c]++;C[size]=c;
    ROW[size]=r;
    U[size]=U[c];D[U[c]]=size;
    D[size]=c;U[c]=size;
    if(H[r]==-1) H[r]=L[size]=R[size]=size;
    else
    {
        L[size]=L[H[r]];R[L[H[r]]]=size;
        R[size]=H[r];L[H[r]]=size;
    }
    size++;
}
void remove1(int c)
{
    int i,j;
    L[R[c]]=L[c];
    R[L[c]]=R[c];
    for(i=D[c];i!=c;i=D[i])
    {
        for(j=R[i];j!=i;j=R[j])
        {
            S[C[j]]--;
            U[D[j]]=U[j];
            D[U[j]]=D[j];
        }
    }
}
void resume1(int c)
{
    int i,j;
    for(i=U[c];i!=c;i=U[i])
    {
        for(j=L[i];j!=i;j=L[j])
        {
            S[C[j]]++;
            U[D[j]]=D[U[j]]=j;
        }
    }
    R[L[c]]=L[R[c]]=c;
}
int dfs(int step)
{
    int i,j,k,temp,nextx,nexty;
    int ans1,ans2,ans;
    if(step==82) return 1;
    i=(step-1)/9+1;
    j=step-(i-1)*9;
    if(vis[i][j]) 
    {
        return dfs(step+1);
    }
    else
    {
        for(k=0;k<=1;k++)
        {
            nextx=i+dir[k][0];
            nexty=j+dir[k][1];
            if(nextx<=9 && nexty<=9 && vis[nextx][nexty]==0)
            {
                ans1=map[i][j];
                ans2=map[nextx][nexty];
                if(ans1>ans2)
                {
                    temp=ans1;
                    ans1=ans2;
                    ans2=temp;
                }
                if(visit[ans1][ans2]==0)
                {
                    vis[i][j]=1;
                    vis[nextx][nexty]=1;
                    visit[ans1][ans2]=1;
                    ans=dfs(step+1);
                    if(ans==1) return 1;
                    vis[i][j]=0;
                    vis[nextx][nexty]=0;
                    visit[ans1][ans2]=0;
                }
            }
        }
        return 0;
    }
}
int check()
{
    int i;
    for(i=0;i<81;i++)
        map[hash1[OK[i]]][hash2[OK[i]]]=hash[OK[i]];
    return dfs(1);
}
int Dance(int k)
{
    int i,j,Min,c;
    if(!R[0])
    {
        if(check()) return 1;
        return 0;
    }
    for(Min=RR,i=R[0];i;i=R[i])
        if(Min>S[i]) Min=S[i],c=i;
    remove1(c);
    for(i=D[c];i!=c;i=D[i])
    {
        for(j=R[i];j!=i;j=R[j])
            remove1(C[j]);
        OK[k]=ROW[i];
        if(Dance(k+1)) return 1;
        for(j=L[i];j!=i;j=L[j])
            resume1(C[j]);
    }
    resume1(c);
    return 0;
}
int main()
{
    int i,j,n,k,t=0,r;
    int ans1,ans2,ii,jj,temp;
    char ch1[5],ch2[5];
    while(scanf("%d",&n)!=EOF && n)
    {
        t++;
        for(i=0;i<=324;i++)
        {
            S[i]=0;
            D[i]=U[i]=i;
            L[i+1]=i;
            R[i]=i+1;
        }R[324]=0;
        size=325;
        r=0;
        memset(map,0,sizeof(map));
        memset(H,-1,sizeof(H));
        memset(visit,0,sizeof(visit));
        memset(vis,0,sizeof(vis));
        for(i=0;i<n;i++)
        {
            scanf("%d%s%d%s",&ans1,ch1,&ans2,ch2);
            ii=ch1[0]-'A'+1;
            jj=ch1[1]-'0';
            map[ii][jj]=ans1;
            vis[ii][jj]=1;

            ii=ch2[0]-'A'+1;
            jj=ch2[1]-'0';
            map[ii][jj]=ans2;
            vis[ii][jj]=1;

            if(ans1>ans2)
            {
                temp=ans1;
                ans1=ans2;
                ans2=temp;
            }
            visit[ans1][ans2]=1;
        }
        for(i=1;i<=9;i++)
        {
            scanf("%s",ch1);
            ii=ch1[0]-'A'+1;
            jj=ch1[1]-'0';
            map[ii][jj]=i;
            vis[ii][jj]=1;
        }
        for(i=1;i<=9;i++)
        {
            for(j=1;j<=9;j++)
            {
                if(map[i][j]!=0)
                {
                    r++;
                    hash1[r]=i;
                    hash2[r]=j;
                    hash[r]=map[i][j];
                    Link(r,(i-1)*9+j);
                    Link(r,81+(i-1)*9+map[i][j]);
                    Link(r,162+(j-1)*9+map[i][j]);
                    Link(r,243+((i-1)/3)*27 + ((j-1)/3)*9+map[i][j]);
                }
                else
                {
                    for(k=1;k<=9;k++)
                    {
                        r++;
                        hash1[r]=i;
                        hash2[r]=j;
                        hash[r]=k;
                        Link(r,(i-1)*9+j);
                        Link(r,81+(i-1)*9+k);
                        Link(r,162+(j-1)*9+k);
                        Link(r,243+(i-1)/3*27 + (j-1)/3*9 + k);
                    }
                }
            }
        }///
        Dance(0);
        /*for(i=0;i<81;i++)
        {
            map[hash1[OK[i]]][hash2[OK[i]]]=hash[OK[i]];
        }*/
        printf("Puzzle %d\n",t);
        for(i=1;i<=9;i++)
        {
            for(j=1;j<=9;j++)
                printf("%d",map[i][j]);
            printf("\n");
        }
    }
    return 0;
}