dance links X

dlx：主要解决精确覆盖，多重覆盖问题。基于十字链表的快速删除与恢复。

学习博客：https://www.cnblogs.com/grenet/p/3145800.html

P4929 【模板】舞蹈链（DLX）

建链表时行的顺序是颠倒的（u,d每次更新的表头），不过并不影响

#include<bits/stdc++.h>
#define N 100100
using namespace std;
int n,m,ans[N];
int d[N],u[N],L[N],R[N],s[N],cnt=0,h[N],row[N],col[N];
void init(){
	for(int i=0;i<=m;i++){
		L[i]=i-1;R[i]=i+1;
		u[i]=d[i]=i;
	}
	L[0]=m;R[m]=0;
	memset(h,-1,sizeof(h));
	memset(s,0,sizeof(s));
	cnt=m+1;
}
void push(int r,int c){
	s[c]++;
	row[cnt]=r;col[cnt]=c;
	u[cnt]=c;
	d[cnt]=d[c];
	u[d[c]]=cnt;
	d[c]=cnt;
	if(h[r]==-1)h[r]=R[cnt]=L[cnt]=cnt;
	else{
		R[cnt]=h[r];
		L[cnt]=L[h[r]];
		R[L[h[r]]]=cnt;
		L[h[r]]=cnt;
	}
	cnt++;
}
void remove(int c){
	R[L[c]]=R[c];L[R[c]]=L[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[col[j]]--;
		}
	}
}
void resume(int c){
	for(int i=u[c];i!=c;i=u[i]){
		for(int j=L[i];j!=i;j=L[j]){
			u[d[j]]=j;
			d[u[j]]=j;
			s[col[j]]++;
		}
	}
	R[L[c]]=c;
	L[R[c]]=c;
}
int dance(int dep){
	if(R[0]==0){
		for(int i=1;i<dep;i++)printf("%d ",ans[i]);
		return 1;
	}
	int c=R[0];
	for(int i=R[0];i!=0;i=R[i])
	if(s[i]<s[c])c=i;
	remove(c);
	for(int i=d[c];i!=c;i=d[i]){
		ans[dep]=row[i];
		for(int j=R[i];j!=i;j=R[j])remove(col[j]);
		if(dance(dep+1))return 1;
		for(int j=L[i];j!=i;j=L[j])resume(col[j]);
	}
	resume(c);
	return 0;
}
int main(){
	scanf("%d%d",&n,&m);
	init();
	int x;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++){
		scanf("%d",&x);
		if(x)push(i,j);
	}
	if(!dance(1))printf("No Solution!");
	return 0;
}

Radar HDU - 2295 （二分+重复覆盖）

二分答案，然后用重复覆盖判断，重复覆盖要注意剪枝

get()是剩余矩阵最少需要的次数

#include<bits/stdc++.h>
#define N 100100
using namespace std;
int n,m,ans[N];
int d[N],u[N],L[N],R[N],s[N],cnt=0,h[N],row[N],col[N],k;
double eps=1e-8;
struct nd{
	double x,y;
}a[N],b[N];
void init(){
	for(int i=0;i<=m;i++){
		L[i]=i-1;R[i]=i+1;
		u[i]=d[i]=i;
	}
	L[0]=m;R[m]=0;
	memset(h,-1,sizeof(h));
	memset(s,0,sizeof(s));
	cnt=m+1;
}
void push(int r,int c){
	s[c]++;
	row[cnt]=r;col[cnt]=c;
	u[cnt]=c;
	d[cnt]=d[c];
	u[d[c]]=cnt;
	d[c]=cnt;
	if(h[r]==-1)h[r]=R[cnt]=L[cnt]=cnt;
	else{
		R[cnt]=h[r];
		L[cnt]=L[h[r]];
		R[L[h[r]]]=cnt;
		L[h[r]]=cnt;
	}
	cnt++;
}
void remove(int c){
	for(int i=d[c];i!=c;i=d[i]){
		R[L[i]]=R[i];
		L[R[i]]=L[i];
	}
}
void resume(int c){
	for(int i=u[c];i!=c;i=u[i]){
		R[L[i]]=i;
		L[R[i]]=i;
	}
}
int vis[55];
int get(){
	memset(vis,0,sizeof(vis));
	int ret=0;
	for(int i=R[0];i!=0;i=R[i]){
		if(!vis[i]){
			ret++;
			for(int j=d[i];j!=i;j=d[j])
			for(int l=R[j];l!=j;l=R[l])
			vis[col[l]]=1;
		}
	}
	return ret;
}
int dance(int dep){
	if(R[0]==0){
		return dep-1<=k;
	}
	if(get()+dep-1>k)return 0;
	int c=R[0];
	for(int i=R[0];i!=0;i=R[i])
	if(s[i]<s[c])c=i;
	for(int i=d[c];i!=c;i=d[i]){
		remove(i);
		for(int j=R[i];j!=i;j=R[j]){
			remove(j);
		}
		if(dance(dep+1))return 1;
		for(int j=L[i];j!=i;j=L[j])resume(j);
		resume(i);
	}
	return 0;
}
double dist(double x,double y,double xx,double yy){
	return sqrt((xx-x)*(xx-x)+(yy-y)*(yy-y));
}
int chk(double mid){
	init();
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			if(mid-dist(a[j].x,a[j].y,b[i].x,b[i].y)>eps)push(i,j);
		}
	}
	if(dance(1))return 1;
	return 0;
}
int main(){
	int qt;
	scanf("%d",&qt);
	while(qt--){
		scanf("%d%d%d",&m,&n,&k);
		for(int i=1;i<=m;i++)scanf("%lf%lf",&a[i].x,&a[i].y);
		for(int i=1;i<=n;i++)scanf("%lf%lf",&b[i].x,&b[i].y);
		double l=0,r=1e4,ans;
		while(r-l>eps){
			double mid=(l+r)/2.0;
			if(chk(mid))r=mid,ans=mid;
			else l=mid;
		}
		printf("%.6f\n",ans);
	}
	return 0;
}

