結果

問題 No.2237 Xor Sum Hoge
コンテスト
ユーザー hitonanode
提出日時 2023-03-11 12:59:39
言語 C++23
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 1,031 ms / 10,000 ms
コード長 1,141 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,865 ms
コンパイル使用メモリ 142,504 KB
実行使用メモリ 7,532 KB
最終ジャッジ日時 2024-09-18 06:19:29
合計ジャッジ時間 20,714 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <vector>
using namespace std;
#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)

#include <atcoder/modint>
#include <atcoder/convolution>
using mint = atcoder::modint998244353;

int main() {
    int N;
    long long B, C;
    cin >> N >> B >> C;

    const int sz = N * 2;
    vector<mint> carry_rev(sz);
    carry_rev.back() = 1;

    vector<mint> trans_even(N + 1), trans_odd(N + 1);
    mint num = 1, den = 1;
    REP(i, N + 1) {
        (i % 2 ? trans_odd : trans_even).at(i) = num / den;
        num *= N - i, den *= i + 1;
    }

    IREP(d, 60) {
        IREP(i, sz) carry_rev.at(sz - 1 - i) = i % 2 ? 0 : carry_rev.at(sz - 1 - i / 2);

        if ((B >> d) & 1) {
            carry_rev.push_back(0);
            carry_rev.erase(carry_rev.begin());
        }

        carry_rev = atcoder::convolution(carry_rev, (C >> d) % 2 ? trans_odd : trans_even);
        carry_rev.resize(sz);
    }

    cout << carry_rev.back().val() << '\n';
}
0