题目链接:[kuangbin带你飞]专题三 Dancing Links D - 神龙的难题

题意

Description

这是个剑与魔法的世界.英雄和魔物同在,动荡和安定并存.但总的来说,库尔特王国是个安宁的国家,人民安居乐业,魔物也比较少.但是.总有一些魔物不时会进入城市附近,干扰人民的生活.就要有一些人出来守护居民们不被魔物侵害.魔法使艾米莉就是这样的一个人.她骑着她的坐骑,神龙米格拉一起消灭干扰人类生存的魔物,维护王国的安定.艾米莉希望能够在损伤最小的前提下完成任务.每次战斗前,她都用时间停止魔法停住时间,然后米格拉他就可以发出火球烧死敌人.米格拉想知道,他如何以最快的速度消灭敌人,减轻艾米莉的负担.

Input

数据有多组,你要处理到EOF为止.每组数据第一行有两个数,n,m,(1<=n,m<=15)表示这次任务的地区范围. 然后接下来有n行,每行m个整数,如为1表示该点有怪物,为0表示该点无怪物.然后接下一行有两个整数,n1,m1 (n1<=n,m1<=m)分别表示米格拉一次能攻击的行,列数(行列不能互换),假设米格拉一单位时间能发出一个火球,所有怪物都可一击必杀.

Output

输出一行,一个整数,表示米格拉消灭所有魔物的最短时间.

思路

精确覆盖是每次选取一行后,就要对这行的所有列诛九族。而本题是重复覆盖,所以删除时,只删除相关的列即可,不诛连。
本题需要进行剪枝,剪枝函数f()代码中已说明。
ps:a的第一道重复覆盖问题,因为数据范围被T了那么多次,欲哭无泪啊,把代码扫了一遍又一遍,就是想不到竟然是数据范围的原因。细心是王道啊。。。

代码

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <vector>using namespace std;const int N = 16*16;
const int MAX = N*N;int G[N][N];//储存地图
int P[MAX];//储存点编号
int U[MAX], D[MAX], L[MAX], R[MAX];//数组模拟链表指向(上下左右)
int C[MAX], M[MAX];//节点所在列与行
int S[N];//储存每列的元素数量
int H[N];//行头指针
int ans;void init(int n, int m)
{for(int i=0; i<=m; i++){L[i+1] = i;R[i] = i+1;U[i] = D[i] = i;S[i] = 0;}for(int i=1; i<=n; i++)H[i] = -1;L[0] = m;R[m] = 0;
}void link(int row, int col, int id)//将节点加入链表
{C[id] = col; M[id] = row;//记录行列U[id] = U[col]; D[U[col]] = id;//上下连接D[id] = col; U[col] = id;if(H[row] == -1)//左右连接(使用表头方便头插)H[row] = L[id] = R[id] = id;else{L[id] = L[H[row]]; R[L[H[row]]] = id;L[H[row]] = id; R[id] = H[row];}S[col]++;
}void remove(int col)//因为是重复覆盖,所以一次循环,只删除该列,不删除相关的
{for(int i=D[col]; i!=col; i=D[i]){L[R[i]] = L[i];R[L[i]] = R[i];}
}void resume(int col)
{for(int i=U[col]; i!=col; i=U[i]){L[R[i]] = i;R[L[i]] = i;}
}bool v[MAX];int f()//剪枝函数
{int num = 0;for(int i=R[0]; i!=0; i=R[i])//0代表未消灭v[i] = 0;for(int i=R[0]; i!=0; i=R[i]){if(!v[i]){v[i] = 1;num++;//(假设攻击范围最大化)能够杀掉它的所有攻击的全部范围进行消灭,类似以它为圆心画圆一样for(int j=U[i]; j!=i; j=U[j])for(int k=R[j]; k!=j; k=R[k])v[C[k]] = 1;}}return num;
}void dance(int k)
{if(k+f() >= ans)return;if(!R[0]){ans = min(ans, k);return;}int col = R[0];for(int i=R[0]; i!=0; i=R[i])if(S[i] < S[col])col = i;for(int i=D[col]; i!=col; i=D[i]){remove(i);for(int j=R[i]; j!=i; j=R[j])remove(j);dance(k+1);for(int j=L[i]; j!=i; j=L[j])resume(j);resume(i);}
}int main()
{int n, m;while(~scanf("%d%d", &n, &m)){int num = 0;for(int i=1; i<=n; i++)for(int j=1; j<=m; j++){scanf("%d", &G[i][j]);if(G[i][j])P[j+(i-1)*m] = ++num;}int r, c;scanf("%d%d", &r, &c);init((n-r+1)*(m-c+1), num);int id = num+1;for(int i=1; i<=n-r+1; i++){for(int j=1; j<=m-c+1; j++)for(int x=i; x<=i+r-1; x++)for(int y=j; y<=j+c-1; y++){if(G[x][y])link(j+(i-1)*(m-c+1), P[y+(x-1)*m], id++);}}ans = 0x3f3f3f3f;dance(0);printf("%d\n", ans);}return 0;
}

