Problem 1016 咒文卷轴 优先队列+前缀和+rmq

题目链接：

题目

Problem 1016 咒文卷轴

Time Limit: 3000 mSec

Memory Limit : 131072 KB

问题描述

小Y 是一个魔法师，有一天他获得了一卷神秘而古老的咒文卷轴，其由N个咒文构成，每一个咒文都有一个威力值ai, 现在小Y可以借助该卷轴释放强力的魔法，一个魔法必须由编号连续的一段的、长度大于等于L且小于等于R的咒文构成，该魔法的威力为构成该魔法的每一个咒文的威力的总和，现在小Y想要释放K个不同的魔法，问最大能产生多大的威力值，若两个魔法是相同的，则其咒文编号的集合应该完全相同。

输入

第一行包含四个整数N,K,L,R,分别表示咒文个数，需要释放的魔法个数，以及魔法包含的咒文术的下限与上限。

接下来n行，每行一个数字ai，表示按编号从小到大每个咒文的威力值。

N<=100000

K<=100000

-1000<=ai<=1000,1<=L<=R<=N并且保证一定能释放k种不同的法术。

输出

仅一行，表示K个不同魔法威力值之和的最大值。要换行。

样例

input

4 3 2 3

3

2

-6

8

output

11

题解

考虑一个简单的版本，0 <= ai <= 1000，那么对于i开始的魔法，结尾为[l, r]之间的魔法肯定是递增的，所以只需要用一个优先队列，先把n个位置的值都丢进去，然后每次取出最大的[i, r]，然后在把[i, r - 1]丢回优先队列，直到r < l就不在丢进去，这样贪心能保证一定是从大到小取。

那么对于ai可以为负数的情况，前缀和就不满足递增了，然而对于一段区间[l, r]的和为两个前缀和相减，sum[r] - sum[l - 1]那么在左端点l固定的情况下，sum[r]越大越好，所以对于刚才那个问题等于是要取出优先队列中，区间[l, r]能取到前缀和值最大的位置mid，然后把区间在分成两部分[l, mid - 1], [mid + 1, r]丢回优先队列中，那么要维护一个区间最大值是个rmq问题，可以倍增st表O(nlogn)预处理一下，每次询问就是O(1)的，优先队列中的优先级就是区间最大值大的优先。这样就保证了取到的前k个法术一定是最大的。

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
#define X first
#define Y second
#define mp make_pair
using namespace std;
typedef __int64 LL;
const int maxn=1e5+10;
int arr[maxn];
LL sumv[maxn];
int n,k,l,r;
struct Node{
	int i,l,r,pos;
	Node(int i,int l,int r,int pos):i(i),l(l),r(r),pos(pos){}
	bool operator < (const Node& tmp) const {
		return sumv[pos]-sumv[i-1]<sumv[tmp.pos]-sumv[tmp.i-1];
	}
};
LL dp[maxn][20];
int pos[maxn][20];
void rmq(){
	for(int i=0;i<=n;i++) dp[i][0]=sumv[i],pos[i][0]=i;
	for(int j=1;(1<<j)<=n+1;j++){
		for(int i=0;i+(1<<j)-1<=n;i++){
			if(dp[i][j-1]<dp[i+(1<<(j-1))][j-1]){
				dp[i][j]=dp[i+(1<<(j-1))][j-1];
				pos[i][j]=pos[i+(1<<(j-1))][j-1];
			}
			else{
				dp[i][j]=dp[i][j-1];
				pos[i][j]=pos[i][j-1];
			}
		}
	}
}
int Max(int l,int r){
	if(l>r) return 10086;
	int k=0;
	while((1<<k)<=r-l+1) k++; k--;
	if(dp[l][k]<dp[r-(1<<k)+1][k]){
		return pos[r-(1<<k)+1][k];
	}
	return pos[l][k];
}
int main(){
	while(scanf("%d%d%d%d",&n,&k,&l,&r)==4){
		sumv[0]=0;
		for(int i=1;i<=n;i++){
			scanf("%d",&arr[i]);
			sumv[i]=sumv[i-1]+arr[i];
		}
		rmq();
		LL ans=0;
		priority_queue<Node> pq;
		for(int i=1;i+r-1<=n;i++){
			pq.push(Node(i,i+l-1,i+r-1,Max(i+l-1,i+r-1)));
		}
		for(int i=n+2-r;i+l-1<=n;i++){
			pq.push(Node(i,i+l-1,n,Max(i+l-1,n)));
		}
		while(k--){
			while(pq.top().l>pq.top().r) pq.pop();
			Node nd=pq.top(); pq.pop();
			ans+=sumv[nd.pos]-sumv[nd.i-1];
			pq.push(Node(nd.i,nd.l,nd.pos-1,Max(nd.l,nd.pos-1)));
			pq.push(Node(nd.i,nd.pos+1,nd.r,Max(nd.pos+1,nd.r)));
		}
		printf("%I64d\n",ans);
	}
	return 0;
}