版权声明:本文为一名蒟蒻的原创文章,大神转载的话顺便说个出处呗。 https://blog.csdn.net/cgh_Andy/article/details/78629633
似乎还没有人写题解,所以访问量是不是应该多点?
首先,这题暴力是能过的(滑稽
然后其实看数据范围得算法可知,用分块来搞就好了
维护的是块内最大和,左起最大和,右起最大和,这个东西用凸包维护就好了
而且斜率还是单调的(加的数>0
接着你发现块内重构 是要
本来想看看能不能再快一点的。。结果越改越慢
如果有人知道怎么样更快 请告诉我一声
可能我用结构体有点多?(代码还是挺能看的吧。。
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=5e4+3,bk=40;
const LL inf=1e15;
char B[1<<14],*S=B,*T=B;
#define gc (S==T&&(T=(S=B)+fread(B,1,1<<14,stdin),S==T)?-1:*S++)
inline int read(){
int x=0,f=1; char ch=gc;
while(ch<'0' || ch>'9'){if(ch=='-')f=-1; ch=gc;}
while(ch>='0' && ch<='9'){x=(x<<1)+(x<<3)+ch-'0'; ch=gc;}
return x*f;
}
int bl[N],a[N];
struct P{double x,y;};
double xl(P x,P y){
return (x.y-y.y)/(x.x-y.x);
}
struct tb{
int tp,o; P s[bk+2];
void ins(P u){
while(tp>1 && xl(u,s[tp-1])>=xl(s[tp],s[tp-1]) )--tp;
s[++tp]=u;
}
void move(int u){
while(o<tp && xl(s[o],s[o+1])>=u)++o;
}
};
struct block{
int l,r,k; LL num;
tb s[3];
void bulid(){
num=0;
for(int i=l;i<=r;++i) a[i]+=k,num+=a[i];
int n=r-l+1; k=0;
s[0].tp=s[1].tp=s[2].tp=0;
s[0].o =s[1].o =s[2].o=1;
for(int i=1;i<=n;++i){
P u=(P){i,-inf}; LL o=0;
for(int j=1;j<i;++j) o+=a[l+j-1];
for(int j=i;j<=n;++j){
o+=a[l+j-1];
if(j!=i)o-=a[l+j-1-i];
u.y=max(u.y,(double)o);
}
s[0].ins(u);
}
P u=(P){0,0};
for(int i=1;i<=n;++i){
u.y+=a[l+i-1],u.x=i;
s[1].ins(u);
}
u=(P){0,0};
for(int i=n;i;--i){
u.y+=a[l+i-1],u.x=n-i+1;
s[2].ins(u);
}
for(int i=0;i<3;++i) s[i].move(0);
}
void add(int u){
num+=(LL)u*(r-l+1),k+=u;
for(int i=0;i<3;++i) s[i].move(-k);
}
LL get(int u){
P o=s[u].s[ s[u].o ];
return o.y+(LL)o.x*k;
}
}b[1302];
int main(){
int n=read(),q=read(),i;
for(i=1;i<=n;++i)
a[i]=read(),bl[i]=(i-1)/bk+1,b[bl[i]].r=i;
for(i=1;i<=bl[n];++i)
b[i].k=0,b[i].l=(i-1)*bk+1,b[i].bulid();
while(q--){
char ch=gc; while(ch!='A' && ch!='Q')ch=gc;
if(ch=='A'){
int l=read(),r=read(),x=read();
for(i=l;bl[i]==bl[l] && i<=r;++i) a[i]+=x;
b[bl[l]].bulid();
for(i=bl[l]+1;i<bl[r];++i) b[i].add(x);
if(bl[l]<bl[r]){
for(i=r;bl[i]==bl[r];--i) a[i]+=x;
b[bl[r]].bulid();
}
}
else{
int l=read(),r=read();
LL ans=0,s=0;
for(i=l;bl[i]==bl[l] && i<=r;++i){
s+=a[i]+b[bl[i]].k;
s=s<0?0:s; ans=s>ans?s:ans;
}
for(i=bl[l]+1;i<bl[r];++i){
ans=max(ans,b[i].get(0));
ans=max(ans,s+b[i].get(1));
s=max(s+b[i].num,b[i].get(2));
s=s<0?0:s;
}
if(bl[l]<bl[r])
for(i=b[bl[r]].l;i<=r;++i){
s+=a[i]+b[bl[r]].k;
s=s<0?0:s; ans=s>ans?s:ans;
}
printf("%lld\n",ans);
}
}
return 0;
}