SP1043 GSS1 - Can you answer these queries I

二分抄代码

于 2020-12-04 00:20:39 发布

阅读量232

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/liufengwei1/article/details/110603216

本文提供了一道来自洛谷平台的动态规划题目解析及C++实现代码。该题涉及区间更新、查询最大值等操作，并通过树状数组和矩阵乘法进行优化，适用于竞赛编程爱好者和技术进阶学习。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

https://www.luogu.com.cn/problem/SP1043

学习自洛谷日报https://www.luogu.com.cn/blog/gkxx-is-here/what-the-hell-is-

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

const int maxl=1e6+10;
const ll inf=1ll<<60;

int n,q,cas;ll ans;
ll a[maxl];
struct matrix
{
	ll a[3][3];
	matrix()
	{
		for(int i=0;i<3;i++)
			for(int j=0;j<3;j++)
				a[i][j]=-inf; 
	}
	matrix operator * (const matrix &b)const
	{
		matrix ret;
		for(int i=0;i<3;i++)
			for(int j=0;j<3;j++)
				for(int k=0;k<3;k++)
					ret.a[i][k]=max(ret.a[i][k],a[i][j]+b.a[j][k]);
		return ret;
	}
};
struct node
{
	int l,r;
	matrix a;
}tr[maxl*4];

inline void upda(matrix &d,ll x)
{
	d.a[0][0]=d.a[0][2]=x;d.a[0][1]=-inf;
	d.a[1][0]=d.a[1][2]=x;d.a[1][1]=0;
	d.a[2][0]=d.a[2][1]=-inf;d.a[2][2]=0;
}

inline void build(int k,int l,int r)
{
	tr[k].l=l;tr[k].r=r;
	if(l==r)
	{
		upda(tr[k].a,a[l]);
		return;
	}
	int mid=(l+r)>>1;
	build(k<<1,l,mid);
	build(k<<1|1,mid+1,r);
	tr[k].a=tr[k<<1].a*tr[k<<1|1].a;
}

inline void prework()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%lld",&a[i]);
	build(1,1,n);
}

inline void upd(int k,int l,int x)
{
	if(tr[k].l==tr[k].r)
	{
		upda(tr[k].a,x);
		return;
	}
	int mid=(tr[k].l+tr[k].r)>>1;
	if(l<=mid)
		upd(k<<1,l,x);
	else 
		upd(k<<1|1,l,x);
	tr[k].a=tr[k<<1].a*tr[k<<1|1].a;
}

inline matrix qry(int k,int l,int r)
{
	if(tr[k].l==l && tr[k].r==r)
		return tr[k].a;
	int mid=(tr[k].l+tr[k].r)>>1;
	if(r<=mid)
		return qry(k<<1,l,r);
	else if(l>mid)
		return qry(k<<1|1,l,r);
	else
		return qry(k<<1,l,mid)*qry(k<<1|1,mid+1,r); 
}

inline void mainwork()
{
	scanf("%d",&q);
	for(int i=1;i<=q;i++)
	{
		int op=1,x,y;
		scanf("%d%d",&x,&y);
		if(op==0)
		{
			upd(1,x,y);
			a[x]=y;
		}
		else
		{
			matrix tmp=qry(1,x,y);
			ans=max(tmp.a[1][0],tmp.a[1][2]);
			printf("%lld\n",ans);
		}
	}
}

inline void print()
{
	
}

int main()
{
	int t=1;
	//scanf("%d",&t);
	for(cas=1;cas<=t;cas++)
	{
		prework();
		mainwork();
		print();
	}
	return 0;
}

ddp