\(\mathcal{Description}\)

  Link.

  给定一棵 \(n\) 个结点的带权树,\(m\) 次单点点权修改,求出每次修改后的带权最大独立集。

  \(n,m\le10^5\)。

\(\mathcal{Solution}\)

  不考虑修改,显然 DP。令 \(f(u,0/1)\) 表示选 / 不选结点 \(u\),\(u\) 子树内的带权最大独立集。那么:

\[\begin{cases}f(u,0)=\sum_v\max\{f(v,0),f(v,1)\}\\f(u,1)=\sum_vf(v,0)+a_u\end{cases}
\]

  引入修改,我们自然想用数据结构维护转移。那么就需要进行树链剖分(不一定是重链剖分,这里以复杂度更小的 LCT 为例)。假设 \(u\) 的实儿子是 \(s_u\),我们单独维护 \(s_u\) 的贡献,而把其它儿子一起考虑。设 \(g(u,0/1)\) 表示选 / 不选结点 \(u\),结点 \(u\) 及其虚儿子的子树们的带权最大独立集。有:

\[\begin{cases}g(u,0)=\sum_{v\not=s_u}\max\{f(v,0),f(v,1)\}\\g(u,1)=\sum_{v\not=s_u}f(v,0)+a_u\end{cases}
\]

  用它来表示 \(f\):

\[\begin{cases}f(u,0)=g(u,0)+\max\{f(s_u,0),f(s_u,1)\}\\f(u,1)=g(u,1)+f(s_u,0)\end{cases}
\]

  写成矩乘:

\[\begin{bmatrix}f(u,0)\\f(u,1)\end{bmatrix}=\begin{bmatrix}g(u,0)&g(u,0)\\g(u,1)&-\infty\end{bmatrix}\begin{bmatrix}f(s_u,0)\\f(s_u,1)\end{bmatrix}
\]

  令 \(G_u=\begin{bmatrix}g(u,0)&g(u,0)\\g(u,1)&-\infty\end{bmatrix}\),每次修改,仅会影响 \(\mathcal O(\log n)\) 个链头,也即是 \(\mathcal O(\log n)\) 个虚实交替的位置的 \(G\) 需要修改,就可以维护了。

  考虑到本题的 LCT 并没有“动”,其实就不需要打翻转标记。仅在 access 操作时虚实儿子会发生变化(ch[x][1] 变虚,y 变实),那么直接在当前结点的矩阵除去 y 的贡献,加上 ch[x][1] 的贡献即可。初始时,只初始化 fa[] 的信息,这样所有边都是虚边,\(f(u,0/1)=g(u,0/1)\),做一遍无修改的 DP 把初始的矩阵挂到树上就行啦。

\(\mathcal{Code}\)

