博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[xsy1294]sub
阅读量:6566 次
发布时间:2019-06-24

本文共 2378 字,大约阅读时间需要 7 分钟。

给出一棵$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
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

转载于:https://www.cnblogs.com/jefflyy/p/8681242.html

你可能感兴趣的文章
RequireJS进阶(二)
查看>>
我设计的网站的分布式架构
查看>>
linux extract rar files
查看>>
Knockout.Js官网学习(监控属性Observables)
查看>>
ORA-12514: TNS: 监听程序当前无法识别连接描述符中请求的服务解决
查看>>
azure之MSSQL服务性能测试
查看>>
Android BitmapFactory.Options
查看>>
前端构建:Less入了个门
查看>>
phonegap(cordova) 自己定义插件代码篇(三)----支付宝支付工具整合
查看>>
linux 批量进行:解压缩某一类压缩文件类型的文件
查看>>
激活modelsim se 10.4 时运行patch_dll.bat不能生成TXT
查看>>
17秋 软件工程 Alpha 事后诸葛亮会议
查看>>
线性空间
查看>>
疑似checkpoint堵塞数据库连接
查看>>
Node.js中针对中文的查找和替换无效的解决方法
查看>>
理解指针的关键
查看>>
如何查看Ubuntu下已安装包版本号
查看>>
我的那些年(2)~我毕业了
查看>>
VS2017 配置ImageMagick
查看>>
Hive任务优化--控制hive任务中的map数和reduce数
查看>>