結果

問題 No.3486 Draw a Rainbow
コンテスト
ユーザー 👑 AngrySadEight
提出日時 2025-12-16 11:03:51
言語 C++17
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 1,040 ms / 4,000 ms
コード長 2,129 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,173 ms
コンパイル使用メモリ 96,256 KB
実行使用メモリ 20,404 KB
最終ジャッジ日時 2026-03-27 20:51:12
合計ジャッジ時間 11,974 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <vector>
using namespace std;
using ll = long long;

const ll mod = 998244353;

vector<ll> fzt(vector<ll> &vec, ll N){
    vector<ll> ret(1 << N);
    for (ll i = 0; i < (1 << N); i++){
        ret[i] = vec[i];
    }
    for (ll i = 0; i < N; i++){
        for (ll j = 0; j < (1 << N); j++){
            if (!((j >> i) & 1)){
                ret[j + (1 << i)] = (ret[j + (1 << i)] + ret[j]) % mod;
            }
        }
    }
    return ret;
}

vector<ll> fmt(vector<ll> &vec, ll N){
    vector<ll> ret(1 << N);
    for (ll i = 0; i < (1 << N); i++){
        ret[i] = vec[i];
    }
    for (ll i = 0; i < N; i++){
        for (ll j = 0; j < (1 << N); j++){
            if ((j >> i) & 1){
                ret[j] = (ret[j] + mod - ret[j - (1 << i)]) % mod;
            }
        }
    }
    return ret;
}

vector<ll> or_convolution(vector<ll> &f, vector<ll> &g, ll N){
    vector<ll> fw = fzt(f, N);
    vector<ll> gw = fzt(g, N);
    vector<ll> hw(1 << N);
    for (ll i = 0; i < (1 << N); i++){
        hw[i] = (fw[i] * gw[i]) % mod;
    }
    vector<ll> h = fmt(hw, N);
    return h;
}

int main(){
    ll N, M, L;
    cin >> N >> M >> L;

    vector<ll> A(N);
    vector<ll> B(N);
    for (ll i = 0; i < N; i++){
        cin >> A[i] >> B[i];
        A[i]--;
        B[i] -= (M + 1);
    }

    vector<ll> pow_2(N + 1);
    pow_2[0] = 1;
    for (ll i = 0; i < N; i++){
        pow_2[i + 1] = (pow_2[i] * 2LL) % mod;
    }
    vector<vector<ll>> cnts(M, vector<ll>(L, 0));
    for (ll i = 0; i < N; i++){
        cnts[A[i]][B[i]]++;
    }

    vector<ll> ans_vec(1 << L, 0);
    ans_vec[0] = 1;
    for (ll i = 0; i < M; i++){
        vector<ll> g_vec(1 << L, 0);
        for (ll j = 1; j < (1 << L); j++){
            ll sub = 1;
            for (ll k = 0; k < L; k++){
                if ((j >> k) & 1){
                    sub = (sub * (pow_2[cnts[i][k]] + mod - 1)) % mod;
                }
            }
            g_vec[j] = sub;
        }
        vector<ll> ans_vec_new = or_convolution(ans_vec, g_vec, L);
        ans_vec = ans_vec_new;
    }
    cout << ans_vec[(1 << L) - 1] << endl;
}
0