#include <cstdio>
#include <cstring>
#include <assert.h> const int MAXN = 1e5, NINF = 0xc0c0c0c0;
int n, m, ecnt, a[MAXN + 5], head[MAXN + 5], f[MAXN + 5][2];
int fa[MAXN + 5], ch[MAXN + 5][2]; inline int max_ ( const int a, const int b ) { return a < b ? b : a; } inline char fgc () {
static char buf[1 << 17], *p = buf, *q = buf;
return p == q && ( q = buf + fread ( p = buf, 1, 1 << 17, stdin ), p == q )
? EOF : *p ++;
} inline int rint () {
int x = 0, f = 1; char s = fgc ();
for ( ; s < '0' || '9' < s; s = fgc () ) f = s == '-' ? -f : f;
for ( ; '0' <= s && s <= '9'; s = fgc () ) x = x * 10 + ( s ^ '0' );
return x * f;
} inline void wint ( int x ) {
if ( x < 0 ) putchar ( '-' ), x = -x;
if ( 9 < x ) wint ( x / 10 );
putchar ( x % 10 ^ '0' );
} struct Edge { int to, nxt; } graph[MAXN * 2 + 5]; struct Matrix {
int n, m, mat[2][2];
Matrix () {}
Matrix ( const int tn, const int tm ): n ( tn ), m ( tm ), mat {} {}
inline int* operator [] ( const int key ) { return mat[key]; }
inline Matrix operator * ( Matrix t ) {
assert ( m == t.n );
Matrix ret ( n, t.m ); memset ( ret.mat, 0xc0, sizeof ret.mat );
for ( int i = 0; i < n; ++ i ) {
for ( int k = 0; k < m; ++ k ) {
for ( int j = 0; j < t.m; ++ j ) {
ret[i][j] = max_ ( ret[i][j], mat[i][k] + t[k][j] );
}
}
}
return ret;
}
} G[MAXN + 5], S[MAXN + 5]; inline void link ( const int s, const int t ) {
graph[++ ecnt] = { t, head[s] };
head[s] = ecnt;
} inline bool nroot ( const int x ) { return ch[fa[x]][0] == x || ch[fa[x]][1] == x; } inline void pushup ( const int x ) {
S[x] = G[x];
if ( ch[x][0] ) S[x] = S[ch[x][0]] * S[x];
if ( ch[x][1] ) S[x] = S[x] * S[ch[x][1]];
} inline void rotate ( const int x ) {
int y = fa[x], z = fa[y], k = ch[y][1] == x;
fa[x] = z; if ( nroot ( y ) ) ch[z][ch[z][1] == y] = x;
ch[y][k] = ch[x][k ^ 1]; if ( ch[x][k ^ 1] ) fa[ch[x][k ^ 1]] = y;
pushup ( ch[fa[y] = x][k ^ 1] = y ), pushup ( x );
} inline void splay ( const int x ) {
for ( int y; nroot ( x ); rotate ( x ) ) {
if ( nroot ( y = fa[x] ) ) {
rotate ( x ^ y ^ ch[y][0] ^ ch[fa[y]][0] ? x : y );
}
}
} inline void access ( int x ) {
for ( int y = 0; x; x = fa[y = x] ) {
splay ( x );
if ( ch[x][1] ) {
G[x][0][0] += max_ ( S[ch[x][1]][0][0], S[ch[x][1]][1][0] );
G[x][1][0] += S[ch[x][1]][0][0];
}
if ( y ) {
G[x][0][0] -= max_ ( S[y][0][0], S[y][1][0] );
G[x][1][0] -= S[y][0][0];
}
G[x][0][1] = G[x][0][0];
ch[x][1] = y, pushup ( x );
}
} inline void initDP ( const int u, const int fath ) {
fa[u] = fath, f[u][1] = a[u], G[u] = Matrix ( 2, 2 );
for ( int i = head[u], v; i; i = graph[i].nxt ) {
if ( ( v = graph[i].to ) ^ fath ) {
initDP ( v, u );
f[u][0] += max_ ( f[v][0], f[v][1] );
f[u][1] += f[v][0];
}
}
G[u][0][0] = G[u][0][1] = f[u][0];
G[u][1][0] = f[u][1], G[u][1][1] = NINF;
S[u] = G[u];
} int main () {
n = rint (), m = rint ();
for ( int i = 1; i <= n; ++ i ) a[i] = rint ();
for ( int i = 1, u, v; i < n; ++ i ) {
u = rint (), v = rint ();
link ( u, v ), link ( v, u );
}
initDP ( 1, 0 );
for ( int x, y; m --; ) {
x = rint (), y = rint ();
access ( x ), splay ( x );
G[x][1][0] += y - a[x], a[x] = y;
pushup ( x ), splay ( 1 );
wint ( max_ ( S[1][0][0], S[1][1][0] ) ), putchar ( '\n' );
}
return 0;
}

