传送门

Description

Solution

时隔多年,补上了这题的长链剖分写法

感觉比点分治要好写的多

我们假设\(pos\)是当前点的\(dfn\),它距离所在链的底端的边的数量是\(len\),距离是\(siz\)

那么我们要求得\(g[pos...pos+len]\)

其中\(g[pos+i]+siz\)表示的是当前点往下长度为\(i\)的最长链的大小

初始情况下,\(g[pos]=-siz[pos]\)

为什么要这么写呢?

因为转移重儿子的时候,我们直接把数组右移了一位,这样子定义使得原先的值仍然有意义

考虑如何转移?

对于一个轻儿子,\(dfn\)为\(pv\)

那么\(当前点,轻儿子g[pos+j+1]\leftarrow siz[pv]+g[pv+j]+w(当前点,轻儿子)-siz\)

考虑如何更新答案

对于\(lca\)为当前点的链

  1. 首先计算一个端点就是当前点的:\((g[pos+i])_{min}+siz\)

  2. 然后是两个端点分别处于不同子树的

    枚举一个轻儿子子树内的链长:\(当前点,轻儿子(g[pos+i])_{min}+siz+(siz[pv]+g[pv+j])+(w(当前点,轻儿子))\)

求min的部分用线段树优化

Code 

#include<bits/stdc++.h>
#define ll long long
#define reg register
#define ri reg int i
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*f;
}
const int MN=1e5+5;
const double eps=1e-4,inf=0x3f3f3f3f;
double siz[MN],g[MN],ans,mid;
struct edge{int to,w,nex;}e[MN<<1];int hr[MN],en;
inline void ins(int x,int y,int w)
{
e[++en]=(edge){y,w,hr[x]};hr[x]=en;
e[++en]=(edge){x,w,hr[y]};hr[y]=en;
}
int dfn[MN],len[MN],mx[MN],w[MN],ind,n,L,R;
struct SegmentTree
{
#define ls (x<<1)
#define rs (x<<1|1)
#define Mid ((l+r)>>1)
double t[MN<<2];
inline void clear(){for(ri=0;i<(MN<<2);++i)t[i]=-inf;}
void Modify(int x,int l,int r,int a,double val)
{
if(l==r) return(void)(t[x]=max(t[x],val));
if(a<=Mid) Modify(ls,l,Mid,a,val);
else Modify(rs,Mid+1,r,a,val);
t[x]=max(t[ls],t[rs]);
}
double Query(int x,int l,int r,int a,int b)
{
if(l==a&&r==b) return t[x];if(a>b)return -inf;
if(b<=Mid) return Query(ls,l,Mid,a,b);
if(a>Mid) return Query(rs,Mid+1,r,a,b);
return max(Query(ls,l,Mid,a,Mid),Query(rs,Mid+1,r,Mid+1,b));
}
}T;
void dfs(int x,int f=0)
{
for(ri=hr[x];i;i=e[i].nex)if(e[i].to^f)
{
dfs(e[i].to,x);
if(len[e[i].to]>=len[mx[x]]) mx[x]=e[i].to,w[x]=e[i].w;
if(len[e[i].to]+1>len[x]) len[x]=len[e[i].to]+1;
}
}
void solve(int x,int f=0)
{
if(!dfn[x]) dfn[x]=++ind;
reg int i,j,pos=dfn[x];
if(mx[x]) solve(mx[x],x),siz[pos]=siz[pos+1]+w[x]-mid;
T.Modify(1,1,n,pos,g[pos]=-siz[pos]);
if(L<=len[x])
{
double tmp=T.Query(1,1,n,pos+L,pos+min(len[x],R));
ans=max(ans,tmp+siz[pos]);
}
for(i=hr[x];i;i=e[i].nex)if(e[i].to!=f&&e[i].to!=mx[x])
{
solve(e[i].to,x);reg int pv=dfn[e[i].to];
for(j=0;j<=len[e[i].to];++j)
{
double tmp=T.Query(1,1,n,pos+max(0,L-j-1),pos+min(len[x],R-j-1));
ans=max(ans,tmp+siz[pv]+siz[pos]+g[pv+j]+e[i].w-mid);
}
for(j=0;j<=len[e[i].to];++j)
{
double tmp=siz[pv]+g[pv+j]+e[i].w-mid-siz[pos];
if(tmp>g[pos+j+1]) T.Modify(1,1,n,pos+j+1,g[pos+j+1]=tmp);
}
}
}
bool check(){T.clear();ans=-inf;solve(1);return ans>=eps;}
int main()
{
n=read();L=read();R=read();
reg int i,x,y;
for(i=1;i<n;++i) x=read(),y=read(),ins(x,y,read());
dfs(1);double l=0.,r=1e6;
for(i=1;i<=40;++i)
{
if(l+eps>=r) break;
mid=(l+r)/2.;
if(check()) l=mid;else r=mid;
}
printf("%.3lf",l);
return 0;
}

