舞蹈鏈—— 16 × 16 數獨模板(ZOJ
阿新 • • 發佈:2018-12-21
原題的資料出錯了,下面給出正確的資料:
Sample Input
–A----C-----O-I
-J–A-B-P-CGF-H-
–D--F-I-E----P-
-G-EL-H----M-J–
----E----C–G—
-I–K-GA-B—E-J
D-GP–J-F----A–
-E—C-B–DP–O-
E–F-M–D--L-K-A
-C--------O-I-L-
H-P-C–F-A–B—
—G-OD—J----H
K—J----H-A-P-L
–B--P–E--K–A-
-H–B--K–FI-C–
–F—C–D--H-N-
Sample Output
FPAHMJECNLBDKOGI
OJMIANBDPKCGFLHE
LNDKGFOIJEAHMBPC
BGCELKHPOFIMAJDN
MFHBELPOACKJGNID
CILNKDGAHBMOPEFJ
DOGPIHJMFNLECAKB
JEKAFCNBGIDPLHOM
EBOFPMIJDGHLNKCA
NCJDHBAEKMOFIGLP
HMPLCGKFIAENBDJO
AKIGNODLBPJCEFMH
KDEMJIFNCHGAOPBL
GLBCDPMHEONKJIAF
PHNOBALKMJFIDCEG
IAFJOECGLDPBHMNK
本題出自 :
#include<bits/stdc++.h> using namespace std; typedef long long ll; const int MAXN = 16*16*16+10; const int MAXM = 16*16*4+10; const int MAXNODE = MAXN*4+MAXM; struct DancingLink { int n,s,ansd;//列數 節點總數 int S[MAXM],A[MAXN],H[MAXN];//S[]該列節點總數 A[]答案 H[]行首指標 int L[MAXNODE],R[MAXNODE],U[MAXNODE],D[MAXNODE]; //L[],R[],U[],D[] 上下左右 int X[MAXNODE],C[MAXNODE];//X[] C[] 行列編號 void init(int n) { //初始化 this->n=n; for(int i=0; i<=n; i++) U[i]=i,D[i]=i,L[i]=i-1,R[i]=i+1; R[n]=0,L[0]=n; s=n+1; memset(S,0,sizeof(S)); memset(H,-1,sizeof(H)); } void DelCol(int c) { //刪除列 L[R[c]]=L[c]; R[L[c]]=R[c]; for(int i=D[c]; i!=c; i=D[i]) for(int j=R[i]; j!=i; j=R[j]) U[D[j]]=U[j],D[U[j]]=D[j],--S[C[j]]; } void ResCol(int c) { //恢復列 for(int i=U[c]; i!=c; i=U[i]) for(int j=L[i]; j!=i; j=L[j]) ++S[C[j]],U[D[j]]=j,D[U[j]]=j; L[R[c]]=c,R[L[c]]=c; } void AddNode(int r,int c) { //新增節點 ++S[c],C[++s]=c,X[s]=r; D[s]=D[c],U[D[c]]=s,U[s]=c,D[c]=s; if(H[r]<0) H[r]=L[s]=R[s]=s;//行首節點 else R[s]=R[H[r]],L[R[H[r]]]=s,L[s]=H[r],R[H[r]]=s; } bool dfs(int d) { //深度,深搜遍歷 if(!R[0]) { ansd=d; return true; } int c=R[0]; for(int i=R[0]; i; i=R[i]) if(S[i]<S[c]) c=i; DelCol(c); for(int i=D[c]; i!=c; i=D[i]) { A[d]=X[i]; for(int j=R[i]; j!=i; j=R[j]) DelCol(C[j]); if(dfs(d+1)) return true; for(int j=L[i]; j!=i; j=L[j]) ResCol(C[j]); } ResCol(c); return false; } } dlx; int encode(int a,int b,int c) { return a*256+b*16+c+1; } void decode(int code,int &a,int &b,int &c) { code--; c=code%16; code/=16; b=code%16; code/=16; a=code; } int main() { char s[20][20]; bool first=true; while(scanf("%s",s[0])!=EOF) { if(first)first=false; else printf("\n"); for(int i=1; i<16; i++) scanf("%s",s[i]); dlx.init(16*16*4); for(int r=0; r<16; r++) for(int c=0; c<16; c++) for(int v=0; v<16; v++) if(s[r][c]=='-'||s[r][c]=='A'+v) { int row=encode(r,c,v); dlx.AddNode(row,encode(0,r,c)); dlx.AddNode(row,encode(1,r,v)); dlx.AddNode(row,encode(2,c,v)); dlx.AddNode(row,encode(3,(r/4)*4+c/4,v)); } dlx.dfs(0); for(int i=0; i<dlx.ansd; i++) { int r,c,v; decode(dlx.A[i],r,c,v); s[r][c]='A'+v; } for(int i=0; i<16; i++) printf("%s\n",s[i]); } }