[cf505E]Mr. Kitayuta vs. Bamboos

二分答案，设$s_{i,j}$表示第$i$天对竹子$j$的操作次数，$h_{i,j}$表示第$i$天结束时竹子$j$的高度，则$h_{i,j}=\max(h_{i-1,j}-ps_{i,j},0)+a_{j}$，合法当且仅当$h_{0,i}=h_{i}$且$h_{m,i}\le ans$

令$h'_{i,j}$表示若$h_{i,j}=h'_{i,j}$，转移后满足$h'_{m,j}\le ans$中最大的，转移：若$h'_{i,j}\ge a_{j}$，则$h'_{i-1,j}=ps_{i,j}+h'_{i,j}-a_{j}$，否则无解

合法当且仅当$h'_{m,i}=ans$、$h'_{0,i}$有解且$h'_{0,i}\ge h_{i}$，有解的必要条件为$h'_{i,j}\ge a_{j}$，因此转移就不需要特判

 这样就可以贪心，初始$s_{i,j}=0$，每一次选择在当前的状态下$\max_{h'_{i,j}<a_{j}}i$最大的$j$并令$s_{now,j}+=1$，用堆维护即可，时间复杂度为$o((mk+n)\log_{2}n\log_{2}ans)$

#include<bits/stdc++.h>
using namespace std;
#define N 100005
#define ll long long
priority_queue<pair<ll,int> >q;
int n,m,k,p,h[N],a[N];
ll f[N];
ll calc(int k){
    if (f[k]/a[k]<m)return f[k]/a[k];
    if (f[k]-1LL*m*a[k]<h[k])return m;
    return m+1;
}
bool pd(ll s){
    for(int i=1;i<=n;i++)f[i]=s;
    while (!q.empty())q.pop();
    for(int i=1;i<=n;i++)q.push(make_pair(-calc(i),i));
    for(int i=1;i<=m;i++){
        for(int j=1;j<=k;j++){
            int x=q.top().second;
            q.pop();
            f[x]+=p;
            q.push(make_pair(-calc(x),x));
        }
        if (-q.top().first<=i)return 0;
    }
    return (-q.top().first>m);
}
int main(){
    scanf("%d%d%d%d",&n,&m,&k,&p);
    for(int i=1;i<=n;i++)scanf("%d%d",&h[i],&a[i]);
    ll l=0,r=1e13;
    for(int i=1;i<=n;i++)l=max(l,1LL*a[i]);
    while (l<r){
        ll mid=(l+r>>1);
        if (pd(mid))r=mid;
        else l=mid+1;
    }
    printf("%lld",l);
}