#562205 (C++17) No.1240 Or Sum of Xor Pair

提出ソース

結果

問題	No.1240 Or Sum of Xor Pair
ユーザー	hitonanode
提出日時	2020-10-03 18:15:33
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	444 ms / 2,000 ms
コード長	5,400 bytes
コンパイル時間	2,196 ms
コンパイル使用メモリ	209,232 KB
実行使用メモリ	20,456 KB
最終ジャッジ日時	2024-07-18 04:45:48
合計ジャッジ時間	17,489 ms
ジャッジサーバーID （参考情報）	judge5 / judge1

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	396 ms 19,624 KB
testcase_01	AC	383 ms 19,444 KB
testcase_02	AC	392 ms 19,560 KB
testcase_03	AC	400 ms 19,548 KB
testcase_04	AC	394 ms 19,752 KB
testcase_05	AC	386 ms 19,512 KB
testcase_06	AC	393 ms 19,576 KB
testcase_07	AC	396 ms 19,556 KB
testcase_08	AC	396 ms 19,524 KB
testcase_09	AC	403 ms 19,496 KB
testcase_10	AC	416 ms 19,536 KB
testcase_11	AC	407 ms 19,544 KB
testcase_12	AC	405 ms 19,608 KB
testcase_13	AC	412 ms 19,680 KB
testcase_14	AC	404 ms 19,596 KB
testcase_15	AC	410 ms 20,380 KB
testcase_16	AC	437 ms 20,292 KB
testcase_17	AC	437 ms 20,244 KB
testcase_18	AC	444 ms 20,244 KB
testcase_19	AC	443 ms 20,360 KB
testcase_20	AC	437 ms 20,224 KB
testcase_21	AC	427 ms 20,364 KB
testcase_22	AC	440 ms 20,236 KB
testcase_23	AC	435 ms 20,308 KB
testcase_24	AC	413 ms 20,456 KB
testcase_25	AC	409 ms 20,284 KB
testcase_26	AC	429 ms 20,444 KB
testcase_27	AC	398 ms 19,476 KB
testcase_28	AC	404 ms 19,504 KB
testcase_29	AC	417 ms 20,408 KB
testcase_30	AC	418 ms 20,032 KB
testcase_31	AC	383 ms 19,852 KB
testcase_32	AC	429 ms 20,316 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
using namespace std;
using lint = long long;
using pint = pair<int, int>;
using plint = pair<lint, lint>;
struct fast_ios {
    fast_ios() { cin.tie(nullptr), ios::sync_with_stdio(false), cout << fixed << setprecision(20); };
} fast_ios_;
#define FOR(i, begin, end) for (int i = (begin), i##_end_ = (end); i < i##_end_; i++)
#define IFOR(i, begin, end) for (int i = (end)-1, i##_begin_ = (begin); i >= i##_begin_; i--)
#define REP(i, n) FOR(i, 0, n)
#define IREP(i, n) IFOR(i, 0, n)
#define ALL(x) (x).begin(), (x).end()
//
template <typename T, typename V>
void ndarray(vector<T>& vec, const V& val, int len) { vec.assign(len, val); }
template <typename T, typename V, typename... Args>
void ndarray(vector<T>& vec, const V& val, int len, Args... args)
{
    vec.resize(len), for_each(begin(vec), end(vec), [&](T& v) { ndarray(v, val, args...); });
}
template <typename T>
bool chmax(T& m, const T q) { return m < q ? (m = q, true) : false; }
template <typename T>
bool chmin(T& m, const T q) { return m > q ? (m = q, true) : false; }
template <typename T1, typename T2>
pair<T1, T2> operator+(const pair<T1, T2>& l, const pair<T1, T2>& r) { return make_pair(l.first + r.first, l.second + r.second); }
template <typename T1, typename T2>
pair<T1, T2> operator-(const pair<T1, T2>& l, const pair<T1, T2>& r) { return make_pair(l.first - r.first, l.second - r.second); }
template <typename T>
vector<T> srtunq(vector<T> vec) { return sort(vec.begin(), vec.end()), vec.erase(unique(vec.begin(), vec.end()), vec.end()), vec; }
template <typename T>
istream& operator>>(istream& is, vector<T>& vec)
{
    return for_each(begin(vec), end(vec), [&](T& v) { is >> v; }), is;
}

