hdu_5029_relief grain(树链剖分)
2016-06-26 21:36
483 查看
题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=5029
题意:给你一个树,然后给你两点,将这两点之间的点涂上颜色,问涂色最多的那个颜色是什么,如果数量相同,就输出编号小的颜色
题解:树链剖分,这个题很巧妙,线段树维护的是一个颜色线段树,然后如果要涂[l,r]的颜色a,就在l处标记a +1,在r+1处标记a -1。然后从左往右依次更新线段树,并query,就能得出答案,我的跑了1400+ms,鸟神的非递归线段树只跑我一半的时间
详解请看:鸟神的这题题解:relief grain
#include<cstdio> #include<cstring> #define F(i,a,b) for(int i=a;i<=b;++i) #define root 0,100001,1 #define ls l,m,rt<<1 #define rs m+1,r,rt<<1|1 #pragma comment(linker, "/STACK:102400000,102400000") const int N=100010; int n,m,x,y,z,idx,ans ,sum[N*4],nxt[10*N],g ,v[N*10],w[N*10],ed,pre ,siz ,dep ,hson ,top ,tid ,ftid ; inline void adg(int x,int y,int z){v[++ed]=y,w[ed]=z,nxt[ed]=g[x],g[x]=ed;} //树链部分 void dfs1(int u,int fa){ pre[u]=fa,siz[u]=1,dep[u]=dep[fa]+1,hson[u]=0; for(int i=g[u];~i;i=nxt[i]){ int vv=v[i]; if(vv!=fa){ dfs1(vv,u); if(siz[vv]>siz[hson[u]])hson[u]=vv; siz[u]+=siz[vv]; } } } void dfs2(int u,int tp){ tid[u]=idx,ftid[idx++]=u,top[u]=tp; if(hson[u])dfs2(hson[u],tp); for(int i=g[u];~i;i=nxt[i]){ int vv=v[i]; if(vv!=pre[u]&&vv!=hson[u])dfs2(vv,vv); } } void Update(int u,int v,int c){ int fu=top[u],fv=top[v]; while(fu!=fv){ if(dep[fu]>=dep[fv]) adg(tid[fu],c,1),adg(tid[u]+1,c,-1),u=pre[fu],fu=top[u]; else adg(tid[fv],c,1),adg(tid[v]+1,c,-1),v=pre[fv],fv=top[v]; } if(dep[u]>dep[v])u=u^v,v=u^v,u=u^v; adg(tid[u],c,1),adg(tid[v]+1,c,-1); } //线段树部分 void update(int x,int v,int l,int r,int rt){ if(l==r){sum[rt]+=v;return;} int m=(l+r)>>1; if(x<=m)update(x,v,ls);else update(x,v,rs); sum[rt]=sum[rt<<1]>sum[rt<<1|1]?sum[rt<<1]:sum[rt<<1|1]; } int query(int l,int r,int rt){ if(l==r)return l; int m=(l+r)>>1; if(sum[rt]==sum[rt<<1])return query(ls); else return query(rs); } int main(){ while(~scanf("%d%d",&n,&m),n+m){ memset(g,-1,sizeof(g)),ed=0; F(i,1,n-1)scanf("%d%d",&x,&y),adg(x,y,0),adg(y,x,0); dfs1(1,0),idx=1,dfs2(1,1); F(i,0,4*N-1)sum[i]=0; memset(g,-1,sizeof(g)),ed=0; F(i,1,m)scanf("%d%d%d",&x,&y,&z),Update(x,y,z); F(u,1,n){ for(int i=g[u];~i;i=nxt[i])update(v[i],w[i],root); ans[ftid[u]]=query(root); } F(i,1,n)printf("%d\n",ans[i]); } return 0; }View Code
相关文章推荐
- 利用Mail Exchanger (MX) Record原理 发送email
- 好用的ProgressBar 和播放动画的Dialog使用
- JTAG无法识别单片机(JTAG Device chain)
- http://blog.csdn.net/zhshulin/article/details/37956105
- 220. Contains Duplicate III
- http://blog.csdn.net/zhshulin/article/details/30779873
- 从硬盘设计思想到RAID改良之道
- 从硬盘设计思想到RAID改良之道
- Sub Main启动与窗体启动
- [217] Contains Duplicate
- Leetcode climbing-stair
- [置顶] ArcEngine Mosaic Dataset 镶嵌数据集总结 (持续补充中......)
- svn svn_uri_basename: Assertion `svn_uri_is_canonical(uri, ((void *)0))' failed.
- constraint the design
- Chain Of Responsibility(职责连)-对象行为型模式
- HDU 1848 Fibonacci again and again(学习sg函数后写出来的第一个题)
- 写一个brainfuck语言的解释器
- hdu 1532 Drainage Ditches(最大流,Edmond Karp)
- lintcode_2 Trailing Zeros
- 人工智能作业homework5------naive bayes/svm文本分类