题意

给定一棵树,设计数据结构支持以下操作

1 u v d  表示将路径 (u,v) 加d(d>=0)

2 u v	表示询问路径 (u,v) 上点权绝对值的和

分析

绝对值之和不好处理,那么我们开两棵线段树,一个存正数,一个存负数.然后对于两棵线段树,都要维护子树sz(有效节点数),sum(有效节点权值之和),lz(加法懒标记).特别的,因为负数可能会加到正数,那么修改一个区间的时候,询问一下这个区间最大的负数加上d有没有变成正数,如果有就暴力从负数线段树中删去这个节点,加入正数线段树中.又题目中满足每一次加的数都是非负的,也就是说一个数最多从负变成正一次,那么时间复杂度均摊是O(nlogn)O(nlogn)O(nlogn)的(每次删除/插入是lognlognlogn的).注意细节就行了…

当然考虑到一个位置要么在正树要么在负树,那么其实可以只写一棵线段树,如果当前区间最大负数+d仍然是负数就打上懒标记返回,否则暴力下传.线段树时间复杂度仍是O(nlogn)O(nlogn)O(nlogn).加上树链剖分就是O(nlog2n)O(nlog^2n)O(nlog2n)啦.

CODE

看着我7000B的代码陷入沉思…为毛打个线段树板题都要220行+

(其实两棵树大多数函数可以直接复制)

