题目描述

输入

第一行包含一个正整数N,表示X国的城市个数. 第二行包含两个正整数L和U,表示政策要求的第一期重建方案中修建道路数的上下限 接下来的N-1行描述重建小组的原有方案,每行三个正整数Ai,Bi,Vi分别表示道路(Ai,Bi),其价值为Vi 其中城市由1..N进行标号

输出

输出最大平均估值,保留三位小数

样例输入

4
2 3
1 2 1
1 3 2
1 4 3

样例输出

2.500

提示

N<=100000,1<=L<=U<=N-1,Vi<=1000000

这题算是长链剖分最经典的一道题了。

首先要找一个最优方案,直接DP显然不好做,那么考虑二分答案然后验证,因为是浮点数,要注意精度问题。

假设当前二分的答案是k,判断答案是否满足时原式也就转化成了$\frac{\sum_{ }^{ }vi}{|S|}>=k$,将分母移到不等式右边得到$\sum_{ }^{ }vi>=k*|S|$

将右边的部分移到左边就变成了$\sum_{ }^{ }vi-k*|S|>=0$

因为vi的个数就是|S|,因此将k*|S|放到Σ里面,判断就变成了$\sum_{ }^{ }(vi-k)>=0$

每次判断只要把每条边的边权和减k,再判断能否有一条路径边权和大于等于0就好了。

怎么判断呢?很容易想到O(n^2)dp,设f[i][j]表示i子树中与i距离为j的链的边权和最大值,枚举另一棵子树找到链长在[L-j,R-j]之内的边权最大值。

O(n^2)dp显然不行,但观察到dp是可合并的以深度为下标的转移方程,因此可以用长链剖分优化成O(n)。

怎么优化呢?

首先对整棵树长链剖分求出树剖序,再把树剖序架到线段树上,因为整棵树是由所有长链组成的,每条长链因为是优先遍历所以在树剖序上是连续的一段。

也就是说树剖序上的每一段都是树剖出的一条长链。那么每个点子树中每个深度的信息就可以都存到这个点往下的长链上。

当做树形DP时,每个点对于重儿子回溯时不做任何操作,直接继承;当轻儿子回溯时枚举轻儿子每个深度的边权和最大值,在长链上找到对应区间求最大值来更新答案。

然后再把这个轻儿子的信息合并到长链上。因为长链上存的是之前所有遍历过的子树合并后的信息,所以相当于每个点子树中有用的信息都在这个点往下的长链上。最后别忘了考虑从上到下的每条直链。

这样DP是O(n)的,再加上线段树的O(log)和二分的O(log)一共是O(nlog2n)。

至于这样DP为什么是O(n)的?

因为每个点对重儿子是直接继承的,而每个点需要被DP当且仅当它是轻儿子时,这时它一定是一个长链的链头,DP的时间复杂度是这个点往下的长链长度,那么DP的总复杂度就是每条长链的链长总和,也就是O(n)。

更加详细的有关长链剖分的讲解参见->长链剖分