Sudoku POJ - 3074（数独问题）

分为四种限制

1.所有格子必须有且仅有一个数，1--81列

2.每行必须有1到9，82--81*2列

3.每列必须有1到9，81*2+1--81*3列

4.每宫必须有1到9，81*3+1--81*4列

对于已经有数的格子，在其对应的列上复制为一。

对于没有数的格子，从1到9 每个数都新加一行

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 11000
using namespace std;
int n,m,cnt,f[N][4],L[N],R[N],s[N],u[N],d[N],row[N],col[N],h[N];
int ans[N],a[100][100];
char str[N];
void init(){
	for(int i=0;i<=m;i++){
		L[i]=i-1;R[i]=i+1;
		u[i]=d[i]=i;
	}
	L[0]=m;R[m]=0;
	memset(h,-1,sizeof(h));
	memset(s,0,sizeof(s));
	cnt=m+1;
//	cout<<cnt<<endl;
}
void push(int r,int c,int x,int y,int z){
	f[r][0]=x;f[r][1]=y;f[r][2]=z;
	s[c]++;
	row[cnt]=r;col[cnt]=c;
	u[cnt]=c;
	d[cnt]=d[c];
	u[d[c]]=cnt;
	d[c]=cnt;
	if(h[r]==-1)h[r]=R[cnt]=L[cnt]=cnt;
	else{
		R[cnt]=h[r];
		L[cnt]=L[h[r]];
		R[L[h[r]]]=cnt;
		L[h[r]]=cnt;
	}
	cnt++;
}
void remove(int c){//删除列及该列所涉及的行 
	R[L[c]]=R[c];L[R[c]]=L[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[col[j]]--;
		}
	}
}
void resume(int c){
	for(int i=u[c];i!=c;i=u[i]){
		for(int j=L[i];j!=i;j=L[j]){
			u[d[j]]=j;
			d[u[j]]=j;
			s[col[j]]++;
		}
	}
	R[L[c]]=c;
	L[R[c]]=c;
}
int dance(int dep){
	if(R[0]==0){
		return 1;
	}
	//if(dep>81)return 0;
	//cout<<dep<<endl;
	int c=R[0];//cout<<c<<endl;
	for(int i=R[0];i!=0;i=R[i])
	if(s[i]<s[c])c=i;
	remove(c);
	
	for(int i=d[c];i!=c;i=d[i]){
		ans[dep]=row[i];
		for(int j=R[i];j!=i;j=R[j])remove(col[j]);
		if(dance(dep+1))return 1;
		for(int j=L[i];j!=i;j=L[j])resume(col[j]);
	}
	resume(c);
	return 0;
}
int main(){
	int t,t1,t2,t3;
	while(scanf("%s",str)){
		if(str[0]=='e')break;
		memset(f,0,sizeof(f));
		m=81*4;n=0; 
		init();
		for(int i=0;i<9;i++)
		for(int j=0;j<9;j++){
			char ch=str[i*9+j];
			if(ch=='.'){
				for(int k=1;k<=9;k++){
					++n;
					t=i*9+j+1;
					t1=i*9+k+81;
					t2=j*9+k+162;
					t3=((i/3)*3+(j/3))*9+k+243;
					push(n,t,i,j,k);
					push(n,t1,i,j,k);
					push(n,t2,i,j,k);
					push(n,t3,i,j,k);
				}
			}
			else{
				++n;
				int x=ch-'0';
				t=i*9+j+1;
				t1=i*9+x+81;
				t2=j*9+x+162;
				t3=((i/3)*3+(j/3))*9+x+243;
				push(n,t,i,j,x);
				push(n,t1,i,j,x);
				push(n,t2,i,j,x);
				push(n,t3,i,j,x);
			}
			//cout<<R[0]<<endl;
		}
	//	cout<<str<<endl;
	//	cout<<"KKK"<<endl;
		dance(1);
	//	cout<<"LLL"<<endl;
		for(int i=1;i<=81;i++){
			int t=ans[i];
			a[f[t][0]][f[t][1]]=f[t][2];
		}
		for(int i=0;i<9;i++)
		for(int j=0;j<9;j++)
		printf("%d",a[i][j]);
		printf("\n");
	}
	return 0;
}

