結果

問題 No.803 Very Limited Xor Subset
ユーザー veqccveqcc
提出日時 2019-03-18 19:28:02
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 2,586 bytes
コンパイル時間 1,003 ms
コンパイル使用メモリ 107,072 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-25 20:14:38
合計ジャッジ時間 2,934 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 2 ms
4,380 KB
testcase_07 AC 2 ms
4,380 KB
testcase_08 AC 1 ms
4,376 KB
testcase_09 AC 1 ms
4,380 KB
testcase_10 AC 1 ms
4,376 KB
testcase_11 AC 2 ms
4,376 KB
testcase_12 AC 1 ms
4,376 KB
testcase_13 AC 2 ms
4,376 KB
testcase_14 AC 2 ms
4,380 KB
testcase_15 AC 2 ms
4,376 KB
testcase_16 AC 2 ms
4,380 KB
testcase_17 AC 2 ms
4,376 KB
testcase_18 AC 2 ms
4,376 KB
testcase_19 AC 2 ms
4,376 KB
testcase_20 AC 2 ms
4,376 KB
testcase_21 AC 2 ms
4,380 KB
testcase_22 AC 2 ms
4,376 KB
testcase_23 AC 2 ms
4,380 KB
testcase_24 AC 2 ms
4,376 KB
testcase_25 AC 2 ms
4,380 KB
testcase_26 AC 2 ms
4,376 KB
testcase_27 AC 2 ms
4,376 KB
testcase_28 AC 2 ms
4,380 KB
testcase_29 AC 2 ms
4,376 KB
testcase_30 AC 2 ms
4,376 KB
testcase_31 AC 2 ms
4,380 KB
testcase_32 AC 2 ms
4,376 KB
testcase_33 AC 2 ms
4,380 KB
testcase_34 AC 2 ms
4,376 KB
testcase_35 AC 2 ms
4,376 KB
testcase_36 AC 2 ms
4,376 KB
testcase_37 AC 2 ms
4,376 KB
testcase_38 AC 2 ms
4,376 KB
testcase_39 AC 2 ms
4,380 KB
testcase_40 AC 2 ms
4,380 KB
testcase_41 AC 2 ms
4,376 KB
testcase_42 AC 2 ms
4,376 KB
testcase_43 AC 2 ms
4,376 KB
testcase_44 AC 2 ms
4,380 KB
testcase_45 AC 2 ms
4,376 KB
testcase_46 AC 2 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <string>
#include <vector>
#include <random>
#include <bitset>
#include <queue>
#include <cmath>
#include <stack>
#include <set>
#include <map>
typedef long long ll;
using namespace std;

const ll MOD = 1000000007LL;
const int MAX_ROW = 600;
const int MAX_COL = 600;

struct BitMatrix {
    int n, m;
    bitset<MAX_COL> val[MAX_ROW];
    BitMatrix(int n_ = 1, int m_ = 1) { n = n_; m = m_; }
    inline bitset<MAX_COL> &operator [] (int i) { return val[i]; }
    inline friend ostream &operator << (ostream &s, BitMatrix M) {
        s << "\n";
        for (int i = 0; i < M.n; i++) {
            for (int j = 0; j < M.m; j++) s << M.val[i][j];
            s << "\n";
        }
        return s;
    }
};

inline BitMatrix operator * (BitMatrix A, BitMatrix B) {
    BitMatrix R(A.n, B.m);
    BitMatrix tB(B.m, B.n);
    for (int i = 0; i < tB.n; i++) for (int j = 0; j < tB.m; j++) tB[i][j] = B[j][i];
    for (int i = 0; i < R.n; i++) for (int j = 0; j < R.m; j++) R[i][j] = (A[i] & tB[j]).any();
    return R;
}

int linear_equation(BitMatrix A, vector<int> b) {
    int rank = 0;
    for (int i = 0; i < A.n; i++) A[i][A.m] = b[i];
    for (int i = 0; i < A.m; i++) {
        int pivot = -1;
        for (int j = rank; j < A.n; j++) if (A[j][i]) { pivot = j; break; }
        if (pivot != -1) {
            swap(A[pivot], A[rank]);
            for (int j = 0; j < A.n; j++) if (j != rank && A[j][i]) A[j] ^= A[rank];
            rank++;
        }
    }
    for (int i = rank; i < A.n; i++) if (A[i][A.m]) return -1;
    return rank;
}

ll modpow(ll a, ll n, ll mod) {
    ll res = 1;
    while (n) {
        if (n & 1) res = res * a % mod;
        a = a * a % mod;
        n >>= 1;
    }
    return res;
}

const int digit = 35;

int main() {
    cin.sync_with_stdio(false);
    cin.tie(0);

    int n, m;
    ll x;
    cin >> n >> m >> x;

    vector <ll> a(n);
    for (int i = 0; i < n; i++) cin >> a[i];

    BitMatrix A(digit + m, n);
    vector <int> b(digit + m);
    for (int d = 0; d < digit; d++) {
        for (int i = 0; i < n; i++) {
            if (a[i] & (1 << d)) A[d][i] = 1;
        }
        if (x & (1 << d)) b[d] = 1;
    }

    for (int d = 0; d < m; d++) {
        int type, l, r;
        cin >> type >> l >> r;
        l--; r--;
        for (int i = l; i <= r; i++) A[d + digit][i] = 1;
        b[d + digit] = type;
    }

    int rank = linear_equation(A, b);
    if (rank == -1) cout << 0 << "\n";
    else cout << modpow(2LL, n - rank, MOD) << "\n";
    return 0;
}
0