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)));elseprintf("%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;}