[DarkBZOJ3636]教义问答手册

最新推荐文章于 2024-08-14 09:58:33 发布

OneInDark

最新推荐文章于 2024-08-14 09:58:33 发布

阅读量339

点赞数

分类专栏： C++ 动态规划文章标签：动态规划

本文链接：https://blog.csdn.net/qq_42101694/article/details/113828810

版权

C++ 同时被 2 个专栏收录

547 篇文章 3 订阅

订阅专栏

动态规划

220 篇文章 2 订阅

订阅专栏

题目

传送门 to DarkBZOJ

题目概要
给出长度为 $n$ 的序列 $a$ ，有 $q$ 次询问：在区间 $[l, r]$ 中选出若干个长度为 $L$ 的子区间，两两无交，使被覆盖的元素的和最大。

数据范围与提示
$n=q\in\{10^3,10^5\}$ 但 $L\le 50$ 。

思路

分治，处理经过中线 $m$ 的询问（即满足 $l\leqslant m\leqslant r$ 的询问）。枚举中间选了区间 $[i,i+L)\;(m-L<i\leqslant m)$ 后，预处理所有区间 $[x, i)$ 和 $[i + L, x]$ 的答案——即在该范围内随意选若干个长度为 $L$ 的子区间的最大求和。显然可以 $\mathcal O(n)$ 递推。然后可以更新一次询问。

不经过中线的询问则划分到两边去。一共 $\mathcal O(\log n)$ 层（即询问的最多被处理次数），而每一层还有额外复杂度 $\mathcal O(nL)$ ，所以总复杂度
$\mathcal O(nL\log n+q\log n)$

本质是猫树。~~虽然它一点也不像树啊~~。

代码

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
typedef long long int_;
inline int readint(){
	int x; scanf("%d",&x); return x;
}

const int MaxN = 200005;
int_ s[MaxN], dp[MaxN]; int L;
void leftDP(int x,int l){
	for(int i=x; i>x-L+1&&i>=l; --i)
		dp[i] = 0;
	if(x-L+1 < l) return ;
	dp[x-L+1] = max(0ll,s[x]-s[x-L]);
	for(int i=x-L; i>=l; --i)
		dp[i] = max(dp[i+1],
			dp[i+L]+s[i+L-1]-s[i-1]);
}
void rightDP(int x,int r){
	for(int i=x; i<x+L-1; ++i)
		dp[i] = 0;
	dp[x+L-1] = max(0ll,s[x+L-1]-s[x-1]);
	for(int i=x+L; i<=r; ++i)
		dp[i] = max(dp[i-1],
			dp[i-L]+s[i]-s[i-L]);
}

struct Query{
	int l, r, id;
};
Query q[MaxN], xyx[MaxN];
int_ ans[MaxN];

void solve(int l,int r,int ql,int qr){
	if(r-l+1 < L || ql > qr) return ;
	int mid = (l+r)>>1;
	for(int j=max(mid-L+2,l); j<=mid; ++j){
		leftDP(j-1,l), rightDP(j+L,r);
		dp[j] = dp[j+L-1] = 0;
		for(int i=ql; i<=qr; ++i)
			if(q[i].l <= j && j+L-1 <= q[i].r)
				ans[q[i].id] = max(ans[q[i].id],
					dp[q[i].l]+dp[q[i].r]
					+s[j+L-1]-s[j-1]);
	}
	leftDP(mid,l), rightDP(mid+1,r);
	for(int i=ql; i<=qr; ++i)
		if(q[i].l <= mid && mid < q[i].r)
			ans[q[i].id] = max(ans[q[i].id],
				dp[q[i].l]+dp[q[i].r]);
	int pl = ql, pr = qr;
	for(int i=ql; i<=qr; ++i)
		if(q[i].r <= mid)
			xyx[pl ++] = q[i];
		else if(q[i].l > mid)
			xyx[pr --] = q[i];
	for(int i=ql; i<=qr; ++i)
		q[i] = xyx[i];
	solve(l,mid,ql,pl-1);
	solve(mid+1,r,pr+1,qr);
}

int main(){
	int n = readint(); L = readint();
	for(int i=1; i<=n; ++i)
		s[i] = s[i-1]+readint();
	int Q = readint();
	for(int i=1; i<=Q; ++i){
		q[i].l = readint();
		q[i].r = readint();
		q[i].id = i;
	}
	solve(1,n,1,Q);
	for(int i=1; i<=Q; ++i)
		printf("%lld\n",ans[i]);
	return 0;
}