#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
char cb[1<<15],*cs=cb,*ct=cb;
#define getc() (cs==ct&&(ct=(cs=cb)+fread(cb,1,1<<15,stdin),cs==ct)?0:*cs++)
template<class T>inline void read(T &res) {
char ch; int flg = 1; for(;!isdigit(ch=getc());)if(ch=='-')flg=-flg;
for(res=ch-'0';isdigit(ch=getc());res=res*10+ch-'0'); res*=flg;
}
const int MAXN = 100005;
const LL INF = 1e16;
int n, m, cnt, tmr, w[MAXN], a[MAXN], fir[MAXN], fa[MAXN], dfn[MAXN], top[MAXN], sz[MAXN], son[MAXN], dep[MAXN];
struct edge { int to, nxt; }e[MAXN<<1];
inline void Add(int u, int v) { e[cnt] = (edge){ v, fir[u] }, fir[u] = cnt++; }
void dfs(int x, int ff) { //树剖
dep[x] = dep[fa[x]=ff] + (sz[x]=1);
for(int v, i = fir[x]; ~i; i = e[i].nxt)
if((v=e[i].to) != ff) {
dfs(v, x), sz[x] += sz[v];
if(sz[v] > sz[son[x]]) son[x] = v;
}
}
void dfs2(int x, int tp) {
top[x] = tp; w[dfn[x] = ++tmr] = a[x];
if(son[x]) dfs2(son[x], tp);
for(int v, i = fir[x]; ~i; i = e[i].nxt)
if((v=e[i].to) != fa[x] && v != son[x])
dfs2(v, v);
} struct node { //写个结构体存最大值以及坐标
LL val; int id;
node(){}
node(LL v, int i):val(v), id(i){}
inline node operator +(const node &o)const {
if(val < o.val) return o;
else return *this;
}
}; namespace PositiveTree { //正线段树
LL sum[MAXN<<2], lz[MAXN<<2]; int sz[MAXN<<2];
inline void update(int i) {
sum[i] = sum[i<<1] + sum[i<<1|1];
sz[i] = sz[i<<1] + sz[i<<1|1];
}
inline void pushdown(int i) {
if(lz[i]) {
lz[i<<1] += lz[i]; sum[i<<1] += 1ll * sz[i<<1] * lz[i];
lz[i<<1|1] += lz[i]; sum[i<<1|1] += 1ll * sz[i<<1|1] * lz[i];
lz[i] = 0;
}
}
void build(int i, int l, int r) {
if(l == r) {
sum[i] = w[l] >= 0 ? w[l] : 0;
sz[i] = w[l] >= 0 ? 1 : 0;
return;
}
int mid = (l + r) >> 1;
build(i<<1, l, mid);
build(i<<1|1, mid+1, r);
update(i);
}
void insert(int i, int l, int r, int x, int y, int val) {
if(x == l && y == r) {
sum[i] += 1ll * sz[i] * val;
lz[i] += val;
return;
}
int mid = (l + r) >> 1;
pushdown(i);
if(y <= mid) insert(i<<1, l, mid, x, y, val);
else if(x > mid) insert(i<<1|1, mid+1, r, x, y, val);
else insert(i<<1, l, mid, x, mid, val), insert(i<<1|1, mid+1, r, mid+1, y, val);
update(i);
}
void modify(int i, int l, int r, int x, int val) {
if(l == r) {
sum[i] = val, sz[i] = 1; return;
}
int mid = (l + r) >> 1;
pushdown(i);
if(x <= mid) modify(i<<1, l, mid, x, val);
else modify(i<<1|1, mid+1, r, x, val);
update(i);
}
LL querysum(int i, int l, int r, int x, int y) {
if(x == l && y == r) return sum[i];
int mid = (l + r) >> 1;
pushdown(i);
LL res;
if(y <= mid) res = querysum(i<<1, l, mid, x, y);
else if(x > mid) res = querysum(i<<1|1, mid+1, r, x, y);
else res = querysum(i<<1, l, mid, x, mid) + querysum(i<<1|1, mid+1, r, mid+1, y);
update(i);
return res;
}
} namespace NegativeTree { //负线段树
node mx[MAXN<<2];
LL sum[MAXN<<2], lz[MAXN<<2]; int sz[MAXN<<2];
inline void update(int i) {
sum[i] = sum[i<<1] + sum[i<<1|1];
mx[i] = mx[i<<1] + mx[i<<1|1];
sz[i] = sz[i<<1] + sz[i<<1|1];
}
inline void pushdown(int i) {
if(lz[i]) {
lz[i<<1] += lz[i]; sum[i<<1] += 1ll * sz[i<<1] * lz[i]; mx[i<<1].val += lz[i];
lz[i<<1|1] += lz[i]; sum[i<<1|1] += 1ll * sz[i<<1|1] * lz[i]; mx[i<<1|1].val += lz[i];
lz[i] = 0;
}
}
void build(int i, int l, int r) {
if(l == r) {
mx[i] = node(w[l] < 0 ? w[l] : -INF, l);
sz[i] = w[l] < 0 ? 1 : 0;
sum[i] = w[l] < 0 ? w[l] : 0; return;
}
int mid = (l + r) >> 1;
build(i<<1, l, mid);
build(i<<1|1, mid+1, r);
update(i);
}
void insert(int i, int l, int r, int x, int y, int val) {
if(x == l && y == r) {
sum[i] += 1ll * sz[i] * val;
lz[i] += val;
mx[i].val += val;
return;
}
int mid = (l + r) >> 1;
pushdown(i);
if(y <= mid) insert(i<<1, l, mid, x, y, val);
else if(x > mid) insert(i<<1|1, mid+1, r, x, y, val);
else insert(i<<1, l, mid, x, mid, val), insert(i<<1|1, mid+1, r, mid+1, y, val);
update(i);
}
void modify(int i, int l, int r, int x) {
if(l == r) {
sum[i] = sz[i] = 0, mx[i].val = -INF; return;
}
int mid = (l + r) >> 1;
pushdown(i);
if(x <= mid) modify(i<<1, l, mid, x);
else modify(i<<1|1, mid+1, r, x);
update(i);
}
node querymx(int i, int l, int r, int x, int y) {
if(x == l && y == r) return mx[i];
int mid = (l + r) >> 1;
pushdown(i);
node res;
if(y <= mid) res = querymx(i<<1, l, mid, x, y);
else if(x > mid) res = querymx(i<<1|1, mid+1, r, x, y);
else res = querymx(i<<1, l, mid, x, mid) + querymx(i<<1|1, mid+1, r, mid+1, y);
update(i);
return res;
}
LL querysum(int i, int l, int r, int x, int y) {
if(x == l && y == r) return sum[i];
int mid = (l + r) >> 1;
pushdown(i);
LL res;
if(y <= mid) res = querysum(i<<1, l, mid, x, y);
else if(x > mid) res = querysum(i<<1|1, mid+1, r, x, y);
else res = querysum(i<<1, l, mid, x, mid) + querysum(i<<1|1, mid+1, r, mid+1, y);
update(i);
return res;
}
void ADD(int l, int r, int val) {
if(r < l) return;
node res = querymx(1, 1, n, l, r);
if(res.val + val >= 0) {
modify(1, 1, n, res.id); //暴力在负线段树中删
PositiveTree::modify(1, 1, n, res.id, res.val + val); //暴力在正线段树中插入
ADD(l, res.id-1, val);
ADD(res.id+1, r, val);
}
else insert(1, 1, n, l, r, val);
}
} inline void Plus(int u, int v, int d) {
while(top[u]!=top[v]) {
if(dep[top[u]] < dep[top[v]]) swap(u, v);
PositiveTree::insert(1, 1, n, dfn[top[u]], dfn[u], d);
NegativeTree::ADD(dfn[top[u]], dfn[u], d);
u = fa[top[u]];
}
if(dep[u] < dep[v]) swap(u, v);
PositiveTree::insert(1, 1, n, dfn[v], dfn[u], d);
NegativeTree::ADD(dfn[v], dfn[u], d);
} inline LL Query(int u, int v) {
LL res = 0;
while(top[u]!=top[v]) {
if(dep[top[u]] < dep[top[v]]) swap(u, v);
res += - NegativeTree::querysum(1, 1, n, dfn[top[u]], dfn[u]) + PositiveTree::querysum(1, 1, n, dfn[top[u]], dfn[u]);
u = fa[top[u]];
}
if(dep[u] < dep[v]) swap(u, v);
return res - NegativeTree::querysum(1, 1, n, dfn[v], dfn[u]) + PositiveTree::querysum(1, 1, n, dfn[v], dfn[u]);
} int main() {
memset(fir, -1, sizeof fir);
read(n), read(m);
for(int i = 1; i <= n; ++i) read(a[i]);
for(int x, y, i = 1; i < n; ++i)
read(x), read(y), Add(x, y), Add(y, x);
dfs(1, 0), dfs2(1, 1);
NegativeTree::build(1, 1, n);
PositiveTree::build(1, 1, n);
int op, u, v, d;
while(m--) {
read(op), read(u), read(v);
if(op == 1) read(d), Plus(u, v, d);
else printf("%lld\n", Query(u, v));
}
}

