結果

問題 No.2977 Kth Xor Pair
ユーザー UrtuseaUrtusea
提出日時 2024-12-19 13:12:56
言語 C++23(gcc13)
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 2,077 ms / 3,000 ms
コード長 2,845 bytes
コンパイル時間 5,692 ms
コンパイル使用メモリ 275,684 KB
実行使用メモリ 37,196 KB
最終ジャッジ日時 2024-12-19 13:13:49
合計ジャッジ時間 43,108 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 3 ms
6,816 KB
testcase_02 AC 4 ms
6,820 KB
testcase_03 AC 4 ms
6,816 KB
testcase_04 AC 4 ms
6,820 KB
testcase_05 AC 4 ms
6,816 KB
testcase_06 AC 4 ms
6,820 KB
testcase_07 AC 793 ms
22,412 KB
testcase_08 AC 1,925 ms
34,172 KB
testcase_09 AC 1,602 ms
32,996 KB
testcase_10 AC 1,640 ms
33,304 KB
testcase_11 AC 1,748 ms
34,136 KB
testcase_12 AC 1,532 ms
34,176 KB
testcase_13 AC 1,720 ms
35,076 KB
testcase_14 AC 853 ms
22,460 KB
testcase_15 AC 1,662 ms
32,416 KB
testcase_16 AC 2,077 ms
35,884 KB
testcase_17 AC 1,726 ms
35,732 KB
testcase_18 AC 1,759 ms
36,808 KB
testcase_19 AC 1,797 ms
35,804 KB
testcase_20 AC 1,677 ms
36,140 KB
testcase_21 AC 1,709 ms
36,156 KB
testcase_22 AC 1,650 ms
36,208 KB
testcase_23 AC 1,591 ms
37,196 KB
testcase_24 AC 1,633 ms
35,864 KB
testcase_25 AC 1,868 ms
36,172 KB
testcase_26 AC 1,646 ms
35,748 KB
testcase_27 AC 499 ms
6,820 KB
testcase_28 AC 337 ms
6,816 KB
testcase_29 AC 351 ms
6,820 KB
testcase_30 AC 513 ms
6,816 KB
testcase_31 AC 431 ms
6,816 KB
testcase_32 AC 275 ms
6,816 KB
testcase_33 AC 293 ms
6,820 KB
testcase_34 AC 269 ms
6,816 KB
testcase_35 AC 280 ms
6,816 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#define int  int32_t
#define Int  int64_t
#define uint uint32_t
#define uInt uint64_t

template <int N, int B, int Build_Mode>
struct Bit_Trie {
    using T = std::conditional_t<(B <= 31), uint, uInt>;
    int idx;
    int cnt[N * B];
    int nxt[N * B][2];

    Bit_Trie() = default;

    void clear() {
        std::memset(cnt, 0, sizeof(cnt[0]) * (idx + 1));
        std::memset(nxt, 0, sizeof(nxt[0]) * (idx + 1));
        idx = 0;
    }

    int get_next(int &p) {
        return p ? p : p = ++idx;
    }

    void insert(const T x, int u = 0) {
        if constexpr (Build_Mode == 1) {
            for (int i = B - 1; i >= 0; i--)
                cnt[u = get_next(nxt[u][x >> i & 1])]++;
        } else {
            for (int i = 0; i < B; i++)
                cnt[u = get_next(nxt[u][x >> i & 1])]++;
        }
    }

    void erase(const T x, int u = 0) {
        if constexpr (Build_Mode == 1) {
            for (int i = B - 1; i >= 0; i--)
                cnt[u = get_next(nxt[u][x >> i & 1])]--;
        } else {
            for (int i = 0; i < B; i++)
                cnt[u = get_next(nxt[u][x >> i & 1])]--;
        }
    }

    int count(const T x, int u = 0) {
        if constexpr (Build_Mode == 1) {
            for (int i = B - 1; i >= 0; i--)
                if ((u = nxt[u][x >> i & 1]) == 0) return 0;
        } else {
            for (int i = 0; i < B; i++)
                if ((u = nxt[u][x >> i & 1]) == 0) return 0;
        }
        return cnt[u];
    }
};

Bit_Trie<200000, 31, 1> bit;

void init() {
    
}

void code() {
    Int n, k;
    std::cin >> n >> k;

    std::vector<int> a(n + 1);
    for (int i = 1; i <= n; i++) {
        std::cin >> a[i];
        bit.insert(a[i]);
    }

    const auto query = [](int x, int lim) -> Int {
        Int res = 0;
        auto cnt = bit.cnt;
        auto nxt = bit.nxt;
        
        int u = 0;
        for (int i = 30; i >= 0; i--) {
            int v1 = x >> i & 1;
            int v2 = lim >> i & 1;
            int v = v1 ^ v2;

            if (v2 == 1)
                res += cnt[nxt[u][v ^ 1]];
            if ((u = nxt[u][v]) == 0)
                break;
        }

        return res + cnt[u];
    };

    const auto check = [n, k, &query, &a](Int x) -> bool {
        Int cnt = 0;
        for (int i = 1; i <= n; i++) {
            cnt += query(a[i], x);
            cnt -= (a[i] ^ a[i]) <= x;
        }
        return cnt / 2 >= k;
    };

    Int l = -1, r = (1LL << 30) - 1;
    while (l + 1 < r) {
        int mid = (l + r) / 2;
        if (check(mid))
            r = mid;
        else
            l = mid;
    }

    std::cout << r << '\n';
}

int main() {
    std::cin.tie(nullptr)->sync_with_stdio(false);

    // for (int i = 1, n = (std::cin >> n, n); i <= n; i++)
    {
        init();
        code();
    }

    return 0;
}
0