BZOJ2595 [Wc2008]游览计划 【状压dp + 最短路】
题目链接
题解
著名的斯坦纳树问题
设\(f[i][j][s]\)表示点\((i,j)\)与景点联通状况为\(s\)的最小志愿者数
设\(val[i][j]\)为\((i,j)\)需要的志愿者数
有两种转移
一种是自己转移
\]
一种是由周围转移过来
\]
第一种\(O(3^{K})\)枚举子集,第二种就是最短路
纪念一下BZOJ500题,截个图,,数字挺整的
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<vector>
#include<map>
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define mp(a,b) make_pair<int,int>(a,b)
#define cls(s) memset(s,0,sizeof(s))
#define cp pair<int,int>
#define LL long long int
using namespace std;
const int maxn = 11,maxm = 1 << 10,INF = 1000000000;
inline int read(){
int out = 0,flag = 1; char c = getchar();
while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();}
while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}
return out * flag;
}
int n,m,K,maxv;
int f[maxn][maxn][maxm],val[maxn][maxn],id[maxn][maxn],vis[maxn][maxn];
int head,tail,X[4] = {0,0,-1,1},Y[4] = {-1,1,0,0};
int S[maxn][maxn];
cp q[10000];
struct tri{
int x,y,s;
};
vector<tri> pre[maxn][maxn][maxm];
void spfa(int s){
head = 0; tail = -1;
REP(i,n) REP(j,m) q[++tail] = mp(i,j),vis[i][j] = true;
cp u; int nx,ny;
while (head <= tail){
u = q[head++];
vis[u.first][u.second] = false;
for (int k = 0; k < 4; k++){
nx = u.first + X[k];
ny = u.second + Y[k];
if (nx < 1 || ny < 1 || nx > n || ny > m) continue;
if (f[nx][ny][s] > f[u.first][u.second][s] + val[nx][ny]){
f[nx][ny][s] = f[u.first][u.second][s] + val[nx][ny];
pre[nx][ny][s].clear();
pre[nx][ny][s].push_back((tri){u.first,u.second,s});
if (!vis[nx][ny]) q[++tail] = mp(nx,ny);
}
}
}
}
void dfs(int x,int y,int s){
S[x][y] = true;
tri u;
for (unsigned int i = 0; i < pre[x][y][s].size(); i++){
u = pre[x][y][s][i];
dfs(u.x,u.y,u.s);
}
}
void work(){
maxv = (1 << K) - 1;
memset(f,0x3f3f3f3f,sizeof(f));
REP(i,n) REP(j,m){
f[i][j][0] = 0;
if (!val[i][j]) f[i][j][1 << id[i][j] - 1] = 0;
}
for (int s = 0; s <= maxv; s++){
REP(i,n) REP(j,m){
for (int e = s; e; e = (e - 1) & s){
if (f[i][j][s] > f[i][j][e] + f[i][j][s ^ e] - val[i][j]){
f[i][j][s] = f[i][j][e] + f[i][j][s ^ e] - val[i][j];
pre[i][j][s].clear();
pre[i][j][s].push_back((tri){i,j,e});
pre[i][j][s].push_back((tri){i,j,s ^ e});
}
}
}
spfa(s);
}
int ans = INF,x,y;
REP(i,n) REP(j,m) if (ans > f[i][j][maxv]) ans = f[i][j][maxv],x = i,y = j;
printf("%d\n",ans);
dfs(x,y,maxv);
REP(i,n){
REP(j,m){
if (id[i][j]) putchar('x');
else if (S[i][j]) putchar('o');
else putchar('_');
}
puts("");
}
}
int main(){
n = read(); m = read();
REP(i,n) REP(j,m){
val[i][j] = read();
if (!val[i][j]) id[i][j] = ++K;
}
work();
return 0;
}
BZOJ2595 [Wc2008]游览计划 【状压dp + 最短路】的更多相关文章
- [WC2008]游览计划(状压dp)
题面太鬼畜不粘了. 题意就是给一张n*m的网格图,每个点有点权,有k个关键点,让你把这k个关键点连成一个联通快的最小代价. 题解 这题nmk都非常小,解法肯定是状压,比较一般的解法插头dp,但不太好写 ...
- luogu4294 [WC2008]游览计划(状压DP/斯坦纳树)
link 题目大意:给定一个网格图,有些点是关键点,选择格点有代价,求把所有关键点联通的最小代价 斯坦纳树模板题 斯坦纳树问题:给定一个图结构,有一些点是关键点,求把这些关键点联通的最小代价e 斯坦纳 ...
- [WC2008]游览计划 状压DP,斯坦纳树
---题面--- 题解: 这是一道斯坦纳树的题,用状压+spfa来解决 什么是斯坦纳树? 一开始还以为是数据结构来着,其实跟最小生成树很像,大致就是最小生成树只能在各个点之间直接相连,而斯坦纳树则允许 ...
- BZOJ2595 Wc2008 游览计划 【斯坦纳树】【状压DP】*
BZOJ2595 Wc2008 游览计划 Description Input 第一行有两个整数,N和 M,描述方块的数目. 接下来 N行, 每行有 M 个非负整数, 如果该整数为 0, 则该方块为一个 ...
- 旅游(CSUST省赛选拔赛2+状压dp+最短路)
题目链接:http://csustacm.com:4803/problem/1016 题目: 思路:状压dp+最短路,比赛的时候有想到状压dp,但是最短路部分写挫了,然后就卡死了,对不起出题人~dis ...
- BZOJ2595 WC2008游览计划(斯坦纳树)
斯坦纳树板子题. 考虑状压dp,设f[i][j][S]表示当前在点(i,j)考虑转移,其所在的联通块包含的关键点集(至少)为S的答案. 转移时首先枚举子集,有f[i][j][S]=min{f[i][j ...
- 斯坦纳树 [bzoj2595][wc2008]游览计划 题解
话说挺早就写过斯坦纳树了,不过当时没怎么总结,也不是很理解……现在来个小结吧~ 斯坦纳树就是包含给定点的最小生成树(个人理解权值应当为正). 一般来讲,给定点的数目应该很小吧...于是我们可以用状压D ...
- [bzoj2595][WC2008]游览计划/[bzoj5180][Baltic2016]Cities_斯坦纳树
游览计划 bzoj-2595 wc-2008 题目大意:题目链接.题目连接. 注释:略. 想法:裸题求斯坦纳树. 斯坦纳树有两种转移方式,设$f[s][i]$表示联通状态为$s$,以$i$为根的最小代 ...
- Codeforces 375C - Circling Round Treasures(状压 dp+最短路转移)
题面传送门 注意到这题中宝藏 \(+\) 炸弹个数最多只有 \(8\) 个,故考虑状压,设 \(dp[x][y][S]\) 表示当前坐标为 \((x,y)\),有且仅有 \(S\) 当中的物品被包围在 ...
随机推荐
- halcon学习相关资料(转载)
https://blog.csdn.net/maweifei/article/details/78162581 论坛.培训 halcon学习网:http://www.ihalcon.com/ 鸟叔机器 ...
- JVM自动内存管理机制--读这篇就GO了
之前看过JVM的相关知识,当时没有留下任何学习成果物,有些遗憾.这次重新复习了下,并通过博客来做下笔记(只能记录一部分,因为写博客真的很花时间),也给其他同行一些知识分享. Java自动内存管理机制包 ...
- 谈谈你对Java异常处理机制的理解
先谈谈我的理解:异常处理机制可以说是让我们编写的程序运行起来更加的健壮,无论是在程序调试.运行期间发生的异常情况的捕获,都提供的有效的补救动作,任何业务逻辑都会存在异常情况,这时只需要记录这些异常情况 ...
- phpcms v9如何给父级单页栏目添加内容
对于phpcms单页的调用相信大家都应该没问题,那么如果我们在后台添加的单页有二层甚至更多的时候,这样在管理内容上是没有给父级栏目添加内容这一功能的!那么我们该怎么实现这个功能并调用呢? 首先我们要修 ...
- Echarts简单图表
一.实现要点 常用可视化图表库 Echarts HighCharts D3.js neo4j (NOSQL) 1.导入js库 <script type="text/javascript ...
- Delphi 实现照片抽奖-原创
有单位年会要用照片抽奖,上网搜了几个都不满意,且居然还要收费.自己写一个算了.只是有一点不爽,Delphi 7 在 Windows 7 64位下有问题,不能双击 dpr 文件直接打开项目! 关于性能: ...
- Matlab中 .' 的作用。
Syntax B = A.' B = transpose(A) Description B = A.' returns the nonconjugate transpose of A, that ...
- NO.3:自学python之路------集合、文件操作、函数
引言 本来计划每周完成一篇Python的自学博客,由于上一篇到这一篇遇到了过年.开学等杂事,导致托更到现在.现在又是一个新的学期,春天也越来越近了(冷到感冒).好了,闲话就说这么多.开始本周的自学Py ...
- Elasticsearch 统计代码例子
aggs avg 平均数 最近15分钟的平均访问时间,upstream_time_ms是每次访问时间,单位毫秒 { "query": { "filtered": ...
- Leftmost Digit(数学)
Description Given a positive integer N, you should output the leftmost digit of N^N. Input The inp ...