Squiggly Sudoku HDU - 4069 （数独问题）

与普通数独类似，只是宫并不是规整的，先预处理出九个宫，然后与普通数独一样

注意：在统计出第一个答案时要及时复制下来，因为可能会被下面的解破坏

Divisibility HDU - 3335 （dlx算法流程的理解）

题目是选取一个数然后删去一些数（冲突，不能再选），dlx多重覆中是先从一列中选取一行，然后在把该行中涉及到的列删去，正好与题目要求相符，每次选的行所对应的数就是最终集合中算选的数。

有冲突的两个数在矩阵中为一，建图然后多重覆盖

get()是剩余矩阵最多需要次数

#include<bits/stdc++.h>
#define N 10010
#define LL long long
using namespace std;
int n,m,ans;
int d[N],u[N],L[N],R[N],s[N],cnt=0,h[N],row[N],col[N];
LL a[N];
void init(){
	for(int i=0;i<=m;i++){
		L[i]=i-1;R[i]=i+1;
		u[i]=d[i]=i;
	}
	L[0]=m;R[m]=0;
	memset(h,-1,sizeof(h));
	memset(s,0,sizeof(s));
	cnt=m+1;
}
void push(int r,int c){
	s[c]++;
	row[cnt]=r;col[cnt]=c;
	u[cnt]=c;
	d[cnt]=d[c];
	u[d[c]]=cnt;
	d[c]=cnt;
	if(h[r]==-1)h[r]=R[cnt]=L[cnt]=cnt;
	else{
		R[cnt]=h[r];
		L[cnt]=L[h[r]];
		R[L[h[r]]]=cnt;
		L[h[r]]=cnt;
	}
	cnt++;
}
void remove(int c){//只删除列 
	for(int i=d[c];i!=c;i=d[i]){
		R[L[i]]=R[i];
		L[R[i]]=L[i];
	}
}
void resume(int c){
	for(int i=u[c];i!=c;i=u[i]){
		R[L[i]]=i;
		L[R[i]]=i;
	}
}
int vis[55];
int get(){
//	memset(vis,0,sizeof(vis));
	int ret=0;
	for(int i=R[0];i!=0;i=R[i]){
		ret++;
	}
	return ret;
}
void dance(int dep){
	if(R[0]==0){
		ans=max(ans,dep-1);
		return;
	}
	if(get()+dep-1<=ans)return;
	int c=R[0];
	for(int i=R[0];i!=0;i=R[i])
	if(s[i]<s[c])c=i;
	for(int i=d[c];i!=c;i=d[i]){
		remove(i);
		for(int j=R[i];j!=i;j=R[j])remove(j);
		dance(dep+1);
		for(int j=L[i];j!=i;j=L[j])resume(j);
		resume(i);
	}
}
int main(){
	int qt;
	scanf("%d",&qt);
	while(qt--){
		scanf("%d",&n);
		m=n;
		ans=0;
		init();
		for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
		for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		if(a[i]%a[j]==0||a[j]%a[i]==0)push(i,j);
	//	cout<<"KK"<<endl;
		dance(1);
		printf("%d\n",ans);
	}
	return 0;
}

A simple math problem. HDU - 4979

对于从n个数中取r个数的所有集合为列，从n个数中取m个作为行，包含对应的r集合为一，建图然后多重覆盖。

直接写会超时，直接打表