博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
洛谷P4513 小白逛公园
阅读量:7198 次
发布时间:2019-06-29

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

区间最大子段和模板题。。

维护四个数组:prefix, suffix, sum, tree

假设当前访问节点为cur

  • prefix[cur]=max(prefix[lson],sum[lson]+preifx[rson])
  • suffix[cur]=max(suffix[rson],sum[rson]+suffix[lson])
  • sum[cur]=sum[lson]+sum[rson]
  • tree[cur]=max(tree[lson], tree[rson], suffix[lson]+suffix[rson])

在query的时候要注意合并,其实和pushup差不多,需要返回一个节点

// luogu-judger-enable-o2#include 
#define INF 0x3f3f3f3fusing namespace std;typedef long long ll;inline int lowbit(int x){ return x & (-x); }inline int read(){ int X = 0, w = 0; char ch = 0; while(!isdigit(ch)) { w |= ch == '-'; ch = getchar(); } while(isdigit(ch)) X = (X << 3) + (X << 1) + (ch ^ 48), ch = getchar(); return w ? -X : X;}inline int gcd(int a, int b){ return a % b ? gcd(b, a % b) : b; }inline int lcm(int a, int b){ return a / gcd(a, b) * b; }template
inline T max(T x, T y, T z){ return max(max(x, y), z); }template
inline T min(T x, T y, T z){ return min(min(x, y), z); }template
inline A fpow(A x, B p, C yql){ A ans = 1; for(; p; p >>= 1, x = 1LL * x * x % yql)if(p & 1)ans = 1LL * x * ans % yql; return ans;}const int N = 500005;struct Node{ int prefix, suffix, tree, sum; Node(int p, int s, int t, int u): prefix(p), suffix(s), tree(t), sum(u){}};int n, m, a[N];int prefix[N<<2], suffix[N<<2], sum[N<<2], tree[N<<2];void push_up(int treeIndex){ int lson = treeIndex << 1, rson = treeIndex << 1 | 1; prefix[treeIndex] = max(prefix[lson], sum[lson] + prefix[rson]); suffix[treeIndex] = max(suffix[rson], sum[rson] + suffix[lson]); sum[treeIndex] = sum[lson] + sum[rson]; tree[treeIndex] = max(tree[lson], tree[rson], suffix[lson] + prefix[rson]);}void buildTree(int treeIndex, int l, int r){ if(l == r){ tree[treeIndex] = prefix[treeIndex] = suffix[treeIndex] = sum[treeIndex] = a[l]; return; } int mid = (l + r) >> 1; buildTree(treeIndex << 1, l, mid); buildTree(treeIndex << 1 | 1, mid + 1, r); push_up(treeIndex);}void modify(int treeIndex, int l, int r, int k, int e){ if(l == r){ tree[treeIndex] = prefix[treeIndex] = suffix[treeIndex] = sum[treeIndex] = e; return; } int mid = (l + r) >> 1; if(k <= mid) modify(treeIndex << 1, l, mid, k, e); else modify(treeIndex << 1 | 1, mid + 1, r, k, e); push_up(treeIndex);}Node query(int treeIndex, int l, int r, int queryL, int queryR){ if(l == queryL && r == queryR){ return Node(prefix[treeIndex], suffix[treeIndex], tree[treeIndex], sum[treeIndex]); } int mid = (l + r) >> 1; if(queryL > mid) return query(treeIndex << 1 | 1, mid + 1, r, queryL, queryR); else if(queryR <= mid) return query(treeIndex << 1, l, mid, queryL, queryR); Node lr = query(treeIndex << 1, l, mid, queryL, mid); Node rr = query(treeIndex << 1 | 1, mid + 1, r, mid + 1, queryR); Node ret = Node(max(lr.prefix, lr.sum + rr.prefix), max(rr.suffix, rr.sum + lr.suffix), max(lr.tree, rr.tree, lr.suffix + rr.prefix), lr.sum + rr.sum); return ret;}int main(){ n = read(), m = read(); for(int i = 1; i <= n; i ++) a[i] = read(); buildTree(1, 1, n); while(m --){ int opt = read(), p = read(), q = read(); if(opt == 1){ if(p > q) swap(p, q); printf("%d\n", query(1, 1, n, p, q).tree); } else if(opt == 2){ modify(1, 1, n, p, q); } } return 0;}

转载于:https://www.cnblogs.com/onionQAQ/p/10507383.html

你可能感兴趣的文章
在RHEL6上配置基于ftp的YUM
查看>>
【行为型】- 状态模式
查看>>
AIX5.3 hacmp配置
查看>>
Linux关闭休眠和屏保模式
查看>>
我的友情链接
查看>>
OpenJudge 4146:数字方格 java穷举法
查看>>
搭建linuxFTP的步骤
查看>>
跟我一起学docker(八)--Dockerfile
查看>>
做SEO,选择大于努力
查看>>
python web开发-flask调试模式
查看>>
3分钟看懂linux磁盘划分
查看>>
Linux-LAMP 默认页,虚拟主机
查看>>
Node.js 特点
查看>>
Linux-日志管理
查看>>
MySQL数据库系统
查看>>
Android Studio 提示帮助文档 一直显示:fetching documentation
查看>>
新 Terraform 提供商: F5 Networks, Nutanix, 腾讯云, Helm
查看>>
SpringBoot框架简介及搭建
查看>>
拯救 Java Code Style 强迫症
查看>>
PDF文档怎样在线合并?
查看>>