結果

問題 No.2169 To Arithmetic
コンテスト
ユーザー emthrm
提出日時 2022-12-21 03:12:18
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 3,160 bytes
コンパイル時間 2,431 ms
コンパイル使用メモリ 219,700 KB
最終ジャッジ日時 2025-02-09 17:48:36
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 5 WA * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define _USE_MATH_DEFINES
#include <bits/stdc++.h>
using namespace std;
#define FOR(i,m,n) for(int i=(m);i<(n);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(v) (v).begin(),(v).end()
using ll = long long;
constexpr int INF = 0x3f3f3f3f;
constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL;
constexpr double EPS = 1e-8;
constexpr int MOD = 998244353;
// constexpr int MOD = 1000000007;
constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1};
constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1};
constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1};
template <typename T, typename U>
inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; }
template <typename T, typename U>
inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; }
struct IOSetup {
  IOSetup() {
    std::cin.tie(nullptr);
    std::ios_base::sync_with_stdio(false);
    std::cout << fixed << setprecision(20);
  }
} iosetup;

vector<ll> solve(vector<ll> a, const vector<int>& ds) {
  const int m = ds.size();
  if (m == 0) return {};
  const int n = a.size();
  for (int i = n - 1; i >= 0; --i) {
    a[i] -= a.front();
  }

  vector<ll> num(m + 1, 0), bit(m + 1, 0);
  FOR(i, 1, n) {
    const ll d = a[i] - a[i - 1];
    const int l = distance(ds.begin(), lower_bound(ALL(ds), d));
    ++num[l];
    bit[l] += d;
  }
  REP(i, m) num[i + 1] += num[i];
  REP(i, m) bit[i + 1] += bit[i];

  vector<ll> ans(m, 0);
  REP(i, m) ans[i] = 1LL * num[i] * ds[i] - bit[i];
  vector<int> inc{0};
  FOR(i, 1, n) {
    if (a[inc.back()] < a[i]) inc.emplace_back(i);
  }
  map<int, vector<pair<int, int>>> mp;
  const auto add = [&](const int l_ind, const int r_ind) -> void {
    const int l = inc[l_ind], r = inc[r_ind];
    mp[(a[r] - a[l]) % (r - l) == 0 ? (a[r] - a[l]) / (r - l) + 1 : (a[r] - a[l] + r - l - 1) / (r - l)].emplace_back(l_ind, r_ind);
  };
  FOR(i, 1, inc.size()) add(i - 1, i);
  set<int> indices;
  REP(i, inc.size()) indices.emplace(i);
  REP(i, m) {
    while (!mp.empty() && mp.begin()->first <= ds[i]) {
      const vector<pair<int, int>> record = mp.begin()->second;
      mp.erase(mp.begin());
      for (const auto& [l_ind, r_ind] : record) {
        const auto it = indices.find(l_ind);
        if (it == indices.end()) continue;
        auto it2 = next(it);
        if (it2 == indices.end() || *it2 != r_ind) continue;
        it2 = indices.erase(it2);
        if (it2 != indices.end()) add(*it, *it2);
      }
    }
    ans[i] += a[*indices.rbegin()] - 1LL * ds[i] * *indices.rbegin();
  }
  return ans;
}

int main() {
  int n, q; cin >> n >> q;
  vector<ll> a(n); REP(i, n) cin >> a[i];
  vector<int> d(q); REP(i, q) cin >> d[i];
  vector<int> neg, pos;
  for (const int d_i : set<int>(ALL(d))) {
    if (d_i < 0) {
      neg.emplace_back(-d_i);
    } else {
      pos.emplace_back(d_i);
    }
  }
  reverse(ALL(neg));
  const vector<ll> ans_pos = solve(a, pos);
  reverse(ALL(a));
  const vector<ll> ans_neg = solve(a, neg);
  REP(i, q) {
    if (d[i] < 0) {
      cout << ans_neg[distance(neg.begin(), lower_bound(ALL(neg), -d[i]))] << '\n';
    } else {
      cout << ans_pos[distance(pos.begin(), lower_bound(ALL(pos), d[i]))] << '\n';
    }
  }
  return 0;
}
0