BZOJ 4127: Abs (树链剖分 线段树求区间绝对值之和 带区间加法)的更多相关文章

  1. BZOJ.4034 [HAOI2015]树上操作 ( 点权树链剖分 线段树 )

    BZOJ.4034 [HAOI2015]树上操作 ( 点权树链剖分 线段树 ) 题意分析 有一棵点数为 N 的树,以点 1 为根,且树点有边权.然后有 M 个 操作,分为三种: 操作 1 :把某个节点 ...

  2. BZOJ.1036 [ZJOI2008]树的统计Count ( 点权树链剖分 线段树维护和与最值)

    BZOJ.1036 [ZJOI2008]树的统计Count (树链剖分 线段树维护和与最值) 题意分析 (题目图片来自于 这里) 第一道树链剖分的题目,谈一下自己的理解. 树链剖分能解决的问题是,题目 ...

  3. BZOJ 3672[NOI2014]购票(树链剖分+线段树维护凸包+斜率优化) + BZOJ 2402 陶陶的难题II (树链剖分+线段树维护凸包+分数规划+斜率优化)

    前言 刚开始看着两道题感觉头皮发麻,后来看看题解,发现挺好理解,只是代码有点长. BZOJ 3672[NOI2014]购票 中文题面,题意略: BZOJ 3672[NOI2014]购票 设f(i)f( ...

  4. bzoj 4196 [Noi2015]软件包管理器 (树链剖分+线段树)

    4196: [Noi2015]软件包管理器 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 2852  Solved: 1668[Submit][Sta ...

  5. 【bzoj4127】Abs 树链剖分+线段树

    题目描述 给定一棵树,设计数据结构支持以下操作 1 u v d 表示将路径 (u,v) 加d 2 u v 表示询问路径 (u,v) 上点权绝对值的和 输入 第一行两个整数n和m,表示结点个数和操作数 ...

  6. bzoj 2157: 旅游【树链剖分+线段树】

    裸的树链剖分+线段树 但是要注意一个地方--我WA了好几次才发现取完相反数之后max值和min值是要交换的-- #include<iostream> #include<cstdio& ...

  7. BZOJ 3589 动态树 (树链剖分+线段树)

    前言 众所周知,90%90\%90%的题目与解法毫无关系. 题意 有一棵有根树,两种操作.一种是子树内每一个点的权值加上一个同一个数,另一种是查询多条路径的并的点权之和. 分析 很容易看出是树链剖分+ ...

  8. BZOJ4127Abs——树链剖分+线段树

    题目描述 给定一棵树,设计数据结构支持以下操作 1 u v d 表示将路径 (u,v) 加d 2 u v 表示询问路径 (u,v) 上点权绝对值的和 输入 第一行两个整数n和m,表示结点个数和操作数 ...

  9. 【BZOJ-2325】道馆之战 树链剖分 + 线段树

    2325: [ZJOI2011]道馆之战 Time Limit: 40 Sec  Memory Limit: 256 MBSubmit: 1153  Solved: 421[Submit][Statu ...

  10. 【BZOJ2243】[SDOI2011]染色 树链剖分+线段树

    [BZOJ2243][SDOI2011]染色 Description 给定一棵有n个节点的无根树和m个操作,操作有2类: 1.将节点a到节点b路径上所有点都染成颜色c: 2.询问节点a到节点b路径上的 ...

随机推荐

  1. springboot整合filter

    新建一个项目 新建Firstfilter类 Firstfliter.java package com.example.filter; import java.io.IOException; impor ...

  2. 小菜鸟之HTML第二课

    JavaScript 运行在浏览器上的一种基于对象和事件的驱动的脚本语言 基于对象(windows – document location histroy 便于调用对象属性和方法 事件驱动 键盘和鼠标 ...

  3. TCP的组包、半包、粘包与分包

    一.概念 1)组包.简单的说就是tcp协议把过大的数据包分成了几个小的包传输,接收方要把同一组的数据包重新组合成一个完整的数据包. 2)半包.指接受方没有接受到一个完整的包,只接受了部分,这种情况主要 ...

  4. 题解 P3627 【[APIO2009]抢掠计划】

    咕了四个小时整整一晚上 P3627 [APIO2009] 抢掠计划(https://www.luogu.org/problemnew/show/P3627) 不难看出答案即为该有向图的最长链长度(允许 ...

  5. 使用JWT的ASP.NET CORE令牌身份验证和授权(无Cookie)——第1部分

    原文:使用JWT的ASP.NET CORE令牌身份验证和授权(无Cookie)--第1部分 原文链接:https://www.codeproject.com/Articles/5160941/ASP- ...

  6. 用<![CDATA[]]>将xml转义为 纯文本

    被<![CDATA[]]>这个标记所包含的内容将表示为纯文本,比如<![CDATA[<]]>表示文本内容"<". 此标记用于xml文档中,我们先 ...

  7. sqlserver关于发布订阅replication_subscription的总结

    (转载)sqlserver关于发布订阅replication_subscription的总结 来自 “ ITPUB博客 ” ,原文地址:http://blog.itpub.net/30126024/v ...

  8. 浅谈WEB中的高并发

    转载:https://www.cnblogs.com/guan-520/p/9575848.html 何谓高并发 高并发指的是:在同时或极短时间内,有大量的请求到达服务端,每个请求都需要服务端耗费资源 ...

  9. linux命令详解——top

    简介 TOP是一个动态显示过程,即可以通过用户按键来不断刷新当前状态.如果在前台执行该命令,它将独占前台,直到用户终止该程序为止.比较准确的说,top命令提供了实时的对系统处理器的状态监视.它将显示系 ...

  10. go语言在Windows系统下编译成linux系统可执行文件

    package main import ( "fmt" "os" "os/exec" ) //filepath: 要编译的文件的路径 fun ...