洛咕 P4556 [Vani有约会]雨天的尾巴
终于把考试题清完了。。。又复活了。。。
树上差分,合并用线段树合并,但是空间会炸。
某大佬:lca和fa[lca]减得时候一定已经存在这个节点了,所以放进vector里,合并完之后减掉就好了。。。
玄学优化就过了。。
// luogu-judger-enable-o2
#include<bits/stdc++.h>
#define il inline
#define vd void
#define N 100000
typedef long long ll;
il int gi(){
int x=0,f=1;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-')f=-1;
ch=getchar();
}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return x*f;
}
int fir[100010],dis[200010],nxt[200010],id,siz[100010],son[100010],fa[100010],top[100010],dep[100010];
il vd link(int a,int b){nxt[++id]=fir[a],fir[a]=id,dis[id]=b;}
typedef struct node* point;
struct node{
int mx,qwq;
point ls,rs;
node(){mx=0,qwq=0,ls=NULL,rs=NULL;}
};
il vd dfs(int x){
siz[x]=1;
for(int i=fir[x];i;i=nxt[i]){
if(fa[x]==dis[i])continue;
dep[dis[i]]=dep[x]+1;fa[dis[i]]=x;
dfs(dis[i]);siz[x]+=siz[dis[i]];
if(siz[dis[i]]>siz[son[x]])son[x]=dis[i];
}
}
il vd dfs2(int x,int tp){
top[x]=tp;
if(son[x])dfs2(son[x],tp);
for(int i=fir[x];i;i=nxt[i])if(dis[i]!=fa[x]&&dis[i]!=son[x])dfs2(dis[i],dis[i]);
}
il int lca(int a,int b){
while(top[a]^top[b])
if(dep[top[a]]>dep[top[b]])a=fa[top[a]];
else b=fa[top[b]];
return dep[a]<dep[b]?a:b;
}
#define mid ((l+r)>>1)
il vd upd(const point&x){
if(x->ls==NULL)x->mx=x->rs->mx,x->qwq=x->rs->qwq;
else if(x->rs==NULL)x->mx=x->ls->mx,x->qwq=x->ls->qwq;
else if(x->ls->mx>=x->rs->mx)x->mx=x->ls->mx,x->qwq=x->ls->qwq;
else x->mx=x->rs->mx,x->qwq=x->rs->qwq;
}
il vd update(point&x,int l,int r,const int&p,const int&s){
if(x==NULL)x=new node;
if(l==r){x->mx+=s,x->qwq=l;return;}
if(p<=mid)update(x->ls,l,mid,p,s);
else update(x->rs,mid+1,r,p,s);
upd(x);
}
point rt[100010];
int ans[100010];
il vd merge(point&x,point y,int l=1,int r=N){
if(x==NULL){x=y;return;}
else if(y==NULL)return;
if(x->ls==NULL&&x->rs==NULL){x->mx+=y->mx;return;}
merge(x->ls,y->ls,l,mid),merge(x->rs,y->rs,mid+1,r);
upd(x);
}
std::vector<int>S[100010];
il vd Merge(int x){
for(int i=fir[x];i;i=nxt[i]){
if(fa[x]==dis[i])continue;
Merge(dis[i]);
merge(rt[x],rt[dis[i]]);
}
for(int i=0;i<S[x].size();++i)update(rt[x],1,N,S[x][i],-1);
if(rt[x]!=NULL&&rt[x]->mx)ans[x]=rt[x]->qwq;
}
int main(){
#ifndef ONLINE_JUDGE
freopen("4556.in","r",stdin);
freopen("4556.out","w",stdout);
#endif
int n=gi(),m=gi(),u,v,w,l;
for(int i=1;i<n;++i)u=gi(),v=gi(),link(u,v),link(v,u);
dfs(1);dfs2(1,1);
while(m--){
u=gi(),v=gi(),w=gi();l=lca(u,v);
update(rt[u],1,N,w,1);update(rt[v],1,N,w,1);
S[l].push_back(w);
S[fa[l]].push_back(w);
}
Merge(1);
for(int i=1;i<=n;++i)printf("%d\n",ans[i]);
return 0;
}
洛咕 P4556 [Vani有约会]雨天的尾巴的更多相关文章
- 洛谷 P4556 [Vani有约会]雨天的尾巴 解题报告
P4556 [Vani有约会]雨天的尾巴 题目背景 深绘里一直很讨厌雨天. 灼热的天气穿透了前半个夏天,后来一场大雨和随之而来的洪水,浇灭了一切. 虽然深绘里家乡的小村落对洪水有着顽固的抵抗力,但也倒 ...
- 洛谷P4556 [Vani有约会]雨天的尾巴(线段树合并)
题目背景 深绘里一直很讨厌雨天. 灼热的天气穿透了前半个夏天,后来一场大雨和随之而来的洪水,浇灭了一切. 虽然深绘里家乡的小村落对洪水有着顽固的抵抗力,但也倒了几座老房子,几棵老树被连根拔起,以及田地 ...
- 2018.08.28 洛谷P4556 [Vani有约会]雨天的尾巴(树上差分+线段树合并)
传送门 要求维护每个点上出现次数最多的颜色. 对于每次修改,我们用树上差分的思想,然后线段树合并统计答案就行了. 注意颜色很大需要离散化. 代码: #include<bits/stdc++.h& ...
- P4556 [Vani有约会]雨天的尾巴(线段树合并+lca)
P4556 [Vani有约会]雨天的尾巴 每个操作拆成4个进行树上差分,动态开点线段树维护每个点的操作. 离线处理完向上合并就好了 luogu倍增lca被卡了5分.....于是用rmq维护.... 常 ...
- P4556 [Vani有约会]雨天的尾巴 (线段树合并)
P4556 [Vani有约会]雨天的尾巴 题意: 首先村落里的一共有n座房屋,并形成一个树状结构.然后救济粮分m次发放,每次选择两个房屋(x,y),然后对于x到y的路径上(含x和y)每座房子里发放一袋 ...
- [题解] P4556 [Vani有约会]雨天的尾巴
[题解] P4556 [Vani有约会]雨天的尾巴 ·题目大意 给定一棵树,有m次修改操作,每次修改 \(( x\) \(y\) \(z )\) 表示 \((x,y)\) 之间的路径上数值 \(z\) ...
- P4556 [Vani有约会]雨天的尾巴(线段树合并)
传送门 一道线段树合并 首先不难看出树上差分 我们把每一次修改拆成四个,在\(u,v\)分别放上一个,在\(lca\)和\(fa[lca]\)各减去一个,那么只要统计一下子树里的总数即可 然而问题就在 ...
- 洛谷4556 [Vani有约会]雨天的尾巴
原题链接 每个点开一个权值线段树,然后用树上差分的方法修改,最后自底向上暴力线段树合并即可. 不过空间较大,会\(MLE\),写个内存池就可以了. #include<cstdio> #in ...
- P4556 [Vani有约会]雨天的尾巴
目录 思路 优化 过程中的问题/疑问 错误 代码 思路 每个节点维护一课线段树(当然是动态开点) 线段树的作用是统计这个节点有多少种粮食型号,以及最多的粮食型号 然后树上差分,u和v点 +1,lca( ...
随机推荐
- 转:只能选择GridView中的一个CheckBox(单选CheckBox)
方法1: protected void GridView2_RowDataBound(object sender, GridViewRowEventArgs e){CheckBox cbx = e.R ...
- [翻译] SAMCoreImageView
SAMCoreImageView https://github.com/soffes/SAMCoreImageView Render a CIImage in an OpenGL thingy so ...
- spider-抓取页面内容
# -*- coding: UTF-8 -*- from HTMLParser import HTMLParser import sys,urllib2,string,re,json reload(s ...
- python 下字符串格式时间比较
python 下有多个有关时间的模块,分别是time.datetime.calendar,今天重点讨论下time写法. 其中time模块,主要有以下方法: ltime=time.time() 获取当前 ...
- 获取索引--------用range()和len()
a = ['Google', 'Baidu', 'Runoob', 'Taobao', 'QQ'] for i in range(len(a)): print(i+1,a[i])
- Istio 1.1尝鲜记
近几天Istio1.1的发布引起了技术界巨大的反响,为了让更多技术爱好者能够亲自体验Istio1.1,公司的技术大佬赶出了这篇尝鲜教程,其中包括环境.安装.可能遇到的问题及解决方式等,希望对大家有所帮 ...
- setnx redis
使用锁 1)setnx(lockkey, 当前时间+过期超时时间) ,如果返回1,则获取锁成功:如果返回0则没有获取到锁,转向2.2.)get(lockkey)获取值oldExpireTime ,并将 ...
- 转载】JQuery中如何传递参数如click(),change()等具体实现
转载地址:http://www.jb51.net/article/36249.htm 有个需求让两个select中option相互转换,这个作业就是给几个按钮添加click()事件接下来为大家介绍下如 ...
- Java虚拟机17:互斥同步、锁优化及synchronized和volatile
互斥同步 互斥同步(Mutual Exclusion & Synchronization)是常见的一种并发正确性保证手段.同步是指子啊多个线程并发访问共享数据时,保证共享数据在同一时刻只能被一 ...
- 3192: [JLOI2013]删除物品
3192: [JLOI2013]删除物品 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 1366 Solved: 794 [Submit][Statu ...