[SPOJ1043,1557,1716,2713,2916,4487,6779,19543]GSS八题系列
[SPOJ1043]GSS1
求区间最大字段和，不带修改，线段树维护总和sum，最大子段和gss，最大左子段和lgss，最大右子段和rgss，即可。
更新答案：
sum=ls->sum+rs->sum;
lgss=max(ls->lgss,ls->sum+rs->lgss);
rgss=max(rs->rgss,rs->sum+ls->rgss);
gss=max(max(ls->gss,rs->gss),ls->rgss+rs->lgss);代码：
#include[spoj1557]GSS2
求最大子段和，重复的数字只能算一次。
这种重复只算一次的题目的思想就是离线给询问排序，预处理a[i]上一次出现的位置pre[i]，把一个区间压缩到一个点，转化成区间加操作。
这题就是先把询问按右端点排序。对序列从左到右操作，操作到第i个点时，保证线段树前i个点(x)的值为去重后的∑(a[j]){x<=j<=i}，并解决右端点为i的询问，其实我们只要给[pre[i]+1,i]区间加上i就可以维护。那么这题就转化为线段树区间加和求历史最值。
维护max，oldmax，add，oldadd分别表示当前最大值，历史最大值，当前加法标记，历史最大加法标记。
下发标记：
if(add!=0||oadd!=0)
{
chkmax(ls->oadd,ls->add+oadd);
chkmax(ls->omx,ls->mx+oadd);
ls->mx+=add; ls->add+=add;
chkmax(rs->oadd,rs->add+oadd);
chkmax(rs->omx,rs->mx+oadd);
rs->mx+=add; rs->add+=add;
add=oadd=0;
}至于这一句话chkmax(ls->omx,ls->mx+oadd)为什么是用当前最大值+历史最大标记来更新历史最大值，因为历史最大值的意义表示加法的操作做了某一个前缀（因为后面的值<0而变小于是舍弃）的结果，历史最大标记的意义也是表示从上一次标记清零开始最大的加法操作前缀，如果要更新历史最大值，就要求加法操作一定要做完至上一次标记清零（也就是还未更新当前最大值），在加上这个最大前缀，才能保证这样做后还是一个前缀。
代码：
#includey;}
bool cmp2(node a,node b){return a.idint l,r;
ll omx,mx,add,oadd;
tree *ls,*rs;
tree()
{
ls=rs=NULL;
l=r=omx=mx=add=oadd=0;
}
void cal(ll c)
{
add+=c;
mx+=c;
chkmax(omx,mx);
chkmax(oadd,add);
}
void pushdown()
{
if(add!=0||oadd!=0)
{
chkmax(ls->oadd,ls->add+oadd);
chkmax(ls->omx,ls->mx+oadd);
ls->mx+=add; ls->add+=add;
chkmax(rs->oadd,rs->add+oadd);
chkmax(rs->omx,rs->mx+oadd);
rs->mx+=add; rs->add+=add;
add=oadd=0;
}
}
void update()
{
mx=max(ls->mx,rs->mx);
omx=max(ls->omx,rs->omx);
}
void build(int lx,int rx)
{
l=lx;r=rx;
if(l==r) {omx=mx=0; return;}
int mid=(l+r)>>1;
(ls=new tree)->build(lx,mid);
(rs=new tree)->build(mid+1,rx);
update();
}
void modify(int lx,int rx,ll c)
{
if(lx==l&&rx==r) {cal(c); return;}
pushdown();
int mid=(l+r)>>1;
if(rx<=mid) ls->modify(lx,rx,c);
else if(lx>mid) rs->modify(lx,rx,c);
else {ls->modify(lx,mid,c);rs->modify(mid+1,rx,c);}
update();
}
ll query(int lx,int rx)
{
if(lx==l&&rx==r) return omx;
pushdown();
int mid=(l+r)>>1;
if(rx<=mid) return ls->query(lx,rx);
else if(lx>mid) return rs->query(lx,rx);
else return max(ls->query(lx,mid),rs->query(mid+1,rx));
update();
}
}*xtr;
int main()
{
n=read();
for(int i=1;i<=n;i++)
a[i]=read();
(xtr=new tree)->build(1,n);
m=read();
for(int i=1;i<=m;i++)
{
q[i].x=read(),q[i].y=read();q[i].id=i;
}
memset(visit,0,sizeof(visit));
memset(pre,0,sizeof(pre));
for(int i=1;i<=n;i++)
{
pre[i]=visit[a[i]+maxn];
visit[a[i]+maxn]=i;
}
sort(q+1,q+m+1,cmp1);
int pnt=1;
for(int i=1;i<=n;i++)
{
xtr->modify(pre[i]+1,i,a[i]);
while(pnt<=m&&q[pnt].y==i) {q[pnt].ans=xtr->query(q[pnt].x,i);pnt++;}
}
sort(q+1,q+m+1,cmp2);
for(int i=1;i<=m;i++)
printf("%lld\n",q[i].ans);
return 0;
}
[SPOJ1716]GSS3
带单点修改的最大子段和，类似GSS1。
代码：
#include[SPOJ2713]GSS4
区间开方，势能分析线段树，同BZOJ3211，代码都是搬过来的。
代码：
type
tree=^treenode;
treenode=record
l,r:longint;
max,min:int64;
sum:int64;
ls,rs:tree;
end;
var
n,i,q,opt,l,r,t,ca:longint;
delta:array[0..100100]of int64;
xtr:tree;
function max(x,y:int64):int64;
begin
if xthen exit(y)
else exit(x);
end;
function min(x,y:int64):int64;
begin
if x>y then exit(y)
else exit(x);
end;
procedure updata(x:tree);
begin
x^.max:=max(x^.ls^.max,x^.rs^.max);
x^.min:=min(x^.ls^.min,x^.rs^.min);
x^.sum:=x^.ls^.sum+x^.rs^.sum;
end;
procedure build(x:tree;l,r:longint);
var
mid:longint;
begin
x^.l:=l;
x^.r:=r;
if l=r then
begin
x^.max:=delta[l];
x^.min:=delta[l];
x^.sum:=delta[l];
x^.ls:=nil;
x^.rs:=nil;
exit;
end;
mid:=(l+r)div 2;
new(x^.ls);new(x^.rs);
build(x^.ls,l,mid);
build(x^.rs,mid+1,r);
updata(x);
end;
function ask(x:tree;l,r:longint):int64;
var
mid:longint;
begin
if (x^.l=l)and(x^.r=r) then exit(x^.sum);
mid:=(x^.l+x^.r)div 2;
if r<=mid then exit(ask(x^.ls,l,r));
if l>mid then exit(ask(x^.rs,l,r));
exit(ask(x^.ls,l,mid)+ask(x^.rs,mid+1,r));
end;
procedure change(x:tree;l,r:longint);
var
mid:longint;
begin
if (x^.max<=1)and(x^.min>=0) then exit;
if x^.l=x^.r then begin x^.sum:=trunc(sqrt(x^.sum)); x^.max:=x^.sum; x^.min:=x^.sum; exit; end;
mid:=(x^.l+x^.r)div 2;
if r<=mid then change(x^.ls,l,r)
else if l>mid then change(x^.rs,l,r)
else
begin
change(x^.ls,l,mid);
change(x^.rs,mid+1,r);
end;
updata(x);
end;
begin
ca:=0;
while(not eof) do
begin
inc(ca);
writeln('Case #',ca,':');
readln(n);
for i:=1 to n do
read(delta[i]);
new(xtr);
build(xtr,1,n);
readln(q);
for i:=1 to q do
begin
readln(opt,l,r);
if l>r then begin t:=l;l:=r;r:=t; end;
if opt=1 then writeln(ask(xtr,l,r))
else change(xtr,l,r);
end;
writeln;
end;
end. [spoj2916]GSS5
分类讨论：
if(y1"%d\n",xtr->query(y1+1,x2-1).sum+xtr->query(x1,y1).rgss+xtr->query(x2,y2).lgss);
else
{
int l1,r1,l2,r2,s;
r1=xtr->query(x1,x2-1).rgss;
r2=xtr->query(x1,y1).rgss;
l1=xtr->query(y1+1,y2).lgss;
l2=xtr->query(x2,y2).lgss;
s=xtr->query(x2,y1).gss;
printf("%d\n",max(s,max(r1+l2,r2+l1)));
} 代码：
#include"%d\n",xtr->query(y1+1,x2-1).sum+xtr->query(x1,y1).rgss+xtr->query(x2,y2).lgss);
else
{
int l1,r1,l2,r2,s;
r1=xtr->query(x1,x2-1).rgss;
r2=xtr->query(x1,y1).rgss;
l1=xtr->query(y1+1,y2).lgss;
l2=xtr->query(x2,y2).lgss;
s=xtr->query(x2,y1).gss;
printf("%d\n",max(s,max(r1+l2,r2+l1)));
}
}
}
return 0;
} [SPOJ4487]
带插入，删除，修改的最大子段和，用splay即可。
读操作符的时候不要用getchar(),卡了两个多小时。。。
因为不同于线段树，splay上的一个点也表示序列的一个位置（子树表示一个区间），所以更新答案有所不同：
sum=l->sum+r->sum+t;
num=l->num+r->num+1;
lgss=max(l->lgss,l->sum+t+max(0,r->lgss));
rgss=max(r->rgss,r->sum+t+max(0,l->rgss));
gss=max(max(l->gss,r->gss),max(0,l->rgss)+t+max(0,r->lgss));优化时间复杂度：O(n)建树（为了方便序列左右端再建两个没用的点，避免加点的时候找不到前驱或者后继），replace不要先delete再insert。
代码：
#includenew tree)->build(lx,mid-1);l->f=this;}
if(rx>mid){(r=new tree)->build(mid+1,rx);r->f=this;}
update();
}
}*str;
tree *find(int k,tree *p)
{
int tmp=p->l->num+1;
if(k==tmp) return p;
if(k>tmp) return find(k-tmp,p->r);
else return find(k,p->l);
}
void ins(int pl,int c,tree *&s)
{
tree *v=new tree;
find(pl-1,s)->splay(s,NUL);
find(pl,s)->splay(s,s);
v->init(c);s->r->l=v;v->f=s->r;
v->splay(s,NUL);
}
void del(int pl,tree *&s)
{
find(pl,s)->splay(s,NUL);
find(pl-1,s)->splay(s,s);
s->r->f=s->l; s->l->r=s->r; s->l->f=NUL; s=s->l;
}
void rpl(int pl,int c,tree *&s)
{
tree *p=find(pl,s);
p->init(c);
p->splay(s,NUL);
}
int query(int lx,int rx,tree *&s)
{
find(lx-1,s)->splay(s,NUL);
find(rx+1,s)->splay(s,s);
return s->r->l->gss;
}
int main()
{
n=read();
for(int i=1;i<=n;i++)
a[i]=read();
NUL=new tree;
(str=new tree)->build(0,n+1);
q=read();
for(int i=1;i<=q;i++)
{
char opt;
scanf("%s",&opt);
int x=read(),y;
if(opt=='I') {y=read();ins(x+1,y,str);}
if(opt=='D') {del(x+1,str);}
if(opt=='R') {y=read();rpl(x+1,y,str);}
if(opt=='Q') {y=read();printf("%d\n",query(x+1,y+1,str));}
}
return 0;
}