#include<set>
#include<map>
#include<cmath>
#include<queue>
#include<stack>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int n;
int L,R;
int tot;
int cnt;
int x,y,z;
double ans;
int s[100010];
int d[100010];
int f[100010];
int v[200010];
int mn[100010];
int to[200010];
int son[100010];
int num[100010];
int head[100010];
int next[200010];
double mx[800010];
double val[200010];
double dis[100010];
double dep[100010];
double INF=1e15;
double eps=1e-4;
void add(int x,int y,int z)
{
tot++;
next[tot]=head[x];
head[x]=tot;
to[tot]=y;
v[tot]=z;
}
void dfs(int x)
{
d[x]=d[f[x]]+1;
mn[x]=d[x];
for(int i=head[x];i;i=next[i])
{
if(to[i]!=f[x])
{
f[to[i]]=x;
dfs(to[i]);
mn[x]=max(mn[x],mn[to[i]]);
if(mn[to[i]]>mn[son[x]])
{
son[x]=to[i];
}
}
}
}
void dfs2(int x)
{
s[x]=++cnt;
if(son[x])
{
dfs2(son[x]);
}
for(int i=head[x];i;i=next[i])
{
if(to[i]!=f[x]&&to[i]!=son[x])
{
dfs2(to[i]);
}
}
}
void build(int rt,int l,int r)
{
mx[rt]=-INF;
if(l==r)
{
num[l]=rt;
return ;
}
int mid=(l+r)>>1;
build(rt<<1,l,mid);
build(rt<<1|1,mid+1,r);
}
void change(int rt,int l,int r,int k,double x)
{
if(l==r)
{
mx[rt]=max(mx[rt],x);
return ;
}
int mid=(l+r)>>1;
if(k<=mid)
{
change(rt<<1,l,mid,k,x);
}
else
{
change(rt<<1|1,mid+1,r,k,x);
}
mx[rt]=max(mx[rt<<1],mx[rt<<1|1]);
}
double query(int rt,int l,int r,int L,int R)
{
if(L>R)
{
return -INF;
}
if(L<=l&&r<=R)
{
return mx[rt];
}
int mid=(l+r)>>1;
double res=-INF;
if(L<=mid)
{
res=max(res,query(rt<<1,l,mid,L,R));
}
if(R>mid)
{
res=max(res,query(rt<<1|1,mid+1,r,L,R));
}
return res;
}
void tree_dp(int x)
{
change(1,1,n,s[x],dis[x]);
for(int i=head[x];i;i=next[i])
{
if(son[x]==to[i])
{
dis[son[x]]=dis[x]+val[i];
tree_dp(son[x]);
}
}
for(int i=head[x];i;i=next[i])
{
if(to[i]!=f[x]&&to[i]!=son[x])
{
dis[to[i]]=dis[x]+val[i];
tree_dp(to[i]);
for(int j=1;j<=mn[to[i]]-d[x];j++)
{
dep[j]=mx[num[s[to[i]]+j-1]];
if(j<=R)
{
ans=max(ans,query(1,1,n,max(s[x],s[x]+L-j),min(s[x]+mn[x]-d[x],s[x]+R-j))+dep[j]-2*dis[x]);
}
}
for(int j=1;j<=mn[to[i]]-d[x];j++)
{
change(1,1,n,s[x]+j,dep[j]);
}
}
}
ans=max(ans,query(1,1,n,s[x]+L,min(s[x]+mn[x]-d[x],s[x]+R))-dis[x]);
}
int main()
{
scanf("%d%d%d",&n,&L,&R);
for(int i=1;i<n;i++)
{
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);
add(y,x,z);
}
dfs(1);
dfs2(1);
double l=0;
double r=1000000;
while(r-l>eps)
{
double mid=(l+r)/2;
for(int i=1;i<=tot;i++)
{
val[i]=v[i]-mid;
}
ans=-INF;
dis[1]=0;
build(1,1,n);
tree_dp(1);
if(ans<0)
{
r=mid;
}
else
{
l=mid;
}
}
printf("%.3lf",l);
}

