bzoj2595 / P4294 [WC2008]游览计划
斯坦纳树
斯坦纳树,是一种神奇的树。它支持在一个连通图上求包含若干个选定点的最小生成树。
前置算法:spfa+状压dp+dfs(大雾)
我们设$f[o][P]$为第$o$个点上状态为$P$的最小代价,其中状态使用二进制存储已经连接了多少个选定点。
初始化:显然对于每个选定点,$f[o][1<<k]=0$,$k$为该选定点在所有选定点中的编号。其他为$inf$
蓝后就是将状态从小到大枚举进行递推
$for(P=0;P<(1<<k);++P)$
对于每层递推,枚举所有点$1$~$o$;
我们先考虑这个点连接了2个不同连通块(链)的状况(并为spfa做好准备)
于是我们枚举$P$的真子集进行递推
$for(j=(P-1)\&P;j;j=(j-1)\&P)$
枚举真子集↑
$f[o][P]=min(f[o][P],f[o][j]+f[o][P$^$j]-val[o]);$
注意该式适用于计算点权,减去$val[o]$是去掉重复点权。如果计算边权需作修改。
但是这样显然远远不够。于是我们用$spfa$通过类似$dp$的形式处理好剩下的所有状况。
对于前面的$f[o][P]$,任何$f[o][P]<inf$都应作为每层spfa的起点(显然spfa也是每层执行)
在共$K$个给定点中,随意找一个给定点作为树根$rt$。
$ans=f[rt][(1<<K)-1]$
对于本题$(extra)$:输出一种方案
我们在每次更新$f[o][P]$时
用$fa[o][P]$记下$f[o][P]$从哪个状态推导而来
最后从树根$rt$用$dfs$倒推回去,更新答案即可。
void dfs(pii u,int o){//u:坐标
par G=fa[F(u)][o];
if(!G.se) return;
use[u.fi][u.se]=;
if(G.fi==u) dfs(u,o^G.se);//两个联通块合并而来,则要从两条路递推回去。
dfs(G.fi,G.se);
}
$code:$
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
typedef pair<int,int> pii;
typedef pair<pii,int> par;
#define mp make_pair
#define fi first
#define se second
int d1[]={,,,-};
int d2[]={,,-,};
int n,m,tot,k,inf,f[][],id[][],a[][];
bool in[][],use[][];
par fa[][]; pii rt;
queue <pii> h;
inline int F(pii x){return id[x.fi][x.se];}
void spfa(int o){
while(!h.empty()){
pii x=h.front(); h.pop(); in[x.fi][x.se]=;
for(int i=;i<;++i){
int r1=x.fi+d1[i],r2=x.se+d2[i],to=id[r1][r2];
if(r1<||r1>=n||r2<||r2>=m) continue;
if(f[to][o]>f[F(x)][o]+a[r1][r2]){
f[to][o]=f[F(x)][o]+a[r1][r2];
fa[to][o]=mp(x,o);
if(!in[r1][r2]) h.push(mp(r1,r2));
}
}
}
}
void dfs(pii u,int o){//逆推找可行方案
par G=fa[F(u)][o];
if(!G.se) return;
use[u.fi][u.se]=;//选择点u
if(G.fi==u) dfs(u,o^G.se);
dfs(G.fi,G.se);
}
int main(){
memset(f,,sizeof(f)); inf=f[][];
scanf("%d%d",&n,&m);
for(int i=;i<n;++i)
for(int j=;j<m;++j){
scanf("%d",&a[i][j]);
id[i][j]=++tot;
if(!a[i][j]) rt=mp(i,j),f[tot][<<k]=,++k;
}
for(int P=;P<(<<k);spfa(P),++P)//每层递推的最后来一次spfa
for(int x=;x<n;++x)
for(int y=;y<m;++y){
int o=id[x][y];
for(int j=(P-)&P;j;j=(j-)&P)//枚举真子集
if(f[o][P]>f[o][j]+f[o][P^j]-a[x][y]){
f[o][P]=f[o][j]+f[o][P^j]-a[x][y];
fa[o][P]=mp(mp(x,y),j);
}
if(f[o][P]<inf) h.push(mp(x,y)),in[x][y]=;//可行点都能spfa
}
printf("%d\n",f[F(rt)][(<<k)-]); dfs(rt,(<<k)-);
for(int i=;i<n;++i,printf("\n"))
for(int j=;j<m;++j){
if(!a[i][j]) printf("x");
else if(use[i][j]) printf("o");
else printf("_");
}
return ;
}
bzoj2595 / P4294 [WC2008]游览计划的更多相关文章
- 【BZOJ2595】[Wc2008]游览计划 斯坦纳树
[BZOJ2595][Wc2008]游览计划 Description Input 第一行有两个整数,N和 M,描述方块的数目. 接下来 N行, 每行有 M 个非负整数, 如果该整数为 0, 则该方块为 ...
- 【BZOJ2595】 [Wc2008]游览计划
BZOJ2595 [Wc2008]游览计划 Solution 考虑这是一个最小费用连通性的问题,既然大家都说这是什么斯坦纳树那就是的吧... 所以我们肯定可以这样设一个dp状态: \(dp_{i,j, ...
- luogu P4294 [WC2008]游览计划
LINK:游览计划 斯坦纳树例题. 斯坦纳树是这样一类问题:带权无向图上有K个关键点 求出包含这K个点的最小生成树. 也就是说 求最小生成树 但是 并不是整张图 仅限于K个点. 可以发现我们利用克鲁斯 ...
- 洛谷 P4294 [WC2008]游览计划
题目链接 不是很会呢,但似乎抄了题解后有点明白了 sol:状态DP显然,其实是要构建一棵最小生成树一样的东西,我自己的理解(可能不是很对哦希望多多指教)f[x][y][zt]就是到x,y这个点,状态为 ...
- BZOJ2595:[Wc2008]游览计划——题解(插头dp)
http://www.lydsy.com/JudgeOnline/problem.php?id=2595 Description Input 第一行有两个整数,N和 M,描述方块的数目. 接下来 N行 ...
- P4294 [WC2008]游览计划
传送门 斯坦纳树 给一个联通图,求 $k$ 个关键点联通的最小生成树权值 设 $f[o][i]$ 表示当前关键点选择状态为 $o$ ,以点 $i$ 为根的树的最小权值 初始 $f[1<<( ...
- P4294 [WC2008]游览计划 (斯坦纳树)
题目链接 差不多是斯坦纳树裸题,不过边权化成了点权,这样在合并两棵子树时需要去掉根结点的权值,防止重复. 题目还要求输出解,只要在转移时记录下路径,然后dfs一遍就好了. #include<bi ...
- BZOJ2595 Wc2008 游览计划 【斯坦纳树】【状压DP】*
BZOJ2595 Wc2008 游览计划 Description Input 第一行有两个整数,N和 M,描述方块的数目. 接下来 N行, 每行有 M 个非负整数, 如果该整数为 0, 则该方块为一个 ...
- BZOJ_2595_[Wc2008]游览计划_斯坦纳树
BZOJ_2595_[Wc2008]游览计划_斯坦纳树 题意: 分析: 斯坦纳树裸题,有几个需要注意的地方 给出矩阵,不用自己建图,但枚举子集转移时会算两遍,需要减去当前点的权值 方案记录比较麻烦,两 ...
随机推荐
- git push 报错:missing Change-Id in commit message footer
使用gerrit后,提交代码会出现如下截图问题: 临时解决: step1:把上面红色的那条gitidir复制下来执行下: step2:执行下面的命令会添加change_id git commit -- ...
- [Java in NetBeans] Lesson 08. If: conditional statement
这个课程的参考视频和图片来自youtube. 主要学到的知识点有: 1. If-else statement if (x > 5) { System.out.println("Inpu ...
- SQL Server2008 R2 数据库镜像实施手册(双机)SQL Server2014同样适用
这篇文章主要介绍了SQL Server2008 R2 数据库镜像实施手册(双机)SQL Server2014同样适用,需要的朋友可以参考下 一.配置主备机 1. 服务器基本信息 主机名称为:HOST_ ...
- go https ajax
这个很好用啊,估计大有用武之地 你会喜欢 //https-ajax.go package main import ( "fmt" "io" "net/ ...
- CentOS7 开放服务端口
CentOS 7 默认是firewall防火墙 如果你想让一个web服务可以被其它机子访问,就得开放这个服务的端口,不然就会被拦截 1. 开放端口命令 firewall-cmd --add-port= ...
- ThreadLocal源码实现。
今天在FileInputStream源码中看到有ThreadLocal,之前一直没有理解过这个类,现在进行补充. ThreadLocal即为线程局部变量,它和同步机制处理的是不同的问题域,同步机制是为 ...
- flask orm
#查询 #查第一行 session.query(User.id,User.userName,User.password).first() #查所有行 sessi ...
- clear/reset select2,重置select2,恢复默认
4.0 version //方法一$('#yourButton').on('click', function() { $('#yourfirstSelect2').val(null).trigger( ...
- 基于MySQL提供的Yum repository安装MySQL5.6
基于MySQL提供的Yum repository安装MySQL5.6 1:下载YUM包 1.1切换到root用户: 1)su - root 2)wget http://dev.mysql.com/ge ...
- schame定义及用处
一.schame详解 http://www.cnblogs.com/Neo-ds/p/4790413.html 1.先明确一点,SQL Server中模式(schema)这个概念是在2005的版本里才 ...