热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

BZOJ4712洪水

永恒的一氧化碳大爷在讨论版里发了个单log的做法,还没来得及看……先写个链剖再说考虑没有修改的情况,我们可以树DP,f[x]表示把x的

永恒的一氧化碳大爷在讨论版里发了个单log的做法,还没来得及看……

先写个链剖再说

考虑没有修改的情况,我们可以树DP,f[x]表示把x的子树截断的最小代价,v[x]表示截断每个点的代价,然后我们求出每个点的s[x],代表x的所有儿子节点的f之和,f[x]就等于min(v[x],s[x])

考虑修改&#xff0c;一次修改操作后&#xff0c;v[x]会增加&#xff0c;f[x]可能增加&#xff0c;可能不变&#xff0c;如果f[x]增加了d&#xff0c;那么从x到根的一段路径上的点的s值和f值都会增加d&#xff0c;直到第一个有v[x]-s[x]<&#61;d的点&#xff0c;即s[x]在加d之后会变得比v[x]大&#xff0c;或原来v[x]就小于等于s[x]&#xff0c;则f[x]不再等于s[x]&#43;d&#xff0c;而改为v[x]&#xff0c;则增加量会改变一次&#xff0c;我们继续对顶上的路径重复上述过程&#xff0c;直到增加量为0或者到根

对于一次处理增加操作&#xff0c;我们可以用链剖加线段树&#xff0c;支持区间加和区间查询v[x]-s[x]的最小值&#xff0c;在线段树上爬来找到第一个v[x]-s[x]<&#61;d的点

