转载请注明出处，谢谢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;
}