Description

有一个N*N的迷宫,其中有一些宝藏,现在,小A要从入口(1,1)出发,到达出口(N,N),每次,小A只能从当前的格子走到上下左右四个格子,为了不空手而归,小A决定要拿到所以的宝藏。请问,他最少要走多少步,才能拿到宝藏?

Input Format

第一行:一个整数N,表示迷宫的大小。

第二到第(N+1)行,每行有N个字符,代表迷宫的情况,其中‘0’代表路径,‘1’代表墙面,‘2’代表宝藏

Output Format

一行,一个整数,表示最少步数,或者输出No Solution,表示无法取到所以宝藏

Sample Input

3

011

020

110

Sample Output

4

Hint

100%,N<=30,宝藏个数<=13

Solution

首先,宝藏个数<=13,自然想到状压DP, 用二进制数表示状态,‘1’代表这个宝藏去过了,‘0’代表没有,

那么可以用F[i][S]表示状态为S时,且位于第i个宝藏的最短步数。

那么容易得到方程F[i][S]=min{F[k][S^(1<<(i-1))]+dis[i][k]},其中1<=k<=宝藏个数,dis[i][k]表示第i个宝藏到第k个宝藏的最少步数

那么最后答案就为min{F[i][(1<<m)-1]+toEnd[i]},其中toEnd[i]表示从宝藏i到点(N ,N )的最少步数.

dis数组和toEnd数组可以用bfs预处理得到,细节也不少,当无法到达终点或有至少一个宝藏拿不到就是无解的情况,而且要特判没有宝藏时直接输出到终点的距离

