【题解】LuoGu2048:[NOI2010]超级钢琴

原题传送门
这道题目的id昭示着它的非同凡响
我可以把这题和别的几题结合起来,产生一种利用堆的节点最优思想
序列合并
最小函数值
以上是两道和这题有共同思想的

本题是要求出前 k k k大的 s u m j − s u m i − 1 sum_j-sum_{i-1} sumjsumi1的和
对于每个节点 i i i,首先寻出 [ i + l − 1 , i + r − 1 ] [i+l-1,i+r-1] [i+l1,i+r1]区间内的 p o s pos pos,使得 s u m p o s − s u m i − 1 > = s u m j − s u m i − 1 ( i + l − 1 < = j < = i + r − 1 ) sum_{pos}-sum_{i-1}>=sum_j-sum_{i-1}(i+l-1<=j<=i+r-1) sumpossumi1>=sumjsumi1(i+l1<=j<=i+r1)
每个 i i i p o s pos pos s t st st表求出
那么每个节点 i i i目前最优的有了,把 n n n个节点最优的放到堆里面,每次取出堆里最大的,假设取出节点 x x x,对应最优解为 y y y,节点 x x x对应的区间为 [ l , r ] [l,r] [l,r]

  • 因为从堆里取出来的,满足全局最优,所以统计到答案里面
  • 又因为这是对于节点 x x x的最优,然后我们要找到最优之外次优的,把区间 [ l , r ] [l,r] [l,r]分成两半
    [ l , y − 1 ] [ y + 1 , r ] [l,y-1][y+1,r] [l,y1][y+1,r],分别求出最优位置放到堆里即可

Code:

#include <bits/stdc++.h>
#define maxn 500010
using namespace std;
int sum[maxn], power[25], st[maxn][25], id[maxn][25], n, k, l, r;
struct node{
	int id, l, r, pos;
	bool operator < (const node &x) const{ return sum[x.pos] - sum[x.id - 1] > sum[pos] - sum[id - 1]; }
};
priority_queue <node> q;

inline int read(){
	int s = 0, w = 1;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') w = -1;
	for (; isdigit(c); c = getchar()) s = (s << 1) + (s << 3) + (c ^ 48);
	return s * w;
}

int calc(int l, int r){
	int ans = 0, x = -1e9;
	for (int i = 20; i >= 0; --i)
		if (l + power[i] - 1 <= r){
			if (st[l][i] > x) x = st[l][i], ans = id[l][i];
			l += power[i];
		}
	return ans;
}

int main(){
	n = read(), k = read(), l = read(), r = read();
	for (int i = 1; i <= n; ++i) st[i][0] = sum[i] = sum[i - 1] + read(), id[i][0] = i;
	power[0] = 1;
	for (int i = 1; i <= 20; ++i) power[i] = power[i - 1] << 1; 
	for (int j = 1; j <= 20; ++j)
		for (int i = 1; i + power[j] - 1 <= n; ++i){
			if (st[i][j - 1] > st[i + power[j - 1]][j - 1]) st[i][j] = st[i][j - 1], id[i][j] = id[i][j - 1];
			else st[i][j] = st[i + power[j - 1]][j - 1], id[i][j] = id[i + power[j - 1]][j - 1];
		}
	for (int i = 1; i <= n - l + 1; ++i){
		int L = i + l - 1, R = min(n, i + r - 1), x = calc(L, R);
		q.push((node){i, L, R, x});
	}
	long long ans = 0;
	while (k--){
		node tmp = q.top(); q.pop();
		ans += sum[tmp.pos] - sum[tmp.id - 1];
		if (tmp.pos > tmp.l) q.push((node){tmp.id, tmp.l, tmp.pos - 1, calc(tmp.l, tmp.pos - 1)});
		if (tmp.pos < tmp.r) q.push((node){tmp.id, tmp.pos + 1, tmp.r, calc(tmp.pos + 1, tmp.r)});
	}
	printf("%lld\n", ans);
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值