codesforces 671D Roads in Yusland
Mayor of Yusland just won the lottery and decided to spent money on something good for town. For example, repair all the roads in the town.
Yusland consists of n intersections connected by n - 1 bidirectional roads. One can travel from any intersection to any other intersection using only these roads.
There is only one road repairing company in town, named "RC company". Company's center is located at the intersection 1. RC company doesn't repair roads you tell them. Instead, they have workers at some intersections, who can repair only some specific paths. The i-th worker can be paid ci coins and then he repairs all roads on a path from ui to some vi that lies on the path from ui to intersection 1.
Mayor asks you to choose the cheapest way to hire some subset of workers in order to repair all the roads in Yusland. It's allowed that some roads will be repaired more than once.
If it's impossible to repair all roads print - 1.
The first line of the input contains two integers n and m (1 ≤ n, m ≤ 300 000) — the number of cities in Yusland and the number of workers respectively.
Then follow n−1 line, each of them contains two integers xi and yi (1 ≤ xi, yi ≤ n) — indices of intersections connected by the i-th road.
Last m lines provide the description of workers, each line containing three integers ui, vi and ci (1 ≤ ui, vi ≤ n, 1 ≤ ci ≤ 109). This means that the i-th worker can repair all roads on the path from vi to ui for ci coins. It's guaranteed that vi lies on the path from ui to 1. Note that vi and ui may coincide.
If it's impossible to repair all roads then print - 1. Otherwise print a single integer — minimum cost required to repair all roads using "RC company" workers.
6 5
1 2
1 3
3 4
4 5
4 6
2 1 2
3 1 4
4 1 3
5 3 1
6 3 2
8
In the first sample, we should choose workers with indices 1, 3, 4 and 5, some roads will be repaired more than once but it is OK. The cost will be equal to 2 + 3 + 1 + 2 = 8 coins.
设f[i]为控制了i点的子树的费用
显然f是存在后效性的,因为当前选的工人会影响祖先的状态
最暴力的方法f[x][i]表示工人最高控制到x,控制了i的子树,不过显然超时
想办法消除后效性
我们考虑把一些工人的属性叠加到一个工人上,那么我们最后查询时只要查询某些能维修完整棵树的最小值
因为我们已经的到了这几棵子树的最优答案,并且子树之间互不影响,那么我们考虑把控制其他子树费用叠加到一个子树工人上
假设儿子有p1,p2,p3,tot=∑f[p]
p1的儿子的工人加上tot-f[p1],其他同理
也就是说,每一个工人都有控制其他所有子树的花费
也就是说,祖先x要用到p儿子这个工人时,就可以直接算出控制p儿子及x的费用
而且没有后效性,因为工人存了包含此工人的最优解,而且工人一直到vi被消掉之前,
都可以覆盖当前点,所以把工人都计算,肯定能算出最优解
因为工人在到vi之前都需要
而且对于计算一个节点的f,我们要给开始的工人赋值,给结束的工人消值,且给子节点所有工人加数
所以我们用dfs序维护以工人为下标的线段树,这样给子节点的工人加数等于线段树的区间加法
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
struct Node
{
int next,to;
}edge[],edgeS[],edgeT[];
int num,numS,numT,head[],headS[],headT[];
int L[],R[],cnt,dfn[],n,m;
ll c[],lazy[],val[],inf=2e15,f[];
void add(int u,int v)
{
num++;
edge[num].next=head[u];
head[u]=num;
edge[num].to=v;
}
void addS(int u,int v)
{
numS++;
edgeS[numS].next=headS[u];
headS[u]=numS;
edgeS[numS].to=v;
}
void addT(int u,int v)
{
numT++;而且没有后效性,因为
edgeT[numT].next=headT[u];
headT[u]=numT;
edgeT[numT].to=v;
}
void pushup(int rt)
{
c[rt]=min(c[rt*],c[rt*+]);
}
void pushdown(int rt)
{
if (lazy[rt])
{
c[rt*]+=lazy[rt];
c[rt*+]+=lazy[rt];
lazy[rt*]+=lazy[rt];
lazy[rt*+]+=lazy[rt];
lazy[rt]=;
}
}
void build(int rt,int l,int r)
{
c[rt]=inf;
if (l==r)
return;
int mid=(l+r)/;
build(rt*,l,mid);
build(rt*+,mid+,r);
}
void prework(int x,int pa)
{int i;
L[x]=cnt+;
for (i=headS[x];i;i=edgeS[i].next)
{
int v=edgeS[i].to;
dfn[v]=++cnt;
}
for (i=head[x];i;i=edge[i].next)
{
int v=edge[i].to;
if (v!=pa)
prework(v,x);
}
R[x]=cnt;
}
void update(int rt,int l,int r,int x,ll d)
{
if (l==r)
{
c[rt]=min(inf,d);
return;
}
pushdown(rt);
int mid=(l+r)/;
if (x<=mid) update(rt*,l,mid,x,d);
else update(rt*+,mid+,r,x,d);
pushup(rt);
}
void change(int rt,int l,int r,int L,int R,ll d)
{
if (L>R) return;
if (l>=L&&r<=R)
{
c[rt]+=d;
c[rt]=min(c[rt],inf);
lazy[rt]+=d;
lazy[rt]=min(lazy[rt],inf);
return;
}
pushdown(rt);
int mid=(l+r)/;
if (L<=mid) change(rt*,l,mid,L,R,d);
if (R>mid) change(rt*+,mid+,r,L,R,d);
pushup(rt);
}
ll query(int rt,int l,int r,int L,int R)
{
if (L>R) return inf;
if (l>=L&&r<=R)
{
return c[rt];
}
pushdown(rt);
int mid=(l+r)/;
ll s=inf;
if (L<=mid) s=min(s,query(rt*,l,mid,L,R));
if (R>mid) s=min(s,query(rt*+,mid+,r,L,R));
pushup(rt);
return s;
}
void dfs(int x,int pa)
{int i;
ll tot=;
for (i=head[x];i;i=edge[i].next)
{
int v=edge[i].to;
if (v==pa) continue;
dfs(v,x);
tot+=f[v];
tot=min(tot,inf);
}
if (x==)
{
f[]=tot;
return;
}
for (i=headS[x];i;i=edgeS[i].next)
{
int v=edgeS[i].to;
update(,,m,dfn[v],tot+val[v]);
}
for (i=headT[x];i;i=edgeT[i].next)
{
int v=edgeT[i].to;
update(,,m,dfn[v],inf);
}
for (i=head[x];i;i=edge[i].next)
{
int v=edge[i].to;
if (v==pa) continue;
change(,,m,L[v],R[v],tot-f[v]);
}
f[x]=query(,,m,L[x],R[x]);
}
int main()
{int i,u,v;
cin>>n>>m;
for (i=;i<=n-;i++)
{
scanf("%d%d",&u,&v);
add(u,v);add(v,u);
}
for (i=;i<=m;i++)
{
scanf("%d%d%I64d",&u,&v,&val[i]);
addS(u,i);addT(v,i);
}
prework(,);
build(,,m);
dfs(,);
if (f[]!=inf) cout<<f[];
else cout<<-;
}
codesforces 671D Roads in Yusland的更多相关文章
- Codeforces 671D Roads in Yusland [树形DP,线段树合并]
洛谷 Codeforces 这是一个非正解,被正解暴踩,但它还是过了. 思路 首先很容易想到DP. 设\(dp_{x,i}\)表示\(x\)子树全部被覆盖,而且向上恰好延伸到\(dep=i\)的位置, ...
- codeforces 671D Roads in Yusland & hdu 5293 Tree chain problem
dp dp优化 dfs序 线段树 算是一个套路.可以处理在树上取链的问题.
- Codeforces 671D. Roads in Yusland(树形DP+线段树)
调了半天居然还能是线段树写错了,药丸 这题大概是类似一个树形DP的东西.设$dp[i]$为修完i这棵子树的最小代价,假设当前点为$x$,但是转移的时候我们不知道子节点到底有没有一条越过$x$的路.如果 ...
- CF 671D Roads in Yusland
弄完之后点进去一看,竟然是div1的D题……最近真是天天被题虐哭 推荐这一篇博客 https://www.cnblogs.com/Sakits/p/8085598.html 感觉讲清楚了,也是基本照着 ...
- Codeforces 671 D. Roads in Yusland
题目描述 Mayor of Yusland just won the lottery and decided to spent money on something good for town. Fo ...
- 【CF671D】Roads in Yusland(贪心,左偏树)
[CF671D]Roads in Yusland(贪心,左偏树) 题面 洛谷 CF 题解 无解的情况随便怎么搞搞提前处理掉. 通过严密(大雾)地推导后,发现问题可以转化成这个问题: 给定一棵树,每条边 ...
- [Codeforces671D]Roads in Yusland
[Codeforces671D]Roads in Yusland Tags:题解 题意 luogu 给定以1为根的一棵树,有\(m\)条直上直下的有代价的链,求选一些链把所有边覆盖的最小代价.若无解输 ...
- 【CF617D】Roads in Yusland
[CF617D]Roads in Yusland 题面 蒯的洛谷的 题解 我们现在已经转化好了题目了,戳这里 那么我们考虑怎么求这个东西,我们先判断一下是否所有的边都能被覆盖,不行的话输出\(-1\) ...
- 【CodeForces】671 D. Roads in Yusland
[题目]D. Roads in Yusland [题意]给定n个点的树,m条从下往上的链,每条链代价ci,求最少代价使得链覆盖所有边.n,m<=3*10^5,ci<=10^9,time=4 ...
随机推荐
- Beta 第六天
今天遇到的困难: github服务器响应很慢 推图的API接口相应较慢,超过了初始设定的最大延迟时间,导致了无法正确返回图片 ListView滑动删除Demo出现了某些Bug,这些Bug可能导致了某些 ...
- 20162328蔡文琛week05
学号 20162328 <程序设计与数据结构>第X周学习总结 教材学习内容总结 面向对象程序设计的核心是类的定义,它代表定义了状态和行为的对象. 变量的作用域依赖于变量声明的位置,作用域决 ...
- TCP和UDP的最完整的区别
TCP UDP TCP与UDP基本区别 1.基于连接与无连接 2.TCP要求系统资源较多,UDP较少: 3.UDP程序结构较简单 4.流模式(TCP)与数据报模式(UDP); ...
- Swift -欢迎界面1页, 延长启动图片的显示时间(LaunchImage)
转自:http://www.hangge.com/blog/cache/detail_1238.html http://www.hangge.com/blog/cache/detail_672.htm ...
- verilog学习笔记(1)_两个小module
第一个小module-ex_module module ex_module( input wire sclk,//声明模块的时候input变量一定是wire变量 input wire rst_n,// ...
- 乐动力APP案例
第一部分 调研, 评测 下载软件并使用起来,描述最简单直观的个人第一次上手体验. 这款软件的主界面功能还是比较完善,里面有多个关于运动相关的数据,还有一些推荐健身教程,记录功能也十分不错,其中最难理解 ...
- js判断flash文件是否加载完毕
轮询判断加载进度 img的加载完成有onload方法,一直不知道该怎么判断swf文件是否加载完毕了? 在应用中使用了轮询判断加载进度值PercentLoaded是否达到100,经测试,可以达到效果. ...
- ESP8266 wifi 模块配置,Wechat+APP控制实现
首先刷入安信可的AiCloud 2.0 SDK文件,AiCloud 2.0具体信息参见AiCloud 1.0 和AiCloud 2.0对比 APP见如下二维码下载. 1.安信可AiCloud 2.0 ...
- CentOS 7 PHP-redis扩展安装,浏览器不显示数据及redis无法储存数据常见问题解决办法
首先使用php -m 可以查看到自己安装了那些扩展. 1.使用wget下载redis压缩包 wget https://github.com/phpredis/phpredis/archive/deve ...
- “认证发布”和“获取展示”,如何在 SharePoint 中正确使用 RSS Feed。
在我们进行的日常工作中,是由一部分信息需要 Share 给其他人或者组织的.SharePoint 虽然支持在某个 Site Collection 中互通信息,但是跨 Site Collection 的 ...