poj3074 Sudoku【DLX】【数独】【精确覆盖问题】

题目链接：http://poj.org/problem?id=3074
大佬题解链接：https://www.cnblogs.com/ylfdrib/archive/2010/10/06/1844559.html
AC代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<queue>
#include<map>
#include<set>
#include<cstring>
#define INF 0x3fffffff
using namespace std;
char a[100];
#define NN 330
#define MM 740
int M,N;
int cntc[NN];
int L[NN*MM],R[NN*MM],U[NN*MM],D[NN*MM],C[NN*MM];
int head;
int adj[MM][NN];
int ans[10][10];
void Remove(int c)
{
     L[R[c]]=L[c];
     R[L[c]]=R[c];
     int i,j;
     for(i=D[c];i!=c;i=D[i])
     {
          for(j=R[i];j!=i;j=R[j])
          {
              U[D[j]]=U[j];
              D[U[j]]=D[j];
              cntc[C[j]]--;
          }
     }
}
void resume(int c)
{
    L[R[c]]=c;
    R[L[c]]=c;
    int i,j;
    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;
              cntc[C[j]]++;
          }
     }
}
int dfs()
{
     if(R[head]==head)
          return 1;
     int minn=INF;
     int c,i,j;
     for(i=R[head];i!=head;i=R[i])
     {
          if(cntc[i]<minn)
          {
               c=i;
               minn=cntc[i];
          }
     }
     Remove(c);
     for(i=D[c];i!=c;i=D[i])
     {
          int r=(i-1)/N;
          int num=(r+8)/9;
          int key=r%9;
          if(key==0)
                 key=9;
          int x=(num+8)/9;
          int y=num%9;
          if(y==0)
                    y=9;
          ans[x][y]=key;
          for(j=R[i];j!=i;j=R[j])
              Remove(C[j]);
          if(dfs())
                    return 1;
          for(j=L[i];j!=i;j=L[j])
              resume(C[j]);
     }
     resume(c);
     return 0;
}
int build()
{
     int i,j,now,pre,first;
     head=0;
     for(j=head;j<N;j++)
     {
          R[j]=j+1;
          L[j+1]=j;
     }
     L[head]=j;
     R[j]=head;
     for(j=1;j<=N;j++)
     {
          pre=j;
          cntc[j]=0;
          for(i =1;i<=M;i++)
          {
                if(adj[i][j])
                {
                    now=i*N+j;
                    C[now]=j;
                    cntc[j]++;
                    D[pre]=now;
                    U[now]=pre;
                    pre=now;
                }
          }
          now=j;
          D[pre]=now;
          U[now]=pre;
          if(cntc[j]==0)
                    return 0;
     }
     for(i=1;i<=M;i++)
     {
          pre=first=-1;
          for(j=1;j<=N;j++)
          {
              if(adj[i][j])
              {
                    now=i*N+j;
                    if(pre!=-1)
                    {
                        R[pre]=now;
                        L[now]=pre;
                    }
                    else
                              first=now;
                    pre=now;
              }
          }
          if(first!=-1)
          {
               now=first;
               R[pre]=now;
               L[now]=pre;
          }
     }
     return 1;
}
int main()
{
     int i,j,k;
     while(~scanf("%s",&a))
     {
             if(strcmp(a,"end")==0)
                    break;
             memset(adj,0,sizeof(adj));
             for(int i=1;i<=9;i++)
             {
                    for(j=1;j<=9;j++)
                    {
                         int t=9*(i-1)+j;
                         if(a[t-1]=='.')
                         {
                              for(k=1;k<=9;k++)
                              {
                                    adj[9*(t-1)+k][t]=1;
                                    adj[9*(t-1)+k][81+(i-1)*9+k]=1;
                                    adj[9*(t-1)+k][162+(j-1)*9+k]=1;
                                    adj[9*(t-1)+k][243+((i-1)/3*3+(j+2)/3-1)*9+k]=1;
                              }
                         }
                         else
                         {
                              k=a[t-1]-'0';
                              adj[9*(t-1)+k][t]=1;
                              adj[9*(t-1)+k][81+(i-1)*9+k]=1;
                              adj[9*(t-1)+k][162+(j-1)*9+k]=1;
                              adj[9*(t-1)+k][243+((i-1)/3*3+(j+2)/3-1)*9+k]=1;
                         }
                    }
             }
             M=729;
             N=324;
             build();
             dfs();
             for(i=1;i<=9;i++)
             {
                    for(j=1;j<=9;j++)
                    {
                         cout<<ans[i][j];
                    }
             }
             cout<<endl;
     }
     return 0;
}

我一定可以的！！！