BZOJ5279: [Usaco2018 Open]Disruption
题目大意:给你一棵n个节点的树,这n条边称为原边,另给出m条带权值的额外边,求删去每条原边后通过给出的m额外条边变回一棵树的最小价值。
题解:
看完题面以为是Tarjan连通性之类的题目,冷静分析后想到是树链剖分,自己真是Too young too simple。
首先将这棵树进行树链剖分,对于每条额外边x-y,可以作为原树上x-y的路径上的任意一条边删去时的答案,所以路径更新最小值即可。
树链剖分+线段树维护区间最小值,边权转点权的技巧直接把这条边的权值赋到儿子节点上,查找更新时不找LCA即可。。。
输出要求是按原边的顺序,这里用了一个小技巧,大家手模感性理解一下就好了。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<string>
#include<map>
#include<iostream>
#include<queue>
using namespace std;
#define isNum(a) (a>='0'&&a<='9')
#define SP putchar(' ')
#define EL putchar('\n')
#define File(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout)
template<class T1>void read(T1 &r_e_a_d);
template<class T1>void write(T1 w_r_i_t_e);
const int N=50005;
int n,m,len=1,x,y,z,a[N],head[N],idn[N],idw[N];
struct EDGE{
int to,next,id;
}edge[N<<1];
void add(int x,int y,int d){
++len;
edge[len].to=y;edge[len].next=head[x];edge[len].id=d;head[x]=len;
}
int dep[N],son[N],siz[N],fa[N];
void dfs1(int u,int father){
dep[u]=dep[father]+1;fa[u]=father;siz[u]=1;
for (register int i=head[u];i;i=edge[i].next){
int v=edge[i].to;
if (v!=father){
dfs1(v,u);siz[u]+=siz[v];
if (son[u]==-1||siz[v]>siz[son[u]]) son[u]=v;
idn[v]=edge[i].id;
}
}
}
int tot,top[N],dfn[N],seg[N];
void dfs2(int u,int tp){
top[u]=tp;
dfn[u]=++tot;
seg[dfn[u]]=u;
if (son[u]==-1) return ;
dfs2(son[u],tp);
for (register int i=head[u];i;i=edge[i].next){
int v=edge[i].to;
if (v!=fa[u]&&v!=son[u]){
dfs2(v,v);
}
}
}
int tree[50005<<2],lazy[50005<<2];
void build(int k,int l,int r){
tree[k]=lazy[k]=1<<30;
if (l==r) return ;
int mid=l+r>>1;
build(k<<1,l,mid);build(k<<1|1,mid+1,r);
}
void pushdown(int k){
tree[k<<1]=min(tree[k<<1],lazy[k]);
tree[k<<1|1]=min(tree[k<<1|1],lazy[k]);
lazy[k<<1]=min(lazy[k],lazy[k<<1]);
lazy[k<<1|1]=min(lazy[k],lazy[k<<1|1]);
lazy[k]=1<<30;
}
void modify(int k,int l,int r,int x,int y,int z){
if (x<=l&&r<=y){
tree[k]=min(tree[k],z);
lazy[k]=min(lazy[k],z);
return ;
}
if (lazy[k]!=1<<30) pushdown(k);
int mid=l+r>>1;
if (x<=mid) modify(k<<1,l,mid,x,y,z);
if (mid<y) modify(k<<1|1,mid+1,r,x,y,z);
tree[k]=min(tree[k<<1],tree[k<<1|1]);
}
int query(int k,int l,int r,int x,int y){
if (x<=l&&r<=y){
return tree[k];
}
if (lazy[k]!=1<<30) pushdown(k);
int qwq=1<<30,mid=l+r>>1;
if (x<=mid) qwq=min(qwq,query(k<<1,l,mid,x,y));
if (mid<y) qwq=min(qwq,query(k<<1|1,mid+1,r,x,y));
return qwq;
}
void solve(int x,int y,int z){
while (top[x]!=top[y]){
if (dep[top[x]]<dep[top[y]]) x^=y^=x^=y;
modify(1,1,n,dfn[top[x]],dfn[x],z);
x=fa[top[x]];
}
if (dep[x]<dep[y]) x^=y^=x^=y;
modify(1,1,n,dfn[y]+1,dfn[x],z);
}
int main(){
memset (son,-1,sizeof (son));
read(n);read(m);
for (register int i=1;i<n;++i){
read(x);read(y);
add(x,y,i);add(y,x,i);
}
dfs1(1,0);
dfs2(1,1);
build(1,1,n);
for (register int i=1;i<=m;++i){
read(x);read(y);read(z);
solve(x,y,z);
}
for (register int i=2;i<=n;++i){
idw[idn[i]]=i;
}
for (register int i=1;i<n;++i){
int ans=query(1,1,n,dfn[idw[i]],dfn[idw[i]]);
if (ans==1<<30) printf ("-1\n");
else write(ans),EL;
}
return 0;
}
template<class T1>void read(T1 &r_e_a_d){
T1 k=0;
char ch=getchar();
int flag=1;
while(!isNum(ch)){
if(ch=='-'){
flag=-1;
}
ch=getchar();
}
while(isNum(ch)){
k=((k<<1)+(k<<3)+ch-'0');
ch=getchar();
}
r_e_a_d=flag*k;
}
template<class T1>void write(T1 w_r_i_t_e){
if(w_r_i_t_e<0){
putchar('-');
write(-w_r_i_t_e);
}else{
if(w_r_i_t_e<10){
putchar(w_r_i_t_e+'0');
}else{
write(w_r_i_t_e/10);
putchar((w_r_i_t_e%10)+'0');
}
}
}
BZOJ5279: [Usaco2018 Open]Disruption的更多相关文章
- [Usaco2018 Open]Disruption
Description Farmer John自豪于他所经营的交通发达的的农场.这个农场是由N块牧场(2≤N≤50,000)组成的,N-1条双向道路将它们连接起来,每一条道路的都为一单位长度.Farm ...
- 【BZOJ5188】 [Usaco2018 Jan]MooTube
BZOJ5188 [Usaco2018 Jan]MooTube 突然发现BZOJ没有题目,放题面. 题意翻译 题面描述 在业余时间,Farmer John创建了一个新的视频共享服务,他将其命名为Moo ...
- BZOJ5196: [Usaco2018 Feb]Taming the Herd(DP暴力)
5196: [Usaco2018 Feb]Taming the Herd Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 78 Solved: 71[ ...
- BZOJ5194: [Usaco2018 Feb]Snow Boots(排序&set)(可线段树优化)
5194: [Usaco2018 Feb]Snow Boots Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 102 Solved: 79[Subm ...
- BZOJ5280: [Usaco2018 Open]Milking Order(二分+拓扑)
5280: [Usaco2018 Open]Milking Order Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 123 Solved: 62[ ...
- BZOJ5281: [Usaco2018 Open]Talent Show(01分数规划&DP)
5281: [Usaco2018 Open]Talent Show Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 166 Solved: 124[S ...
- BZOJ5484: [Usaco2018 Dec]Sort It Out
5484: [Usaco2018 Dec]Sort It Out https://www.lydsy.com/JudgeOnline/problem.php?id=5484 Sol. 考虑没有在被喊叫 ...
- BZOJ5487: [Usaco2018 Dec]Cowpatibility
Description 研究证明,有一个因素在两头奶牛能否作为朋友和谐共处这方面比其他任何因素都来得重要--她们是不是喜欢同 一种口味的冰激凌!Farmer John的N头奶牛(2≤N≤50,000) ...
- [并查集+LCA USACO18OPEN ] Disruption
https://www.luogu.org/problemnew/show/P4374 一看这道题就是一个妙题,然后题解什么树链剖分...珂朵莉树... 还不如并查集来的实在!我们知道并查集本来就是路 ...
随机推荐
- PHP 简易聊天室 利用redis的订阅发布功能
demo:http://www.200ok.fun:8083/api/chat/list 前言:这个种方式太耗redis连接数,每次订阅都会新起一个进程,仅供练手使用,切勿用于生产环境. 原理:1.P ...
- python遍历文件(替换)
#!/usr/local/bin/python # -*- coding: UTF-8 -*- #coding:gbk import re import os w_str="" x ...
- 学习懈怠的时候,可以运行Qt自带的Demo,或者Delphi控件自带的Demo,或者Cantu书带的源码,运行一下Boost的例子(搞C++不学习Boost/Poco/Folly绝对是一大损失,有需要使用库要第一时间想到)(在六大的痛苦经历说明,我的理论性确实不强,更适合做实践)
这样学还不用动脑子,而且熟悉控件也需要时间,而且慢慢就找到感觉了,就可以精神抖擞的恢复斗志干活了.或者Cantu书带的源码. 并且可以使用Mac SSD运行Qt的Demo,这样运行速度快一点. 此外, ...
- Elasticsearch学习笔记(十四)relevance score相关性评分的计算(1)
一.多shard场景下relevance score不准确问题 1.问题描述: 多个shard下,如果每个shard包含指定搜索条件的document数量不均匀的情况下, ...
- Vue keep-alive如何实现只缓存部分页面
prop: include: 字符串或正则表达式.只有匹配的组件会被缓存. exclude: 字符串或正则表达式.任何匹配的组件都不会被缓存. 在2.1.0版本Vue中 常见用法: // 组件 exp ...
- 循环结构 : while .. for
# ###循环结构 : while .. for ''' while 循环 可以提高代码的效率,减少代码的冗余 while 条件表达式: code1 code2 如果条件表达式成立,返回True , ...
- vi命令复制粘贴
2.复制粘贴 yy :复制当前行 p :粘贴到光标所在行的下一行
- windows下安装配置postgreSQL
1.下载 postgresql-10.4-1-windows-x64.exe 进行安装 2.环境配置(1)文本使用的IDE是VS2010,我们需要配置包含目录(include).库目录(lib).链接 ...
- C# Json解析Json = "{\"EX_RETURN\":[{\"MATNR\":\"test\"}] }";
string jtext = "{\"jiangsu\":[{\"wuxi\":\"无锡\"},{\"suzhou\&q ...
- memcached----------linux下安装memcached,以及php的memcached扩展。
1.通过wget http://www.memcached.org/files/memcached-1.4.24.tar.gz下载最新源码2.解压tar -xf memcached-1.4.24.ta ...