标签：tmp rw int Sum Maximum mid lw Subsequence mx

正题

题目链接:https://www.luogu.com.cn/problem/CF280D

题目大意

一个长度为\(n\)的序列，\(m\)次操作

1. 修改一个数
2. 询问一个区间中选出\(k\)段不交子段使得和最大

\(1\leq n\leq 10^5,1\leq m\leq 2\times 10^5,1\leq k\leq 20\)

解题思路

考虑模拟一下费用流，发现费用流的每次增广就是不停找到一个最大的子段取反。

用线段树维护最大子段，然后取反的话维护一个正的一个反的。

每次暴力做\(k\)次再倒流回去就好了。

时间复杂度：\(O(nk\log n)\)

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stack>
#define mp(x,y) make_pair(x,y)
using namespace std;
const int N=1e5+10;
struct node{
	int lw,rw,w;
	int mx,l,r,L,R;
}w[N<<2],v[N<<2];
int n,m,lazy[N<<2];
stack<pair<int,int> >s;
node Merge(node L,node R){
	node tmp;tmp.w=L.w+R.w;
	if(L.mx>R.mx)tmp.mx=L.mx,tmp.l=L.l,tmp.r=L.r;
	else tmp.mx=R.mx,tmp.l=R.l,tmp.r=R.r;
	if(L.rw+R.lw>tmp.mx)tmp.mx=L.rw+R.lw,tmp.l=L.R,tmp.r=R.L;
	if(L.lw>=L.w+R.lw)tmp.lw=L.lw,tmp.L=L.L;
	else tmp.lw=L.w+R.lw,tmp.L=R.L;
	if(R.rw>=R.w+L.rw)tmp.rw=R.rw,tmp.R=R.R;
	else tmp.rw=R.w+L.rw,tmp.R=L.R;
	return tmp;
}
void Downdata(int x){
	if(!lazy[x])return;
	lazy[x*2]^=1;swap(w[x*2],v[x*2]);
	lazy[x*2+1]^=1;swap(w[x*2+1],v[x*2+1]);
	lazy[x]=0;return;
}
void Change(int x,int L,int R,int l,int r){
	if(L==l&&R==r){swap(w[x],v[x]);lazy[x]^=1;return;}
	int mid=(L+R)>>1;Downdata(x);
	if(r<=mid)Change(x*2,L,mid,l,r);
	else if(l>mid)Change(x*2+1,mid+1,R,l,r);
	else Change(x*2,L,mid,l,mid),Change(x*2+1,mid+1,R,mid+1,r);
	w[x]=Merge(w[x*2],w[x*2+1]);
	v[x]=Merge(v[x*2],v[x*2+1]);
	return;
}
void Updata(int x,int L,int R,int pos,int val){
	if(L==R){
		w[x].w=w[x].mx=w[x].lw=w[x].rw=val;w[x].l=w[x].r=w[x].L=w[x].R=pos;
		v[x].w=v[x].mx=v[x].lw=v[x].rw=-val;v[x].l=v[x].r=v[x].L=v[x].R=pos;
		return;
	}
	int mid=(L+R)>>1;Downdata(x);
	if(pos<=mid)Updata(x*2,L,mid,pos,val);
	else Updata(x*2+1,mid+1,R,pos,val);
	w[x]=Merge(w[x*2],w[x*2+1]);
	v[x]=Merge(v[x*2],v[x*2+1]);
	return;
}
node Ask(int x,int L,int R,int l,int r){
	if(L==l&&R==r)return w[x];
	int mid=(L+R)>>1;Downdata(x);
	if(r<=mid)return Ask(x*2,L,mid,l,r);
	if(l>mid)return Ask(x*2+1,mid+1,R,l,r);
	return Merge(Ask(x*2,L,mid,l,mid),Ask(x*2+1,mid+1,R,mid+1,r));
}
int main()
{
	scanf("%d",&n);
	for(int i=1,x;i<=n;i++)
		scanf("%d",&x),Updata(1,1,n,i,x);
	scanf("%d",&m);
	while(m--){
		int op;scanf("%d",&op);
		if(op==0){
			int x,w;
			scanf("%d%d",&x,&w);
			Updata(1,1,n,x,w);
			node tmp=Ask(1,1,n,4,9);
			tmp.w++;tmp.w--;
		}
		else{
			int l,r,k,ans=0;
			scanf("%d%d%d",&l,&r,&k);
			while(k){
				node tmp=Ask(1,1,n,l,r);
				if(tmp.mx>0){
					ans+=tmp.mx;k--;
					s.push(mp(tmp.l,tmp.r));
					Change(1,1,n,tmp.l,tmp.r);
				}
				else break;
			}
			printf("%d\n",ans);
			while(!s.empty())
				Change(1,1,n,s.top().first,s.top().second),s.pop();
		}
	}
	return 0;
}

标签：tmp,rw,int,Sum,Maximum,mid,lw,Subsequence,mx
来源： https://www.cnblogs.com/QuantAsk/p/15328345.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。