题意:给你一个01矩阵,让你选择尽可能少的行数,使得这些行的并集能够覆盖到所有列。

DLX算法求解重复覆盖问题模板,使用估价函数进行剪枝。

#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
const int maxn=55+5;
const int maxr=55+5;
const int maxnode=55*55+5;
//ÐбàºÅ´Ó1¿ªÊ¼£¬ÁбàºÅΪ1~n£¬½áµã0ÊDZíÍ·½áµã£»½áµã1~nÊǸ÷Áж¥²¿µÄÐéÄâ½áµã
struct DLX{
int n,sz;//ÁÐÊý£¬½áµã×ÜÊý
int S[maxn];//¸÷ÁнáµãÊý
int row[maxnode],col[maxnode];//¸÷½áµãËùÔÚµÄÐÐÁбàºÅ
int L[maxnode],R[maxnode],U[maxnode],D[maxnode];
int ansd,ans[maxr];//½â
void init(int n){//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;
sz=n+1; ansd=1000000007;
memset(S,0,sizeof(S));
}
void addRow(int r,vector<int> columns){
int first=sz;
for(int i=0;i<columns.size();++i){
int c=columns[i];
L[sz]=sz-1;
R[sz]=sz+1;
D[sz]=c;
U[sz]=U[c];
D[U[c]]=sz;
U[c]=sz;
row[sz]=r;
col[sz]=c;
++S[c];
++sz;
}
R[sz-1]=first;
L[first]=sz-1;
}
//˳×ÅÁ´±íA£¬±éÀú³ýsÍâµÄÆäËûÔªËØ
#define FOR(i,A,s) for(int i=A[s];i!=s;i=A[i])
void remove(int c){
L[R[c]]=L[c];
R[L[c]]=R[c];
FOR(i,D,c){
FOR(j,R,i){
U[D[j]]=U[j];
D[U[j]]=D[j];
--S[col[j]];
}
}
}
void restore(int c){
FOR(i,U,c){
FOR(j,L,i){
++S[col[j]];
U[D[j]]=j;
D[U[j]]=j;
}
}
L[R[c]]=c;
R[L[c]]=c;
}
bool dfs(int d){
// printf("%d",d);
if(R[0]==0){//ÕÒµ½½â
ansd=d;//¼Ç¼½âµÄ³¤¶È
return 1;
}
//ÕÒ½áµãÊý×îСµÄÁÐc
int c=R[0];//µÚÒ»¸öδɾ³ýµÄÁÐ
FOR(i,R,0){
if(S[i]<S[c]){
c=i;
}
}
remove(c);//ɾ³ýµÚcÁÐ
FOR(i,D,c){//ÓýáµãiËùÔÚÐи²¸ÇµÚcÁÐ
ans[d]=row[i];
FOR(j,R,i){
remove(col[j]);//ɾ³ý½áµãiËùÔÚÐÐÄܸ²¸ÇµÄËùÓÐÆäËûÁÐ
}
if(dfs(d+1)){
return 1;
}
FOR(j,L,i){
restore(col[j]);//»Ö¸´½áµãiËùÔÚÐÐÄܸ²¸ÇµÄÆäËûËùÓÐÁÐ
}
}
restore(c);//»Ö¸´µÚcÁÐ
return 0;
}
bool solve(vector<int>& v){
v.clear();
if(!dfs(0)){
return 0;
}
for(int i=0;i<ansd;++i){
v.push_back(ans[i]);
}
return 1;
} void remov2(int c){
FOR(i,D,c){
L[R[i]]=L[i];
R[L[i]]=R[i];
}
}
void restor2(int c){
FOR(i,U,c){
L[R[i]]=R[L[i]]=i;
}
}
int f(){//¹À¼Ûº¯Êý
bool vis[maxn];
int res=0;
FOR(i,R,0){
vis[i]=1;
}
FOR(i,R,0){
if(vis[i]){
++res;
vis[i]=0;
FOR(j,D,i){
FOR(k,R,j){
vis[col[k]]=0;
}
}
}
}
return res;
}
void df2(int d){
if(d+f()>=ansd){
return;
}
if(R[0]==0){//ÕÒµ½½â
ansd=min(ansd,d);//¼Ç¼½âµÄ³¤¶È
return;
}
//ÕÒ½áµãÊý×îСµÄÁÐc
int c=R[0];//µÚÒ»¸öδɾ³ýµÄÁÐ
FOR(i,R,0){
if(S[i]<S[c]){
c=i;
}
}
FOR(i,D,c){
remov2(i);
FOR(j,R,i){
remov2(j);
}
df2(d+1);
FOR(j,L,i){
restor2(j);
}
restor2(i);
}
} }dlx;
int n,m;
bool a[60][60];
int main(){
//freopen("hdu3498.in","r",stdin);
int x,y;
while(scanf("%d%d",&n,&m)!=EOF){
dlx.init(n);
memset(a,0,sizeof(a));
for(int i=1;i<=n;++i){
a[i][i]=1;
}
for(int i=1;i<=m;++i){
scanf("%d%d",&x,&y);
a[x][y]=a[y][x]=1;
}
for(int i=1;i<=n;++i){
vector<int> columns;
for(int j=1;j<=n;++j){
if(a[i][j]){
columns.push_back(j);
}
}
if(!columns.empty()){
dlx.addRow(i,columns);
}
}
dlx.df2(0);
printf("%d\n",dlx.ansd);
}
return 0;
}