FZU 1686 神龙的难题(重复覆盖问题舞蹈链)相关推荐

  1. FZU 1686 神龙的难题(DLX反复覆盖)

    FZU 1686 神龙的难题 题目链接 题意:中文题 思路:每个1看成列,每个位置作为左上角的矩阵看成行.dlx反复覆盖就可以 代码: #include <cstdio> #include ...

  2. FZU - 1686 神龙的难题 (舞蹈链,可重覆盖)

    这是个剑与魔法的世界.英雄和魔物同在,动荡和安定并存.但总的来说,库尔特王国是个安宁的国家,人民安居乐业,魔物也比较少.但是.总有一些魔物不时会进入城市附近,干扰人民的生活.就要有一些人出来守护居民们 ...

  3. FZU Problem 1686 神龙的难题(舞蹈链)

    Problem 1686 神龙的难题 Accept: 860    Submit: 2548 Time Limit: 1000 mSec    Memory Limit : 32768 KB  Pro ...

  4. dancing links(舞蹈链)——求解精准覆盖及重复覆盖问题

    以下转自:https://blog.csdn.net/the_star_is_at/article/details/53425736 问题描述: 给定一个n*m的矩阵,有些位置为1,有些位置为0.如果 ...

  5. fzu 1686(DLX 重复点覆盖)

    可以重复的点覆盖, 写法和一般的DLX 很相似,但是在找到选择一行的时候,只删去这一行上所有点的列.   还要注意的是启发函数起了很大的剪枝作用. Problem 1686 神龙的难题 Accept: ...

  6. HDU 3498 whosyourdaddy (可重复覆盖舞蹈链)

    题意: 给出N个点,M个边构成图.每选取一个点都可以覆盖其相邻点,问要覆盖所有点最少选几个点 思路: 每个点都作为一个点集,进行爆搜覆盖,舞蹈链模板题. 代码: #include <bits/s ...

  7. fzu-神龙的难题 舞蹈链之可重复覆盖

    这是个剑与魔法的世界.英雄和魔物同在,动荡和安定并存.但总的来说,库尔特王国是个安宁的国家,人民安居乐业,魔物也比较少.但是.总有一些魔物不时会进入城市附近,干扰人民的生活.就要有一些人出来守护居民们 ...

  8. hdu 2295 Radar DLX 重复覆盖问题

    http://acm.hdu.edu.cn/showproblem.php?pid=229 题意: 一个国家有n个城市,m个雷达,我们同时操作的雷达数最多为k,给出城市与国家的坐标,求小于等于k的操作 ...

  9. hdu 2298 Radar 重复覆盖

    #include <bits/stdc++.h> using namespace std; struct DLX {const static int maxn=64,maxm=64,max ...

最新文章

  1. 解决linux下中文文件名显示乱码问题
  2. Hadoop参数汇总
  3. 苹果电脑怎么设置佳博标签打印机_热销斑马标签打印机恢复出厂设置的方法
  4. lda主题评论文本python_利用python做LDA文本分析,该从哪里入手呢?
  5. 肝!用 Python 加密文件
  6. final const java_Java中的final关键字 与 C#中的const, readonly关键字
  7. Mac终端给命令设置别名alias的办法
  8. linux命令chown和chmod什么区别
  9. java 获取本机的IP和hostname
  10. 3级城市选择(数据库版)
  11. Git的17条基本用法
  12. 5. JavaScript RegExp 类型
  13. Matlab画图常用的指令是啥,matlab画图常用命令
  14. 服务器2012怎么安装无线网卡驱动,如何安装usb无线网卡驱动,教您如何安装电脑usb驱动...
  15. java ftpclient prompt,一个JAVA FTP客户端制作的收获
  16. Linux终端解压压缩文件,如何从终端解压zip文件?
  17. @Cacheable(sync=true) only allows a single cache on
  18. ospf及lsa总结
  19. 在html标签中写alert,文本不会写入页面,但会显示在alert()消息中。 JavaScript和HTML...
  20. 西门子PLC 冷启动、暖启动、热启动的一种通俗解释

热门文章

  1. Vue2 模板template的四种写法
  2. 工作记录:EasyExcel对于不定长Excel表头的处理,写excel
  3. pandas错位计算
  4. 试卷上出现的爆笑答案
  5. 【经典算法学习-排序篇】直接选择排序
  6. 欧美开源法案频出,你准备好了吗?
  7. 太原铁警严厉打击倒卖车票 查获假票60张
  8. 闪电网络预备知识一:闪电网络概念以及为什么使用闪电网络
  9. 基于YOLOv5的疲劳驾驶检测系统(Python+清新界面+数据集)
  10. webbuilder8开发经验(二)