fzu2082 过路费 树链剖分

2024-06-14 09:08

本文主要是介绍fzu2082 过路费 树链剖分,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题意:中文题。。

思路:树链剖分后,在线段树上单点更新,区间查询。详见代码:

/*********************************************************file name: fzu2082.cppauthor : kereocreate time:  2015年01月21日 星期三 09时21分42秒
*********************************************************/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<set>
#include<map>
#include<vector>
#include<stack>
#include<cmath>
#include<string>
#include<algorithm>
using namespace std;
typedef long long ll;
const int sigma_size=26;
const int N=100+50;
const int MAXN=50000+50;
const int inf=0x3fffffff;
const double eps=1e-8;
const int mod=100000000+7;
#define L(x) (x<<1)
#define R(x) (x<<1|1)
#define PII pair<int, int>
#define mk(x,y) make_pair((x),(y))
int n,m,edge_cnt,cnt;
int head[MAXN],sz[MAXN],son[MAXN],fa[MAXN],top[MAXN],dep[MAXN],pos[MAXN];
struct Edge{int u,v,w,next;
}edge[MAXN<<1];
struct node{int l,r;ll sum;
}segtree[MAXN<<2];
void init(){edge_cnt=cnt=0;memset(head,-1,sizeof(head));
}
void addedge(int u,int v,int w){edge[edge_cnt].u=u; edge[edge_cnt].v=v;  edge[edge_cnt].w=w; edge[edge_cnt].next=head[u]; head[u]=edge_cnt++;
}
void dfs1(int u,int pre,int depth){sz[u]=1; son[u]=0; dep[u]=depth; fa[u]=pre;for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].v;if(v == pre)continue;dfs1(v,u,depth+1);sz[u]+=sz[v];if(sz[son[u]]<sz[v])son[u]=v;}
}
void dfs2(int u,int tp){pos[u]=cnt++; top[u]=tp;if(son[u]!=0)dfs2(son[u],top[u]);for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].v;if(v == fa[u] || v == son[u])continue;dfs2(v,v);}
}
void build(int rt,int l,int r){segtree[rt].l=l; segtree[rt].r=r; segtree[rt].sum=0;if(l == r)return ;int mid=(l+r)>>1;build(L(rt),l,mid); build(R(rt),mid+1,r);
}
void push_up(int rt){segtree[rt].sum=segtree[L(rt)].sum+segtree[R(rt)].sum;
}
void update(int rt,int p,int x){if(segtree[rt].l == segtree[rt].r){segtree[rt].sum=x;return ;}int mid=(segtree[rt].l+segtree[rt].r)>>1;if(p<=mid)update(L(rt),p,x);else update(R(rt),p,x);push_up(rt);
}
ll query(int rt,int l,int r){if(segtree[rt].l == l && segtree[rt].r == r)return segtree[rt].sum;int mid=(segtree[rt].l+segtree[rt].r)>>1;if(r<=mid)return query(L(rt),l,r);else if(l>mid)return query(R(rt),l,r);else return query(L(rt),l,mid)+query(R(rt),mid+1,r);
}
ll solve(int u,int v){ll ans=0;while(top[u]!=top[v]){if(dep[top[u]]<dep[top[v]])swap(u,v);ans+=query(1,pos[top[u]],pos[u]);u=fa[top[u]];}if(dep[u]>dep[v])swap(u,v);if(u!=v)ans+=query(1,pos[u]+1,pos[v]);return ans;
}
int main(){while(~scanf("%d%d",&n,&m)){init();for(int i=1;i<n;i++){int u,v,w;scanf("%d%d%d",&u,&v,&w);addedge(u,v,w); addedge(v,u,w);}dfs1(1,1,1); dfs2(1,1); build(1,1,n-1);for(int i=0;i<edge_cnt;i+=2){int u=edge[i].u,v=edge[i].v,w=edge[i].w;if(dep[u]<dep[v])swap(u,v);update(1,pos[u],w);}for(int i=0;i<m;i++){int id,u,v;scanf("%d",&id);if(id == 0){int x,w;scanf("%d%d",&x,&w);x=(x-1)*2;u=edge[x].u; v=edge[x].v;if(dep[u]<dep[v])swap(u,v);update(1,pos[u],w);}else{scanf("%d%d",&u,&v);printf("%lld\n",solve(u,v));}}}return 0;
}