由于修改操作的增加量是非负的&#xff0c;并且每次只更改一个点&#xff0c;所以至多有n&#43;m次s[x]变得比v[x]大的时刻&#xff0c;处理一个这样的时刻是O(log^2 n)的&#xff0c;所以总复杂度O((n&#43;m) log^2 n)

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
#define MAXN 200010
#define MAXM 1010
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
#define ll long long
struct vec{ll to;ll fro;
};
vec mp[MAXN*2];
ll tai[MAXN],cnt;
ll n,m;
ll v[MAXN],s[MAXN],f[MAXN];
ll fa[MAXN],son[MAXN],tp[MAXN],siz[MAXN],dfn[MAXN],ndf[MAXN],tim;
ll mn[MAXN<<2],ch[MAXN<<2];
inline void be(ll x,ll y){mp[&#43;&#43;cnt].to&#61;y;mp[cnt].fro&#61;tai[x];tai[x]&#61;cnt;
}
inline void bde(ll x,ll y){be(x,y);be(y,x);
}
void dfs(ll x){ll i,y;siz[x]&#61;1;bool flag&#61;0;s[x]&#61;INF;for(i&#61;tai[x];i;i&#61;mp[i].fro){y&#61;mp[i].to;if(!siz[y]){fa[y]&#61;x;dfs(y);if(!flag){flag&#61;1;s[x]&#61;0;}s[x]&#43;&#61;f[y];siz[x]&#43;&#61;siz[y];if(siz[son[x]]}
void dfs2(ll x,ll z){ll i,y;ndf[dfn[x]&#61;&#43;&#43;tim]&#61;x;tp[x]&#61;z;if(son[x]){dfs2(son[x],z);for(i&#61;tai[x];i;i&#61;mp[i].fro){y&#61;mp[i].to;if(!dfn[y]){dfs2(y,y);}}}
}
inline void ud(ll x){mn[x]&#61;min(mn[x<<1],mn[x<<1|1]);
}
inline void toch(ll x,ll y){ch[x]&#43;&#61;y;mn[x]-&#61;y;
}
inline void pd(ll x){if(ch[x]){toch(x<<1,ch[x]);toch(x<<1|1,ch[x]);ch[x]&#61;0;}
}
void build(ll x,ll y,ll z){if(y&#61;&#61;z){mn[x]&#61;v[ndf[y]]-s[ndf[y]];return ;}ll mid&#61;y&#43;z>>1;build(x<<1,y,mid);build(x<<1|1,mid&#43;1,z);ud(x);
}
ll change(ll x,ll y,ll z,ll l,ll r,ll cv){if(y&#61;&#61;z){toch(x,cv);if(mn[x]&#43;cv<&#61;cv){return ndf[y];}else{return 0;}}pd(x);ll t&#61;0;ll mid&#61;y&#43;z>>1;if(y&#61;&#61;l&&z&#61;&#61;r){if(mn[x]>cv){toch(x,cv);return 0;}}if(r<&#61;mid){t&#61;change(x<<1,y,mid,l,r,cv);}else if(l>mid){t&#61;change(x<<1|1,mid&#43;1,z,l,r,cv);}else{t&#61;change(x<<1|1,mid&#43;1,z,mid&#43;1,r,cv);if(!t){t&#61;change(x<<1,y,mid,l,mid,cv);}}ud(x);return t;
}
ll ask(ll x,ll y,ll z,ll p){if(y&#61;&#61;z){return mn[x];}pd(x);ll mid&#61;y&#43;z>>1;if(p<&#61;mid){return ask(x<<1,y,mid,p);}else{return ask(x<<1|1,mid&#43;1,z,p);}
}
void tochange(ll x,ll y){if(y<&#61;0||!x){return ;}while(x){ll t&#61;change(1,1,n,dfn[tp[x]],dfn[x],y);if(!t){x&#61;fa[tp[x]];}else{tochange(fa[t],ask(1,1,n,dfn[t])&#43;y);return ;}}
}
int main(){ll i,x,y;char o[2];scanf("%lld",&n);for(i&#61;1;i<&#61;n;i&#43;&#43;){scanf("%lld",&v[i]);}for(i&#61;1;i}/*
5
5 7 4 1 5
2 1
3 2
4 3
5 3
5
C 3 10
C 1 6
Q 1*/





推荐阅读
  • 本文主要解析了Open judge C16H问题中涉及到的Magical Balls的快速幂和逆元算法,并给出了问题的解析和解决方法。详细介绍了问题的背景和规则,并给出了相应的算法解析和实现步骤。通过本文的解析,读者可以更好地理解和解决Open judge C16H问题中的Magical Balls部分。 ... [详细]
  • 本文讨论了使用差分约束系统求解House Man跳跃问题的思路与方法。给定一组不同高度,要求从最低点跳跃到最高点,每次跳跃的距离不超过D,并且不能改变给定的顺序。通过建立差分约束系统,将问题转化为图的建立和查询距离的问题。文章详细介绍了建立约束条件的方法,并使用SPFA算法判环并输出结果。同时还讨论了建边方向和跳跃顺序的关系。 ... [详细]
  • CF:3D City Model(小思维)问题解析和代码实现
    本文通过解析CF:3D City Model问题,介绍了问题的背景和要求,并给出了相应的代码实现。该问题涉及到在一个矩形的网格上建造城市的情景,每个网格单元可以作为建筑的基础,建筑由多个立方体叠加而成。文章详细讲解了问题的解决思路,并给出了相应的代码实现供读者参考。 ... [详细]
  • 开发笔记:实验7的文件读写操作
    本文介绍了使用C++的ofstream和ifstream类进行文件读写操作的方法,包括创建文件、写入文件和读取文件的过程。同时还介绍了如何判断文件是否成功打开和关闭文件的方法。通过本文的学习,读者可以了解如何在C++中进行文件读写操作。 ... [详细]
  • 本文讨论了一个数列求和问题,该数列按照一定规律生成。通过观察数列的规律,我们可以得出求解该问题的算法。具体算法为计算前n项i*f[i]的和,其中f[i]表示数列中有i个数字。根据参考的思路,我们可以将算法的时间复杂度控制在O(n),即计算到5e5即可满足1e9的要求。 ... [详细]
  • 李逍遥寻找仙药的迷阵之旅
    本文讲述了少年李逍遥为了救治婶婶的病情,前往仙灵岛寻找仙药的故事。他需要穿越一个由M×N个方格组成的迷阵,有些方格内有怪物,有些方格是安全的。李逍遥需要避开有怪物的方格,并经过最少的方格,找到仙药。在寻找的过程中,他还会遇到神秘人物。本文提供了一个迷阵样例及李逍遥找到仙药的路线。 ... [详细]
  • 本文由编程笔记#小编整理,主要介绍了关于数论相关的知识,包括数论的算法和百度百科的链接。文章还介绍了欧几里得算法、辗转相除法、gcd、lcm和扩展欧几里得算法的使用方法。此外,文章还提到了数论在求解不定方程、模线性方程和乘法逆元方面的应用。摘要长度:184字。 ... [详细]
  • 本文介绍了Codeforces Round #321 (Div. 2)比赛中的问题Kefa and Dishes,通过状压和spfa算法解决了这个问题。给定一个有向图,求在不超过m步的情况下,能获得的最大权值和。点不能重复走。文章详细介绍了问题的题意、解题思路和代码实现。 ... [详细]
  • 本文介绍了设计师伊振华受邀参与沈阳市智慧城市运行管理中心项目的整体设计,并以数字赋能和创新驱动高质量发展的理念,建设了集成、智慧、高效的一体化城市综合管理平台,促进了城市的数字化转型。该中心被称为当代城市的智能心脏,为沈阳市的智慧城市建设做出了重要贡献。 ... [详细]
  • 本文介绍了一种划分和计数油田地块的方法。根据给定的条件,通过遍历和DFS算法,将符合条件的地块标记为不符合条件的地块,并进行计数。同时,还介绍了如何判断点是否在给定范围内的方法。 ... [详细]
  • 本文介绍了解决二叉树层序创建问题的方法。通过使用队列结构体和二叉树结构体,实现了入队和出队操作,并提供了判断队列是否为空的函数。详细介绍了解决该问题的步骤和流程。 ... [详细]
  • 动态规划算法的基本步骤及最长递增子序列问题详解
    本文详细介绍了动态规划算法的基本步骤,包括划分阶段、选择状态、决策和状态转移方程,并以最长递增子序列问题为例进行了详细解析。动态规划算法的有效性依赖于问题本身所具有的最优子结构性质和子问题重叠性质。通过将子问题的解保存在一个表中,在以后尽可能多地利用这些子问题的解,从而提高算法的效率。 ... [详细]
  • This article discusses the efficiency of using char str[] and char *str and whether there is any reason to prefer one over the other. It explains the difference between the two and provides an example to illustrate their usage. ... [详细]
  • 先看官方文档TheJavaTutorialshavebeenwrittenforJDK8.Examplesandpracticesdescribedinthispagedontta ... [详细]
  • Android工程师面试准备及设计模式使用场景
    本文介绍了Android工程师面试准备的经验,包括面试流程和重点准备内容。同时,还介绍了建造者模式的使用场景,以及在Android开发中的具体应用。 ... [详细]
author-avatar
mobiledu2502886443
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有