結果

問題 No.961 Vibrant Fillumination
ユーザー PachicobuePachicobue
提出日時 2019-12-09 02:59:52
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,550 ms / 5,000 ms
コード長 4,002 bytes
コンパイル時間 3,089 ms
コンパイル使用メモリ 225,688 KB
実行使用メモリ 37,476 KB
最終ジャッジ日時 2023-09-01 22:43:25
合計ジャッジ時間 71,562 ms
ジャッジサーバーID
（参考情報） 		judge16 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
5,360 KB
testcase_01 AC 4 ms
5,364 KB
testcase_02 AC 1,333 ms
35,936 KB
testcase_03 AC 1,330 ms
36,128 KB
testcase_04 AC 1,382 ms
36,148 KB
testcase_05 AC 1,406 ms
36,692 KB
testcase_06 AC 1,326 ms
36,464 KB
testcase_07 AC 1,321 ms
35,944 KB
testcase_08 AC 1,326 ms
36,728 KB
testcase_09 AC 1,335 ms
36,076 KB
testcase_10 AC 1,470 ms
36,024 KB
testcase_11 AC 1,448 ms
36,432 KB
testcase_12 AC 1,438 ms
35,808 KB
testcase_13 AC 1,443 ms
37,088 KB
testcase_14 AC 1,489 ms
36,084 KB
testcase_15 AC 1,487 ms
36,212 KB
testcase_16 AC 1,523 ms
36,580 KB
testcase_17 AC 1,486 ms
37,152 KB
testcase_18 AC 1,507 ms
35,996 KB
testcase_19 AC 1,454 ms
35,732 KB
testcase_20 AC 1,550 ms
36,432 KB
testcase_21 AC 1,482 ms
36,376 KB
testcase_22 AC 1,452 ms
36,048 KB
testcase_23 AC 1,426 ms
35,992 KB
testcase_24 AC 1,467 ms
35,972 KB
testcase_25 AC 1,537 ms
36,248 KB
testcase_26 AC 1,213 ms
36,432 KB
testcase_27 AC 816 ms
35,988 KB
testcase_28 AC 1,188 ms
36,564 KB
testcase_29 AC 789 ms
36,432 KB
testcase_30 AC 589 ms
36,052 KB
testcase_31 AC 378 ms
36,164 KB
testcase_32 AC 492 ms
36,456 KB
testcase_33 AC 4 ms
5,416 KB
testcase_34 AC 1,361 ms
35,964 KB
testcase_35 AC 1,332 ms
35,972 KB
testcase_36 AC 1,298 ms
37,212 KB
testcase_37 AC 1,321 ms
35,760 KB
testcase_38 AC 1,298 ms
36,040 KB
testcase_39 AC 1,369 ms
36,940 KB
testcase_40 AC 1,336 ms
36,368 KB
testcase_41 AC 1,389 ms
36,040 KB
testcase_42 AC 1,347 ms
36,012 KB
testcase_43 AC 1,387 ms
36,276 KB
testcase_44 AC 1,367 ms
36,212 KB
testcase_45 AC 1,397 ms
36,644 KB
testcase_46 AC 1,379 ms
36,440 KB
testcase_47 AC 1,391 ms
37,476 KB
testcase_48 AC 1,370 ms
36,048 KB
testcase_49 AC 1,362 ms
37,156 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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 = 300000;                                             // 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;
    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), e(N + 1, 0);
    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]];
        e[i]     = E[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])^(bi%2) : 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;
}
0