Solution -「洛谷 P4719」「模板」"动态 DP" & 动态树分治的更多相关文章

  1. 【洛谷】P3919 【模板】可持久化线段树(主席树)

    题目 传送门:QWQ 分析 主席树的模板,囤着 代码 #include <bits/stdc++.h> using namespace std; ; ], rs[N*], root[N*] ...

  2. 洛谷:P2292 [HNOI2004]L语言(DP+Trie树)

    P2292 [HNOI2004]L语言 题目链接:https://www.luogu.org/problemnew/show/P2292 题目描述 标点符号的出现晚于文字的出现,所以以前的语言都是没有 ...

  3. 【洛谷5298】[PKUWC2018] Minimax(树形DP+线段树合并)

    点此看题面 大致题意: 有一棵树,给出每个叶节点的点权(互不相同),非叶节点\(x\)至多有两个子节点,且其点权有\(p_x\)的概率是子节点点权较大值,有\(1-p_x\)的概率是子节点点权较小值. ...

  4. 【洛谷4585】[FJOI2015] 火星商店问题(线段树分治)

    点此看题面 大致题意: 有\(n\)家店,每个商品有一个标价.每天,都可能有某家商店进货,也可能有某人去购物.一个人在购物时,会于编号在区间\([L_i,R_i]\)的商店里挑选一件进货\(d_i\) ...

  5. 「区间DP」「洛谷P1043」数字游戏

    「洛谷P1043」数字游戏 日后再写 代码 /*#!/bin/sh dir=$GEDIT_CURRENT_DOCUMENT_DIR name=$GEDIT_CURRENT_DOCUMENT_NAME ...

  6. 「 洛谷 」P2768 珍珠项链

    珍珠项链 题目限制 内存限制:125.00MB 时间限制:1.00s 标准输入输出 题目知识点 动态规划 \(dp\) 矩阵 矩阵乘法 矩阵加速 矩阵快速幂 题目来源 「 洛谷 」P2768 珍珠项链 ...

  7. 「 洛谷 」P4539 [SCOI2006]zh_tree

    小兔的话 推荐 小兔的CSDN [SCOI2006]zh_tree 题目限制 内存限制:250.00MB 时间限制:1.00s 标准输入输出 题目知识点 思维 动态规划 \(dp\) 区间\(dp\) ...

  8. 「 洛谷 」P2151 [SDOI2009]HH去散步

    小兔的话 欢迎大家在评论区留言哦~ HH去散步 题目限制 内存限制:125.00MB 时间限制:1.00s 标准输入 标准输出 题目知识点 动态规划 \(dp\) 矩阵 矩阵乘法 矩阵加速 矩阵快速幂 ...

  9. 「BZOJ2733」「洛谷3224」「HNOI2012」永无乡【线段树合并】

    题目链接 [洛谷] 题解 很明显是要用线段树合并的. 对于当前的每一个连通块都建立一个权值线段树. 权值线段树处理操作中的\(k\)大的问题. 如果需要合并,那么就线段树暴力合并,时间复杂度是\(nl ...

随机推荐

  1. nvm切换node版本出现乱码 exit status 1:

    nvm切换nodejs版本出现exit status 1:乱码 跟着网上的教程一步一步做,还是出现问题.浪费一下午的时间 最后发现却因为我没用CMD管理员权限运行 扑街 解决方法: 用管理员身份运行就 ...

  2. SQL高级优化(三)之存储引擎

    一.MySQL数据库引擎简介 1. ISAM(indexed Sequential Access Method) ​ ISAM 是一个定义明确且历经时间考验的数据表格管理方法,它在设计之时就考虑到数据 ...

  3. Spark本地环境实现wordCount单词计数

    注:图片如果损坏,点击文章链接:https://www.toutiao.com/i6814778610788860424/ 编写类似MapReduce的案例-单词统计WordCount 要统计的文件为 ...

  4. Jeesite富文本编辑框ckeditor显示错误

    Jeesite富文本编辑框ckeditor显示错误 原文链接:https://www.toutiao.com/i6485135618190869005/ Jeesite中Control都会继承一个Ba ...

  5. Linux防火墙--IPtables企业级配置策略思路

    一.防火墙简介 防火墙定义:是通过有机结合各类用于安全管理与筛选的软件和硬件设备,帮助计算机网络于其内.外网之间构建一道相对隔绝的保护屏障,以保护用户资料与信息安全性的一种技术. 防火墙发展应用:最早 ...

  6. 【Java常用类】StringBuffer、StringBuilder

    Stringbuffer.StringBuilder String.StringBuffer.StringBuilder三者的异同? String:不可变的字符序列:底层使用char[]存储 Stri ...

  7. spring 事务传播性

    一.什么是事务传播性 大白话讲就是,方法之间互相调用的时候,事务如何传播,比如A()调用B(),B()的事务是和A()共用一个事务(失败一起提交)? 还是新事务(两者事务互不影响)?,还是说B()不需 ...

  8. 集合框架-Map集合-HashMap存储自定义对象

    1 package cn.itcast.p6.hashmap.demo; 2 3 import java.util.HashMap; 4 import java.util.Iterator; 5 im ...

  9. WPS修改批注部分的字体颜色?

    今天遇到一个问题,就是复制文档的时候有几块红色字体想改成黑色,怎么也改不成功,通过修改字体颜色无效,通过百度找到了解决方法记录一下. 解决方法 审阅--显示标记--点击插入和删除(去掉前面的对钩即可) ...

  10. SaltStack 的基本概念与工作原理 架构设计

    随着云计算技术的快速普及与发展,越来越多的企业开始学习和搭建自己的云平台代替传统的 IT 交付模式,企业的 IT 环境也随之越来越复杂,常规的运维方法与技术已经无法满足现在云环境中系统的配置与变更.基 ...