洛谷 P6021 洪水
题意
- 给定一棵有 \(n\) 个结点的树,点有点权;一共有 \(m\) 次操作,每次操作包括以下两种:
- 在一个点的子树中删去一些结点,使得该子树中所有叶结点与该子树的根结点不连通,并且使删去的点的点权和最小,求出这个最小值。
- 修改一个点的点权。
- \(n \le 2 \times 10^5\) ,\(m\) 与 \(n\) 同级 (题目中没有说明QAQ)
题目分析
显然是动态 DP。
据说这题珂以线段树二分但是我太蒻只会暴力 DDP
虽然题目中没有说明,但是这棵树以 \(1\) 号结点为根。
如果没有修改操作,这是一道简单的树形动规。
首先考虑没有修改操作的情况。
记 \(w_u\) 表示点 \(u\) 的权值, \(f_u\) 表示以 \(u\) 为根的子树的最小答案。
那么显然有
\]
方程的意义为:要么删去结点 \(u\) ,要么不删去结点 \(u\) ,在结点 \(u\) 的孩子中删去结点,两者取最小。
特殊的,对于叶结点,\(f_u = w_u\) 。
然后考虑修改点权的操作。
注意到每次修改点权只会对该结点到根的路径上的 \(f\) 值产生影响。
所以每次修改时暴力从该结点一直修改到 \(1\) 号结点。
最坏时间复杂度 $ \Theta (mn)$ ,稳稳地 TLE 。
接下来就需要请出主角:动态 DP 了。
动态 DP
动态 DP 用于解决一类被套上修改点权操作的基础 DP 题目。
其主要思想就是利用线段树维护每个结点的 DP 值。
但是由于大多数状态转移方程不满足结合律,所以无法直接使用线段树维护。
但是我们知道矩阵乘法运算是满足结合律的。
所以我们就需要把状态转移方程用矩阵乘法的形式表达出来。
模板题相信大家都过了。
不会动态 DP 的同学可以先去写一下这道模板题,里面的题解介绍的很详细。
构造矩阵
终于回到这题了。
让我们再来玩一玩刚才的状态转移方程。
\]
由于这是一棵树,所以先对它进行重链剖分,然后 按照套路 把轻儿子的答案的总和拎出来;也就是说,记 \(g_u\) 表示以 \(u\) 为根的子树中轻儿子的答案总和,则可以将原方程转化为:
\]
然后考虑把这个方程转化为矩阵乘法的形式。
在这里,我们需要重新定义一下矩阵乘法:对于矩阵 \(A\) ,\(B\) ,若矩阵 \(C\) 满足 \(C=A*B\) ,则有:
\]
至于这个运算为什么满足结合律,请读者自行验证。
那么我们相当于要找到一个转移矩阵 \(U\) ,使得 :
\]
注意新定义的矩阵乘法不满足交换律。说的好像原来的矩阵乘法满足交换律一样。
首先 \(U\) 是一个 \(2 \times2\) 的矩阵。
根据新定义的矩阵乘法,可以得出:
\(f_u=\min(f_{son_u}+U_{1,1}, 0+U_{1,2})\)
\(0 = \min(f_{son_u}+U_{2,1}, 0+U_{2,2})\)
再结合原方程,即可得出
* \begin{pmatrix} g_u & w_u \\ \infty & 0 \end{pmatrix}
= \begin{pmatrix} f_{u} \\ 0 \end{pmatrix}\]
然后用线段树维护转移矩阵就可以了。
在实现时我使用了一个矩阵数组来记录转移矩阵,听说这样可以减常数。
另外由于我太蒻了,写的树剖常数大,所以这份代码需要吸氧才能通过。
Code
#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
typedef long long ll;
const int MAXN = 200010;
const ll INF = 1e17+1;
int n, m;;
struct edge{
int ne, to;
}e[MAXN<<1];
int fir[MAXN], num = 0;
inline void join(int a, int b)
{
e[++num].ne = fir[a];
fir[a] = num;
e[num].to = b;
}
struct mat{
ll ele[2][2];
mat(int type = 0)
{
for(int i=0;i<2;i++)
for(int j=0;j<2;j++)
ele[i][j] = INF;
if(type == 1) {
for(int i=0;i<2;i++)
ele[i][i] = 0;
}
}
ll& operator()(const int ix, const int iy){return ele[ix][iy];}
// 重定义的矩阵乘法
inline friend mat operator*(mat mx, mat my)
{
mat res(0);
for(int i=0;i<2;i++)
for(int k=0;k<2;k++)
for(int j=0;j<2;j++)
res(i,j)=min(res(i,j), mx(i,k)+my(k,j));
return res;
}
};
ll val[MAXN], f[MAXN];
mat g[MAXN];
int dep[MAXN], pa[MAXN], siz[MAXN], son[MAXN], top[MAXN], dfn[MAXN], rev[MAXN], end[MAXN], cnt = 0;
void dfs1(int u, int fa)
{
siz[u] = 1; dep[u] = dep[fa] + 1; pa[u] = fa;
for(int i=fir[u];i;i=e[i].ne)
{
int v = e[i].to;
if(v == fa) continue;
dfs1(v, u);
siz[u] += siz[v];
if(siz[son[u]] < siz[v]) son[u] = v;
}
}
void dfs2(int u, int t)
{
dfn[u] = ++cnt; rev[cnt] = u; top[u] = t; end[t] = max(end[t], cnt);
g[u](0, 0) = 0; g[u](0, 1) = val[u];
g[u](1, 0) = INF; g[u](1, 1) = 0;
f[u] = 0;
if(son[u]) {
dfs2(son[u], t);
f[u] += f[son[u]];
}
else g[u](0, 0) = INF; //为了保证叶结点转移的合法
for(int i=fir[u];i;i=e[i].ne)
{
int v = e[i].to;
if(v == pa[u] || v == son[u]) continue;
dfs2(v, v);
g[u](0, 0) += f[v]; //g 数组的转移不包括重儿子
f[u] += f[v];
}
if(!son[u]) f[u] = val[u]; // 特判叶结点
else f[u] = min(f[u], val[u]); // 记得两者取最小
}
struct {
int l, r;
mat val;
}t[MAXN<<2];
inline void pushUp(int k)
{
t[k].val = t[k<<1].val*t[k<<1|1].val;
}
void build(int l, int r, int k)
{
t[k].l = l; t[k].r = r;
if(l == r) {
t[k].val = g[rev[l]];
return ;
}
int mid = t[k].l+t[k].r>>1;
build(l, mid, k<<1);
build(mid+1, r, k<<1|1);
pushUp(k);
}
void update(int x, int k)
{
if(t[k].l == t[k].r) {
t[k].val = g[rev[x]]; //由于有在树外记录转移矩阵,直接赋值即可
return ;
}
int mid = t[k].l+t[k].r>>1;
if(x <= mid) update(x, k<<1);
else update(x, k<<1|1);
pushUp(k);
}
mat query(int x, int y, int k)
{
if(t[k].l == x && t[k].r == y) return t[k].val;
int mid = t[k].l+t[k].r>>1;
if(y <= mid) return query(x, y, k<<1);
else if(x >= mid+1) return query(x, y, k<<1|1);
else return query(x, mid, k<<1)*query(mid+1, y, k<<1|1);
}
inline void updatePath(int x, ll z)
{
val[x] += z;
g[x](0, 1) = val[x];
mat bef, aft;
while(x)
{
bef = query(dfn[top[x]], end[top[x]], 1);
// 记录更新之前的矩阵
update(dfn[x], 1);
aft = query(dfn[top[x]], end[top[x]], 1);
x = pa[top[x]];
g[x](0, 0) += min(aft(0, 0), aft(0, 1)) - min(bef(0, 0), bef(0, 1));
// 更新转移矩阵,为下一条链的修改做准备,注意这里需要先减去原先的值
}
}
inline ll querySubtree(int x)
{
mat res = query(dfn[x], end[x], 1); // 查询时只需要从当前点到链尾即可
return min(res(0, 0), res(0, 1));
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;++i)
scanf("%lld",&val[i]);
for(int i=1;i<n;i++)
{
int a, b;
scanf("%d%d",&a,&b);
join(a, b);
join(b, a);
}
dfs1(1, 0);
dfs2(1, 1);
for(int i=1;i<=n;i++)
end[i] = end[top[i]];
build(1, n, 1);
scanf("%d",&m);
while(m--)
{
char opt; int x; ll z;
cin>>opt;
if(opt == 'C') {
scanf("%d%lld",&x,&z);
updatePath(x, z);
}
else {
scanf("%d",&x);
printf("%lld\n",querySubtree(x));
}
}
return 0;
}
洛谷 P6021 洪水的更多相关文章
- Solution -「洛谷 P6021」洪水
\(\mathcal{Description}\) Link. 给定一棵 \(n\) 个点的带点权树,删除 \(u\) 点的代价是该点点权 \(a_u\).\(m\) 次操作: 修改单点点权. ...
- [洛谷3457][POI2007]POW-The Flood
洛谷题目链接:[POI2007]POW-The Flood 题意翻译 Description 你手头有一张该市的地图.这张地图是边长为 m∗n 的矩形,被划分为m∗n个1∗1的小正方形.对于每个小正方 ...
- 洛谷1640 bzoj1854游戏 匈牙利就是又短又快
bzoj炸了,靠离线版题目做了两道(过过样例什么的还是轻松的)但是交不了,正巧洛谷有个"大牛分站",就转回洛谷做题了 水题先行,一道傻逼匈牙利 其实本来的思路是搜索然后发现写出来类 ...
- 洛谷P1352 codevs1380 没有上司的舞会——S.B.S.
没有上司的舞会 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题目描述 Description Ural大学有N个职员,编号为1~N.他们有 ...
- 洛谷P1108 低价购买[DP | LIS方案数]
题目描述 “低价购买”这条建议是在奶牛股票市场取得成功的一半规则.要想被认为是伟大的投资者,你必须遵循以下的问题建议:“低价购买:再低价购买”.每次你购买一支股票,你必须用低于你上次购买它的价格购买它 ...
- 洛谷 P2701 [USACO5.3]巨大的牛棚Big Barn Label:二维数组前缀和 你够了 这次我用DP
题目背景 (USACO 5.3.4) 题目描述 农夫约翰想要在他的正方形农场上建造一座正方形大牛棚.他讨厌在他的农场中砍树,想找一个能够让他在空旷无树的地方修建牛棚的地方.我们假定,他的农场划分成 N ...
- 洛谷P1710 地铁涨价
P1710 地铁涨价 51通过 339提交 题目提供者洛谷OnlineJudge 标签O2优化云端评测2 难度提高+/省选- 提交 讨论 题解 最新讨论 求教:为什么只有40分 数组大小一定要开够 ...
- 洛谷P1371 NOI元丹
P1371 NOI元丹 71通过 394提交 题目提供者洛谷OnlineJudge 标签云端评测 难度普及/提高- 提交 讨论 题解 最新讨论 我觉得不需要讨论O long long 不够 没有取 ...
- 洛谷P1538迎春舞会之数字舞蹈
题目背景 HNSDFZ的同学们为了庆祝春节,准备排练一场舞会. 题目描述 在越来越讲究合作的时代,人们注意的更多的不是个人物的舞姿,而是集体的排列. 为了配合每年的倒计时,同学们决定排出——“数字舞蹈 ...
- 洛谷八月月赛Round1凄惨记
个人背景: 上午9:30放学,然后因为学校举办读书工程跟同学去书城选书,中午回来开始打比赛,下午又回老家,中间抽出一点时间调代码,回家已经8:50了 也许是7月月赛时“连蒙带骗”AK的太幸运然而因同学 ...
随机推荐
- java中的字符串数组
本文主要讲述java中的字符串数组 字符串数组的声明有如下几种形式: // 第一种方式:new // 注意在String的后面[]中不需要添加字符串数组的长度.否则报错. String[] arr_1 ...
- Django项目启动 AttributeError: ‘str‘ object has no attribute ‘decode‘ 问题
Watching for file changes with StatReloader Performing system checks... System check identified no i ...
- 自研ORM框架实现工作单元模式
平时我们使用事务,需要显示的Try Catch 并且开启事务 提交事务 异常回滚事务 三步骤,使用工作单元后则只需要Commit. 1.接口定义 1 /// <summary> 2 /// ...
- 有向图的拓扑排序——DFS
在有向图的拓扑排序--BFS这篇文章中,介绍了有向图的拓扑排序的定义以及使用广度优先搜索(BFS)对有向图进行拓扑排序的方法,这里再介绍另一种方法:深度优先搜索(DFS). 算法 考虑下面这张图: 首 ...
- 分享自己亲测过的Visualstudio 2019中开发Typescript时,设置自动编译生成js,无需手工运行命令生成的方法。
步骤1)右键web项目,添加 tsconfig.json文件. 步骤2)确保配置如下,编译版本可自行设置,这里主要关注编译目标目录和自动编译设置: { "compileOnSave" ...
- python之路53 ajax补充返回序列化数据,多对多创建三种方式,django内置序列化组件(drf前身),批量操作数据,自定义分页器,form组件
ajax补充说明 主要是针对回调函数args接收到的响应数据 1.后端request.is_ajax() 用于判断当前请求是否由ajax发出 2.后端返回的三板斧都会被args接收不再影响整个浏览器页 ...
- Coolify系列-解决WARNING: IPv4 forwarding is disabled. Networking will not work.以及开启防火墙端口
背景 我在windows电脑安装了一个VM,使用VM开启了Linux服务器,运行docker,然后遇到了这个报错. 解决 首先:在宿主机上执行 echo "net.ipv4.ip_forwa ...
- 异常概念&异常体系-异常分类
异常概念&异常体系 异常,就是不正常的意思.在生活中:医生说,你的身体某个部分有异常,该部位和正常相比有点不同,该部位的功能将受影响,在程序中的意思就是: 异常:指的是程序在执行过程中,出现的 ...
- ubuntu下yaml-cpp安装与使用
安装 从GitHub上下载源码编译安装:git clone https://github.com/jbeder/yaml-cpp.git: 进入源码目录并创建一个 build 目录:cd yaml-c ...
- 数据库系列:MySQL慢查询分析和性能优化
1 背景 我们的业务服务随着功能规模扩大,用户量扩增,流量的不断的增长,经常会遇到一个问题,就是数据存储服务响应变慢. 导致数据库服务变慢的诱因很多,而RD最重要的工作之一就是找到问题并解决问题. 下 ...