正题
很容易发现这个答案具有单调性,所以选择花费1个log的代价取二分。
顺着考虑发现很难,我们倒着考虑,二分到的值为x。
我们就假设一开始全部的高度都为x,每天每根竹子缩短,每天的末尾可以选择k根竹子把它拉高p。
问存不存在一种方案使得每根竹子在任意时刻高度不为负数,而且最后高度,利用放缩法容易证明这是等价的。
接下来的问题就很好办了,假设一根竹子高度为A,那么他最多能坚持的敲击,用一个堆来维护能坚持最小敲击的竹子是哪一根,假设它坚持不到第i天,显然不可行,否则就把它拉高,然后重新计算它能坚持多少次敲击。
假设它可以满足当前的高度已经可以满足最后的条件了,就直接把它弹出,这样只要看看堆里面是否有元素就可以了。
这样贪心正确是正确的,首先给可以坚持最小敲击的竹子拔高是正确的,因为不可能去拔高一个可以坚持更多天数的竹子,其次,如果满足天数的条件,一根竹子被拔高的次数也就确定下来了。
#include<bits/stdc++.h>
using namespace std;
const int N=100010;
int n,m,k;
struct node{
int x;
long long d;
bool operator<(const node q)const{
return d>q.d;
}
}X;
priority_queue<node> q;
long long p;
long long h[N],a[N];
int f[N];
bool check(long long x){
memset(f,0,sizeof(f));
while(!q.empty()) q.pop();
for(int i=1;i<=n;i++) if(x-a[i]*m<h[i]) q.push((node){i,x/a[i]});
for(int i=1;i<=m;i++)
for(int j=1;j<=k;j++) if(!q.empty()){
X=q.top();q.pop();
if(X.d<i) return false;
f[X.x]++;
if(x-a[X.x]*m+f[X.x]*p<h[X.x]) q.push((node){X.x,(x+f[X.x]*p)/a[X.x]});
}
return q.empty();
}
int main(){
scanf("%d %d %d %lld",&n,&m,&k,&p);
for(int i=1;i<=n;i++) scanf("%lld %lld",&h[i],&a[i]);
long long l=0,r=1e13,ans=0;
while(l<=r){
long long mid=(l+r)/2;
if(check(mid)) r=(ans=mid)-1;
else l=mid+1;
}
printf("%lld\n",ans);
}