結果

問題 No.2627 Unnatural Pitch
ユーザー zeta7532zeta7532
提出日時 2024-01-15 02:02:55
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 121 ms / 4,000 ms
コード長 1,552 bytes
コンパイル時間 2,219 ms
コンパイル使用メモリ 211,112 KB
実行使用メモリ 7,936 KB
最終ジャッジ日時 2024-09-28 02:10:05
合計ジャッジ時間 5,518 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
using namespace std;
using ll = long long;
const ll mod = 998244353;
#define fi first
#define se second
#define rep(i,n) for(ll i=0;i<n;i++)
#define all(x) x.begin(),x.end()
#define faster ios::sync_with_stdio(false);cin.tie(nullptr)

ll f(ll x,ll y){
    if(x<=0) return 0;
    return (x+y-1)/y;
}

ll g(ll x,ll y,ll z){
    ll ok=(z-x)/y-100,ng=(z-x)/y+100;
    while(ng-ok>1){
        ll mid=(ok+ng)/2;
        if(x+mid*y<z) ok=mid;
        else ng=mid;
    }
    return x+ng*y;
}

int main() {
    ll N,K,L,U;
    cin >> N >> K >> L >> U;
    vector<ll> A(N);
    rep(i,N) cin >> A[i];
    ll ok=-1,ng=1e12+500;
    while(ng-ok>1){
        ll x=(ok+ng)/2;
        ll a=0,b=0;
        rep(i,N){
            if(A[i]<x) a++;
            if(x+U-L<A[i]) b++;
        }
        if(a<=b) ok=x;
        else ng=x;
    }
    ll now=0;
    ll l=ok-K+1;
    ll u=l+U-L;
    rep(i,N){
        now+=max((ll)0,f(l-A[i],K));
        now+=max((ll)0,f(A[i]-u,K));
    }
    ll ans=now;
    vector<ll> cum(K*2-1,0);
    rep(i,N){
        ll xl=(A[i]+1)%K;
        ll minl=g(xl,K,max(A[i],ok-K+2));
        while(minl<=ok+K) cum[minl-(ok-K+2)]++,minl+=K;
        ll xu=(A[i]+L+K-U%K)%K;
        ll minu=g(xu,K,ok-K+2);
        while(minu<=ok+K){
            if(minu+U-L<=A[i]) cum[minu-(ok-K+2)]--;
            minu+=K;
        }    
    }
    rep(i,K*2-1){
        now+=cum[i];
        ans=min(ans,now);
    }
    cout << ans << endl;
    return 0;
}
0