题目链接:https://www.spoj.pl/problems/GSS2/
题目大意:求子序列的最大子序列和,但是相同值只加一次.
思路:
可以看成是Necklace的升级版,不同的地方在于求的是区间最大子区间,而不是区间和,所以线段树维护的值当然也不太一样,其余都一样了.
那么不一样在神马地方捏?
首先结点中保存的值不再是某位置的值而是从某位置起加到当前位置(sum[])以及从某位置加起到当前位置的最大和(psum[]).为什么呢?
易知一个询问[x,y]的答案为
ret[x,y]=max{
max(A[x],A[x]+A[x+1],...,A[x]+A[x+1]+...+A[y]) ,
max (A[x+1],A[x+1]+A[x+2],...,A[x+1]+A[x+2]+...+A[y])
,...,
max(A[i],A[i]+A[i+1],...,A[i]+A[i+1]+...+A[y])
,..,
max(A[y]) }
仔细看看,肯定能看出来的.
如果这个原因看出来了之后,那么就知道了这是一道区间更新的线段树,所以cov[]标记也是必要的了,当然psum也需要一个标记pcov[],cov[]标记是显然的,那pcov[]呢?
显然是为了对pmax[]向下更新用的.
比如
8
-1 -2 -3 4 -4 -5 -6 -7
1
3 5
这么一组数据,那么如果pmax[]没有向下更新,那么答案是0而不是4.
最后唯一要注意的是向下更新时要注意pcov[]的连续性,这样才能保证psum[]的连续性.
代码:
#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <ctype.h>
#include <math.h>
#include <time.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;
//#define ull unsigned __int64
//#define ll __int64
#define ull unsigned long long
#define ll long long
#define son1 New(p.xl,xm,p.yl,ym),(rt<<2)-2
#define son2 New(p.xl,xm,min(ym+1,p.yr),p.yr),(rt<<2)-1
#define son3 New(min(xm+1,p.xr),p.xr,p.yl,ym),rt<<2
#define son4 New(min(xm+1,p.xr),p.xr,min(ym+1,p.yr),p.yr),rt<<2|1
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define middle (l+r)>>1
#define MOD 1000000007
#define esp (1e-8)
const int INF=0x3F3F3F3F;
const double DINF=10000.00;
//const double pi=acos(-1.0);
const int N=100010;
int n,m,tot;
struct node{
int x,y,id;
void write(int _id){scanf("%d%d",&x,&y),id=_id;}
bool operator < (const node &tmp)const{
return y < tmp.y;
}
}q[N];
ll A[N],ret[N];
ll pmax[N<<2],nmax[N<<2],pcov[N<<2],ncov[N<<2];
int pos[N<<1];
ll max(ll x,ll y){return x>y? x:y;}
void pushUp(int rt){
int ls=rt<<1,rs=ls|1;
nmax[rt]=max(nmax[ls],nmax[rs]);
pmax[rt]=max(pmax[ls],pmax[rs]);
}
void pushDown(int rt){
int ls=rt<<1,rs=ls|1;
pcov[ls]=max(pcov[ls],ncov[ls]+pcov[rt]);
pmax[ls]=max(pmax[ls],nmax[ls]+pcov[rt]);
ncov[ls]+=ncov[rt],nmax[ls]+=ncov[rt];
pcov[rs]=max(pcov[rs],ncov[rs]+pcov[rt]);
pmax[rs]=max(pmax[rs],nmax[rs]+pcov[rt]);
ncov[rs]+=ncov[rt],nmax[rs]+=ncov[rt];
ncov[rt]=pcov[rt]=0;
}
void update(int l,int r,int rt,int L,int R,ll c){
if(L<=l && r<=R){
pcov[rt]=max(pcov[rt],ncov[rt]+=c);
pmax[rt]=max(pmax[rt],nmax[rt]+c);
nmax[rt]+=c;
return;
}
pushDown(rt);
int mid=middle;
if(L<=mid) update(lson,L,R,c);
if(mid<R) update(rson,L,R,c);
pushUp(rt);
}
ll query(int l,int r,int rt,int L,int R){
if(L<=l && r<=R) return pmax[rt];
pushDown(rt);
int mid=middle;
ll ret=0;
if(L<=mid) ret=max(ret,query(lson,L,R));
if(mid<R) ret=max(ret,query(rson,L,R));
return ret;
}
void init(){
memset(pmax,0,sizeof(pmax));
memset(nmax,0,sizeof(nmax));
memset(pcov,0,sizeof(pcov));
memset(ncov,0,sizeof(ncov));
memset(pos,0,sizeof(pos));
int i;
for(i=1;i<=n;i++) scanf("%lld",&A[i]);
scanf("%d",&m);
for(i=1;i<=m;i++) q[i].write(i);
sort(q+1,q+m+1);
}
void sof(){
int i,j;
for(i=1,j=1;i<=n;i++){
update(1,n,1,pos[A[i]+N]+1,i,A[i]);
pos[A[i]+N]=i;
while(q[j].y==i){
ret[q[j].id]=query(1,n,1,q[j].x,q[j].y);
if(++j>m) goto end;
}
}end:
for(i=1;i<=m;i++) printf("%lld\n",ret[i]);
}
int main(){
//freopen("1.in","r",stdin);
//freopen("me.out","w",stdout);
//int T,cas;scanf("%d",&T);for(cas=1;cas<=T;cas++){
scanf("%d",&n);init(),sof();
return 0;
}