poj2763
//Accepted 11676 KB 2344 ms
/*
source:poj2763
time :2015.5.29
by :songt
*/
/*题解:
树链剖分
基于边权,路径查询
wind第一次在s节点,假如她走到了p节点,那么下次开始的时候他就在p节点
*/
#include <cstdio>
#include <cstring>
;
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=;
}
int s[imax_n];
int cnt_s;
int vis[imax_n];
void dfs1(int u)
{
fa[u]=;
deep[u]=;
cnt_s=;
s[cnt_s++]=u;
memset(vis,,sizeof(vis));
while (cnt_s)
{
];
if (vis[x])
{
cnt_s--;
if (x!=u)
{
num[fa[x]]+=num[x];
|| num[fa[x]]<num[x])
son[fa[x]]=x;
}
continue;
}
vis[x]=;
num[x]=;
;i=next[i])
{
int v=edge[i].v;
if (v!=fa[x])
{
fa[v]=x;
deep[v]=deep[x]+;
s[cnt_s++]=v;
}
}
}
}
void dfs2(int u)
{
top[u]=u;
memset(vis,,sizeof(vis));
cnt_s=;
s[cnt_s++]=u;
while (cnt_s)
{
];
if (vis[x])
{
cnt_s--;
continue ;
}
vis[x]=;
p[x]=pos++;
fp[p[x]]=x;
) continue ;
;i=next[i])
{
int v=edge[i].v;
if (v!=fa[x] && v!=son[x])
{
top[v]=v;
s[cnt_s++]=v;
}
}
top[son[x]]=top[x];
s[cnt_s++]=son[x];
}
}
struct Tree
{
int l,r;
long long sum;
}f[imax_n*];
void build(int t,int l,int r)
{
f[t].l=l;
f[t].r=r;
f[t].sum=;
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;
return ;
}
;
*t,k,value);
*t+,k,value);
f[t].sum=f[*t].sum+f[*t+].sum;
}
long long query(int t,int l,int r)
{
if (f[t].l==l && f[t].r==r)
{
return f[t].sum;
}
;
*t,l,r);
else
{
*t+,l,r);
*t,l,mid)+query(*t+,mid+,r);
}
}
void swap(int &a,int &b)
{
int t=a;
a=b;
b=t;
}
long long find(int u,int v)
{
int f1=top[u],f2=top[v];
;
while (f1!=f2)
{
if (deep[f1]<deep[f2])
{
swap(f1,f2);
swap(u,v);
}
sum+=query(,p[f1],p[u]);
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]);
return sum;
}
];
int n,m;
int res;
int main()
{
)
//scanf("%d%d%d",&n,&m,&res);
{
init();
;i<n-;i++)
{
scanf(],&e[i][],&e[i][]);
addEdge(e[i][],e[i][]);
addEdge(e[i][],e[i][]);
}
dfs1();
dfs2();
build(,,pos-);
;i<n-;i++)
{
]]<deep[e[i][]]) swap(e[i][],e[i][]);
update(,p[e[i][]],e[i][]);
}
int kind;
int u,v;
;i<m;i++)
{
scanf("%d",&kind);
)
{
scanf("%d",&u);
printf("%lld\n",find(u,res));
res=u;
}
else
{
scanf("%d%d",&u,&v);
update(,p[e[u-][]],v);
}
}
}
;
}
poj2763的更多相关文章
- 【lct】poj2763 Housewife Wind
题意:给你一棵树,边带权,支持两种操作:修改某条边的权值:查询两点之间的最短路. lct主要实现单点修改和路径和. 修改x结点的值只需将x Splay到其所在辅助树的根,然后修改其值,再maintai ...
- poj2763(树链剖分 - 边权)
poj2763 题意 给定一个树形图,某人原来在 s 点,每条边(路)有通过的时间花费,有两种操作:1. 查询某人到 u 点花费的时间 2. 更新某条路的时间花费. 分析 权值在边上,可以把它们 &q ...
- POJ2763 Housewife Wind 树链剖分 边权
POJ2763 Housewife Wind 树链剖分 边权 传送门:http://poj.org/problem?id=2763 题意: n个点的,n-1条边,有边权 修改单边边权 询问 输出 当前 ...
- POJ2763 Housewife Wind
Time Limit: 4000MS Memory Limit: 65536K Total Submissions: 9701 Accepted: 2661 Description Aft ...
- poj2763 树链剖分(线段树)
注意这里都是把边放到线段树中,所以lca的时候,要注意如果top[x]==top[y] && x==y 的时候已经完成了. 仔细想想边和点的不同之处!!! #include<ma ...
- POJ2763 Housewife Wind(DFS序)
题目:单边修改,树链查询. 这题是边权,不是点权,不过也可以看作是点权. 然后其实就和BZOJ2819一样. #include<cstdio> #include<cstring> ...
- poj2763(树链剖分)
题目链接:http://poj.org/problem?id=2763 题意:定一棵带边权的树,要求支持两种操作:1)询问树中某两点间的距离. 2)修改某条边的权值. 分析:树链剖分,边权修改,路径求 ...
- 树链剖分——边权poj2763
边权操作起来也和点权一样,只要把边的权值映射到点上即可,要注意的地方是向上爬的过程中和点权不太一样,还有个特判(WA了几次..) 完整代码 #include<cstring> #inclu ...
- poj2763树链剖分边权+区间和
自己写的比原来的板子常数小了不少嘻嘻,边权处理起来比点权要复杂一下 由于根节点没有被映射的边,其被访问到的顺序是0,直接排除在线段树外 #include<iostream> #includ ...
随机推荐
- ES6(一)ECMAscript6介绍
nvm-windows Node.js是JavaScript语言的服务器运行环境,对ES6的支持度比浏览器更高.通过Node,可以体验更多ES6的特性.建议使用版本管理工具nvm,来安装Node,因为 ...
- chrome 跨域设置
右击chrome快捷方式,在启动指令后面添加--disable-web-security,然后保存.如下: "C:\Program Files (x86)\Google\Chrome\App ...
- Git常用命令速查表
- Python常用内置函数总结
一.数学相关 1.绝对值:abs(-1)2.最大最小值:max([1,2,3]).min([1,2,3])3.序列长度:len('abc').len([1,2,3]).len((1,2,3))4.取模 ...
- (转载)数据库出现ORA-00283/ORA-01610的问题
在这里需要感谢棉花糖给予无私帮助,真的谢谢他!http://blog.itpub.net/67668/viewspace-353270/处理过程可以参照http://www.itpub.net/vie ...
- 浅谈对Js闭包的理解
理解Js的闭包,首先让我们先看几个概念 执行环境(executive environment)每个函数都有自己的执行环境,匿名函数默认为全局环境. 作用域链(scope chain)子函数继承父函数, ...
- C++11 智能指针
C++ 11标准库引入了几种智能指针 unique_ptr shared_ptr weak_ptr C++内存管理机制是当一个变量或对象从作用域过期的时候就会从内存中将他干掉.但是如果变量只是一个指针 ...
- .net中excel遇到的一些问题
Excel.Application app; 错误 1 未能找到类型或命名空间名称“Excel”(是否缺少 using 指令或程序集引用?) 解决方法: 添加引用 Interop.Excel.dll ...
- [原创]WinForm分页控件制作
先简单说一下思路: 1.做一个分页控件的导航类,即记录总页数.当前页.每页记录数,下一页.上一页.跳转等操作的页数变更. class PageNavigation{/// <summary> ...
- bzoj 3211: 花神游历各国
#include<cstdio> #include<cmath> #include<iostream> #define M 100006 using namespa ...