結果

問題 No.1240 Or Sum of Xor Pair
ユーザー carrot46carrot46
提出日時 2020-09-28 11:30:09
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 496 ms / 2,000 ms
コード長 1,906 bytes
コンパイル時間 1,965 ms
コンパイル使用メモリ 171,672 KB
実行使用メモリ 15,212 KB
最終ジャッジ日時 2024-07-01 20:21:40
合計ジャッジ時間 17,850 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 399 ms
13,524 KB
testcase_01 AC 400 ms
13,652 KB
testcase_02 AC 397 ms
13,524 KB
testcase_03 AC 406 ms
13,552 KB
testcase_04 AC 403 ms
13,552 KB
testcase_05 AC 403 ms
13,548 KB
testcase_06 AC 404 ms
13,548 KB
testcase_07 AC 405 ms
13,544 KB
testcase_08 AC 403 ms
13,540 KB
testcase_09 AC 400 ms
13,544 KB
testcase_10 AC 431 ms
13,680 KB
testcase_11 AC 411 ms
13,840 KB
testcase_12 AC 416 ms
13,832 KB
testcase_13 AC 419 ms
13,784 KB
testcase_14 AC 421 ms
13,836 KB
testcase_15 AC 484 ms
15,084 KB
testcase_16 AC 476 ms
15,084 KB
testcase_17 AC 483 ms
14,960 KB
testcase_18 AC 478 ms
15,084 KB
testcase_19 AC 474 ms
15,080 KB
testcase_20 AC 474 ms
15,084 KB
testcase_21 AC 481 ms
15,080 KB
testcase_22 AC 475 ms
15,084 KB
testcase_23 AC 496 ms
15,212 KB
testcase_24 AC 477 ms
15,084 KB
testcase_25 AC 476 ms
15,084 KB
testcase_26 AC 470 ms
15,084 KB
testcase_27 AC 399 ms
13,516 KB
testcase_28 AC 399 ms
13,520 KB
testcase_29 AC 467 ms
15,080 KB
testcase_30 AC 429 ms
14,304 KB
testcase_31 AC 431 ms
14,304 KB
testcase_32 AC 458 ms
15,084 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
//#include <chrono>
//#pragma GCC optimize("Ofast")
using namespace std;
#define reps(i,s,n) for(int i = s; i < n; i++)
#define rep(i,n) reps(i,0,n)
#define Rreps(i,n,e) for(int i = n - 1; i >= e; --i)
#define Rrep(i,n) Rreps(i,n,0)
#define ALL(a) a.begin(), a.end()
#define fi first
#define se second

using ll = long long;
using vec = vector<ll>;
using mat = vector<vec>;

ll N,M,H,W,Q,K,A,B;
string S;
typedef pair<ll, ll> P;
const ll INF = (1LL<<58);

template <class T> void adamard(vector<T> &a, bool inv){
    int n = a.size();
    for(int i = 1; i < n; i <<= 1){
        int m = n - i;
        for(int j = 0; j < m; ++j){
            j += j&i;
            T temp = a[j|i];
            a[j|i] = a[j] - temp;
            a[j] += temp;
            if(inv){
                a[j] /= 2;
                a[j|i] /= 2;
            }
        }
    }
}

template <class T> vector<T> xor_convolution(vector<T> &a, vector<T> &b){
    //can_divide ? " / 2" : " * 2.inv()"
    int n(1);
    while(n < max((int)a.size(), (int)b.size())) n <<= 1;
    vector<T> a_cpy(n), b_cpy(n);
    copy(ALL(a), a_cpy.begin());
    copy(ALL(b), b_cpy.begin());
    adamard(a_cpy, false);
    adamard(b_cpy, false);
    rep(i, n) a_cpy[i] *= b_cpy[i];
    adamard(a_cpy, true);
    return a_cpy;
}

int main() {
    const int max_A = 1<<18;
    cin>>N>>K;
    vec a(N), cnt(max_A, 0), cnt_cnv;
    rep(i, N) {cin>>a[i]; ++cnt[a[i]];}
    cnt_cnv = xor_convolution(cnt, cnt);
    ll res(-accumulate(ALL(a), 0LL));
    for(int n = 1; n < max_A; n <<= 1){
        vec cnt_without_bitk(max_A, 0);
        int r = max_A - n;
        rep(i, r){
            i += i & n;
            cnt_without_bitk[i] = cnt[i];
        }
        vec cnt_without_bitk_cnv = xor_convolution(cnt_without_bitk, cnt_without_bitk);
        rep(i, K) res += (cnt_cnv[i] - cnt_without_bitk_cnv[i]) * n;
    }
    cout<<res / 2<<endl;
}
0