結果

問題 No.1240 Or Sum of Xor Pair
ユーザー MisterMister
提出日時 2020-09-26 14:52:40
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 258 ms / 2,000 ms
コード長 2,024 bytes
コンパイル時間 711 ms
コンパイル使用メモリ 77,304 KB
実行使用メモリ 6,252 KB
最終ジャッジ日時 2024-06-29 02:16:05
合計ジャッジ時間 9,774 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 209 ms
5,340 KB
testcase_01 AC 209 ms
5,468 KB
testcase_02 AC 209 ms
5,376 KB
testcase_03 AC 212 ms
5,376 KB
testcase_04 AC 211 ms
5,492 KB
testcase_05 AC 212 ms
5,484 KB
testcase_06 AC 212 ms
5,376 KB
testcase_07 AC 210 ms
5,488 KB
testcase_08 AC 210 ms
5,396 KB
testcase_09 AC 211 ms
5,488 KB
testcase_10 AC 218 ms
5,552 KB
testcase_11 AC 217 ms
5,444 KB
testcase_12 AC 221 ms
5,476 KB
testcase_13 AC 221 ms
5,476 KB
testcase_14 AC 221 ms
5,496 KB
testcase_15 AC 255 ms
6,056 KB
testcase_16 AC 253 ms
5,996 KB
testcase_17 AC 255 ms
5,996 KB
testcase_18 AC 255 ms
6,000 KB
testcase_19 AC 256 ms
6,096 KB
testcase_20 AC 255 ms
6,000 KB
testcase_21 AC 255 ms
6,128 KB
testcase_22 AC 256 ms
6,000 KB
testcase_23 AC 258 ms
6,132 KB
testcase_24 AC 255 ms
6,252 KB
testcase_25 AC 254 ms
6,000 KB
testcase_26 AC 253 ms
6,128 KB
testcase_27 AC 209 ms
5,376 KB
testcase_28 AC 208 ms
5,468 KB
testcase_29 AC 238 ms
5,996 KB
testcase_30 AC 221 ms
5,736 KB
testcase_31 AC 224 ms
5,732 KB
testcase_32 AC 234 ms
6,124 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <numeric>
#include <vector>

using lint = long long;
constexpr int K = 18;

// 自身とのxor畳み込み
void xor_convolution(std::vector<lint>& xs) {
    for (int i = 1; i < (1 << K); i <<= 1) {
        for (int b = 0; b < (1 << K); ++b) {
            if (b & i) continue;
            auto l = xs[b],
                 r = xs[b | i];
            xs[b] = l + r;
            xs[b | i] = l - r;
        }
    }

    for (auto& x : xs) x *= x;

    for (int i = 1; i < (1 << K); i <<= 1) {
        for (int b = 0; b < (1 << K); ++b) {
            if (b & i) continue;
            auto l = xs[b],
                 r = xs[b | i];
            xs[b] = l + r;
            xs[b | i] = l - r;
        }
    }

    for (auto& x : xs) x >>= K;
}

void solve() {
    int n, a;
    std::cin >> n >> a;

    std::vector<int> xs(n);
    for (auto& x : xs) std::cin >> x;

    lint tot;  // 全体の個数
    {
        int sz = 0;
        std::vector<lint> cnt(1 << K, 0);
        for (auto x : xs) {
            ++cnt[x];
            ++sz;
        }

        xor_convolution(cnt);
        cnt[0] -= sz;                 // 0にi=jなるものが含まれることに注意
        for (auto& c : cnt) c >>= 1;  // i>jを除外
        tot = std::accumulate(cnt.begin(), cnt.begin() + a, 0LL);
    }

    lint ans = 0;
    for (int k = 1; k < (1 << K); k <<= 1) {
        // k桁目が0であるもののみの個数
        lint num;
        {
            int sz = 0;
            std::vector<lint> cnt(1 << K, 0);
            for (auto x : xs) {
                if (x & k) continue;
                ++cnt[x];
                ++sz;
            }

            xor_convolution(cnt);
            cnt[0] -= sz;
            for (auto& c : cnt) c >>= 1;
            num = std::accumulate(cnt.begin(), cnt.begin() + a, 0LL);
        }

        ans += k * (tot - num);
    }

    std::cout << ans << "\n";
}

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

    solve();

    return 0;
}
0