题目描述

www.lydsy.com/JudgeOnline/upload/201805/day1(1).pdf

题解

首先来看这个我们要最大化的东西。

deep[u]+deep[v]-deep[lca(u,v)]-deep[lca(u',v')]

后面的那个东西看起来不太合群,我们可以把前后拆开。

deep[u]+deep[v]-deep[lca(u,v)]

我们发现这其实就是u到根的链和v到根的链的并。

然后它还等于(deep[u]+deep[v]+dis[u][v])/2

因为deep数组我们可以直接求出,所以我们就把一颗有根树上的问题放到了无根树上,也就是可以去掉lca的影响了。

然后考虑枚举第二颗树的LCA,那么一组合法的点应当在这个点的两颗不同的子树中。

然后对第一棵树边分,发现这颗边分树也是一颗二叉树,每个叶子结点代表原树上的一个点。

于是这题的做法来了,我们在dfs第二颗树的时候,像线段树合并一样合并边分树,因为不管叶子的情况下,每个节点都代表一条边,每条边连接着两个点。

这样我们对这个点记一个lans和rans分别代表左端点的最优答案和右端点的最优答案。

维护答案的形式为deep[x]+dis(x,edge)

合并的时候顺带计算答案。

注意,要考虑u和v重合的情况。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define inf 1e18
#define N 740009
using namespace std;
typedef long long ll;
ll lans[N*],rans[N*],ans,dis[N],deep[N][],nowdeep,val[N*];
int atp,size[N],sum,nowroot,root,ha,finaldep[N];
int fa[N*],ls[N*],rs[N*],tr[N*][],dian,T[N],n,id[N*];
bool jin[N<<];
inline ll rd(){
ll x=;char c=getchar();bool f=;
while(!isdigit(c)){if(c=='-')f=;c=getchar();}
while(isdigit(c)){x=(x<<)+(x<<)+(c^);c=getchar();}
return f?-x:x;
}
struct tu{
int head[N],tot;
struct edge{int n,to;ll l;}e[N<<];
void clear(){memset(head,,sizeof(head));tot=;}
inline void add(int u,int v,ll l){e[++tot].n=head[u];e[tot].to=v;head[u]=tot;e[tot].l=l;}
}E[];
struct node{int n,to;ll l;}e[N<<];
int head[N<<],tot=;
inline void add(int u,int v,ll l){
e[++tot].n=head[u];e[tot].to=v;head[u]=tot;e[tot].l=l;
e[++tot].n=head[v];e[tot].to=u;head[v]=tot;e[tot].l=l;
}
void dfs1(int u,int fa){
int now=;
for(int i=E[].head[u];i;i=E[].e[i].n)if(E[].e[i].to!=fa){
int v=E[].e[i].to;dis[v]=dis[u]+E[].e[i].l;
if(!now){add(u,v,E[].e[i].l);now=u;}
else{++atp;add(now,atp,);add(atp,v,E[].e[i].l);now=atp;}
dfs1(v,u);
}
}
void getroot(int u,int fa){
size[u]=;
for(int i=head[u];i;i=e[i].n)if(e[i].to!=fa&&!jin[i]){
int v=e[i].to;
getroot(v,u);
size[u]+=size[v];
if(max(size[v],sum-size[v])<nowroot){root=i;nowroot=max(size[v],sum-size[v]);ha=size[v];}
}
}
void getdeep(int u,int fa,int dep){
for(int i=head[u];i;i=e[i].n)if(e[i].to!=fa&&!jin[i]){
int v=e[i].to;
deep[v][dep]=deep[u][dep]+e[i].l;
getdeep(v,u,dep);
}
}
int solve(int u,int s,int dep){
if(s==){finaldep[u]=dep;return u;}
int now=++atp;
root=atp+;nowroot=atp;sum=s;
getroot(u,);
jin[root]=jin[root^]=;
int x=e[root].to,y=e[root^].to,xs=ha,ys=s-ha;val[atp]=e[root].l;
getdeep(x,y,dep);getdeep(y,x,dep);
fa[ls[now]=solve(x,xs,dep+)]=now;fa[rs[now]=solve(y,ys,dep+)]=now;
return now;
}
inline int ins(int x){
int now=x,pre=;
for(int i=finaldep[x];i;--i){
++dian;id[dian]=fa[now];lans[dian]=rans[dian]=-inf;
if(ls[fa[now]]==now)lans[dian]=dis[x]+deep[x][i-],tr[dian][]=pre;
if(rs[fa[now]]==now)rans[dian]=dis[x]+deep[x][i-],tr[dian][]=pre;
pre=dian;now=fa[now];
}
return dian;
}
inline int merge(int x,int y){
if(!x||!y)return x^y;
ans=max(ans,(lans[x]+rans[y]+val[id[x]])/-nowdeep);
ans=max(ans,(rans[x]+lans[y]+val[id[x]])/-nowdeep);
lans[x]=max(lans[x],lans[y]);rans[x]=max(rans[x],rans[y]);
tr[x][]=merge(tr[x][],tr[y][]);tr[x][]=merge(tr[x][],tr[y][]);
return x;
}
void dfs2(int u,int fa,ll d){
T[u]=ins(u);
ans=max(ans,dis[u]-d);
for(int i=E[].head[u];i;i=E[].e[i].n)if(E[].e[i].to!=fa){
int v=E[].e[i].to;
dfs2(v,u,d+E[].e[i].l);nowdeep=d;
T[u]=merge(T[u],T[v]);
}
}
int main(){
n=rd();int u,v;ll w;ans=-inf;
for(int i=;i<n;++i){u=rd();v=rd();w=rd();E[].add(u,v,w);E[].add(v,u,w);}
for(int i=;i<n;++i){u=rd();v=rd();w=rd();E[].add(u,v,w);E[].add(v,u,w);}
atp=n;dfs1(,);
solve(,atp,);
dfs2(,,);
printf("%lld\n",ans);
return ;
}

