[xsy1294]sub
给出一棵$N$个节点的无根树,节点$i$有权值$v_i$。现在有$M$次操作,操作有如下两种:
$1\ x\ y$ 将节点$x$的权值$v_x$修改为$y$
$2$ 选择一个联通块(也可以不选择),使得点权和最大。输出这个点权和
树剖==
考虑先做一次DP,$f_x$表示在以$x$为根的子树中,选择$x$的情况下的最大联通块和,那么$f_x=v_x+\sum\limits_{u\in son_x}\max(f_u,0)$($v_x$表示点$x$的权值)
然后树剖,每一个点的权值设为$f_x-f_{heavy_x}$($heavy_x$表示点$x$的重儿子),然后可以发现,对于一条重链所代表的区间,求最大子段和就是这条链上所有点的最大答案
为了避免不同重链之间互相影响,我们可以在重链之间加上一个空节点,权值为$-\infty$
考虑修改,修改一个点只会影响它沿重链往上跳时经过的(重链顶端节点的父亲节点),所以我们只需要修改重链顶端的父亲节点就好了
#include<stdio.h>
typedef long long ll;
const ll inf=1000000000ll;
ll max(ll a,ll b){return a>b?a:b;}
struct zt{
ll s,ms,ls,rs;
}f[400010];
zt merge(zt l,zt r){
zt c;
c.s=l.s+r.s;
c.ls=max(l.ls,l.s+r.ls);
c.rs=max(r.rs,r.s+l.rs);
c.ms=max(max(l.ms,r.ms),l.rs+r.ls);
return c;
}
void pushup(int x){f[x]=merge(f[x<<1],f[x<<1|1]);}
ll p[200010];
void build(int l,int r,int x){
if(l==r){
f[x].s=p[l];
f[x].ms=f[x].ls=f[x].rs=max(p[l],0);
return;
}
int mid=(l+r)>>1;
build(l,mid,x<<1);
build(mid+1,r,x<<1|1);
pushup(x);
}
void modify(int p,ll v,int l,int r,int x){
if(l==r){
f[x].s+=v;
f[x].ms=f[x].ls=f[x].rs=max(f[x].s,0);
return;
}
int mid=(l+r)>>1;
if(p<=mid)
modify(p,v,l,mid,x<<1);
else
modify(p,v,mid+1,r,x<<1|1);
pushup(x);
}
zt query(int L,int R,int l,int r,int x){
if(L<=l&&r<=R)return f[x];
int mid=(l+r)>>1;
if(R<=mid)return query(L,R,l,mid,x<<1);
if(L>mid)return query(L,R,mid+1,r,x<<1|1);
return merge(query(L,R,l,mid,x<<1),query(L,R,mid+1,r,x<<1|1));
}
int h[100010],nex[200010],to[200010],v[100010],fa[100010],siz[100010],son[100010],bl[100010],pos[100010],M;
ll dp[100010];
void add(int a,int b){
M++;
to[M]=b;
nex[M]=h[a];
h[a]=M;
}
void dfs(int x){
int i,mx=0,k=0;
siz[x]=1;
dp[x]=v[x];
for(i=h[x];i;i=nex[i]){
if(to[i]!=fa[x]){
fa[to[i]]=x;
dfs(to[i]);
dp[x]+=max(dp[to[i]],0);
siz[x]+=siz[to[i]];
if(siz[to[i]]>mx){
mx=siz[to[i]];
k=to[i];
}
}
}
son[x]=k;
}
void dfs(int x,int chain){
pos[x]=++M;
bl[x]=chain;
if(son[x]){
dp[x]-=max(dp[son[x]],0);
dfs(son[x],chain);
}
p[pos[x]]=dp[x];
for(int i=h[x];i;i=nex[i]){
if(to[i]!=fa[x]&&to[i]!=son[x]){
M++;
p[M]=-inf;
dfs(to[i],to[i]);
}
}
}
void modify(int x,ll d){
d-=v[x];
v[x]+=d;
ll t1,t2;
while(x){
t1=query(pos[bl[x]],M,1,M,1).ls;
modify(pos[x],d,1,M,1);
t2=query(pos[bl[x]],M,1,M,1).ls;
d=t2-t1;
x=fa[bl[x]];
}
}
int main(){
int n,m,i,x,y;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)scanf("%d",v+i);
for(i=1;i<n;i++){
scanf("%d%d",&x,&y);
add(x,y);
add(y,x);
}
dfs(1);
M=0;
dfs(1,1);
build(1,M,1);
while(m--){
scanf("%d",&i);
if(i==1){
scanf("%d%d",&x,&y);
modify(x,y);
}else
printf("%d\n",f[1].ms);
}
}
[xsy1294]sub的更多相关文章
- 【XSY1294】sub 树链剖分
题目描述 给你一棵\(n\)个点的无根树,节点\(i\)有权值\(v_i\).现在有\(m\)次操作,操作有如下两种: \(1~x~y\):把\(v_x\)改成\(y\). \(2\):选择一个连通块 ...
随机推荐
- [ZJOI2007]棋盘制作 (单调栈)
[ZJOI2007]棋盘制作 题目描述 国际象棋是世界上最古老的博弈游戏之一,和中国的围棋.象棋以及日本的将棋同享盛名.据说国际象棋起源于易经的思想,棋盘是一个8 \times 88×8大小的黑白相间 ...
- tomcat内存配置(二)
Tomcat本身不能直接在计算机上运行,需要依赖于硬件基础之上的操作系统和一个Java虚拟机.Tomcat的内存溢出本质就是JVM内存溢出,所以在本文开始时,应该先对JavaJVM有关内存方面的知识进 ...
- source改变当前路径
转摘自:http://hi.baidu.com/homappy/item/90e416525d2faf958c12edb7 Shell 脚本执行有三种方法 bash 脚本名 sh 脚本名 chmod ...
- JavaScript 知识点之escape()与decodeURI()
解释:escape() 函数可对字符串进行编码,这样就可以在所有的计算机上读取该字符串. 语法:escape(string)参数描述string必需.要被转义或编码的字符串. 返回值:已编码的 str ...
- ES6学习笔记(一)——Promise
Promise 是 ES6 提供的一种异步编程的解决方案: 将异步操作以同步操作的流程表达出来,避免了层层嵌套的回调函数(解决异步函数回调地狱的问题).Promise 对象保存着异步操作的结果. 首先 ...
- Spring 中 AbstractExcelView 支持根据模板生成Excel文件. 通过设置 view 的 URL 属性指定模板的路径
注意:1. 模板需放在 WEB-INF 目录下2. 指定模板路径时不需要添加扩展名, Spring将自动添加 .xls 到URL 属性中.3. 在指定URL前需先设置 view 的 Applicat ...
- sql数据库的链接方式
今天看见了一个数据库的链接方法,给转载了,记得我刚刚学DAO的时候老是要记载这些东西,所以就上博客园上面看了看,就转过来了... MySQL: String Driver="com.mysq ...
- finally return 执行顺序问题
网上有很多人探讨Java中异常捕获机制try...catch...finally块中的finally语句是不是一定会被执行?很多人都说不是,当然他们的回答是正确的,经过我试验,至少有两种情况下fina ...
- SDK登录cognos
通过SDK登录cognos 一种是拼xml,如这里的实现https://github.com/cosysoft/cognos-tools/blob/master/src/com/ibm/cognos/ ...
- bzoj 1878 SDOI2009树状数组 离线操作
本来想写2120的,结果想起来了这个 我们先对于询问左端点排序,用树状数组存区间字母个数,对于每种字母, 第一次出现的位置记录为1,剩下的记录为0,然后记录下,每种颜色 后面第一个和他相同颜色的位置 ...