SPOJ 1557 Can you answer these queries II（GSS2 线段树）

转载请注明出处，谢谢/article/2566293.html
by---cxlove

题目：求出区间的最大子段和，而且相同的数只能计算一次

http://www.spoj.pl/problems/GSS2/

嘛~~~这个题还是很不错的，不会，哭~~~~喵

在线的处理不了，只能离线处理

按询问区间的右端点排序，然后依次枚举a[i]，pos[a[i]]表示对于数字a[i]，上一个位置在哪

那么只需要将pos[a[i]]与i之间的部分全部加上a[i]

当查询的右区间为当前更新的位置的时候，查询便是历史区间最值。

但是实现起来并不简单~~~

维护4个值：old_cover 历史最大更新 old_max 历史最大值

now_cover 当前的更新值 now_max 当前区间的最大值

具体的看代码吧，各种纠结，尤其是Push_Down部分

#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
#include<algorithm>
#include<set>
#define inf (1ull<<63)-1
#define N 100005
#define maxn 100005
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define pb(a) push_back(a)
#define mem(a,b) memset(a,b,sizeof(a))
#define eps 1e-9
#define zero(a) fabs(a)<eps
#define LL long long
#define ULL unsigned long long
#define lson (step<<1)
#define rson (step<<1|1)
#define MOD 1000000007
#define mp(a,b) make_pair(a,b)
using namespace std;
struct Node{
	int left,right;
	int old_max,now_max,old_cover,now_cover;
}L[N*4];
struct Que{
	int left,right,id;
	bool operator<(const Que t)const{
		return right<t.right;
	}
}Q
;
int n,q,a
,pos[N*2+5],ans
;
void Bulid(int step,int l,int r){
	L[step].left=l;
	L[step].right=r;
	L[step].old_max=L[step].now_max=L[step].now_cover=L[step].old_cover=0;
	if(l==r) return;
	int m=(l+r)/2;
	Bulid(lson,l,m);
	Bulid(rson,m+1,r);
}
void Push_Down(int step){
	L[lson].old_cover=max(L[lson].old_cover,L[lson].now_cover+L[step].old_cover);
	L[rson].old_cover=max(L[rson].old_cover,L[rson].now_cover+L[step].old_cover);

	L[lson].old_max=max(L[lson].old_max,L[lson].now_max+L[step].old_cover);
	L[rson].old_max=max(L[rson].old_max,L[rson].now_max+L[step].old_cover);

	L[lson].now_cover+=L[step].now_cover;
	L[rson].now_cover+=L[step].now_cover;

	L[lson].now_max+=L[step].now_cover;
	L[rson].now_max+=L[step].now_cover;

	L[step].now_cover=L[step].old_cover=0;
}
void Push_Up(int step){
	L[step].now_max=max(L[lson].now_max,L[rson].now_max);
	L[step].old_max=max(L[lson].old_max,L[rson].old_max);
}
void Update(int step,int l,int r,int c){
	if(L[step].left==l&&L[step].right==r){
		L[step].old_cover=max(L[step].old_cover,L[step].now_cover+=c);
		L[step].old_max=max(L[step].old_max,L[step].now_max+=c);
		return;
	}
	Push_Down(step);
	int m=(L[step].left+L[step].right)/2;
	if(r<=m) Update(lson,l,r,c);
	else if(l>m) Update(rson,l,r,c);
	else{
		Update(lson,l,m,c);
		Update(rson,m+1,r,c);
	}
	Push_Up(step);
}
int Query(int step,int l,int r){
	if(L[step].left==l&&L[step].right==r)
		return L[step].old_max;
	Push_Down(step);
	int m=(L[step].left+L[step].right)/2;	
	if(r<=m) return Query(lson,l,r);
	else if(l>m) return Query(rson,l,r);
	else return max(Query(lson,l,m),Query(rson,m+1,r));
}
int main(){
	while(scanf("%d",&n)!=EOF){
		Bulid(1,1,n);
		for(int i=1;i<=n;i++)
			scanf("%d",&a[i]);
		scanf("%d",&q);
		for(int i=0;i<q;i++){
			scanf("%d%d",&Q[i].left,&Q[i].right);
			Q[i].id=i;
		}
		sort(Q,Q+q);
		mem(pos,0);
		for(int i=1,j=0;i<=n;i++){
			Update(1,pos[a[i]+N]+1,i,a[i]);	
			pos[a[i]+N]=i;
			while(j<q&&Q[j].right==i){
				ans[Q[j].id]=Query(1,Q[j].left,Q[j].right);
				j++;
			}
			if(j>=q) break;
		}
		for(int i=0;i<q;i++) printf("%d\n",ans[i]);
	}
	return 0;
}