[CTSC2018]暴力写挂的更多相关文章

  1. [CTSC2018]暴力写挂——边分树合并

    [CTSC2018]暴力写挂 题面不错 给定两棵树,两点“距离”定义为:二者深度相加,减去两棵树上的LCA的深度(深度指到根节点的距离) 求最大的距离. 解决多棵树的问题就是降维了. 经典的做法是边分 ...

  2. [LOJ#2553][CTSC2018]暴力写挂

    [LOJ#2553][CTSC2018]暴力写挂 试题描述 temporaryDO 是一个很菜的 OIer .在 4 月,他在省队选拔赛的考场上见到了<林克卡特树>一题,其中 \(k = ...

  3. BZOJ5341: [Ctsc2018]暴力写挂

    BZOJ5341: [Ctsc2018]暴力写挂 https://lydsy.com/JudgeOnline/problem.php?id=5341 分析: 学习边分治. 感觉边分治在多数情况下都能用 ...

  4. BZOJ5341[Ctsc2018]暴力写挂——边分治+虚树+树形DP

    题目链接: CSTC2018暴力写挂 题目大意:给出n个点结构不同的两棵树,边有边权(有负权边及0边),要求找到一个点对(a,b)满足dep(a)+dep(b)-dep(lca)-dep'(lca)最 ...

  5. UOJ400/LOJ2553 CTSC2018 暴力写挂 边分治、虚树

    传送门--UOJ 传送门--LOJ 跟隔壁通道是一个类型的 要求的式子中有两个LCA,不是很方便,因为事实上在这种题目中LCA一般都是枚举的对象-- 第二棵树上的LCA显然是动不了的,因为没有其他的量 ...

  6. bzoj 5341: [Ctsc2018]暴力写挂

    Description Solution 边分治+边分树合并 这个题很多做法都是启发式合并的复杂度的,都有点卡 以前有个套路叫做线段树合并优化启发式合并,消掉一个 \(log\) 这个题思路类似,建出 ...

  7. 并不对劲的bzoj5341:loj2553:uoj400:p4565:[Ctsc2018]暴力写挂

    题目大意 有两棵\(n\)(\(n\leq366666\))个节点的树,\(T\)和\(T'\),有边权 \(dep(i)\)表示在\(T\)中\(i\)号点到\(1\)号点的距离,\(dep'(i) ...

  8. 题解 「CTSC2018暴力写挂」

    题目传送门 题目大意 给出两个大小为 \(n\) 的树,求出: \[\max\{\text{depth}(x)+\text{depth}(y)-\text{depth}(\text{LCA}(x,y) ...

  9. 【CTSC2018】暴力写挂(边分治,虚树)

    [CTSC2018]暴力写挂(边分治,虚树) 题面 UOJ BZOJ 洛谷 题解 发现第二棵树上的\(LCA\)的深度这玩意没法搞,那么枚举在第二棵树上的\(LCA\). 然后剩下的部分就是\(dep ...

随机推荐

  1. Docker Data Center系列(一)- 快速搭建云原生架构的实践环境

    本系列文章演示如何快速搭建一个简单的云原生架构的实践环境. 基于这个基础架构,可以持续部署微服务架构的应用栈,演练敏捷开发过程,提升DevOps实践能力. 1 整体规划 1.1 拓扑架构 1.2 基础 ...

  2. 我现在有个表,里面有100个不同的单词,每个单词对应有大概20个词组,我想通过sql,每个单词随机获取对应的3个词组,请问怎么写可以实现?

    闲来无事刷技术论坛,看到一个这样的问题: 我现在有个表,里面有100个不同的单词,每个单词对应有大概20个词组,我想通过sql,每个单词随机获取对应的3个词组,请问怎么写可以实现? 感觉题材很新颖,角 ...

  3. MySQL常用数值函数

    数值函数: 用来处理很多数值方面的运算,使用数值函数,可以免去很多繁杂的判断求值的过程,能够大大提高用户的工作效率. 1.ABS(x):返回 x 的绝对值 mysql> select abs(- ...

  4. SQLServer之创建INSTEAD OF INSERT,UPDATE,DELETE触发器

    INSTEAD OF触发器工作原理 INSTEAD OF表示并不执行其所定义的操作INSERT,UPDATE ,DELETE,而仅是执行触发器本身,即当对表进行INSERT.UPDATE 或 DELE ...

  5. About Pull Strings 英语走后门议论文

    About pull strings Author : Pleiades_Antares 1. From ancient times to the present, the "going b ...

  6. Java开发学习心得(一):SSM环境搭建

    目录 Java开发学习心得(一):SSM环境搭建 1 SSM框架 1.1 Spring Framework 1.2 Spring MVC Java开发学习心得(一):SSM环境搭建 有一点.NET的开 ...

  7. 英语口语练习系列-C06-购物

    <水调歌头>·苏轼 明月几时有,把酒问青天. 不知天上宫阙,今夕是何年? 我欲乘风归去,又恐琼楼玉宇, 高处不胜寒. 起舞弄清影,何似在人间! 转朱阁,低绮户,照无眠. 不应有恨,何事长向 ...

  8. Kafka配置项unclean.leader.election.enable造成consumer出现offset重置现象

    消费端出现offset重置为latest, earliest现象,类似log: (org.apache.kafka.clients.consumer.internals.Fetcher.handleF ...

  9. JavaScript的函数闭包详细解释

    闭包是指有权访问另一个函数作用域中的变量的函数 一.创建闭包的常见的方式: 就是在一个函数内部创建另一个函数,通过另一个函数访问这个函数的局部变量. //通过闭包可以返回局部变量 function b ...

  10. 5000量子位支持量子编程,D-Wave推出下一代量子计算平台计划

    5000量子位支持量子编程,D-Wave推出下一代量子计算平台计划 近日,全球量子商用化重要参与者 D-Wave 公司又有大动作:推出其5000量子比特量子计算的发展蓝图.D-Wave 下一代量子计算 ...