結果

問題 No.2627 Unnatural Pitch
コンテスト
ユーザー hos.lyric
提出日時 2024-02-09 21:57:36
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 200 ms / 4,000 ms
コード長 3,009 bytes
コンパイル時間 1,110 ms
コンパイル使用メモリ 116,148 KB
実行使用メモリ 8,192 KB
最終ジャッジ日時 2024-09-28 15:14:02
合計ジャッジ時間 3,887 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cassert>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <functional>
#include <iostream>
#include <limits>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>

using namespace std;

using Int = long long;

template <class T1, class T2> ostream &operator<<(ostream &os, const pair<T1, T2> &a) { return os << "(" << a.first << ", " << a.second << ")"; };
template <class T> ostream &operator<<(ostream &os, const vector<T> &as) { const int sz = as.size(); os << "["; for (int i = 0; i < sz; ++i) { if (i >= 256) { os << ", ..."; break; } if (i > 0) { os << ", "; } os << as[i]; } return os << "]"; }
template <class T> void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; }
template <class T> bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; }
template <class T> bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; }
#define COLOR(s) ("\x1b[" s "m")


// floor(a / b)
Int divFloor(Int a, Int b) {
  return a / b - (((a ^ b) < 0 && a % b != 0) ? 1 : 0);
}

// ceil(a / b)
Int divCeil(Int a, Int b) {
  return a / b + (((a ^ b) > 0 && a % b != 0) ? 1 : 0);
}


constexpr Int INF = 100100100100100100LL;

int N;
Int K, L, U;
vector<Int> A;
Int W;

// [b, b + W]
Int calc(Int b) {
  const Int c = b + W;
  Int ret = 0;
  for (int i = 0; i < N; ++i) {
    if (A[i] < b) {
      ret += ((b - A[i]) + K - 1) / K;
    }
    if (c < A[i]) {
      ret += ((A[i] - c) + K - 1) / K;
    }
  }
  return ret;
}

int main() {
  for (; ~scanf("%d%lld%lld%lld", &N, &K, &L, &U); ) {
    A.resize(N);
    for (int i = 0; i < N; ++i) {
      scanf("%lld", &A[i]);
    }
    W = U - L;
    sort(A.begin(), A.end());
    
    Int lo = divFloor(A[0] - W, K) - 1;
    Int hi = divCeil(A[N - 1] + W, K) + 1;
    for (; lo + 1 < hi; ) {
      const Int mid = (lo + hi) / 2;
      ((calc(K * mid) > calc(K * (mid + 1))) ? lo : hi) = mid;
    }
    
    const Int bL = K * (hi - 1);
    const Int bR = K * (hi + 1);
cerr<<"W = "<<W<<", [bL, bR] = ["<<bL<<", "<<bR<<"]"<<endl;
    vector<Int> fs(bR - bL + 1, 0);
    fs[0] += calc(bL);
    for (int i = 0; i < N; ++i) {
      for (Int b = A[i] + max(divCeil(bL - A[i], K), 0LL) * K + 1; b <= bR; b += K) {
        ++fs[b - bL];
      }
      for (Int c = A[i] - max(divCeil(A[i] - (bR + W), K), 0LL) * K; c > bL + W; c -= K) {
        --fs[c - W - bL];
      }
    }
    for (Int b = bL + 1; b <= bR; ++b) {
      fs[b - bL] += fs[b - 1 - bL];
    }
    
    const Int ans = *min_element(fs.begin(), fs.end());
    printf("%lld\n", ans);
#ifdef LOCAL
cerr<<"fs = "<<fs<<endl;
// for(Int b=bL;b<=bR;++b)cerr<<calc(b)<<" ";cerr<<endl;
for(Int b=bL;b<=bR;++b)assert(calc(b)==fs[b-bL]);
#endif
  }
  return 0;
}
0