题意:区间最大字段和
思路:猫树模板题或合并子区间查询线段树。
猫树
#include<bits/stdc++.h>
using namespace std;
const int N=100010;
int n,len,l,m;
int lg[N<<2],a[N],pos[N],p[21][N],s[21][N];//考虑到当前位置的最大子段和,p强制端点,s非强制端点(普通的最大子段和)
void build(int l,int r,int k,int d){
if(l==r){
pos[l]=k;//记录这个位置的编号
return ;
}
int f1,f2;//端点,非端点
int mid=(l+r)>>1;
p[d][mid]=s[d][mid]=f1=f2=a[mid];
f2=max(f2,0);
for(int i=mid-1;i>=l;i--){
f1+=a[i];
f2+=a[i];
p[d][i]=max(p[d][i+1],f1);
s[d][i]=max(s[d][i+1],f2);
f2=max(f2,0);
}
p[d][mid+1]=s[d][mid+1]=f1=f2=a[mid+1];
f2=max(f2,0);
for(int i=mid+2;i<=r;i++){
f1+=a[i];
f2+=a[i];
p[d][i]=max(p[d][i-1],f1);
s[d][i]=max(s[d][i-1],f2);
f2=max(f2,0);
}
build(l,mid,k<<1,d+1);
build(mid+1,r,k<<1|1,d+1);
return;
}
int query(int x,int y){
if(x==y)
return a[x];
int d=lg[pos[x]]-lg[pos[x]^pos[y]];
return max(max(s[d][x],s[d][y]),p[d][x]+p[d][y]);
}
int main(){
int x,y;
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
for(len=2;len<n;len<<=1);
l=len<<1;
for(int i=2;i<=l;i++)
lg[i]=lg[i>>1]+1;
build(1,len,1,1);
scanf("%d",&m);
while(m--){
scanf("%d%d",&x,&y);
printf("%d\n",query(x,y));
}
return 0;
}
线段树
#include<bits/stdc++.h>
using namespace std;
const int maxn=50010;
int n,m,x,y;
struct segmentTree{
int lmax,rmax,sum,dat;
}tr[maxn<<2];
int a[maxn];
void up(int p){
tr[p].sum=tr[p<<1].sum+tr[p<<1|1].sum;
tr[p].lmax=max(tr[p<<1].lmax,tr[p<<1].sum+tr[p<<1|1].lmax);
tr[p].rmax=max(tr[p<<1|1].rmax,tr[p<<1|1].sum+tr[p<<1].rmax);
tr[p].dat=max(max(tr[p<<1].dat,tr[p<<1|1].dat),tr[p<<1].rmax+tr[p<<1|1].lmax);
}
void build(int p,int l,int r){
if(l==r){
tr[p].dat=tr[p].sum=tr[p].lmax=tr[p].rmax=a[l];
return ;
}
int mid=(l+r)>>1;
build(p<<1,l,mid);
build(p<<1|1,mid+1,r);
up(p);
}
segmentTree query(int p,int l,int r,int ql,int qr){
if(ql<=l&&qr>=r)
return tr[p];
int mid=(l+r)>>1;
if(ql>mid)
return query(p<<1|1,mid+1,r,ql,qr);
if(qr<=mid)
return query(p<<1,l,mid,ql,qr);
else{
segmentTree ans,a,b;
a=query(p<<1,l,mid,ql,qr);
b=query(p<<1|1,mid+1,r,ql,qr);
ans.sum=a.sum+b.sum;
ans.dat=max(a.dat,a.rmax+b.lmax);
ans.dat=max(ans.dat,b.dat);
ans.lmax=max(a.lmax,a.sum+b.lmax);
ans.rmax=max(b.rmax,b.sum+a.rmax);
return ans;
}
}
int main(){
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
int x,y;
build(1,1,n);
scanf("%d",&m);
for(int i=1;i<=m;i++){
scanf("%d%d",&x,&y);
printf("%d\n",query(1,1,n,x,y).dat);
}
return 0;
}
dalao说一个算法学不会就写一遍,还学不会就再写两遍,三遍.....