結果
| 問題 |
No.961 Vibrant Fillumination
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2019-12-08 23:29:58 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 2,362 ms / 5,000 ms |
| コード長 | 3,959 bytes |
| コンパイル時間 | 2,704 ms |
| コンパイル使用メモリ | 220,852 KB |
| 最終ジャッジ日時 | 2025-01-08 09:52:13 |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 48 |
ソースコード
#include <bits/stdc++.h>
using ull = unsigned long long;
using usize = std::size_t;
int main()
{
std::cin.tie(nullptr), std::ios::sync_with_stdio(false);
usize N;
std::cin >> N;
std::vector<ull> h(N + 1, 0ULL); // ハッシュ値
for (usize i = 1; i <= N; i++) { std::cin >> h[i]; }
constexpr ull X = 200000; // eps用(i/Xを足してdistinctにする)
std::vector<ull> eps(X + 1, 0); // eps用(i/Xを足してdistinctにする)
std::vector<ull> A(N + 1, 0), B(N + 1, 0), C(N + 1, 0), D(N + 1, 0); // フィルム座標(座圧前)
std::vector<ull> e(N + 1, 0); // 色
std::vector<ull> unzip{0LL};
std::map<ull, usize> zip;
for (usize i = 1; i <= N; i++) {
std::cin >> A[i] >> B[i] >> C[i] >> D[i] >> e[i];
A[i]++, B[i]++, C[i]++, D[i]++;
A[i] = A[i] * X + (eps[A[i]]++);
B[i] = B[i] * X + (eps[B[i]]++);
C[i] = C[i] * X + (eps[C[i]]++);
D[i] = D[i] * X + (eps[D[i]]++);
unzip.push_back(A[i]), unzip.push_back(B[i]), unzip.push_back(C[i]), unzip.push_back(D[i]);
}
usize QN;
std::cin >> QN;
std::vector<usize> P(QN), Q(QN); // クエリ座標(座圧前)
for (usize i = 0; i < QN; i++) {
std::cin >> P[i] >> Q[i];
P[i]++, Q[i]++;
P[i] = P[i] * X + (eps[P[i]]++);
Q[i] = Q[i] * X + (eps[Q[i]]++);
unzip.push_back(P[i]), unzip.push_back(Q[i]);
}
std::sort(unzip.begin(), unzip.end());
unzip.erase(std::unique(unzip.begin(), unzip.end()), unzip.end());
const usize L = unzip.size();
for (usize i = 0; i < L; i++) { zip[unzip[i]] = i; }
std::vector<usize> a(N + 1, 0), b(N + 1, 0), c(N + 1, L), d(N + 1, L);
std::vector<int> lr(L + 1, 0), du(L + 1, 0); // 左右矢印,上下矢印
for (usize i = 1; i <= N; i++) {
a[i] = zip[A[i]], b[i] = zip[B[i]], c[i] = zip[C[i]], d[i] = zip[D[i]];
lr[a[i]] = du[b[i]] = int(i), lr[c[i]] = du[d[i]] = -int(i);
}
std::vector<usize> p(QN), q(QN);
for (usize i = 0; i < QN; i++) { p[i] = zip[P[i]], q[i] = zip[Q[i]]; }
constexpr usize BS = 800;
std::vector<usize> inds(QN);
std::iota(inds.begin(), inds.end(), 0UL);
std::sort(inds.begin(), inds.end(), [&](const usize i, const usize j) {
const usize bi = p[i] / BS, bj = p[j] / BS;
return bi != bj ? bi < bj : q[i] != q[j] ? q[i] < q[j] : p[i] < p[j];
});
std::vector<usize> color_cnt(N + 1, 0);
ull hash = 0;
auto add_film = [&](const usize i) { // フィルムiを追加
const usize c = e[i];
if (color_cnt[c]++ == 0) { hash ^= h[c]; }
};
auto del_film = [&](const usize i) { // フィルムiを削除
const usize c = e[i];
if (--color_cnt[c] == 0) { hash ^= h[c]; }
};
auto add_del_film = [&](const int i) {
if (i > 0) { add_film(usize(i)); }
if (i < 0) { del_film(usize(-i)); }
};
std::vector<ull> answer(QN);
usize x = 0, y = 0;
auto von = [&](const usize x, const usize y) { return b[(usize)std::abs(lr[x])] <= y and y < d[(usize)std::abs(lr[x])]; };
auto hon = [&](const usize x, const usize y) { return a[(usize)std::abs(du[y])] <= x and x < c[(usize)std::abs(du[y])]; };
for (const usize i : inds) {
const usize nx = p[i], ny = q[i]; // 次のクエリ点
for (; x < nx; x++) {
if (von(x + 1, y)) { add_del_film(lr[x + 1]); }
}
for (; x > nx; x--) {
if (von(x, y)) { add_del_film(-lr[x]); }
}
for (; y < ny; y++) {
if (hon(x, y + 1)) { add_del_film(du[y + 1]); }
}
for (; y > ny; y--) {
if (hon(x, y)) { add_del_film(-du[y]); }
}
answer[i] = hash;
}
for (usize i = 0; i < QN; i++) { std::cout << answer[i] << "\n"; }
return 0;
}