k-Maximum Subsequence Sum

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;
}