这篇关于fzu2082 过路费 树链剖分的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1059975

相关文章

树链剖分 + 后缀数组 - E. Misha and LCP on Tree

E. Misha and LCP on Tree  Problem's Link   Mean:  给出一棵树,每个结点上有一个字母。每个询问给出两个路径,问这两个路径的串的最长公共前缀。 analyse: 做法:树链剖分+后缀数组. 记录每条链的串,正反都需要标记,组成一个长串. 然后记录每条链对应的串在大串中的位置,对大串求后缀数组,最后询问就是在一些链

【HDU】5029 Relief grain 树链剖分+离线标记法

传送门:【HDU】5029 Relief grain 题目分析:这题的方法太美妙了! 首先看到这是一颗树,那么很容易想到树链剖分。然后想到可以将查询排个序,然后每一种查询执行完以后更新每个点的最优值。但是这样进行的复杂度太高!尤其是当z给的没有一样的时候尤其如此。 那么我们是否可以找到更加高效的算法? 答案是肯定的! 先简化一下问题,如果这些操作是在线段上进行的,我们怎么求解?

【HDU】5221 Occupation【树链剖分】

传送门:【HDU】5221 Occupation 题目分析: 最直接的想法,用一棵树链剖分维护路径,一棵dfs序线段树维护子树。因为每次最多修改一个点,所以修改的时候我们暴力修改每个点就可以了。 my  code: my~~code: #pragma comment(linker, "/STACK:102400000,102400000")#include <stdio.h>#in

【HDU】5436 Transmigration tree【树链剖分+dp+rmq】

题目链接:【HDU】5436 Transmigration tree #pragma comment(linker, "/STACK:16777216")#include <stdio.h>#include <string.h>#include <vector>#include <algorithm>using namespace std ;typedef long long LL ;

【HDU】5566 Clarke and room【树链剖分+AC自动机】

题目链接:Clarke and room #include <bits/stdc++.h>using namespace std ;typedef long long LL ;#define ls ( o << 1 )#define lson ls , l , m#define rs ( o << 1 | 1 )#define rson rs , m + 1 , r#define rt

Open3D mesh 模型精细化处理--中点剖分

目录 一、概述 1.1原理 1.2实现步骤 二、代码实现 2.1关键函数 输入参数 输出参数 三、实现效果 3.1原始mesh 3.2精细化mesh Open3D点云算法汇总及实战案例汇总的目录地址: Open3D点云算法与点云深度学习案例汇总(长期更新)-CSDN博客 一、概述         在三维模型处理过程中,精细化处理(subdivision)是一个

Dominant Indices【模板】长链剖分

~~~~~      Dominant Indices ~~~~~      总题单链接 对比重链剖分和长链剖分 ~~~~~      什么是重链剖分:根据子树的大小将树拆成多条互不相交的链。 ~~~~~      什么是长链剖分:根据字数的深度将树拆成多条互不相交的链。 ~~~~~      重链剖分中的重儿子:子树大小最大的儿子。 ~~~~~      长链剖分中的重儿子:

hdu5052 Yaoge’s maximum profit 树链剖分

一棵树上,从u走到v,在某点买入,咋之后的某点卖出,求最大利润。 维护正着走和反着走的最大利润。 #include<iostream>#include<cstdio>#include<cmath>#include<algorithm>#include<cstring>#include<vector>#include<set>#include<map>#include<que

poj3237--Tree(树链剖分+线段树)

题目链接:poj--3237 题意很简单,给出n个节点的一棵树,有三种操作: 1、C修改第i条边的值为v 2、N改变节点a到b内边的权值的符号(取反) 3、Q询问节点a到b内权值的最大值 首先树链剖分,将边整合到线段树上,线段树数组cl,因为存在取反操作,所以最大值可能是由最小值取反得到,所以记录最大和最小值,cl[i][0]记录第i段的最大值,cl[i][1]记录最小值,lazy做标记

BZOJ 3531 旅行【树链剖分】

[Sdoi2014] 简单的树链剖分 以每个信仰应该建一个线段树,空间复杂度为 O(5×1010) O(5\times 10^{10}) 因此会爆空间,所以需要动态申请空间。 // whn6325689// Mr.Phoebe// http://blog.csdn.net/u013007900#include <algorithm>#include <