k-Maximum Subsequence Sum
题目描述:
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3638
题解
数据较小时,考虑费用流。为什么我没想到(雾)。
由于本题n较大,所以考虑进行人工增广。
用线段树维护区间最大子串,每次增广将其取反即可。
本题有4倍经验!!!
代码
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<ctime>
#include<algorithm>
#define N 100010
using namespace std;
int n,m,ans,tot,s[N];
struct info{int sum,ml,mr,num,pl,pr,l,r;}q[22];
struct node{int tag;info minn,maxn;}t[N*4];
class seg_tree
{
void change(int x)
{
t[x].tag^=1;swap(t[x].maxn,t[x].minn);
}
void pushdown(int x)
{
if(!t[x].tag)return;
int lc=x<<1,rc=lc+1;t[x].tag=0;
change(lc);change(rc);
}
info merge(info a,info b)
{
info res;res.sum=a.sum+b.sum;
res.ml=a.ml;res.mr=b.mr;
res.pl=a.pl;res.pr=b.pr;
if(a.sum+b.ml>res.ml)
res.ml=a.sum+b.ml,res.pl=b.pl;
if(b.sum+a.mr>res.mr)
res.mr=b.sum+a.mr,res.pr=a.pr;
res.num=a.mr+b.ml;res.l=a.pr;res.r=b.pl;
if(a.num>res.num)
res.num=a.num,res.l=a.l,res.r=a.r;
if(b.num>res.num)
res.num=b.num,res.l=b.l,res.r=b.r;
return res;
}
void update(int x)
{
int lc=x<<1,rc=lc+1;
t[x].maxn=merge(t[lc].maxn,t[rc].maxn);
t[x].minn=merge(t[lc].minn,t[rc].minn);
}
public:
void build(int x,int l,int r)
{
t[x].tag=0;
if(l==r)
{
info res;
res.sum=res.num=res.ml=res.mr=s[l];
res.l=res.r=res.pl=res.pr=l;
t[x].maxn=res;
res.sum=res.num=res.ml=res.mr=-s[l];
t[x].minn=res;
return;
}
int mid=l+r>>1,lc=x<<1,rc=lc+1;
build(lc,l,mid);build(rc,mid+1,r);
update(x);
}
void modify(int x,int l,int r,int des,int val)
{
if(l==r)
{
info res;
res.sum=res.num=res.ml=res.mr=val;
res.l=res.r=res.pl=res.pr=l;
t[x].maxn=res;
res.sum=res.num=res.ml=res.mr=-val;
t[x].minn=res;
return;
}
int mid=l+r>>1,lc=x<<1,rc=lc+1;
pushdown(x);
if(des<=mid)modify(lc,l,mid,des,val);
else modify(rc,mid+1,r,des,val);
update(x);
}
void reverse(int x,int l,int r,int ql,int qr)
{
if(ql<=l&&r<=qr){change(x);return;}
int mid=l+r>>1,lc=x<<1,rc=lc+1;
pushdown(x);
if(ql<=mid)reverse(lc,l,mid,ql,qr);
if(qr>mid)reverse(rc,mid+1,r,ql,qr);
update(x);
}
info qry(int x,int l,int r,int ql,int qr)
{
if(ql<=l&&r<=qr)return t[x].maxn;
int mid=l+r>>1,lc=x<<1,rc=lc+1;
pushdown(x);
if(qr<=mid)return qry(lc,l,mid,ql,qr);
if(ql>mid)return qry(rc,mid+1,r,ql,qr);
return merge(qry(lc,l,mid,ql,qr),qry(rc,mid+1,r,ql,qr));
}
}T;
int main()
{
int tp,x,y,k;
scanf("%d",&n);
for(int i=1;i<=n;i++)scanf("%d",&s[i]);
T.build(1,1,n);
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
scanf("%d%d%d",&tp,&x,&y);
if(tp==1)
{
scanf("%d",&k);ans=0;tot=0;
while(k--)
{
info p=T.qry(1,1,n,x,y);
if(p.num<0)break;
T.reverse(1,1,n,p.l,p.r);
ans+=p.num;q[++tot]=p;
}
for(int i=1;i<=tot;i++)
T.reverse(1,1,n,q[i].l,q[i].r);
printf("%d\n",ans);
}
else T.modify(1,1,n,x,y);
}
return 0;
}