题目链接: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;
}
我一定可以的!!!