当初挂了好几次,主要有以下几点原因:

  1. 记忆化时tmp忘记初始化为INF,导致答案都为-1
  2. 没有考虑没有宝藏的情况(m==0)
  3. 没有考虑到不了终点的情况
  4. 方程原来写成转移到DP(i, S ^ (1 << (i - 1)),应该是'(p-1)'
  5. 读入时临时字符串数组开小,导致RE

Code

#include <cstdio>
#include <queue>
#include <cstring>
#include <algorithm>
using namespace std; const int dx[4] = {0, 0, 1, -1};
const int dy[4] = {1, -1, 0, 0};
struct point {
point(int a, int b) {x = a, y = b;}
int x, y;
};
struct info {
int x, y;
} bag[32];
int n, a[32][32], m, f[32][1 << 16], dis[32][32], toEnd[32], toBag[32]; int d[32][32];
bool vis[32][32];
queue<point> q;
bool flag;
void bfs(int k, int sx, int sy) {
memset(vis, 0, sizeof(vis));
memset(d, -1, sizeof(d));
d[sx][sy] = 0;
while (!q.empty()) q.pop();
q.push(point(sx, sy));
vis[sx][sy] = 1; while (!q.empty()) {
point u = q.front();
int x = u.x, y = u.y;
q.pop(); for (int i = 0; i < 4; ++i) {
int nx = x + dx[i], ny = y + dy[i];
if (nx <= 0 || ny <= 0 || nx > n || ny > n || vis[nx][ny] || a[nx][ny]) continue;
vis[nx][ny] = 1;
d[nx][ny] = d[x][y] + 1;
q.push(point(nx, ny));
}
} for (int i = k + 1; i <= m; ++i) {
int x = bag[i].x, y = bag[i].y;
if (d[x][y] == -1) {
if (k == 0) flag = 1;
continue;
}
if (k == 0) toBag[i] = d[x][y];
else dis[k][i] = dis[i][k] = d[x][y];
}
if (d[n][n] == -1) flag = 1;
toEnd[k] = d[n][n];
} int DP(int p, int S) {
int &tmp = f[p][S];
if (tmp != -1) return tmp; tmp = (int)1e9;
for (int i = 1; i <= m; ++i) {
if (i == p) continue;
if (!((1 << (i - 1))&S)) continue;
tmp = min(tmp, DP(i, S ^ (1 << (p - 1))) + dis[i][p]);
} return tmp;
} int main() {
scanf("%d", &n);
memset(dis, 0x3f, sizeof(dis));
for (int i = 1; i <= n; ++i) {
char s[32];
scanf("%s\n", s);
for (int j = 0; j < n; ++j) {
a[i][j + 1] = s[j] - '0';
if (a[i][j + 1] == 2) {
a[i][j + 1] = 0;
bag[++m].x = i;
bag[m].y = j + 1;
}
}
} bfs(0, 1, 1);
if (flag) {
printf("No Solution\n");
return 0;
}
for (int i = 1; i <= m; ++i)
bfs(i, bag[i].x, bag[i].y); int Ans = (int)1e9;
memset(f, -1, sizeof(f));
for (int i = 1; i <= m; ++i)
f[i][1 << (i - 1)] = toBag[i];
for (int i = 1; i <= m; ++i)
Ans = min(Ans, DP(i, (1 << m) - 1) + toEnd[i]);
if (m == 0) Ans = toEnd[0];
printf("%d\n", Ans);
return 0;
}

Tirp(状压DP)的更多相关文章

  1. BZOJ 1087: [SCOI2005]互不侵犯King [状压DP]

    1087: [SCOI2005]互不侵犯King Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 3336  Solved: 1936[Submit][ ...

  2. nefu1109 游戏争霸赛(状压dp)

    题目链接:http://acm.nefu.edu.cn/JudgeOnline/problemShow.php?problem_id=1109 //我们校赛的一个题,状压dp,还在的人用1表示,被淘汰 ...

  3. poj3311 TSP经典状压dp(Traveling Saleman Problem)

    题目链接:http://poj.org/problem?id=3311 题意:一个人到一些地方送披萨,要求找到一条路径能够遍历每一个城市后返回出发点,并且路径距离最短.最后输出最短距离即可.注意:每一 ...

  4. [NOIP2016]愤怒的小鸟 D2 T3 状压DP

    [NOIP2016]愤怒的小鸟 D2 T3 Description Kiana最近沉迷于一款神奇的游戏无法自拔. 简单来说,这款游戏是在一个平面上进行的. 有一架弹弓位于(0,0)处,每次Kiana可 ...

  5. 【BZOJ2073】[POI2004]PRZ 状压DP

    [BZOJ2073][POI2004]PRZ Description 一只队伍在爬山时碰到了雪崩,他们在逃跑时遇到了一座桥,他们要尽快的过桥. 桥已经很旧了, 所以它不能承受太重的东西. 任何时候队伍 ...

  6. bzoj3380: [Usaco2004 Open]Cave Cows 1 洞穴里的牛之一(spfa+状压DP)

    数据最多14个有宝藏的地方,所以可以想到用状压dp 可以先预处理出每个i到j的路径中最小权值的最大值dis[i][j] 本来想用Floyd写,无奈太弱调不出来..后来改用spfa 然后进行dp,这基本 ...

  7. HDU 1074 Doing Homework (状压dp)

    题意:给你N(<=15)个作业,每个作业有最晚提交时间与需要做的时间,每次只能做一个作业,每个作业超出最晚提交时间一天扣一分 求出扣的最小分数,并输出做作业的顺序.如果有多个最小分数一样的话,则 ...

  8. 【BZOJ1688】[Usaco2005 Open]Disease Manangement 疾病管理 状压DP

    [BZOJ1688][Usaco2005 Open]Disease Manangement 疾病管理 Description Alas! A set of D (1 <= D <= 15) ...

  9. 【BZOJ1725】[Usaco2006 Nov]Corn Fields牧场的安排 状压DP

    [BZOJ1725][Usaco2006 Nov]Corn Fields牧场的安排 Description Farmer John新买了一块长方形的牧场,这块牧场被划分成M列N行(1<=M< ...

随机推荐

  1. 201521123064 《Java程序设计》第1周学习总结

    1. 本章学习总结 1.了解Java与C/C++的区别,Java语言所写程序较为繁琐,C/C++较为简洁. 2.学会使用cmd进行程序的编译及运行. 3.学会利用JDK辅助Java编程. 4.使用ec ...

  2. 搭建连接MySql的三层架构的ASP.NetCore2.0的WebApi

    里我们用三层架构搭建一个连接MySql的ASP.netCore模板的WebApi项目 首先添加WebApi项目(ASP.NetCore版本) 右键解决方案>新建项目> 选择Web>A ...

  3. 移动商城第四篇【Controller配置、添加品牌之文件上传和数据校验】

    Controller层配置 编写SpringMVC的配置文件 springmvc.xml <?xml version="1.0" encoding="UTF-8&q ...

  4. 接口测试入门(4)--接口自动化测试框架 / list和map用法 / 随机选取新闻 (随机数生成) / 接口相关id映射

    一.接口自动化测试框架 为了更好的组织测试方法,测试用例并且持续集成,我们选择了  java+testNG(测试用例组织)+gitlab(代码版本管理)+Jenkins(持续集成工具) 作为一整套的自 ...

  5. HashMap、HashTable、ArrayList、LinkedList、Vector区别

    HashTable和HashMap区别 ①继承不同. public class Hashtable extends Dictionary implements Map public class Has ...

  6. CA认证和颁发吊销证书

    摘要:涉及到网络安全这一块,想必大家都听过CA吧.像百度.淘宝.京东等这些知名网站,每年都要花费一笔money来买CA证书.但其实简单的企业内的CA认证,我们自己就可以实现,今天小编我就讲解一下怎么在 ...

  7. RobotFramework自动化测试框架-移动手机自动化测试Click Element关键字的使用

    Click Element关键字用来模拟点击APP界面上的一个元素,该关键字接收一个参数[ locator ] ,这里的locator指的是界面元素的定位方式. 示例1:使用Click Element ...

  8. Hadoop就是一个别人造好的轮子

    这个想法源自于我看了<Hadoop: The Definitive Guide>的Part I Ch 2中MapReduce的引入和介绍,书中先说了怎么通过原始的办法处理数据,然后引入到如 ...

  9. jq 通过标签属性匹配标签

    有时候会有这样的需求 在一堆的标签中定位到某个标签去获取值 这时候我们就可以通过使用自定义属性将这一堆标签打上不同的标记 通过属性定位标签了 下面是代码: <input class=" ...

  10. SqlServer2008 导入导出txt或Execl数据

    --右键user表所在的数据库,然后任务--导出数据,然后根据提示设置就行 --从txt中导入 EXEC master..xp_cmdshell 'bcp Northwind.dbo.sysusers ...