P3714 [BJOI2017]树的难题 点分治+线段树合并
题目描述
分析
路径问题考虑点分治
对于一个分治中心,我们可以很容易地得到从它开始的一条路径的价值和长度
问题就是如何将不同的路径合并
很显然,对于同一个子树中的所有路径,它们起始的颜色是相同的
因此我们可以将一个节点的所有子结点按照颜色排序
这个可以在建图之前处理好
然后开两个权值线段树,一棵存储与当前节点起始颜色相同的所有路径,另一棵存储不同的所有路径
当节点的颜色改变时,把相同的那一堆合并到另一堆即可
时间复杂度 \(nlog^2n\)
代码
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#define rg register
inline int read(){
rg int x=0,fh=1;
rg char ch=getchar();
while(ch<'0' || ch>'9'){
if(ch=='-') fh=-1;
ch=getchar();
}
while(ch>='0' && ch<='9'){
x=(x<<1)+(x<<3)+(ch^48);
ch=getchar();
}
return x*fh;
}
const int maxn=4e5+5;
int h[maxn],tot=1,n,m,lef,rig,val[maxn];
struct asd{
int to,nxt,val;
}b[maxn];
void ad(rg int aa,rg int bb,rg int cc){
b[tot].to=bb;
b[tot].nxt=h[aa];
b[tot].val=cc;
h[aa]=tot++;
}
int siz[maxn],maxsiz[maxn],rt,totsiz;
bool vis[maxn];
void getroot(rg int now,rg int lat){
siz[now]=1,maxsiz[now]=0;
for(rg int i=h[now];i!=-1;i=b[i].nxt){
rg int u=b[i].to;
if(u==lat || vis[u]) continue;
getroot(u,now);
siz[now]+=siz[u];
maxsiz[now]=std::max(maxsiz[now],siz[u]);
}
maxsiz[now]=std::max(maxsiz[now],totsiz-siz[now]);
if(maxsiz[now]<maxsiz[rt]) rt=now;
}
struct trr{
int lch,rch,val;
trr(){
lch=rch=0;
val=-0x3f3f3f3f;//一定要初始化成无穷小
}
}tr[maxn*40];
int cnt,rt1,rt2,tp,ans=-0x3f3f3f3f;
void push_up(rg int da){
tr[da].val=std::max(tr[tr[da].lch].val,tr[tr[da].rch].val);
}
int xg(rg int da,rg int l,rg int r,rg int wz,rg int val){
if(!da){
da=++cnt;
tr[da].lch=tr[da].rch=0,tr[da].val=-0x3f3f3f3f;
}
if(l==r){
tr[da].val=std::max(tr[da].val,val);
return da;
}
rg int mids=(l+r)>>1;
if(wz<=mids) tr[da].lch=xg(tr[da].lch,l,mids,wz,val);
else tr[da].rch=xg(tr[da].rch,mids+1,r,wz,val);
push_up(da);
return da;
}
int bing(rg int aa,rg int bb,rg int l,rg int r){
if(!aa || !bb) return aa+bb;
if(l==r){
tr[aa].val=std::max(tr[aa].val,tr[bb].val);
tr[bb].val=-0x3f3f3f3f;
return aa;
}
rg int mids=(l+r)>>1;
tr[aa].lch=bing(tr[aa].lch,tr[bb].lch,l,mids);
tr[aa].rch=bing(tr[aa].rch,tr[bb].rch,mids+1,r);
push_up(aa);
return aa;
}
int cx(rg int da,rg int l,rg int r,rg int nl,rg int nr){
if(!da || l>r) return -0x3f3f3f3f;
if(l>=nl && r<=nr) return tr[da].val;
rg int mids=(l+r)>>1,nans=-0x3f3f3f3f;
if(nl<=mids) nans=std::max(nans,cx(tr[da].lch,l,mids,nl,nr));
if(nr>mids) nans=std::max(nans,cx(tr[da].rch,mids+1,r,nl,nr));
return nans;
}
struct jie{
int val,dep;
jie(){}
jie(rg int aa,rg int bb){
val=aa,dep=bb;
}
}sta[maxn];
bool cmp(rg jie aa,rg jie bb){
return aa.val>bb.val;
}
std::vector<jie> g[maxn];
void dfs(rg int now,rg int lat,rg int nval,rg int ndep,rg int latcol){
if(ndep>rig) return;
sta[++tp]=jie(nval,ndep);
for(rg int i=h[now];i!=-1;i=b[i].nxt){
rg int u=b[i].to;
if(u==lat || vis[u]) continue;
dfs(u,now,(latcol==b[i].val)?nval:nval+val[b[i].val],ndep+1,b[i].val);
}
}
void solve(rg int now){
vis[now]=1;
rt1=rt2=cnt=0;
rg int latcol=0,jud=0;
for(rg int i=h[now];i!=-1;i=b[i].nxt){
rg int u=b[i].to;
if(!vis[u]){
tp=jud=0;
dfs(u,now,val[b[i].val],1,b[i].val);
if(b[i].val==latcol) jud=1;
else {
rt1=bing(rt1,rt2,1,n);
rt2=0;
}
for(rg int j=1;j<=tp;j++){
if(jud) ans=std::max(ans,cx(rt2,1,n,std::max(1,lef-sta[j].dep),rig-sta[j].dep)+sta[j].val-val[latcol]);
ans=std::max(ans,cx(rt1,1,n,std::max(1,lef-sta[j].dep),rig-sta[j].dep)+sta[j].val);
if(sta[j].dep>=lef && sta[j].dep<=rig) ans=std::max(ans,sta[j].val);
}
for(rg int j=1;j<=tp;j++) rt2=xg(rt2,1,n,sta[j].dep,sta[j].val);
latcol=b[i].val;
}
}
for(rg int i=h[now];i!=-1;i=b[i].nxt){
rg int u=b[i].to;
if(!vis[u]){
totsiz=siz[u],rt=0;
getroot(u,now);
solve(rt);
}
}
}
int main(){
memset(h,-1,sizeof(h));
n=read(),m=read(),lef=read(),rig=read();
for(rg int i=1;i<=m;i++) val[i]=read();
rg int aa,bb,cc;
for(rg int i=1;i<n;i++){
aa=read(),bb=read(),cc=read();
g[aa].push_back(jie(cc,bb)),g[bb].push_back(jie(cc,aa));
}
for(rg int i=1;i<=n;i++) std::sort(g[i].begin(),g[i].end(),cmp);
for(rg int i=1;i<=n;i++){
for(rg int j=0;j<g[i].size();j++){
ad(i,g[i][j].dep,g[i][j].val);
}
}
maxsiz[0]=0x3f3f3f3f,rt=0,totsiz=n;
getroot(1,0);
solve(rt);
printf("%d\n",ans);
return 0;
}
P3714 [BJOI2017]树的难题 点分治+线段树合并的更多相关文章
- [BJOI2017]树的难题 点分治 线段树
题面 [BJOI2017]树的难题 题解 考虑点分治. 对于每个点,将所有边按照颜色排序. 那么只需要考虑如何合并2条链. 有2种情况. 合并路径的接口处2条路径颜色不同 合并路径的接口处2条路径颜色 ...
- [BJOI2017]树的难题 点分治,线段树合并
[BJOI2017]树的难题 LG传送门 点分治+线段树合并. 我不会写单调队列,所以就写了好写的线段树. 考虑对于每一个分治中心,把出边按颜色排序,这样就能把颜色相同的子树放在一起处理.用一棵动态开 ...
- UVALive 7148 LRIP【树分治+线段树】
题意就是要求一棵树上的最长不下降序列,同时不下降序列的最小值与最大值不超过D. 做法是树分治+线段树,假设树根是x,y是其当前需要处理的子树,对于子树y,需要处理出两个数组MN,MX,MN[i]表示以 ...
- BZOJ4012[HNOI2015]开店——树链剖分+可持久化线段树/动态点分治+vector
题目描述 风见幽香有一个好朋友叫八云紫,她们经常一起看星星看月亮从诗词歌赋谈到 人生哲学.最近她们灵机一动,打算在幻想乡开一家小店来做生意赚点钱.这样的 想法当然非常好啦,但是她们也发现她们面临着一个 ...
- LOJ#6463 AK YOI 树分治+线段树合并
传送门 既然是树上路径统计问题,不难想到要使用树分治,这里以点分治为例 由点分治的性质,每层只需要考虑经过重心的路径 因为需要维护路径长度在一定范围内的最大权值和,所以要用一个数据结构维护一下到根节点 ...
- 【loj6145】「2017 山东三轮集训 Day7」Easy 动态点分治+线段树
题目描述 给你一棵 $n$ 个点的树,边有边权.$m$ 次询问,每次给出 $l$ .$r$ .$x$ ,求 $\text{Min}_{i=l}^r\text{dis}(i,x)$ . $n,m\le ...
- 【BZOJ4372】烁烁的游戏 动态树分治+线段树
[BZOJ4372]烁烁的游戏 Description 背景:烁烁很喜欢爬树,这吓坏了树上的皮皮鼠.题意:给定一颗n个节点的树,边权均为1,初始树上没有皮皮鼠.烁烁他每次会跳到一个节点u,把周围与他距 ...
- 【bzoj4372】烁烁的游戏 动态点分治+线段树
题目描述 给一颗n个节点的树,边权均为1,初始点权均为0,m次操作:Q x:询问x的点权.M x d w:将树上与节点x距离不超过d的节点的点权均加上w. 输入 第一行两个正整数:n,m接下来的n-1 ...
- 【bzoj3730】震波 动态点分治+线段树
题目描述 在一片土地上有N个城市,通过N-1条无向边互相连接,形成一棵树的结构,相邻两个城市的距离为1,其中第i个城市的价值为value[i].不幸的是,这片土地常常发生地震,并且随着时代的发展,城市 ...
随机推荐
- 题解-FJOI2018 领导集团问题
题面 FJOI2018 领导集团问题 给一棵树 \(T(|T|=n)\),每个点有个权值 \(w_i\),从中选出一个子点集 \(P=\{x\in {\rm node}|x\in T\}\),使得 \ ...
- 计算机网络——HTTP知识点整理
1.HTTP简介 HTTP协议(HyperText Transfer Protocol,超文本传输协议)是用于从WWW服务器传输超文本到本地浏览器的传送协议.它可以使浏览器更加高效,使网络传输减少.它 ...
- centos7网卡bond配置--自己另一篇文章的补充
这篇文章是自己另一篇文章的第二种方法的一个完善的补充 https://www.cnblogs.com/zzf0305/p/9588585.html 1 备份网卡配置文件2 使用nmcli命令配置bon ...
- C++ 虚函数表与多态 —— 关键字 override 的用法
override 仅能用于虚函数,他属于C++新特性,是重写覆盖的意思,他的存在仅仅是为了提高代码的可阅读性: 作用: 1. 提示程序的阅读者,这个函数是重写父类的功能. 2. 防止程序员在重写父类的 ...
- 高手查看Linux系统用户命令-测评
一.Linux查看用户命令-测评 查看linux所有的用户 cat /etc/passwd 查看普通用户.系统用户(1-499) root:x:0:0:root:/root:/bin/bash < ...
- UML—20—003
博客班级 https://edu.cnblogs.com/campus/fzzcxy/2018SE1 作业要求 https://edu.cnblogs.com/campus/fzzcxy/2018SE ...
- 为什么类只能用public修饰?
为什么类只能使用public修饰? 首先,类只能使用public修饰是一个伪命题,应该说我们只见到过使用public修饰的类,还有一些类没有访问修饰符,此时访问权限为default.其次,类实际上分为 ...
- 【漏洞复现】Struts2-045分析(CVE-2017-5638)
如果需要大佬写好的脚本,可以直接去github上面搜 struts2 - 045 一个还比较出名的漏洞,因为涉及到利用Gopher协议反弹shell,所以写篇文章来简单学习下这个漏洞. Struts2 ...
- 在matlab 画箭头
[转载]在matlab 画箭头 原文地址:在matlab 画箭头作者:纯情小郎君 完整见链接http://www.mathworks.com/matlabcentral/fx_files/14056/ ...
- Python读取word文档内容
1,利用python读取纯文字的word文档,读取段落和段落里的文字. 先读取段落,代码如下: 1 ''' 2 #利用python读取word文档,先读取段落 3 ''' 4 #导入所需库 5 fro ...