// output
template <typename T, typename V>
ostream& dmpseq(ostream&, const T&, const string&, const string&, const string&);
#if __cplusplus >= 201703L
template <typename... T>
ostream& operator<<(ostream& os, const tuple<T...>& tpl)
{
    return apply([&os](auto&&... args) { ((os << args << ','), ...); }, tpl), os;
}
#endif
//
template <typename T1, typename T2>
ostream& operator<<(ostream& os, const pair<T1, T2>& p) { return os << '(' << p.first << ',' << p.second << ')'; }
template <typename T>
ostream& operator<<(ostream& os, const vector<T>& x) { return dmpseq<vector<T>, T>(os, x, "[", ",", "]"); }
template <typename T>
ostream& operator<<(ostream& os, const deque<T>& x) { return dmpseq<deque<T>, T>(os, x, "deq[", ",", "]"); }
template <typename T>
ostream& operator<<(ostream& os, const set<T>& x) { return dmpseq<set<T>, T>(os, x, "{", ",", "}"); }
template <typename T, typename TH>
ostream& operator<<(ostream& os, const unordered_set<T, TH>& x) { return dmpseq<unordered_set<T, TH>, T>(os, x, "{", ",", "}"); }
template <typename T>
ostream& operator<<(ostream& os, const multiset<T>& x) { return dmpseq<multiset<T>, T>(os, x, "{", ",", "}"); }
template <typename TK, typename T>
ostream& operator<<(ostream& os, const map<TK, T>& x) { return dmpseq<map<TK, T>, pair<TK, T>>(os, x, "{", ",", "}"); }
template <typename TK, typename T, typename TH>
ostream& operator<<(ostream& os, const unordered_map<TK, T, TH>& x) { return dmpseq<unordered_map<TK, T, TH>, pair<TK, T>>(os, x, "{", ",", "}"); }
template <typename T, typename V>
ostream& dmpseq(ostream& os, const T& seq, const string& pre, const string& sp, const string& suf)
{
    return os << pre, for_each(begin(seq), end(seq), [&](V x) { os << x << sp; }), os << suf;
}
template <typename T>
void print(const vector<T>& x) { dmpseq<vector<T>, T>(cout, x, "", " ", "\n"); }
#ifdef HITONANODE_LOCAL
constexpr int D = 9;
#define dbg(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ") " << __FILE__ << endl
#else
constexpr int D = 19;
#define dbg(x) {}
#endif

template <typename T, typename F>
void walsh_hadamard(std::vector<T>& seq, F f)
{
    const int n = seq.size();
    assert(__builtin_popcount(n) == 1);
    for (int w = 1; w < n; w *= 2) {
        for (int i = 0; i < n; i += w * 2) {
            for (int j = 0; j < w; j++) {
                f(seq[i + j], seq[i + j + w]);
            }
        }
    }
}

// xor convolution
// ret[i] = \sum_j x[j] * y[i ^ j]
// if T is integer, ||x||_1 * ||y||_1 < numeric_limits<T>::max()
template <typename T>
std::vector<T> xorconv(std::vector<T> x, std::vector<T> y)
{
    const int n = x.size();
    assert(__builtin_popcount(n) == 1);
    assert(x.size() == y.size());
    auto f = [](T& lo, T& hi) {
        T c = lo + hi;
        hi = lo - hi, lo = c;
    };
    if (x == y) {
        walsh_hadamard(x, f), y = x;
    } else {
        walsh_hadamard(x, f), walsh_hadamard(y, f);
    }
    for (size_t i = 0; i < x.size(); i++) {
        x[i] *= y[i];
    }
    auto finv = [](T& lo, T& hi) {
        T c = lo + hi;
        hi = (lo - hi) / 2, lo = c / 2;
    };
    walsh_hadamard(x, finv);
    return x;
}


int main()
{
    int N, X;
    cin >> N >> X;

    vector<int> A(N);
    cin >> A;
    vector<lint> cnt(1 << D);
    for (auto a : A) cnt[a]++;
    cnt = xorconv(cnt, cnt);
    cnt[0] -= N;
    lint ret = 0;

    REP(d, D)
    {
        vector<lint> cntd(1 << D);
        int nd = 0;
        for (auto a : A) if (((a >> d) & 1) == 0) {
            cntd[a]++, nd++;
        }
        cntd = xorconv(cntd, cntd);
        cntd[0] -= nd;
        REP(i, X)
        {
            ret += (1LL << d) * (cnt[i] - cntd[i]) / 2;
        }
    }
    cout << ret << '\n';
}

yukicoder

結果

テストケース

ソースコード