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 一看这道题就是一个妙题,然后题解什么树链剖分...珂朵莉树... 还不如并查集来的实在!我们知道并查集本来就是路 ...
随机推荐
- css学习_css3伸缩布局 flex布局
1.flex布局 案例一: 案例二: 保证不至于缩放得太小或太大 案例三:flex的值不一定要写成几份,可以写成固定值 案例四: 竖着3等分(父容器按照高度3等分) !!案例 -----用fle ...
- Ultimate Guide to WeChat for Business 2019
Ultimate Guide to WeChat for Business 2019 By Iaroslav Kudritskiy (source :https://rocketbots.io/blo ...
- WPF气泡样式弹窗效果
页面设计需求,做了一个气泡形状的弹出框,效果如下: 设计思路如下: 1. 使用Path绘制气泡的尖尖,将这个放到顶层: 2. 在用border绘制长方形框,将这个放到底层,并且设置Margin值,使得 ...
- 扩展视图之xpath用法
在视图扩展中,需要定位扩展字段需要显示的位置,通过xpath来实现定位 odoo 视图函数 在整个项目文件中,结构并不是十分明显,虽然它也遵循MVC设计,类比django的MTV模式,各个模块区分的十 ...
- try catch的使用场景
- windows下gitbash中使用zip命令
参考: https://ranxing.wordpress.com/2016/12/13/add-zip-into-git-bash-on-windows/
- selenium各种定位方法(转)
selenium使用 Xpath CSS JavaScript jQuery的定位方法 (治疗selenium各种定位不到,点击不了的并发症) 2017年07月28日 22:47:36 阅读数:369 ...
- 关于最小生成树(并查集)prime和kruskal
适合对并查集有一定理解的人. 新手可能看不懂吧.... 并查集简单点说就是将相关的2个数字联系起来 比如 房子 1 2 3 4 5 6 ...
- java String补足
regionMatches()方法: equals 比较内容 == 比较的是地址
- 微信小程序案例大全
微信小程序demo:足球,赛事分析 小程序简易导航 小程序demo:办公审批 小程序Demo:电魔方 小程序demo:借阅伴侣 微信小程序demo:投票 微信小程序demo:健康生活 小程序demo: ...