「JXOI2017」加法
LOJ
显然可以想到二分答案。现在考虑如何\(check\)
容易发现,假设考虑完前\(i-1\)个位置,对于位置\(i\),我们肯定需要选择经过它的右端点尽量大的区间。
这就可以用一个堆来维护一下。
然后我们还需要知道每个位置被覆盖了多少次。这就再写颗线段树就行了。
代码:
#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<queue>
#include<set>
#include<bitset>
#include<map>
using namespace std;
#define rg register
void read(int &x){
char ch;bool ok;
for(ok=0,ch=getchar();!isdigit(ch);ch=getchar())if(ch=='-')ok=1;
for(x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());if(ok)x=-x;
}
const int maxn=3e5+10;
int n,c[maxn],T,m,k,a,b[maxn];
int f[maxn*4],la[maxn*4];
struct oo{int l,r;}s[maxn];
priority_queue<oo>q;
bool operator<(oo a,oo b){return a.r<b.r;}
bool cmp(oo a,oo b){return a.l<b.l;}
void update(int x){f[x]=f[x<<1]+f[x<<1|1];}
void pushdown(int x,int l,int r,int mid){
la[x<<1]+=la[x],la[x<<1|1]+=la[x];
f[x<<1]+=la[x]*(mid-l+1);f[x<<1|1]+=la[x]*(r-mid);
la[x]=0;
}
void change(int x,int l,int r,int a,int b){
if(a<=l&&b>=r)return la[x]++,f[x]+=r-l+1,void();
int mid=(l+r)>>1;if(la[x])pushdown(x,l,r,mid);
if(a<=mid)change(x<<1,l,mid,a,b);
if(b>mid)change(x<<1|1,mid+1,r,a,b);
update(x);
}
int get(int x,int l,int r,int a){
if(l==r)return f[x];
int mid=(l+r)>>1;if(la[x])pushdown(x,l,r,mid);
return a<=mid?get(x<<1,l,mid,a):get(x<<1|1,mid+1,r,a);
}
bool check(int now){
int x=1,t=k;while(q.size())q.pop();
for(rg int i=1;i<=n;i++){
while(x<=m&&s[x].l<=i)q.push(s[x]),x++;
while(get(1,1,n,i)*a+c[i]<now){
if(!q.size()||q.top().r<i)return 0;
change(1,1,n,q.top().l,q.top().r),t--;
if(t<0)return 0;q.pop();
}
}
return 1;
}
void build(int x,int l,int r){
f[x]=la[x]=0;if(l==r)return ;
int mid=(l+r)>>1;
build(x<<1,l,mid),build(x<<1|1,mid+1,r);
}
int main(){
read(T);
while(T--){
read(n),read(m),read(k),read(a);
for(rg int i=1;i<=n;i++)read(c[i]);
for(rg int i=1;i<=m;i++)read(s[i].l),read(s[i].r);
sort(s+1,s+m+1,cmp);int l=1,r=1e9;
while(l<=r){
build(1,1,n);int mid=(l+r)>>1;
if(check(mid))l=mid+1;
else r=mid-1;
}
printf("%d\n",r);
}
}