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的更多相关文章
随机推荐
- [分享] 封装工具ES4配置文件解释
[分享] 封装工具ES4配置文件解释 LiQiang 发表于 2015-2-3 14:41:21 https://www.itsk.com/thread-346132-1-4.html [分享] 封装 ...
- 虚拟机EAL: Error reading from file descriptor
这个是虚拟机安装固有的BUG,代码差异如下: diff --git a/lib/librte_eal/linuxapp/igb_uio/igb_uio.c b/lib/librte_eal/linux ...
- Java源文件编译成功但是运行时加载不到文件
最近系统重装了一些,Java等环境变量都需要重新配置,配置好以后编写了一个Java源文件编译了一下,通过Javac编译源文件,编译成功,但是再通过Java运行时没找到报出找不到加载文件或者加载文件不存 ...
- gdb 常用内容
gdb exegdb exe coregdb -p info m TAB ^関数の先頭 info b ^list the breakpoint set args -a test ^引数設定 show ...
- map erase iterator
错误写法: map<int, int> m; for (map<int, int>::iterator it = m.begin(); it != m.end(); it++) ...
- ASP.NET MVC多表单提交
多表单提交需要写清路径,以便主程序可以找到 方法一:直接写路径, action="~/Home/other1" "~"表示从根目录开始 方法2:@using ...
- 用PHP解析类JSON字符串为数组的实现
题目:把字符串嵌套关系转换成数组,字符串只包含成对中括号.数字和逗号字符串:(1,(1,2,(1,(1,2,(1)),3)),3,(1,(1,2,((1((1,(1,2,(1,2,3),4,5),3) ...
- jQuery 点击显示再次点击隐藏
<html> <head> <script type="text/javascript" src="/jquery/jquery.js&qu ...
- 安装Windows 10后PDF补丁丁等程序界面变得模糊的解决办法
对于使用高分辨率屏幕且屏幕缩放比例在 100%以上的用户,升级到 Windows 10 后将发现许多程序的界面,例如QQ.电脑管家.Windows本身的服务管理程序等等,都变得非常模糊,<PDF ...
- 深入理解iOS开发中的BitCode功能
前言 做iOS开发的朋友们都知道,目前最新的Xcode7,新建项目默认就打开了bitcode设置.而且大部分开发者都被这个突如其来的bitcode功能给坑过导致项目编译失败,而这些因为bitcode而 ...