【DLX算法】hdu3498 whosyourdaddy的更多相关文章

  1. 关于用舞蹈链DLX算法求解数独的解析

    欢迎访问——该文出处-博客园-zhouzhendong 去博客园看该文章--传送门 描述 在做DLX算法题中,经常会做到数独类型的题目,那么,如何求解数独类型的题目?其实,学了数独的构建方法,那么DL ...

  2. 搜索:DLX算法

    精确覆盖问题:在一个0-1矩阵中,选定部分行,使得每一列都有且只有一个1.求解一种选法 舞蹈链(Dance Link),也就是一个循环十字链表,可以快速的删掉和恢复某行某列 结合了舞蹈链的搜索就称作D ...

  3. 【DLX算法】poj2676 Sudoku

    DLX算法求解精确覆盖问题模板.赛场上可以参见白书. #include<cstdio> #include<cstring> #include<vector> usi ...

  4. DLX算法一览

    目录: 1 X思想的了解. 链表的递归与回溯. 具体操作. 优化. 一些应用与应用中的再次优化(例题). 练手题 X思想的了解. 首先了解DLX是什么? DLX是一种多元未饱和型指令集结构,DLX 代 ...

  5. 精确覆盖DLX算法模板另一种写法

    代码 struct DLX { int n,id; int L[maxn],R[maxn],U[maxn],D[maxn]; ]; int H[ms]; ) //传列长 { n=nn; ;i<= ...

  6. 精确覆盖DLX算法模板

    代码 struct DLX { int n,id; int L[maxn],R[maxn],U[maxn],D[maxn]; ]; ) //传列长 { n=nn; ;i<=n;i++) U[i] ...

  7. 舞蹈链(DLX)

    舞蹈链(DLX) Tags:搜索 作业部落 评论地址 一.概述 特别特别感谢这位童鞋His blog 舞蹈链是一种优美的搜索,就像下面这样跳舞- 舞蹈链用于解决精确覆盖或者重复覆盖的问题 你可以想象成 ...

  8. 数独求解 DFS && DLX

    题目:Sudoku 题意:求解数独.从样例和结果来看应该是简单难度的数独 思路:DFS 设置3个数组,row[i][j] 判断第i行是否放了j数字,col[i][j] 判断第i列是否放了j数字.squ ...

  9. POJ2676,HDU4069解决数独的两种实现:DFS、DLX

    搜索实现:解决数独有两种思考策略,一种是枚举当前格能填的数字的种数,这里有一优化策略就是先搜索能填入种数小的格子:另一种是考虑处理某一行(列.宫)时,对于某一个没用过的数字,若该行(列.宫)只有一个可 ...

随机推荐

  1. 【leetcode 简单】 第八题 删除排序数组中的重复项

    给定一个排序数组,你需要在原地删除重复出现的元素,使得每个元素只出现一次,返回移除后数组的新长度. 不要使用额外的数组空间,你必须在原地修改输入数组并在使用 O(1) 额外空间的条件下完成. 示例 1 ...

  2. Sublime之插件的安装(三)

    今天在写js的时候,突然想到一个问题就是能不能快速的对齐的插件,当当当~~~sublime这么神器当然有,那就是:Alignment插件 如果你写的代码是这样的: var a = , b =, ccc ...

  3. 使用vscode实现git同步

    用了git最方便的就是项目同步管理,回到家打开vscode只需要点击一下pull就能全部同步过来.是不是很方便....毕竟之前我都是拿u盘拷贝回家或者存到云盘再下载下来..   我这里之前用的是国内的 ...

  4. javascript经典小游戏代码集合

    http://www.jb51.net/Special/349.htm

  5. 使用JSON语法创建JS对象(重要)

    JS对象的键值可以加单引号或者不加或者加双引号 JSON语法提供了一种更简单的方式来创建对象,可以避免书写函数,也可避免用new关键字,可以直接创建一个JS对象,使用一个花括号,然后将每个属性写成&q ...

  6. 64_r2

    ruby-gnomecanvas2-0.90.4-7.fc26.3.x86_64.rpm 13-Feb-2017 08:00 75794 ruby-gnomecanvas2-devel-0.90.4- ...

  7. javascript反混淆之packed混淆(一)

    javascript反混淆之packed混淆(一) 什么是JavaScript反混淆,在理解这个概念前我们先来看下什么是代码混淆,代码混淆,是将计算机程序的代码,转换成一种功能上等价,但是难于阅读和理 ...

  8. 设计模式之笔记--外观模式(Facade)

    外观模式(Facade) 定义 外观模式(Facade),为子系统中的一组接口提供一个一致的界面,此模式定义了一个高层接口,这个接口使得这一子系统更加容易使用. 类图 描述 Facade:外观类,外观 ...

  9. 设计模式--工厂模式 caffe_layer注册

    来源:http://www.cnblogs.com/zhouqiang/archive/2012/07/20/2601365.html 来源:http://blog.luoyetx.com/2016/ ...

  10. JavaScript变量命名规则:匈牙利命名法

    匈牙利命名法语法 变量名=类型+对象描述 类型指变量的类型 对象描述指对象名字全称或名字的一部分,要求有明确含义,命名要容易记忆容易理解. 提示 虽然JavaScript变量表面上没有类型,但是Jav ...