源代碼-dlx算法問題，求解答，高分啊

dlx算法問題，求解答，高分啊

有個算法面試-面試官問了dancinglink-完全不知道這是什麼啊-後來知道可能是解精確覆蓋問題的-求大家發份代碼-並且能夠在站內解釋一下-才給這麼多分哦，求個數獨問題的dancinglink解法，有源代碼，線上解釋。

最佳回答：

數獨問題的dancinglink解法源代碼如下，有問題隨時聯系我。
dlx是從數據結構角度優化01矩陣精確覆蓋和重復覆蓋的數據結構，它用十字鏈表只存貯矩陣中的非0元，而01矩陣精確覆蓋dfs過程中矩陣會越來越稀疏而且每次恢復現場會浪費大量時間，dlx恰好能解決這兩個問題。本題關鍵是將數獨問題轉化為01矩陣精確覆蓋。數獨轉化為精確覆蓋問題的方法還是參照Knuth的論文，如果讀取到一個格子是空的，那麼加9行，分別表示這個格子填1到9這9個數字，如果讀取到的格子是一個數字，那麼就加一行就可以了，然後列有9*9*4列，前81列表示這一行表示填的是第i行第j列的格子，接下來81列表示第i行填寫k，接下來81列表示第j列填寫k，最後81列表示對應九宮格填寫k。
實例代碼：
//poj 3074
//sep9
#include
#include
#define INT_MAX 2147483647
using namespace std;
const int MAX=1024;
const int col_num=9*9*4;
const int head=0;
const int delta[]={1,82,163,244};
int cnt[MAX],st[MAX];
int left[MAX*MAX],right[MAX*MAX],up[MAX*MAX],down[MAX*MAX];
int row[MAX*MAX],col[MAX*MAX];

int K,M;//k:node's idx M:row's number

struct ANS
{
int r,c,k;
}ans[MAX*MAX];

void init()
{
left[head]=col_num;
right[head]=1;

up[head]=down[head]=head;
for(int i=1;i<=col_num;++i){
left[i]=i-1;
right[i]=(i+1)%(col_num+1);
up[i]=down[i]=i;
cnt[i]=0;
col[i]=i;
row[i]=0;
}
M=0;
K=col_num;

}

int make_col_head(int c)
{
++K;
++cnt[c];
col[K]=c;
row[K]=M;

left[K]=right[K]=K;

up[K]=c;
down[K]=down[c];
up[down[K]]=K;
down[up[K]]=K;
return K;

}

void addcol(int ids,int c)
{
++K;
++cnt[c];
col[K]=c;
row[K]=M;

left[K]=ids;
right[K]=right[ids];
left[right[K]]=K;
right[left[K]]=K;

up[K]=c;
down[K]=down[c];
up[down[K]]=K;
down[up[K]]=K;

}

void addrow(int i,int j,int k)
{
++M;
ans[M].r=i;
ans[M].c=j;
ans[M].k=k+1;

int ids=make_col_head(9*i+j+delta[0]);
addcol(ids,9*i+k+delta[1]);
addcol(ids,9*j+k+delta[2]);
addcol(ids,9*(i/3*3+j/3)+k+delta[3]);

}

void remove(int c)
{
left[right[c]]=left[c];
right[left[c]]=right[c];
for(int i=down[c];i!=c;i=down[i])
for(int j=right[i];j!=i;j=right[j]){
up[down[j]]=up[j];
down[up[j]]=down[j];
--cnt[col[j]];
}

}

void resume(int c)
{
for(int i=up[c];i!=c;i=up[i])
for(int j=left[i];j!=i;j=left[j]){
down[up[j]]=j;
up[down[j]]=j;
++cnt[col[j]];
}
left[right[c]]=c;
right[left[c]]=c;

}

bool dfs(int k)
{
if(right[head]==head){
char s[128];
for(int i=0;i<k;++i)
s[ans[st[i]].r*9+ans[st[i]].c]=ans[st[i]].k+'0';
s[81]='\0';
puts(s);
return true;
}

int s=INT_MAX,c=0;
for(int i=right[head];i!=head;i=right[i]){
if(cnt[i]<s){
s=cnt[i];
c=i;
}
}
remove(c);
for(int i=down[c];i!=c;i=down[i]){
st[k]=row[i];
for(int j=right[i];j!=i;j=right[j])
remove(col[j]);
if(dfs(k+1))
return true;
for(int j=left[i];j!=i;j=left[j])
resume(col[j]);
}
resume(c);
return false;
}

int main()
{
char s[128];
while(scanf("%s",s)==1&&s[0]!='e'){
init();
for(int i=0;i<9;++i)
for(int j=0;j<9;++j)
if(s[i*9+j]=='.'){
for(int k=0;k<9;++k)
addrow(i,j,k);
}
else
addrow(i,j,s[i*9+j]-'1');

dfs(0);
}
return 0;

}