[bzoj2286] [洛谷P2495] [sdoi2015] 消耗战
Description
在一场战争中,战场由 \(n\) 个岛屿和 \(n-1\) 个桥梁组成,保证每两个岛屿间有且仅有一条路径可达。现在,我军已经侦查到敌军的总部在编号为1的岛屿,而且他们已经没有足够多的能源维系战斗,我军胜利在望。已知在其他 \(k\) 个岛屿上有丰富能源,为了防止敌军获取能源,我军的任务是炸毁一些桥梁,使得敌军不能到达任何能源丰富的岛屿。由于不同桥梁的材质和结构不同,所以炸毁不同的桥梁有不同的代价,我军希望在满足目标的同时使得总代价最小。
侦查部门还发现,敌军有一台神秘机器。即使我军切断所有能源之后,他们也可以用那台机器。机器产生的效果不仅仅会修复所有我军炸毁的桥梁,而且会重新随机资源分布(但可以保证的是,资源不会分布到1号岛屿上)。不过侦查部门还发现了这台机器只能够使用 \(m\) 次,所以我们只需要把每次任务完成即可。
Input
第一行一个整数 \(n\),代表岛屿数量。
接下来 \(n-1\) 行,每行三个整数 \(u,v,w\),代表 \(u\) 号岛屿和 \(v\) 号岛屿由一条代价为 \(c\) 的桥梁直接相连,保证 \(1\leq u,v \leq n\) 且$ 1 \leq c \leq 100000$。
第 \(n+1\) 行,一个整数 \(m\),代表敌方机器能使用的次数。
接下来 \(m\) 行,每行一个整数 \(k_i\) ,代表第 \(i\) 次后,有 \(k_i\) 个岛屿资源丰富,接下来 \(k\) 个整数 $ h_1,h_2,…h_k$,表示资源丰富岛屿的编号。
Output
输出有 \(m\) 行,分别代表每次任务的最小代价。
Sample Input
10
1 5 13
1 9 6
2 1 19
2 4 8
2 3 91
5 6 8
7 5 4
7 8 31
10 7 9
3
2 10 6
4 5 7 8 3
3 9 4 6
Sample Output
12
32
22
HINT
对于100%的数据,\(2\leq n \leq 250000,1 \leq m,\sum k_i \leq 500000,1 \leq k_i \leq n-1\)
想法
注意到数据范围中重要的提示 $\sum k_i \leq 500000 $
于是我们对于每次询问把需要用到的点(最多 \(2k\) 个)挑出来建成一棵新树,然后在新树上进行树形dp就行了
说的好容易的样子…实际上“新树”就是“虚树”,它的建树过程是十分巧妙的。
首先把询问点根据原树DFS序排序,显然这些点都要出现在虚树中来,而且为了保证结构不被破坏,另外一些跟他们有关系的点都要加入到虚树中来。我们用一个栈,维护原树上的一条链,自栈底到栈顶,深度由小变大。每次考虑插入询问点进栈。如果插入点的祖先是栈顶元素,那么直接插入即可,因为反正是一条链上的结点。如果不是的话,那么只有可能分居他们的lca的两棵子树中。现在我们就需要分类讨论,如果栈顶元素的下一位的深度比lca深,那么我们需要不断弹出栈顶元素,并且在弹出之前与栈顶下一位连一下边。直到lca深度比栈顶元素深,此时把lca加入栈,把需要插入的点加入栈,继续往下处理。又因为我们是按DFS序做的,这样就可以保证我们开始说的,维护的是树上的一条链。之后再DP就可以了。
——by ljh2000
代码
#include<cstdio>
#include<iostream>
#include<algorithm>
#define INF 1e14
typedef long long ll;
using namespace std;
const int N = 250005;
struct node{
int v,len;
node *next;
}pool[N*2],*h[N],pool2[N*2],*h2[N];
int cnt,cnt2;
void addedge(int u,int v,int l){
node *p=&pool[++cnt],*q=&pool[++cnt];
p->v=v;p->next=h[u];h[u]=p;p->len=l;
q->v=u;q->next=h[v];h[v]=q;q->len=l;
}
void addedge2(int u,int v,int l){
node *p=&pool2[++cnt2];
p->v=v;p->next=h2[u];h2[u]=p;p->len=l;
}
int n,m;
int dfn[N],f[N][20],fmin[N][20],dep[N],tot;
void dfs(int u){
int v;
dfn[u]=++tot;
for(node *p=h[u];p;p=p->next)
if(!dfn[v=p->v]){
f[v][0]=u; fmin[v][0]=p->len;
for(int j=1;j<20;j++){
f[v][j]=f[f[v][j-1]][j-1];
fmin[v][j]=min(fmin[f[v][j-1]][j-1],fmin[v][j-1]);
}
dep[v]=dep[u]+1;
dfs(v);
}
}
int lca(int x,int y){
if(dep[x]<dep[y]) swap(x,y);
for(int i=19;i>=0;i--)
if(dep[f[x][i]]>=dep[y]) x=f[x][i];
if(x==y) return x;
for(int i=19;i>=0;i--)
if(f[x][i]!=f[y][i])
x=f[x][i],y=f[y][i];
return f[x][0];
}
int mindis(int x,int y){
if(dep[x]<dep[y]) swap(x,y);
if(x==y) return 0;
int t=N;
for(int i=19;i>=0;i--)
if(dep[f[x][i]]>=dep[y]){
t=min(t,fmin[x][i]);
x=f[x][i];
}
return t;
}
ll w[N];
int g[N];
void dp(int u){
int v;
w[u]=0;
for(node *p=h2[u];p;p=p->next){
dp(v=p->v);
w[u]+=min(w[v],(ll)p->len);
}
h2[u]=NULL; //clear 保证每次邻接矩阵清空效率为O(k)
if(g[u]) w[u]=INF;
}
int st[N],top;
int K,k[N];
bool cmp(int x,int y) { return dfn[x]<dfn[y]; }
void work(){
cnt2=0; //clear
sort(k,k+K,cmp);
//build tree
top=0; st[top++]=1;
for(int i=0;i<K;i++){
int now=k[i],f=lca(now,st[top-1]);
while(1){
if(top<=1) { if(st[top-1]!=f) st[top++]=f; break; }
if(dep[f]>=dep[st[top-2]]){
if(f!=st[top-1]) addedge2(f,st[top-1],mindis(f,st[top-1]));
if(f==st[top-2]) top--;
else st[top-1]=f;
break;
}
addedge2(st[top-2],st[top-1],mindis(st[top-1],st[top-2]));
top--;
}
if(now!=st[top-1]) st[top++]=now;
}
while(top>1) {
addedge2(st[top-2],st[top-1],mindis(st[top-1],st[top-2]));
top--;
}
dp(1);
printf("%lld\n",w[1]);
}
int main()
{
int u,v,l;
scanf("%d",&n);
for(int i=1;i<n;i++){
scanf("%d%d%d",&u,&v,&l);
addedge(u,v,l);
}
dep[1]=1; dfs(1);
scanf("%d",&m);
while(m--){
scanf("%d",&K);
for(int i=0;i<K;i++) scanf("%d",&k[i]),g[k[i]]=1;
work();
for(int i=0;i<K;i++) g[k[i]]=0;
}
return 0;
}
[bzoj2286] [洛谷P2495] [sdoi2015] 消耗战的更多相关文章
- 洛谷 P2495 [SDOI2011]消耗战(虚树,dp)
题面 洛谷 题解 虚树+dp 关于虚树 了解一下 具体实现 inline void insert(int x) { if (top == 1) {s[++top] = x; return ;} int ...
- 洛谷P2495 [SDOI2011]消耗战(虚树dp)
P2495 [SDOI2011]消耗战 题目链接 题解: 虚树\(dp\)入门题吧.虚树的核心思想其实就是每次只保留关键点,因为关键点的dfs序的相对大小顺序和原来的树中结点dfs序的相对大小顺序都是 ...
- ●洛谷P2495 [SDOI2011]消耗战
题链: https://www.luogu.org/problemnew/show/P2495题解: 虚树入门,树形dp 推荐博客:http://blog.csdn.net/lych_cys/arti ...
- 洛谷P2495 [SDOI2011]消耗战(虚树)
题面 传送门 题解 为啥一直莫名其妙\(90\)分啊--重构了一下代码才\(A\)掉-- 先考虑直接\(dp\)怎么做 树形\(dp\)的时候,记一下断开某个节点的最小值,就是从根节点到它的路径上最短 ...
- [洛谷P2495][SDOI2011]消耗战
题目大意:有一棵$n(n\leqslant2.5\times10^5)$个节点的带边权的树,$m$个询问,每次询问给出$k(\sum\limits_{i=1}^mk_i\leqslant5\times ...
- 洛谷P3324 [SDOI2015]星际战争
题目:洛谷P3324 [SDOI2015]星际战争 思路: 类似<导弹防御塔>,因为题目保证有解,花费时间小于最终答案时一定无法消灭所有敌人,只要花费时间大于等于最终答案都可以消灭所有敌人 ...
- LOJ #2185 / 洛谷 P3329 - [SDOI2015]约数个数和(莫比乌斯函数)
LOJ 题面传送门 / 洛谷题面传送门 题意: 求 \(\sum\limits_{i=1}^n\sum\limits_{j=1}^md(ij)\),\(d(x)\) 为 \(x\) 的约数个数. \( ...
- AC日记——[SDOI2011]消耗战 洛谷 P2495
[SDOI2011]消耗战 思路: 建虚树走树形dp: 代码: #include <bits/stdc++.h> using namespace std; #define INF 1e17 ...
- [虚树模板] 洛谷P2495 消耗战
题意:给定树上k个点,求切断这些点到根路径的最小代价.∑k <= 5e5 解:虚树. 构建虚树大概是这样的:设加入点与栈顶的lca为y,比较y和栈中第二个元素的DFS序大小关系. 代码如下: i ...
随机推荐
- JMETER+JENKINS接口测试持续集成
FIDDER+ANT+JENKINS+JMETER+SVN+tomcat接口测试集成 操作流程: 1.测试人员通过FIDDER过滤抓取接口调用信息,导出成jmx文件.(jmeter支持命令行方式调用j ...
- Linux 内核存取 I/O 和内存空间
一个 PCI 设备实现直至 6 个 I/O 地址区. 每个区由要么内存要么 I/O 区组成. 大部分 设备实现它们的 I/O 寄存器在内存区中, 因为通常它是一个完善的方法(如同在" I/O ...
- Jenkins安装部署与使用
一.Jenkins平台安装部署 Jenkins官网免费获取Jenkins软件,官网地址为:http://mirrors.jenkins-ci.org/下载稳定的Jenkins版本.由于Jenkins是 ...
- apache WEB服务器安装(包括虚拟主机)
一.apache下载编译安装 yum install apr apr-devel apr-util apr-util-devel gcc-c++ wget tar -y cd /usr/src wge ...
- Keras lstm 文本分类示例
#基于IMDB数据集的简单文本分类任务 #一层embedding层+一层lstm层+一层全连接层 #基于Keras 2.1.1 Tensorflow 1.4.0 代码: '''Trains an LS ...
- Struts2 注释类型
Struts 2 应用程序可以使用Java5注释作为替代XML和Java属性配置.这里是清单的不同的类别有关的最重要的注解: 命名空间注释(动作注释): @ Namespace注释允许在Action类 ...
- Iptables-linux服务器做路由转发
https://blog.csdn.net/liang_operations/article/details/80747510 实现内部服务器C可以经过服务器B进行上网. 3.1服务器双网卡,一块配置 ...
- 第二阶段:4.商业需求文档MRD:4.PRD-用例和规则
类似之前的泳道图 可以在下面添加一些描述 有时候用图还是会有一些限制 不能够很好的表达
- 清晰架构(Clean Architecture)的Go微服务: 设计原则
我最近写了一个Go微服务应用程序,这个程序的设计来自三个灵感: 清晰架构"Clean Architecture"¹ and SOLID (面向对象设计)² 设计 原则³ Sprin ...
- 后缀自动机&回文自动机学习笔记
在学了一天其实是边学边摆之后我终于大概$get$后缀自动机了,,,就很感动,于是时隔多年我终于决定再写篇学习笔记辽$QwQ$ $umm$和$FFT$学习笔记一样,这是一篇单纯的$gql$的知识总结博, ...