BZOJ1758[Wc2010]重建计划——分数规划+长链剖分+线段树+二分答案+树形DP的更多相关文章

  1. 2019.01.21 bzoj1758: [Wc2010]重建计划(01分数规划+长链剖分+线段树)

    传送门 长链剖分好题. 题意简述:给一棵树,问边数在[L,R][L,R][L,R]之间的路径权值和与边数之比的最大值. 思路: 用脚指头想都知道要01分数规划. 考虑怎么checkcheckcheck ...

  2. BZOJ.1758.[WC2010]重建计划(分数规划 点分治 单调队列/长链剖分 线段树)

    题目链接 BZOJ 洛谷 点分治 单调队列: 二分答案,然后判断是否存在一条长度在\([L,R]\)的路径满足权值和非负.可以点分治. 对于(距当前根节点)深度为\(d\)的一条路径,可以用其它子树深 ...

  3. BZOJ 1758 / Luogu P4292 [WC2010]重建计划 (分数规划(二分/迭代) + 长链剖分/点分治)

    题意 自己看. 分析 求这个平均值的最大值就是分数规划,二分一下就变成了求一条长度在[L,R]内路径的权值和最大.有淀粉质的做法但是我没写,感觉常数会很大.这道题可以用长链剖分做. 先对树长链剖分. ...

  4. [WC2010]重建计划(长链剖分+线段树+分数规划)

    看到平均值一眼分数规划,二分答案mid,边权变为w[i]-mid,看是否有长度在[L,R]的正权路径.设f[i][j]表示以i为根向下j步最长路径,用长链剖分可以优化到O(1),查询答案线段树即可,复 ...

  5. bzoj 1758 [Wc2010]重建计划 分数规划+树分治单调队列check

    [Wc2010]重建计划 Time Limit: 40 Sec  Memory Limit: 162 MBSubmit: 4345  Solved: 1054[Submit][Status][Disc ...

  6. 2019.03.11 COGS2652 秘术(天文密葬法)(分数规划+长链剖分)

    传送门 题意:nnn个点的树,每个点两个值a,ba,ba,b,问长度为mmm的路径∑ai∑bi\frac{\sum a_i}{\sum b_i}∑bi​∑ai​​的最大值. 思路:一眼要01分数规划, ...

  7. 2018牛客网暑假ACM多校训练赛(第七场)I Tree Subset Diameter 动态规划 长链剖分 线段树

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round7-I.html 题目传送门 -  https://www.n ...

  8. BZOJ.3653.谈笑风生(长链剖分/线段树合并/树状数组)

    BZOJ 洛谷 \(Description\) 给定一棵树,每次询问给定\(p,k\),求满足\(p,a\)都是\(b\)的祖先,且\(p,a\)距离不超过\(k\)的三元组\(p,a,b\)个数. ...

  9. BZOJ.3252.攻略(贪心 长链剖分/线段树)

    题目链接 贪心,每次选价值最大的一条到根的链.比较显然(不选白不选). 考虑如何维护这个过程.一个点的价值选了就没有了,而它只会影响它子树里的点,可以用DFS序+线段树修改.而求最大值也可以用线段树. ...

随机推荐

  1. maven 打war包tomcat服务器乱码问题

    今天用maven3的命令打war包,命令是mvn clean package -Dmaven.test.skip=true,打包后放在tomcat跑起来后发现tomcat的日志出现乱码. 后来在pom ...

  2. AC_Dream 1224 Robbers(贪心)

    题意:n个抢劫犯分别抢到的金钱是k1, k2, k3,...,一共得到的金钱是m, 但是在分钱的时候是按照x1/y, x2/y, x3/y,....的比例进行分配的!这样的话 一些抢劫犯就会觉得不公平 ...

  3. 扫雷游戏制作过程(C#描述):第五节、菜单操作(续)

    前言 这里给出教程原文地址. 该项目已经放在github上托管. 发布版已经分享到百度网盘 菜单操作(续) 接着节前一章节的内容,我们继续完善菜单栏的功能. 我们首先,先完善Rank的选项,我们希望我 ...

  4. javascript执行机制

    文的目的就是要保证你彻底弄懂javascript的执行机制,如果读完本文还不懂,可以揍我. 不论你是javascript新手还是老鸟,不论是面试求职,还是日常开发工作,我们经常会遇到这样的情况:给定的 ...

  5. [PA2014]Muzeum

    [PA2014]Muzeum 题目大意: 有\(n\)件展品和\(m\)个警卫,每件展品有一个坐标\((x_i,y_i)\)和价值\(v_i\),每个警卫的坐标为\((x_i,y_i)\).每个警卫面 ...

  6. thinkphp 检测验证码

    /** * 检测验证码 * @param integer $id 验证码ID * @return boolean 检测结果 */function check_verify($code, $id = 1 ...

  7. jstack 查看线程状态

    使用jstack pid命令可以查看JVM的线程状态,其中值得关注的线程状态有:死锁,Deadlock(重点关注)执行中,Runnable等待资源,Waiting on condition(重点关注) ...

  8. Spring---七大核心模块

    核心容器(Spring Core) 核心容器提供Spring框架的基本功能.Spring以bean的方式组织和管理Java应用中的各个组件及其关系.Spring使用BeanFactory来产生和管理B ...

  9. GB2312汉字区位码、交换码和机内码转换方法 (ZT)

    GB2312汉字区位码.交换码和机内码转换方法 (ZT) 为了适应计算机处理汉字信息的需要,1981年我国颁布了GB2312国家标准.该标准选出6763个常用汉字(其中,一级常用汉字3755个,二级汉 ...

  10. Java学习笔记三十:Java小项目之租车系统

    Java小项目之租车系统 一:项目背景介绍: 根据所学知识,编写一个控制台版的“呱呱租车系统” 功能: 1.展示所有可租车辆: 2.选择车型.租车量: 3.展示租车清单,包含:总金额.总载货量以及其车 ...