P4064 [JXOI2017]加法-题解报告

0x00 前言

氧化钙,调了半天,二分假了。

前置知识:优先队列,贪心,二分,一种可以维护区间和的数据结构(线段树,树状数组…)。

0x01 二分

首先答案肯定是成单调性的,如果答案越大显然越难凑出来,从而我们利用二分判断答案是否可行即可。

0x02 贪心

贪心思路还是需要一点悟性,因为其中涉及到的贪心比较多,但也比较细致。

下面一一来分析分析。

我们假设我们需要判断的值是 x x x

我们去枚举点,从小到大枚举。(正常人会这么想)

我们在选择区间的时候,对于已经大于了 x x x a i a_i ai,我们不需要加区间,直接跳过。

a i < x a_i<x ai<x 时,需要用区间去覆盖。

因为在处理第 i i i 个点时,前 i − 1 i-1 i1 个点必然已经处理过了,所以选择区间的 l ≤ i l\le i li 就好了。

r r r 肯定是越大越好,因为贪心的想,选择的区间有限,覆盖的越远后面肯定加的越多,所以 l l l 合理的情况下 r r r 越大越好。

关于 l l l 我们可以在排序时预处理出来,没有什么好说的。

对于 r r r,我们可以采用优先队列处理。

然后对于区间加的和,用树状数组或线段树处理即可(线段树稍稍麻烦一些,且常数大,初始化复杂,这里用的是树状数组)。

最后,如果这个点加完了所有区间仍然不满足条件或者没有可行区间选择则答案不行。

0x03 复杂度

二分: log ⁡ n \log n logn

二分判断函数:

  • 遍历每一个点: O ( n ) O(n) O(n)

  • 优先队列和树状数组: O ( log ⁡ n ) O(\log n) O(logn)

所以复杂度是 O ( n log ⁡ 3 n ) O({n\log^3 n}) O(nlog3n)

0x04 代码

#include<bits/stdc++.h>
#define inf 0x3f3f3f3f
#define lowbit(x) x&-x
using namespace std;
const int N=2e5+5;
const int M=2e5+5;
int n,m,k,p;
int a[N],tr[N];
struct node{
	int L,R;
	friend bool operator<(node x,node y){
		return x.R<y.R;
	}
}s[M];
bool cmp(node x,node y){
	return x.L==y.L?x.R>y.R:x.L<y.L;
}
void add(int x,int v){
	for(;x<=n;x+=lowbit(x)) tr[x]+=v;
}
int ask(int x){
	int res=0;
	for(;x;x-=lowbit(x)) res+=tr[x];
	return res;
}
bool check(int P){
	memset(tr,0,sizeof(tr));
	priority_queue<node>q;
	int dq=1,u=0;
	for(int i=1;i<=n;i++) add(i,a[i]-a[i-1]);
	for(int i=1;i<=n;i++){
		while(dq<=m && s[dq].L<=i){
			if(s[dq].R>=i) q.push((node){s[dq].L,s[dq].R});
			dq++;
		}
		while(ask(i)<P){
			if(q.empty()) return 0;
			if(u==k) return 0;
			if(q.top().L<=i && q.top().R>=i){
				u++;
				add(q.top().L,p);
				add(q.top().R+1,-p);	
			}
			q.pop();
		}
	}
	return 1;
}
int f(int L,int R){
	if(L==R) return L;
	if(L+1==R){
		if(check(R)) return R;
		return L;
	}
	int mid=(L+R)/2;
	if(check(mid)) return f(mid,R);
	else return f(L,mid-1);
}
int main(){
	int T;
	scanf("%d",&T);
	while(T--){
		scanf("%d %d %d %d",&n,&m,&k,&p);
		int amin=inf;
		for(int i=1;i<=n;i++) scanf("%d",&a[i]),amin=min(amin,a[i]);
		for(int i=1;i<=m;i++) scanf("%d %d",&s[i].L,&s[i].R);
		sort(s+1,s+1+m,cmp);
		printf("%d\n",f(0,1e8));
	}
	return 0;
} 
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值