【POJ3740】Easy Finding DLX（Dancing Links）精确覆盖问题

题意：多组数据。每组数据给你几行数，要求选出当中几行，使得每一列都有且仅有一个1，询问是可不可行，或者说能不能找出来。

题解：1、暴搜。2、DLX（Dancing links）。

本文写的是DLX。

算法參考白书P406或者http://www.cnblogs.com/grenet/p/3145800.html

我说一些仔细的东西，就是删除操作的形状是

    |

——|————

——|————

——|————

被删除的点们之间的联系不用删，能够保留。准确地说它并非删去了这些点，而是删去这个形。

并且恢复时要反着恢复。

首先先确定从哪一列删除。进行一次remove。然后枚举这一列的每一行，对其进行remove，然后dfs，然后再resume。跳出循环时再resume确定列。

附代码，应该非常好看，仅仅要略耐心 一点点点点点。。我保证绝对照网上其它人的可读性强。

哦，对了。并不须要写点、行、列的单独的remove和resume函数，那太傻了，写模板也不须要，我想这应该永远用不上的。

。。。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 400
#define NN 10000
using namespace std;
struct DLX
{
	int U[NN],D[NN],L[NN],R[NN],C[NN];
	int H[NN],T[NN],cnt;
	inline void init(int m)
	{
		cnt=0;
		memset(H,0,sizeof(H));
		for(cnt=0;cnt<=m;cnt++)
		{
			U[cnt]=D[cnt]=C[cnt]=cnt;
			L[cnt]=cnt-1,R[cnt]=cnt+1;
		}
		L[0]=--cnt,R[cnt]=0;
	}
	inline void newnode(int x,int y)
	{
		C[++cnt]=y;T[y]++;

		if(!H[x])H[x]=L[cnt]=R[cnt]=cnt;
		else L[cnt]=H[x],R[cnt]=R[H[x]];
		R[H[x]]=L[R[H[x]]]=cnt,H[x]=cnt;

		U[cnt]=U[y],D[cnt]=y;
		U[y]=D[U[y]]=cnt;
	}
	inline void scan(int n,int m)
	{
		int i,j,k;
		for(i=1;i<=n;i++)for(j=1;j<=m;j++)
		{
			scanf("%d",&k);
			if(k)newnode(i,j);
		}
	}
	inline void remove(int x)
	{
		for(int i=D[x];i!=x;i=D[i])
		{
			for(int j=R[i];j!=i;j=R[j])
			{
				U[D[j]]=U[j];
				D[U[j]]=D[j];
				T[C[j]]--;
			}
		}
		L[R[x]]=L[x];
		R[L[x]]=R[x];
	}
	inline void resume(int x)
	{
		for(int i=U[x];i!=x;i=U[i])
		{
			for(int j=L[i];j!=i;j=L[j])
			{
				U[D[j]]=j;
				D[U[j]]=j;
				T[C[j]]++;
			}
		}
		L[R[x]]=x;
		R[L[x]]=x;
	}
	inline bool dfs()
	{
		if(!R[0])return true;
		int S=R[0],W=T[S],i,j;
		for(i=R[S];i;i=R[i])if(T[i]<W)
		{
			W=T[i];
			S=i;
		}
		remove(S);
		for(i=D[S];i!=S;i=D[i])
		{
			for(j=R[i];j!=i;j=R[j])remove(C[j]);
			if(dfs())return true;
			for(j=L[i];j!=i;j=L[j])resume(C[j]);
		}
		resume(S);
		return false;
	}
}dlx;
int main()
{
//	freopen("test.in","r",stdin);
//	freopen("my.out","w",stdout);
	int n,m;
	while(~scanf("%d%d",&n,&m))
	{
		dlx.init(m);
		dlx.scan(n,m);
		if(dlx.dfs())puts("Yes, I found it");
		else puts("It is impossible");
	}
	return 0;
}

专属那些输出调试者。

	inline void print_line(int x)
	{
		int i=x;
		do{
			printf("%d->",i);
			i=R[i];
		}while(i!=x);
		puts("");
	}
	inline void print_list(int x)
	{
		int i=x;
		do{
			printf("%d->",i);
			i=D[i];
		}while(i!=x);
		puts("");
	}
	inline void print(int x)
	{
		int i=x;
		do{
			print_list(i);
			puts("|");
			i=R[i];
		}while(i!=x);
		puts("");
		puts("");
	}

        inline void remove_point(int x)
	{
		D[U[x]]=D[x];
		U[D[x]]=U[x];
		R[L[x]]=R[x];
		L[R[x]]=L[x];
	}
	inline void resume_point(int x)
	{
		D[U[x]]=x;
		U[D[x]]=x;
		R[L[x]]=x;
		L[R[x]]=x;
	}
	inline void remove_line(int x)
	{
		int i=x;
		do{
			U[D[i]]=U[i];
			D[U[i]]=D[i];
			i=R[i];
		}while(i!=x);
	}
	inline void resume_line(int x)
	{
		int i=x;
		do{
			U[D[i]]=i;
			D[U[i]]=i;
			i=L[i];
		}while(i!=x);
	}
	inline void remove_list(int x)
	{
		int i=x;
		do{
			R[L[i]]=R[i];
			L[R[i]]=L[i];
			i=D[i];
		}while(i!=x);
	}
	inline void resume_list(int x)
	{
		int i=x;
		do{
			R[L[i]]=i;
			L[R[i]]=i;
			i=U[i];
		}while(i!=x);
	}