[UVA 10605]钻石矿道【启发式搜索+状态压缩】

[UVA 10605]钻石矿道

很好的一道题就是了，做了两遍了，要理清楚其中的状态转移，其余的倒是不难处理，这里的dp[ i ][ j ]表示的是我们这个一块，也就是它们被合并成一个块的包含j这个点的最短路径。g[ i ]表示的是状态i下的这个块的出去的最短距离。f[ i ]表示的是这个状态下出去的最短距离，没必要在一个块中。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
#define MP(x, y) make_pair(x, y)
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
const int maxN = 15;
const int dir[4][2] =
{-1, 0,0, -1,0, 1,1, 0
};
int N, M, tot, h[maxN], g[1<<maxN], f[1<<maxN];
char mp[maxN][maxN];
struct node
{int x, y;node(int a=0, int b=0):x(a), y(b) {}
}a[maxN * maxN];
inline int _dis(node e1, node e2) { return abs(e1.x - e2.x) + abs(e1.y - e2.y); }
bool vis[maxN][maxN];
pair<node, int> now;
queue<pair<node, int>> Q;
inline int bfs(int sx, int sy)
{memset(vis, false, sizeof(vis));while(!Q.empty()) Q.pop();Q.push(MP(node(sx, sy), 0));vis[sx][sy] = true;while(!Q.empty()){now = Q.front(); Q.pop();for(int i=0, xx, yy; i<4; i++){xx = now.first.x + dir[i][0]; yy = now.first.y + dir[i][1];if(vis[xx][yy]) continue;vis[xx][yy] = true;if(mp[xx][yy] == '#') return now.second + 1;Q.push(MP(node(xx, yy), now.second + 1));}}return 0;
}
int dis[maxN], dp[1<<maxN][maxN];
inline void pre_did()
{for(int i=0; i<tot; i++) dis[i] = bfs(a[i].x, a[i].y);for(int i=0; i<tot; i++) dp[1<<i][i] = dis[i];for(int i=1; i<(1<<tot); i++){for(int j=0; j<tot; j++){if((1<<j) & i){for(int k=0; k<tot; k++){if(((1<<k) & i) && k != j)  //不能是两个相同的点{dp[i][j] = min(dp[i][j], dp[i^(1<<j)][k] + _dis(a[j], a[k]));}}}}for(int j=0; j<tot; j++) g[i] = min(g[i], dp[i][j]);}for(int i=1; i<(1<<tot); i++){for(int j=1; j<=i; j++){if((i & j) == j){f[i] = min(f[i], f[i^j] + g[j]);}}}printf("%d\n", f[(1<<tot) - 1]);
}
inline void init()
{tot = 0;memset(dp, INF, sizeof(dp));dp[0][0] = 0;memset(g, INF, sizeof(g));memset(f, INF, sizeof(f));g[0] = f[0] = 0;
}
int main()
{
//    freopen("uva10605.in", "r", stdin);
//    freopen("uva10605.out", "w", stdout);int T; scanf("%d", &T);while(T--){scanf("%d%d", &N, &M);init();for(int i=1; i<=N; i++){scanf("%s", mp[i] + 1);for(int j=1; j<=M; j++) if(mp[i][j] == '*') a[tot++] = node(i, j);}pre_did();}return 0;
}

[UVA 10605]钻石矿道【启发式搜索+状态压缩】相关推荐

UVA 11825 状态压缩DP+子集思想
很明显的状态压缩思想了.把全集分组,枚举每个集合的子集,看一个子集是否能覆盖所有的点,若能,则f[s]=max(f[s],f[s^s0]+1).即与差集+1比较. 这种枚举集合的思想还是第一次遇到,果 ...
UVA - 1633 Dyslexic Gollum 状态压缩
问题输入N,K(N<=400,K<=10),求长度为N且不含有长度至少为K的连续回文子串的01字符串有多少个? 分析输入k,要判断k,k+1的情况,因为一个字符串有长度为k的回文子串, ...
[状态压缩DP] COJ 1129 送货到家
第一道状态压缩DP: 这道题要求一个无向图的最小权回路,要求经过所有点,所以可以任选一个点(这里选0)作为起点,以后的状态f[s, i]表示从0出发到i结束的最小权路径,最终求得f[1<< ...
327 玉米田（状态压缩dp）
1. 问题描述: 农夫约翰的土地由 M×N 个小方格组成,现在他要在土地里种植玉米.非常遗憾,部分土地是不育的,无法种植.而且,相邻的土地不能同时种植玉米,也就是说种植玉米的所有方格之间都不会有公共边 ...
TSOJ 好好做题(屑)——递推状态压缩+高精度
文章目录题目描述解题思路喜闻乐见的AC代码补充一:高精度补充二:状态压缩完整注释代码后记题目描述选修程序设计和算法课程的学生人数为 n,任课老师设置了 m 道练习题目(其中:1 &l ...
sdoi2009 [动态规划状态压缩DP] 学校食堂
背景飘逸的EWF组合~ 描述小F的学校在城市的一个偏僻角落,所有学生都只好在学校吃饭.学校有一个食堂,虽然简陋,但食堂大厨总能做出让同学们满意的菜肴.当然,不同的人口味也不一定相同,但每个人的口味 ...
洛谷P1562 还是N皇后（DFS+状态压缩+位运算）
八皇后问题的介绍在此不再赘述,只贴一下经典八皇后问题的实现代码(参考刘汝佳 <算法竞赛入门经典>) void search(int i) {if(i>n){ans++;return; ...
第一章动态规划状态压缩DP
1.基本概述状态压缩dp和状态机一样,都是一种特殊的状态表示方式.状态机用一系列小状态表示某一状态.状态压缩dp用二进制数进行表示.虽然看代码起来时间复杂度比较高,但是很多的情况都给剪枝掉了. 状态 ...
洛谷 P-4045 密码(AC自动机+状态压缩+数位DP+乱搞)
洛谷 P-4045 密码记AC的第一道黑题! 题意:已知一段密码包含了一些字符串,然后求满足条件的密码有多少个,数量小于42时还得全部输出思路: 一开始WA了两个点,不知道WA的什么,索性把读入的 ...

[UVA 10605]钻石矿道【启发式搜索+状态压缩】

[UVA 10605]钻石矿道

[UVA 10605]钻石矿道【启发式搜索+状态压缩】相关推荐

最新文章

热门文章