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 i−1 个点必然已经处理过了,所以选择区间的 l ≤ i l\le i l≤i 就好了。
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;
}