[SPOJ 1716 GSS3] Can you answer these queries III(动态 DP) | 错题本

题目

Can you answer these queries III

分析

本题当然可以用传统线段树做
f i f_i fi 表示以 a i a_i ai 结尾的最大子段和, g i g_i gi 表示前 i i i 个数的最大子段和。则 f i = max ⁡ { f i − 1 + a i , a i } g i = max ⁡ { f i , g i − 1 } = max ⁡ { f i − 1 + a i , a i , g i − 1 } \begin{aligned} f_i &= \max\{f_{i - 1} + a_i, a_i\} \\ g_i &= \max\{f_i, g_{i - 1}\} \\ &=\max\{f_{i - 1} + a_i, a_i, g_{i - 1}\} \end{aligned} figi=max{fi1+ai,ai}=max{fi,gi1}=max{fi1+ai,ai,gi1} 更改一下矩阵乘法的定义( + → max ⁡ , × → + + \to \max, \times \to + +max,×+)不难得到矩阵转移(我习惯列乘行) ( f i g i 0 ) = ( f i − 1 g i − 1 0 ) × ( a i − ∞ a i a i 0 a i − ∞ − ∞ 0 ) \begin{pmatrix} f_i \\ g_i \\ 0 \end{pmatrix} = \begin{pmatrix} f_{i - 1} \\ g_{i - 1} \\ 0 \end{pmatrix} \times \begin{pmatrix} a_i & -\infty & a_i \\ a_i & 0 & a_i \\ -\infty & -\infty & 0 \end{pmatrix} figi0=fi1gi10×aiai0aiai0
线段树维护区间矩阵乘积即可。

错因

  • 算矩阵乘法时初始矩阵应初始化为全 − ∞ -\infty ,而不是全 0 0 0

代码

注意矩阵初始化和计算答案的细节。

#include <algorithm>
#include <cstdio>
#include <cstring>

const int INF = 0x3f3f3f3f;
const int MAXN = 50000;
const int MAXP = 3;

int N, Q;
int A[MAXN + 5];

struct Matrix {
	int n, m;
	int Mat[MAXP + 5][MAXP + 5];

	Matrix(int _n = 0, int _m = 0) {
		n = _n, m = _m;
		for (int i = 1; i <= n; i++)
			for (int j = 1; j <= m; j++)
				Mat[i][j] = -INF;
	}

	int* operator [] (const int &i) {
		return Mat[i];
	}

	Matrix operator * (Matrix other) {
		Matrix ret(n, other.m);
		for (int i = 1; i <= ret.n; i++)
			for (int j = 1; j <= ret.m; j++)
				for (int k = 1; k <= m; k++)
					ret[i][j] = std::max(ret[i][j], Mat[i][k] + other[k][j]);
		return ret;
	}

	void Debug() {
		for (int i = 1; i <= n; i++, puts(""))
			for (int j = 1; j <= m; j++)
				printf("%d ", Mat[i][j]);
	}
};

struct SegmentTree {
	#define lch (u << 1)
	#define rch (u << 1 | 1)

	Matrix Dp[(MAXN << 2) + 5];

	void PushUp(int u) {
		Dp[u] = Dp[lch] * Dp[rch];
	}

	void Build(int u, int lft, int rgt) {
		if (lft == rgt) {
			Dp[u].n = Dp[u].m = 3;
			Dp[u][2][2] = Dp[u][3][3] = 0;
			Dp[u][1][1] = Dp[u][1][3] = A[lft];
			Dp[u][2][1] = Dp[u][2][3] = A[lft];
			Dp[u][1][2] = Dp[u][3][1] = Dp[u][3][2] = -INF;
			return;
		}
		int mid = (lft + rgt) >> 1;
		Build(lch, lft, mid);
		Build(rch, mid + 1, rgt);
		PushUp(u);
	}

	void Modify(int u, int lft, int rgt, int pos, int val) {
		if (lft == rgt) {
			Dp[u][1][1] = Dp[u][1][3] = val;
			Dp[u][2][1] = Dp[u][2][3] = val;
			return;
		}
		int mid = (lft + rgt) >> 1;
		if (pos <= mid)
			Modify(lch, lft, mid, pos, val);
		else
			Modify(rch, mid + 1, rgt, pos, val);
		PushUp(u);
	}

	Matrix Query(int u, int lft, int rgt, int l, int r) {
		if (l <= lft && rgt <= r)
			return Dp[u];
		int mid = (lft + rgt) >> 1;
		if (mid < l)
			return Query(rch, mid + 1, rgt, l, r);
		if (r <= mid)
			return Query(lch, lft, mid, l, r);
		return Query(lch, lft, mid, l, r) * Query(rch, mid + 1, rgt, l, r);
	}

	#undef lch
	#undef rch
}T;

int main() {
	scanf("%d", &N);
	for (int i = 1; i <= N; i++)
		scanf("%d", &A[i]);
	T.Build(1, 1, N);
	scanf("%d", &Q);
	while (Q--) {
		int opt, x, y;
		scanf("%d%d%d", &opt, &x, &y);
		if (opt == 0)
			T.Modify(1, 1, N, x, y);
		else {
			Matrix Ans = T.Query(1, 1, N, x, y);
			printf("%d\n", std::max(Ans[2][1], Ans[2][3]));
		}
	}
	return 0;
}
  • 2
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
洛谷的SPOJ需要注册一个SPOJ账号并进行绑定才能进行交题。您可以按照以下步骤进行注册: 1. 打开洛谷网站(https://www.luogu.com.cn/)并登录您的洛谷账号。 2. 在网站顶部导航栏中找到“题库”选项,将鼠标悬停在上面,然后选择“SPOJ”。 3. 在SPOJ页面上,您会看到一个提示,要求您注册SPOJ账号并进行绑定。点击提示中的链接,将会跳转到SPOJ注册页面。 4. 在SPOJ注册页面上,按照要求填写您的用户名、密码和邮箱等信息,并完成注册。 5. 注册完成后,返回洛谷网站,再次进入SPOJ页面。您会看到一个输入框,要求您输入刚刚注册的SPOJ用户名。输入用户名后,点击“绑定”按钮即可完成绑定。 现在您已经成功注册并绑定了SPOJ账号,可以开始在洛谷的SPOJ题库上刷题了。祝您顺利完成编程练习!\[1\]\[2\] #### 引用[.reference_title] - *1* *3* [(洛谷入门系列,适合洛谷新用户)洛谷功能全解](https://blog.csdn.net/rrc12345/article/details/122500057)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [luogu p7492 序列](https://blog.csdn.net/zhu_yin233/article/details/122051384)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值