fzu2028
//Accepted 7324 KB 203 ms /* source:fzu2028 time :2015.5.29 by :songt */ /*题解: 树链剖分 单点更新,求路径和 */ #include <cstdio> #include <cstring> ; int max(int a,int b) { return a>b?a:b; } void swap(int &a,int &b) { int t=a; a=b; b=t; } struct Edge { int u,v; Edge(){} Edge(int u,int v):u(u),v(v){} } edge[*imax_n]; int head[imax_n]; *imax_n]; int tot; void addEdge(int u,int v) { edge[tot]=Edge(u,v); next[tot]=head[u]; head[u]=tot++; } int fa[imax_n],deep[imax_n],num[imax_n],son[imax_n]; int p[imax_n],fp[imax_n],top[imax_n]; int pos; void init() { memset(head,-,sizeof(head)); memset(next,-,sizeof(next)); tot=; memset(son,-,sizeof(son)); pos=; } void dfs1(int u,int pre,int depth) { deep[u]=depth; fa[u]=pre; num[u]=; ;i=next[i]) { int v=edge[i].v; if (v!=pre) { dfs1(v,u,depth+); num[u]+=num[v]; || num[son[u]]<num[v]) son[u]=v; } } } void dfs2(int u,int sp) { p[u]=pos++; fp[p[u]]=u; top[u]=sp; ) return ; dfs2(son[u],sp); ;i=next[i]) { int v=edge[i].v; if (v!=fa[u] && v!=son[u]) { dfs2(v,v); } } } struct Tree { int l,r; long long sum; int tmax; }f[imax_n*]; void build(int t,int l,int r) { f[t].l=l; f[t].r=r; f[t].sum=; f[t].tmax=; if (l==r) { return ; } ; build(*t,l,mid); build(*t+,mid+,r); } void update(int t,int k,int value) { if (f[t].l==k && f[t].r==k) { f[t].sum=value; f[t].tmax=value; return ; } ; *t,k,value); *t+,k,value); f[t].sum=f[*t].sum+f[*t+].sum; f[t].tmax=max(f[*t].tmax,f[*t+].tmax); } long long query(int t,int l,int r,int &tmax) { if (f[t].l==l && f[t].r==r) { tmax=f[t].tmax; return f[t].sum; } ; *t,l,r,tmax); else { *t+,l,r,tmax); else { int tmax1,tmax2; long long sum1,sum2; sum1=query(*t,l,mid,tmax1); sum2=query(*t+,mid+,r,tmax2); tmax=max(tmax1,tmax2); return sum1+sum2; } } } long long find(int u,int v,int tmax) { int f1=top[u],f2=top[v]; ; ; int t; while (f1!=f2) { if (deep[f1]<deep[f2]) { swap(f1,f2); swap(u,v); } sum+=query(,p[f1],p[u],t); tmp=max(tmp,t); u=fa[f1]; f1=top[u]; } if (u==v) return sum; if (deep[u]>deep[v]) swap(u,v); sum+=query(,p[son[u]],p[v],t); tmp=max(tmp,t); tmax=tmp; return sum; } ]; int n,m; int main() { ) { init(); ;i<n-;i++) { scanf(],&e[i][],&e[i][]); addEdge(e[i][],e[i][]); addEdge(e[i][],e[i][]); } dfs1(,,); dfs2(,); //printf("after dfs\n"); build(,,pos-); //printf("after build\n"); ;i<n-;i++) { ]]<deep[e[i][]]) swap(e[i][],e[i][]); update(,p[e[i][]],e[i][]); } //printf("after insert\n"); int kind; int u,v,c; int tmax; ;i<m;i++) { scanf("%d%d%d",&kind,&u,&v); ) { update(,p[e[u-][]],v); } else { printf("%lld\n",find(u,v,tmax)); } } } ; }
fzu2028的更多相关文章
随机推荐
- 打开hibernate文件报警告
在web工程中新建一个,Hibetnate项目,打开配置文件报警高错误,This project is not a MyEclipse Hiberbate Project. Assuming Hibe ...
- 弹性盒布局(FlexBox)
CSS3引入了一种新的布局模式——Flexbox布局,即伸缩布局盒模型(Flexible Box),用来提供一个更加有效的方式制定.调整和分布一个容器里项目布局,即使它们的大小是未知或者动态的,这里简 ...
- OpenBSD内核之引导PBR
OpenBSD引导的第二部PBR,也是活动分区的一个扇区的代码,由第一步的MBR加载到0x7C00处,manpage里详细的讲解了过程和大致实现 biosboot(8) (http://man.ope ...
- oracle 中数据库完全导入导出:cmd命令行模式(转载)
http://www.3lian.com/edu/2012/12-01/47252.html Oracle数据导入导出imp/exp就相当于oracle数据还原与备份.exp命令可以把数据从远程数据库 ...
- KMP算法解析(转自图灵社区)
KMP算法是一个很精妙的字符串算法,个人认为这个算法十分符合编程美学:十分简洁,而又极难理解.笔者算法学的很烂,所以接触到这个算法的时候也是一头雾水,去网上看各种帖子,发现写着各种KMP算法详解的转载 ...
- 前端面试题2016--CSS
介绍一下标准的CSS的盒子模型?低版本IE的盒子模型有什么不同的? (1)有两种,IE 盒子模型.W3C 盒子模型:(2)盒模型:内容(content).填充(padding).边界(margin). ...
- Javascript设计模式(摘译)
说明: 未完成...更新中.... 一.javascipt设计模式分类 设计模式分类有很多标准,最流行的三种如下 1) creational -- 主要关注对象创建 Creational des ...
- Sql Server Alter语句
1:向表中添加字段 Alter table [表名] add [列名] 类型 2: 删除字段 Alter table [表名] drop column [列名] 3: 修改表中字段类型 (可以修 ...
- 时间控件之赋值问题:datetimebox
1.datetimebox不显示毫秒数: <input class="easyui-datetimebox" name="birthday" data-o ...
- DEVEXPRESS控件使用技巧记录-GRIDCONTROL
1. 存在父子表的时候,只允许父表的一条记录展开,其他记录都收起 方法:Feature Browser -> Master-Detail -> behavior -> AllowOn ...