結果

問題 No.2643 Many Range Sums Problems
ユーザー SSRS
提出日時 2024-02-19 23:13:12
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,597 bytes
コンパイル時間 2,415 ms
コンパイル使用メモリ 210,824 KB
最終ジャッジ日時 2025-02-19 17:39:15
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 27 WA * 7
権限があれば一括ダウンロードができます

ソースコード

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];
};
auto dist = [&](int u, int v){
return s[u] - s[v];
};
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 (dist(i, i + 1) < 0 || dist(i, 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;
}
}
} else {
int u = i, v = i + 1;
while (u != v){
if (d[u] > d[v]){
mn[u] = max(mn[u], -s[i] + s[i + 1] + X[u]);
mx[u] = min(mx[u], -s[i] + s[i + 1] + X[u] + K);
u = R[u];
} else {
mn[v] = max(mn[v], s[i] - s[i + 1] + X[v] - K);
mx[v] = min(mx[v], s[i] - s[i + 1] + X[v]);
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;
}
}
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0