SPOJ 2916 Can you answer these queries V（GSS5 线段树）

转载请注明出处，谢谢http://blog.csdn.net/acm_cxlove/article/details/7854526       by---cxlove

题目：和GSS3类似，只不过查询麻烦一点，要注意细节。区间的起点和终点都在指定区间内。

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

更新部分类似，查询部分也一样。

主要是把每一种查询分清楚。

我是分为三种情况

第一种：两个区间完全没有重叠：

则是x1    y1      x2   y2的形势，那么最终的区间肯定是包括(y1,x2)这段的，然后再对[x1,y1]取右最值，对[x2,y2]取左最值就可以了。

第二种：完全嵌套情况

即x1  x2  y2  y1这种情况，这里有两种情况，一个是，区间完全处于[x2,y2]内，这个肯定没有问题，另外一种是起点在[x1,x2),而终点在[x2,y2]。那么对于[x1,x2)取右最值，对[x2,y2]取左最值就行了。

第三种情况：部分重叠

即 x1 x2  y1  y2这种情况

分为4种：(1) 起点和终点都在重叠部分，也就是[x2,y1]求一次最值就行了

                (2)终点在[x2,y1],起点在[x1,x2),这样就是对于前者求右最值，对于后者取左最值

                (3)起点在[x2,y1],终点在(y1,y2]，对于前者求右最值，对于后者取左最值

                (4) 起点在[x1,x2),终点在(y1,y2]。这样中间的部分[x2,y1]肯定是包含的，再对前者求右最值对于后者求左最值就OK了。

有点乱，不过能够1A，哇哈哈

#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 10005
#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;
int n,a[N],q,sum[N]={0};
struct Node{
	int left,right,lx,rx,mx;
	int all(){
		return sum[right]-sum[left-1];
	}
}L[N*4];
void Push_Up(int step){
	L[step].lx=max(L[lson].lx,L[lson].all()+L[rson].lx);
	L[step].rx=max(L[rson].rx,L[rson].all()+L[lson].rx);
	L[step].mx=max(max(L[lson].mx,L[rson].mx),L[lson].rx+L[rson].lx);
}
void Bulid(int step,int l,int r){
	L[step].left=l;L[step].right=r;
	if(l==r){
		L[step].lx=L[step].rx=L[step].mx=a[l];
		return ;
	}
	int m=(l+r)/2;
	Bulid(lson,l,m);
	Bulid(rson,m+1,r);
	Push_Up(step);
}
int Q_L(int step,int l,int r){
	if(l>r) return 0;
	if(L[step].left==l&&L[step].right==r)
		return L[step].lx;
	int m=(L[step].left+L[step].right)/2;
	if(r<=m) return Q_L(lson,l,r);
	else if(l>m) return Q_L(rson,l,r);
	else return max(Q_L(lson,l,m),max(sum[m]-sum[l-1],sum[m]-sum[l-1]+Q_L(rson,m+1,r)));
}
int Q_R(int step,int l,int r){
	if(l>r) return 0;
	if(L[step].left==l&&L[step].right==r)
		return L[step].rx;
	int m=(L[step].left+L[step].right)/2;
	if(r<=m) return Q_R(lson,l,r);
	else if(l>m) return Q_R(rson,l,r);
	else return max(Q_R(rson,m+1,r),max(sum[r]-sum[m],sum[r]-sum[m]+Q_R(lson,l,m)));
}
int Query(int step,int l,int r){
	if(l>r) return 0;
	if(L[step].left==l&&L[step].right==r)
		return L[step].mx;
	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(max(Query(lson,l,m),Query(rson,m+1,r)),Q_R(lson,l,m)+Q_L(rson,m+1,r));
}
int main(){
	int t;
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		for(int i=1;i<=n;i++) {scanf("%d",&a[i]);sum[i]=sum[i-1]+a[i];}
		Bulid(1,1,n);
		scanf("%d",&q);
		while(q--){
			int x1,y1,x2,y2;
			scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
			if(x2>y1) printf("%d\n",sum[x2-1]-sum[y1]+Q_L(1,x2,y2)+Q_R(1,x1,y1));
			else if(y1>=y2) printf("%d\n",max(Query(1,x2,y2),Q_L(1,x2,y2)+Q_R(1,x1,x2-1)));
			else
				printf("%d\n",max(Query(1,x2,y1),max(Q_R(1,x1,x2-1)+Q_L(1,x2,y1),max(Q_R(1,x2,y1)+Q_L(1,y1+1,y2),Q_R(1,x1,x2-1)+Q_L(1,y1+1,y2)+sum[y1]-sum[x2-1]))));
		}
	}
	return 0;
}