結果

問題 No.2643 Many Range Sums Problems
コンテスト
ユーザー SSRS
提出日時 2024-02-19 23:15:04
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 111 ms / 8,000 ms
コード長 2,472 bytes
コンパイル時間 2,367 ms
コンパイル使用メモリ 211,720 KB
最終ジャッジ日時 2025-02-19 17:40:10
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 34
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
const int LOG = 13;
const long long INF = 1000000000000000000;
int main(){
  int N, K;
  cin >> N >> K;
  vector<int> R(N);
  vector<long long> X(N);
  for (int i = 0; i < N; i++){
    cin >> R[i] >> X[i];
  }
  vector<vector<pair<long long, int>>> E(N + 1);
  for (int i = 0; i < N; i++){
    E[i].push_back(make_pair(X[i], R[i]));
    E[R[i]].push_back(make_pair(-X[i], i));
  }
  vector<int> d(N + 1, 0);
  vector<long long> s(N + 1, 0);
  for (int i = N - 1; i >= 0; i--){
    d[i] = d[R[i]] + 1;
    s[i] = s[R[i]] + X[i];
  }
  vector<vector<int>> pp(LOG, vector<int>(N + 1, -1));
  vector<vector<long long>> sum(LOG, vector<long long>(N + 1, 0));
  for (int i = 0; i < N; i++){
    pp[0][i] = R[i];
    sum[0][i] = X[i];
  }
  for (int i = 0; i < LOG - 1; i++){
    for (int j = 0; j <= N; j++){
      if (pp[i][j] != -1){
        pp[i + 1][j] = pp[i][pp[i][j]];
        sum[i + 1][j] = sum[i][j] + sum[i][pp[i][j]];
      }
    }
  }
  auto lca = [&](int u, int v){
    if (d[u] > d[v]){
      swap(u, v);
    }
    for (int i = LOG - 1; i >= 0; i--){
      if (((d[v] - d[u]) >> i & 1) == 1){
        v = pp[i][v];
      }
    }
    if (u == v){
      return u;
    }
    for (int i = LOG - 1; i >= 0; i--){
      if (pp[i][u] != pp[i][v]){
        u = pp[i][u];
        v = pp[i][v];
      }
    }
    return R[u];
  };
  vector<bool> ok(N, true);
  vector<long long> mn(N, -INF), mx(N, INF);
  for (int i = 0; i < N; i++){
    int w = lca(i, i + 1);
    if (s[i] - s[i + 1] < 0 || s[i] - s[i + 1] > K){
      vector<bool> is_on_path(N, false);
      int u = i, v = i + 1;
      while (u != v){
        if (d[u] > d[v]){
          is_on_path[u] = true;
          u = R[u];
        } else {
          is_on_path[v] = true;
          v = R[v];
        }
      }
      for (int j = 0; j < N; j++){
        if (!is_on_path[j]){
          ok[j] = false;
        }
      }
    }
    int u = i, v = i + 1;
    while (u != v){
      if (d[u] > d[v]){
        mn[u] = max(mn[u], -s[i] + s[i + 1]);
        mx[u] = min(mx[u], -s[i] + s[i + 1] + K);
        u = R[u];
      } else {
        mn[v] = max(mn[v], s[i] - s[i + 1] - K);
        mx[v] = min(mx[v], s[i] - s[i + 1]);
        v = R[v];
      }
    }
  }
  for (int i = 0; i < N; i++){
    if (mn[i] > mx[i]){
      ok[i] = false;
    }
  }
  for (int i = 0; i < N; i++){
    if (ok[i]){
      cout << "Yes" << endl;
    } else {
      cout << "No" << endl;
    }
  }
}
0