问题描述:给定一个环形序列,进行在线操作,每次修改一个元素,输出环上的最大连续子列的和。
出题者的简单解题报告:把环从一个地方,切断拉成一条直线,用线段树记录当前区间的非空最大子列和当前区间的非空最小子列。
如果环上的数都是正整数,答案是:环上数的总和-根结点的非空最小子列;
否则,答案是:max{根结点的非空最大子列,环上数的总和-根结点的非空最小子列},每次问答的复杂度是O(logN)。
对于一个结点,我们要存储以下信息:
int maxSeg[MAXN<<2], minSeg[MAXN<<2]; //结点[a, b]的最大非空连续子列和最小非空连续子列
int lmax[MAXN<<2]; //结点[a, b]包含a的最大非空连续子列
int rmax[MAXN<<2]; //结点[a, b]包含b的最大非空连续子列
int lmin[MAXN<<2]; //结点[a, b]包含a的最小非空连续子列
int rmin[MAXN<<2]; //结点[a, b]包含b的最小非空连续子列
int sum[MAXN<<2]; //结点[a, b]的总和
结点[a, m],[m+1, b]分别为[a, b]的左右儿子
动态规划思想由结点左右儿子的信息得到父节点的信息,设父节点下标为rt
那么:状态转移方程如下:
maxSeg[rt] = max(max(maxSeg[rt<<1], maxSeg[rt<<1|1]), rmax[rt<<1]+lmax[rt<<1|1]); //求父节点的最大非空连续子列
minSeg[rt] = min(min(minSeg[rt<<1], minSeg[rt<<1|1]), rmin[rt<<1]+lmin[rt<<1|1]);//求父节点的最小非空连续子列
lmax[rt] = max(lmax[rt<<1], sum[rt<<1]+lmax[rt<<1|1]);//求父节点的包含a的最大非空连续子列
rmax[rt] = max(rmax[rt<<1|1], sum[rt<<1|1]+rmax[rt<<1]);//求父节点的包含b的最大非空连续子列
lmin[rt] = min(lmin[rt<<1], sum[rt<<1]+lmin[rt<<1|1]); //求父节点的包含a的最小非空连续子列
rmin[rt] = min(rmin[rt<<1|1], sum[rt<<1|1]+rmin[rt<<1]);//求父节点的包含b的最小非空连续子列
#include <iostream> #include <cstdio> #include <cstring> using namespace std; #define MAXN 100005 int maxSeg[MAXN<<2], minSeg[MAXN<<2]; int lmax[MAXN<<2], rmax[MAXN<<2]; int lmin[MAXN<<2], rmin[MAXN<<2]; int sum[MAXN<<2]; void PushUp(int rt) { int l = rt<<1; int r = l+1; sum[rt] = sum[l] + sum[r]; maxSeg[rt] = max(max(maxSeg[l], maxSeg[r]), rmax[l]+lmax[r]); minSeg[rt] = min(min(minSeg[l], minSeg[r]), rmin[l]+lmin[r]); lmax[rt] = max(lmax[l], sum[l]+lmax[r]); rmax[rt] = max(rmax[r], sum[r]+rmax[l]); lmin[rt] = min(lmin[l], sum[l]+lmin[r]); rmin[rt] = min(rmin[r], sum[r]+rmin[l]); } void Build(int l, int r, int rt) { if (l == r) { scanf("%d", &sum[rt]); minSeg[rt] = lmax[rt] = rmax[rt] = lmin[rt] = rmin[rt] = maxSeg[rt] = sum[rt]; return; } int m = (l+r)>>1; Build(l, m, rt<<1); Build(m+1, r, rt<<1|1); PushUp(rt); } void Update(int id, int c, int l, int r, int rt) { if (l == r) { sum[rt] = maxSeg[rt] = minSeg[rt] = c; lmax[rt] = rmax[rt] = lmin[rt] = rmin[rt] = c; return; } int m = (l+r)>>1; if (m >= id) Update(id, c, l, m, rt<<1); else Update(id, c, m+1, r, rt<<1|1); PushUp(rt); } int main() { int n, m, ans; scanf ("%d", &n); Build(1, n, 1); scanf("%d", &m); while (m--) { int a, b; scanf ("%d%d", &a, &b); Update(a, b, 1, n, 1); if (sum[1] == maxSeg[1]) //序列全为非负数的时候 ans = sum[1] - minSeg[1]; else ans = max(maxSeg[1], sum[1]-minSeg[1]); printf ("%d\n", ans); } return 0; }