SPOJ - GSS2

    发现一道数据结构,写好了,WA,改好了,WA,改好了,WA...
    离线再用线段树维护的思路挺好想的,但是这个线段树求得MAX要包含历史所有版本中的,加上有lazy标记比较难维护。
    最开始的思路是一个M1[]表示当前最大值,delta[]表示lazy标记,这部分正常做应该没问题。
    然后再来一个M2[]表示M1[]的历史最大值。感觉没了...其实这样有的节点会出错...
    考虑一个lazy标记了两次的节点,它向下传递信息的时候根本没有考虑第一个标记的情况...被第二次直接覆盖了...
    so...
    解决方法也不难,对delta[]记一个历史最大值Md[x],向下传递的时候一起传就好了。
    更新M2[x]=max(M2[x],M1[x]-delta[x]+Md[x]);
    这样就记录所有标记的情况,不会产生标记覆盖了。
#include <map>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define MAX(x,y) (x>y?x:y)
typedef long long LL;
inline LL read(){
    LL x=0,f=1; char ch=getchar();
    while (ch<'0'||ch>'9') { if (ch=='-') f=-1; ch=getchar(); }
    while (ch>='0'&&ch<='9') { x=x*10+ch-'0'; ch=getchar(); }
    return x*f;
}

const LL MAXN=100005;
map<LL,LL> Q;
LL n,m,a[MAXN],c[MAXN],M1[MAXN<<2],M2[MAXN<<2],Md[MAXN<<2],delta[MAXN<<2];
struct Data{
    LL l,r,o;
    bool operator <(const Data &C)const{ return r<C.r||r==C.r&&l<C.l; }
}B[MAXN];
void Treeadd(LL x,LL l,LL r,LL ll,LL rr,LL k){
    if (l==ll&&r==rr) {
        delta[x]+=k; Md[x]=MAX(Md[x],delta[x]);
        M1[x]+=k; M2[x]=MAX(M2[x],M1[x]-delta[x]+Md[x]);
        return;
    }
    if (delta[x]!=0||Md[x]!=0) {
        delta[x<<1]+=delta[x]; delta[x<<1|1]+=delta[x];
        Md[x<<1]=MAX(Md[x<<1],delta[x<<1]); Md[x<<1|1]=MAX(Md[x<<1|1],delta[x<<1|1]);
        Md[x<<1]=MAX(Md[x<<1],Md[x]+delta[x<<1]-delta[x]);
        Md[x<<1|1]=MAX(Md[x<<1|1],Md[x]+delta[x<<1|1]-delta[x]);
        M1[x<<1]+=delta[x]; M1[x<<1|1]+=delta[x]; delta[x]=Md[x]=0;
        M2[x<<1]=MAX(M2[x<<1],M1[x<<1]-delta[x<<1]+Md[x<<1]);
        M2[x<<1|1]=MAX(M2[x<<1|1],M1[x<<1|1]-delta[x<<1|1]+Md[x<<1|1]);
    }
    LL mid=(l+r)>>1;
    if (ll>mid) Treeadd(x<<1|1,mid+1,r,ll,rr,k);
    else if (mid>=rr) Treeadd(x<<1,l,mid,ll,rr,k);
    else {
        Treeadd(x<<1,l,mid,ll,mid,k);
        Treeadd(x<<1|1,mid+1,r,mid+1,rr,k);
    }
    M1[x]=MAX(M1[x<<1],M1[x<<1|1]);
    M2[x]=MAX(M2[x],M1[x]);
}
LL Treeask(LL x,LL l,LL r,LL ll,LL rr){
    //cout<<x<<' '<<l<<' '<<r<<' '<<ll<<' '<<rr<<endl;
    //cout<<x<<' '<<M1[x]<<' '<<M2[x]<<endl;
    if (l==ll&&r==rr) return M2[x];
    if (delta[x]!=0||Md[x]!=0) {
        delta[x<<1]+=delta[x]; delta[x<<1|1]+=delta[x];
        Md[x<<1]=MAX(Md[x<<1],delta[x<<1]); Md[x<<1|1]=MAX(Md[x<<1|1],delta[x<<1|1]);
        Md[x<<1]=MAX(Md[x<<1],Md[x]+delta[x<<1]-delta[x]);
        Md[x<<1|1]=MAX(Md[x<<1|1],Md[x]+delta[x<<1|1]-delta[x]);
        M1[x<<1]+=delta[x]; M1[x<<1|1]+=delta[x]; delta[x]=Md[x]=0;
        M2[x<<1]=MAX(M2[x<<1],M1[x<<1]-delta[x<<1]+Md[x<<1]);
        M2[x<<1|1]=MAX(M2[x<<1|1],M1[x<<1|1]-delta[x<<1|1]+Md[x<<1|1]);
    }
    LL mid=(l+r)>>1;
    if (ll>mid) return Treeask(x<<1|1,mid+1,r,ll,rr);
    else if (mid>=rr) return Treeask(x<<1,l,mid,ll,rr);
    else return MAX(Treeask(x<<1,l,mid,ll,mid),Treeask(x<<1|1,mid+1,r,mid+1,rr));
}
int main(){
    LL i,j;
    n=read();
    for (i=1;i<=n;++i) a[i]=read();
    m=read();
    for (i=1;i<=m;++i) { B[i].l=read(); B[i].r=read(); B[i].o=i; }
    sort(B+1,B+1+m);
    //for (i=1;i<=m;++i) cout<<B[i].l<<' '<<B[i].r<<' '<<B[i].o<<endl;
    Q.clear();
    memset(M1,0,sizeof(M1));
    memset(M2,0,sizeof(M2));
    memset(delta,0,sizeof(delta));
    for (i=j=1;i<=n;++i) {
        Treeadd(1,1,n,Q[a[i]]+1,i,a[i]); Q[a[i]]=i;
        for (;j<=m&&i==B[j].r;++j) c[B[j].o]=Treeask(1,1,n,B[j].l,B[j].r);
        //cout<<Treeask(1,1,n,2,8)<<"lxsb"<<endl;
    }
    for (i=1;i<=m;++i) cout<<c[i]<<endl;

    return 0;
}

/*
10 
-1 -2 -3 4 4 5 6 -6 -6 100 
7 
1 10 
1 4 
1 5 
2 6 
4 9 
3 9 
2 8

*/
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值