poj2750 线段树+动态规划

问题描述:给定一个环形序列,进行在线操作,每次修改一个元素,输出环上的最大连续子列的和。

出题者的简单解题报告:把环从一个地方,切断拉成一条直线,用线段树记录当前区间的非空最大子列和当前区间的非空最小子列。

如果环上的数都是正整数,答案是:环上数的总和-根结点的非空最小子列;

否则,答案是: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;
}


你可能感兴趣的:(c,存储,Build)