题意:多组数据,每组数据给你几行数,要求选出当中几行。使得每一列都有且仅有一个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);
}

然后数据我就不贴了,直接去看3740discuss吧,有两组。

【POJ3740】Easy Finding DLX(Dancing Links)精确覆盖问题的更多相关文章

  1. HDU 3111 Sudoku ( Dancing Links 精确覆盖模型 )

    推荐两篇学DLX的博文: http://bbs.9ria.com/thread-130295-1-1.html(这篇对DLX的工作过程演示的很详细) http://yzmduncan.iteye.co ...

  2. 【转】Dancing Links精确覆盖问题

    原文链接:http://sqybi.com/works/dlxcn/ (只转载过来一部分,全文请看原文,感觉讲得很好~)正文    精确覆盖问题    解决精确覆盖问题    舞蹈步骤    效率分析 ...

  3. hihoCoder #1321 : 搜索五•数独 (Dancing Links ,精确覆盖)

    hiho一下第102周的题目. 原题地址:http://hihocoder.com/problemset/problem/1321 题意:输入一个9*9数独矩阵,0表示没填的空位,输出这个数独的答案. ...

  4. hust 1017 dancing links 精确覆盖模板题

    最基础的dancing links的精确覆盖题目 #include <iostream> #include <cstring> #include <cstdio> ...

  5. ZOJ 3209 Treasure Map (Dancing Links 精确覆盖 )

    题意 :  给你一个大小为 n * m 的矩形 , 坐标是( 0 , 0 ) ~ ( n , m )  .然后给你 p 个小矩形 . 坐标是( x1 , y1 ) ~ ( x2 , y2 ) , 你选 ...

  6. POJ3074 Sudoku —— Dancing Links 精确覆盖

    题目链接:http://poj.org/problem?id=3074 Sudoku Time Limit: 1000MS   Memory Limit: 65536K Total Submissio ...

  7. HUST1017 Exact cover —— Dancing Links 精确覆盖 模板题

    题目链接:https://vjudge.net/problem/HUST-1017 1017 - Exact cover 时间限制:15秒 内存限制:128兆 自定评测 7673 次提交 3898 次 ...

  8. HDU5046 Airport dancing links 重复覆盖+二分

    这一道题和HDU2295是一样 是一个dancing links重复覆盖解决最小支配集的问题 在给定长度下求一个最小支配集,只要小于k就行 然后就是二分答案,每次求最小支配集 只不过HDU2295是浮 ...

  9. dancing link 精确覆盖 重复覆盖 (DLX)

    申明:因为转载的没有给出转载链接,我就把他的链接附上,请尊重原创: http://www.cnblogs.com/-sunshine/p/3358922.html 如果谁知道原创链接 给一下,请尊重原 ...

随机推荐

  1. 常州模拟赛d4t3 字符串划分

    题目描述 给你一串由小写字母组成的字符串,希望你把它划分成一些小段,使得每一小段字符串中的字母 都不相同,并且希望分的段数尽量少. 然后,把这些小段按字典序排序后输出,中间由一个空格分隔. 例如:字符 ...

  2. Java 面试参考指南 — 同步

    同步 在多线程程序中,同步修饰符用来控制对临界区代码的访问.其中一种方式是用synchronized关键字来保证代码的线程安全性.在Java中,synchronized修饰的代码块或方法不会被多个线程 ...

  3. bzoj3609 [Heoi2014]人人尽说江南好 博弈

    [Heoi2014]人人尽说江南好 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 581  Solved: 420[Submit][Status][D ...

  4. cf701E Connecting Universities

    Treeland is a country in which there are n towns connected by n - 1 two-way road such that it's poss ...

  5. Laravel 5 Form 和 HTML 的使用

    最近在用 laravel 5 做例子,在做到表单的时候,习惯性的使用 Form::open() 结果发现提示错误,没有这个类, 好吧,找了找,发现 在laravel 5 中,把 from 和 html ...

  6. CSV模块

    Python csv模块的使用 1.csv简介 CSV (Comma Separated Values),即逗号分隔值(也称字符分隔值,因为分隔符可以不是逗号),是一种常用的文本 格式,用以存储表格数 ...

  7. LINUX下面NetworkManager和network冲突的问题

    https://blog.csdn.net/ID_EAGLE/article/details/74085409

  8. JavaScript ES6中,export与export default

    自述: 本来是对new Vue()和export default比较懵的,查了一下,发现我理解错了两者的关系,也没意识到export与export default的区别,先简单的记录一下基本概念,后续 ...

  9. Oracle ORA-01033: ORACLE initialization or shutdown in progress

    先说明,我出现此错误的原因是:我手动通过drop语句删除表空间,结果磁盘中文件还存在,然后我手动删除了文件,重启了oracle服务,再去连接oracle时就出现了这个错误. 网上也有“连接Oracle ...

  10. Spring中Bean的定义继承

    以下内容引用自http://wiki.jikexueyuan.com/project/spring/bean-definition-inheritance.html: Bean定义继承 bean定义可 ...