Blog来自PaperCloud,未经允许,请勿转载,TKS!

[WC2010]重建计划(长链剖分版)的更多相关文章

  1. [WC2010]重建计划 长链剖分

    [WC2010]重建计划 LG传送门 又一道长链剖分好题. 这题写点分治的人应该比较多吧,但是我太菜了,只会长链剖分. 如果你还不会长链剖分的基本操作,可以看看我的长链剖分总结. 首先一看求平均值最大 ...

  2. 「WC2010」重建计划(长链剖分/点分治)

    「WC2010」重建计划(长链剖分/点分治) 题目描述 有一棵大小为 \(n\) 的树,给定 \(L, R\) ,要求找到一条长度在 \([L, R]\) 的路径,并且路径上边权的平均值最大 \(1 ...

  3. BZOJ1758[Wc2010]重建计划——分数规划+长链剖分+线段树+二分答案+树形DP

    题目描述 输入 第一行包含一个正整数N,表示X国的城市个数. 第二行包含两个正整数L和U,表示政策要求的第一期重建方案中修建道路数的上下限 接下来的N-1行描述重建小组的原有方案,每行三个正整数Ai, ...

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

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

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

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

  6. WC2010 BZOJ1758 重建计划_长链剖分

    题目大意: 求长度$\in [L,U]$的路径的最大边权和平均值. 题解 首先二分就不用说了,分数规划大家都懂. 这题有非常显然的点分治做法,但还是借着这个题学一波长链剖分. 其长链剖分本身也没啥,就 ...

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

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

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

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

  9. 洛谷 P4292 - [WC2010]重建计划(长链剖分+线段树)

    题面传送门 我!竟!然!独!立!A!C!了!这!道!题!incredible! 首先看到这类最大化某个分式的题目,可以套路地想到分数规划,考虑二分答案 \(mid\) 并检验是否存在合法的 \(S\) ...

随机推荐

  1. python-django中使用事务以及小坑

    django中使用事务 一.导入事务模块 from django.db import transaction 二.对相应的业务进行事务操作 方式一:为整个函数进行事务操作 @transaction.a ...

  2. 【洛谷 P2597】 [ZJOI2012]灾难(LCA)

    题目链接 考虑建一棵树,使一个生物灭绝时他的子树都会灭绝,显然这样答案就是以每个点为根的子树大小-1. 为什么原图不是一棵树,因为一个生物可能会以多个生物为食,所以按拓扑序来建树,把每个遍历到的点的父 ...

  3. C++中头文件与源文件的作用详解

    一.C++ 编译模式 通常,在一个 C++ 程序中,只包含两类文件―― .cpp 文件和 .h 文件.其中,.cpp 文件被称作 C++ 源文件,里面放的都是 C++ 的源代码:而 .h 文件则被称作 ...

  4. 我与SAP成都研究院吴院长的二三事

    这几天Jerry没怎么看手机,今天才注意到,昨天SAP中国研究院公众号上发布了一篇文章:SAP高管说: 体验经济时代下的SAP客户体验.仔细一看,这不是咱SAP成都研究院的吴院长么. 在今年没有发生部 ...

  5. Python学习资料收集

    1.Learn Python the hard way 2.Google's Python lesson 3.Python最佳实践 4.Full Stack Python 5.explore flas ...

  6. netcore里使用jwt做登陆授权

    1 什么是JWT? JWT是一种用于双方之间传递安全信息的简洁的.URL安全的表述性声明规范.JWT作为一个开放的标准(RFC 7519),定义了一种简洁的,自包含的方法用于通信双方之间以Json对象 ...

  7. 漫谈五种IO模型(主讲IO多路复用)

    首先引用levin的回答让我们理清楚五种IO模型 1.阻塞I/O模型 老李去火车站买票,排队三天买到一张退票. 耗费:在车站吃喝拉撒睡 3天,其他事一件没干. 2.非阻塞I/O模型 老李去火车站买票, ...

  8. java加密算法-SHA1

    public class SHAUtil { /*** * SHA加密 生成40位SHA码 * @param 待加密字符串 * @return 返回40位SHA码 */ public static S ...

  9. 史诗级干货-python爬虫之增加CSDN访问量

    史诗级干货-python爬虫之增加CSDN访问量 搜索微信公众号:'AI-ming3526'或者'计算机视觉这件小事' 获取更多算法.机器学习干货 csdn:https://blog.csdn.net ...

  10. 零基础如何学好Python 之int 数字整型类型 定义int()范围大小转换

    本文主题是讲python数字类型python int整型使用方法及技巧.它是不可变数据类型中的一种,它的一些性质和字符串是一样的,注意是整型不是整形哦. Python int有多种数字类型:整型int ...