POJ 3074&&3076 舞蹈链
点击打开链接
题意：完成数独
思路：直接用舞蹈链的不可重复覆盖直接就可以过了，建图时列为81+81+81+81，第一个为81个空都要有元素，然后第二个是9行，有9个数要添，第三个就是9列有9个元素要添，最后一个是9块有9个元素要添，然后建完模版OK
#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=1010;
int L[maxn*400],R[maxn*400],U[maxn*400],D[maxn*400];//节点的上下左右四个方向的链表
int C[maxn*400],H[maxn],cnt[maxn],ans[maxn],Row[maxn*400];//C列H行cnt列链表中元素个数
int n,m,id,len;
void init(int lll){
for(int i=0;i<=lll;i++){
cnt[i]=0;U[i]=D[i]=i;
L[i+1]=i;R[i]=i+1;
}
R[lll]=0;id=lll+1;
memset(H,-1,sizeof(H));
}
void Link(int r,int c){
cnt[c]++;C[id]=c;Row[id]=r;
U[id]=U[c];D[U[c]]=id;
D[id]=c;U[c]=id;
if(H[r]==-1) H[r]=L[id]=R[id]=id;
else{
L[id]=L[H[r]];R[L[H[r]]]=id;
R[id]=H[r];L[H[r]]=id;
}
id++;
}
void Remove(int Size){
L[R[Size]]=L[Size];
R[L[Size]]=R[Size];
for(int i=D[Size];i!=Size;i=D[i]){
for(int j=R[i];j!=i;j=R[j]){
U[D[j]]=U[j];D[U[j]]=D[j];
cnt[C[j]]--;
}
}
}
void Resume(int Size){
for(int i=D[Size];i!=Size;i=D[i]){
for(int j=R[i];j!=i;j=R[j]){
U[D[j]]=j;D[U[j]]=j;
cnt[C[j]]++;
}
}
L[R[Size]]=Size;R[L[Size]]=Size;
}
int Dance(int k){
int c;
if(!R[0]){
len=k;
return 1;
}
for(int Min=inf,i=R[0];i;i=R[i])
if(Min>cnt[i]) Min=cnt[i],c=i;
Remove(c);
for(int i=D[c];i!=c;i=D[i]){
ans[k]=Row[i];
for(int j=R[i];j!=i;j=R[j]) Remove(C[j]);
if(Dance(k+1)) return 1;
for(int j=L[i];j!=i;j=L[j]) Resume(C[j]);
}
Resume(c);
return 0;
}
char str[110];
int X[maxn],Y[maxn],VAL[maxn];
int astr[110];
int main(){
while(scanf("%s",str)!=-1){
if(strcmp(str,"end")==0) break;
init(324);
int kkk=1;
for(int i=1;i<=9;i++){
for(int j=1;j<=9;j++){
if(str[(i-1)*9+j-1]=='.'){
int pos=(i-1)/3*3+j/3;
if(j%3!=0) pos++;
for(int l=1;l<=9;l++){
X[kkk]=i;Y[kkk]=j;VAL[kkk]=l;
Link(kkk,(i-1)*9+j);
Link(kkk,81+(i-1)*9+l);
Link(kkk,162+(j-1)*9+l);
Link(kkk++,243+(pos-1)*9+l);
}
}else{
int pos=(i-1)/3*3+j/3;
if(j%3!=0) pos++;
int l=str[(i-1)*9+j-1]-'0';
X[kkk]=i;Y[kkk]=j;VAL[kkk]=l;
Link(kkk,(i-1)*9+j);
Link(kkk,81+(i-1)*9+l);
Link(kkk,162+(j-1)*9+l);
Link(kkk++,243+(pos-1)*9+l);
}
}
}
int fff=Dance(0);
for(int i=0;i<len;i++){
int lll=(X[ans[i]]-1)*9+Y[ans[i]]-1;
astr[lll]=VAL[ans[i]];
}
for(int i=0;i<81;i++) printf("%d",astr[i]);
printf("\n");
}
return 0;
} 点击打开链接
题意：与3074一样完成数独
思路：改一下数据范围什么的简单处理一下就可以了
#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=5010;
int L[maxn*1010],R[maxn*1010],U[maxn*1010],D[maxn*1010];//节点的上下左右四个方向的链表
int C[maxn*1010],H[maxn],cnt[maxn],ans[maxn],Row[maxn*1010];//C列H行cnt列链表中元素个数
int n,m,id,len;
void init(int lll){
for(int i=0;i<=lll;i++){
cnt[i]=0;U[i]=D[i]=i;
L[i+1]=i;R[i]=i+1;
}
R[lll]=0;id=lll+1;
memset(H,-1,sizeof(H));
}
void Link(int r,int c){
cnt[c]++;C[id]=c;Row[id]=r;
U[id]=U[c];D[U[c]]=id;
D[id]=c;U[c]=id;
if(H[r]==-1) H[r]=L[id]=R[id]=id;
else{
L[id]=L[H[r]];R[L[H[r]]]=id;
R[id]=H[r];L[H[r]]=id;
}
id++;
}
void Remove(int Size){
L[R[Size]]=L[Size];
R[L[Size]]=R[Size];
for(int i=D[Size];i!=Size;i=D[i]){
for(int j=R[i];j!=i;j=R[j]){
U[D[j]]=U[j];D[U[j]]=D[j];
cnt[C[j]]--;
}
}
}
void Resume(int Size){
for(int i=D[Size];i!=Size;i=D[i]){
for(int j=R[i];j!=i;j=R[j]){
U[D[j]]=j;D[U[j]]=j;
cnt[C[j]]++;
}
}
L[R[Size]]=Size;R[L[Size]]=Size;
}
int Dance(int k){
int c;
if(!R[0]){
len=k;
return 1;
}
for(int Min=inf,i=R[0];i;i=R[i])
if(Min>cnt[i]) Min=cnt[i],c=i;
Remove(c);
for(int i=D[c];i!=c;i=D[i]){
ans[k]=Row[i];
for(int j=R[i];j!=i;j=R[j]) Remove(C[j]);
if(Dance(k+1)) return 1;
for(int j=L[i];j!=i;j=L[j]) Resume(C[j]);
}
Resume(c);
return 0;
}
char str[20][20];
int X[maxn],Y[maxn],VAL[maxn];
int astr[110];
int main(){
while(scanf("%s",str[0])!=-1){
for(int i=1;i<16;i++) scanf("%s",str[i]);
init(1024);
int kkk=1;
for(int i=1;i<=16;i++){
for(int j=1;j<=16;j++){
if(str[i-1][j-1]=='-'){
int pos=(i-1)/4*4+j/4;
if(j%4!=0) pos++;
for(int l=1;l<=16;l++){
X[kkk]=i;Y[kkk]=j;VAL[kkk]=l;
Link(kkk,(i-1)*16+j);
Link(kkk,256+(i-1)*16+l);
Link(kkk,512+(j-1)*16+l);
Link(kkk++,768+(pos-1)*16+l);
}
}else{
int pos=(i-1)/4*4+j/4;
if(j%4!=0) pos++;
int l=str[i-1][j-1]-'A'+1;
X[kkk]=i;Y[kkk]=j;VAL[kkk]=l;
Link(kkk,(i-1)*16+j);
Link(kkk,256+(i-1)*16+l);
Link(kkk,512+(j-1)*16+l);
Link(kkk++,768+(pos-1)*16+l);
}
}
}
int fff=Dance(0);
for(int i=0;i<len;i++){
int lll=(X[ans[i]]-1)*16+Y[ans[i]]-1;
astr[lll]=VAL[ans[i]];
}
for(int i=0;i<256;){
int lllll=16;
while(lllll--){
char ch=astr[i++]+'A'-1;
printf("%c",ch);
}
printf("\n");
}
printf("\n");
}
return 0;
}