結果

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

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,820 KB
testcase_02 AC 4 ms
6,816 KB
testcase_03 AC 4 ms
6,816 KB
testcase_04 AC 4 ms
6,820 KB
testcase_05 AC 4 ms
6,820 KB
testcase_06 AC 3 ms
6,816 KB
testcase_07 AC 974 ms
22,460 KB
testcase_08 AC 2,273 ms
35,620 KB
testcase_09 AC 2,001 ms
32,628 KB
testcase_10 AC 1,923 ms
33,248 KB
testcase_11 AC 2,010 ms
33,860 KB
testcase_12 AC 1,823 ms
33,844 KB
testcase_13 AC 2,045 ms
35,012 KB
testcase_14 AC 869 ms
22,348 KB
testcase_15 AC 1,712 ms
30,864 KB
testcase_16 AC 2,213 ms
36,960 KB
testcase_17 AC 2,170 ms
35,504 KB
testcase_18 AC 2,231 ms
35,824 KB
testcase_19 AC 2,255 ms
36,560 KB
testcase_20 AC 2,156 ms
36,520 KB
testcase_21 AC 2,227 ms
36,732 KB
testcase_22 AC 2,159 ms
35,960 KB
testcase_23 AC 2,098 ms
36,904 KB
testcase_24 AC 2,117 ms
36,368 KB
testcase_25 AC 2,437 ms
36,840 KB
testcase_26 AC 2,036 ms
37,204 KB
testcase_27 AC 473 ms
7,260 KB
testcase_28 AC 316 ms
6,820 KB
testcase_29 AC 344 ms
6,816 KB
testcase_30 AC 461 ms
6,816 KB
testcase_31 AC 455 ms
6,820 KB
testcase_32 AC 276 ms
6,820 KB
testcase_33 AC 280 ms
6,816 KB
testcase_34 AC 281 ms
6,820 KB
testcase_35 AC 312 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