結果

問題 No.803 Very Limited Xor Subset
ユーザー heno239
提出日時 2019-03-14 20:01:53
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 3,797 bytes
コンパイル時間 1,230 ms
コンパイル使用メモリ 121,308 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-06-27 06:41:09
合計ジャッジ時間 3,026 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 WA * 3
other AC * 13 WA * 29 RE * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<string>
#include<cstdio>
#include<vector>
#include<cmath>
#include<algorithm>
#include<functional>
#include<iomanip>
#include<queue>
#include<ciso646>
#include<random>
#include<map>
#include<set>
#include<complex>
#include<bitset>
#include<stack>
#include<unordered_map>
#include<utility>
using namespace std;
typedef long long ll;
typedef unsigned int ui;
const ll mod = 1000000007;
typedef long double ld;
const ll INF = mod * mod;
typedef pair<int, int> P;
#define stop char nyaa;cin>>nyaa;
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define Rep(i,sta,n) for(int i=sta;i<n;i++)
#define rep1(i,n) for(int i=1;i<=n;i++)
#define per1(i,n) for(int i=n;i>=1;i--)
#define Rep1(i,sta,n) for(int i=sta;i<=n;i++)
typedef complex<ld> Point;
const ld eps = 1e-8;
const ld pi = acos(-1.0);
typedef pair<ld, ld> LDP;
typedef pair<ll, ll> LP;



typedef vector<int> vec;
typedef pair<vec, int> vect;
typedef vector<vect> mat;
//bit.ver
//rankを返す
int Gaussian_elimination(mat &A) {
	int n = A.size(); int m = A[0].first.size();
	int le = 0;
	rep(j, m) {
		Rep(i, le, n) {
			if (A[i].first[j]) {
				swap(A[le], A[i]);
				rep(i2, n) {
					if (i2 == le)continue;
					if (A[i2].first[j]) {
						rep(k, m) {
							A[i2].first[k] ^= A[le].first[k];
						}
					}
				}
				le++; break;
			}
		}
	}
	return le;
}
int n, m, x;
mat A,B;
bool use[300];
struct zyoken {
	int t, l, r;
};
zyoken z[300];

int mem[300];
void solve() {
	cin >> n >> m >> x;
	A.resize(n);
	rep(i, n) {
		A[i].first.resize(30);
		int a; cin >> a; mem[i] = a;
		per(j, 30) {
			if (a&(1 << j)) {
				A[i].first[j] = 1;
			}
		}
		A[i].second = i;
	}
	int r=Gaussian_elimination(A);
	//cout << r << "hello?" << endl;
	//debug
	/*rep(i, n) {
		rep(j, 30) {
			cout << A[i].first[j];
		}
		cout << endl;
	}
*/
	rep(i, r) {
		rep(j, 30) {
			if (A[i].first[j] && (x&(1 << j))) {
				use[i] = true;
			}
		}
	}
	int s = 0;
	rep(i, r) {
		if (use[i]) {
			rep(j, 30) {
				if (A[i].first[j]) {
					s += (1 << j);
				}
			}
		}
	}
	rep(i, m) {
		int typ, le, ri; cin >> typ >> le >> ri;
		le--; ri--;
		z[i] = { typ,le,ri };
	}
	//cout << s << endl;
	if (s != x) {
		cout << 0 << endl; return;
	}
	bool cnt[300];
	rep(i, r) {
		if (use[i]) {
			cnt[A[i].second] = true;
		}
	}
	B.resize(n-r);
	bitset<300> b;
	rep(i, m) {
		int num = 0;
		Rep1(j, z[i].l, z[i].r) {
			if (cnt[j])num ^= 1;
		}
		if (z[i].t^num)b[i] = 1;
	}

	Rep(i, r, n) {
		B[i - r].first.resize(m);
		fill(use, use + n, false);
		use[A[i].second] = true;
		int xu = mem[A[i].second];

		rep(j, r) {
			rep(k, 30) {
				if (A[j].first[k] && (xu&(1 << k))) {
					use[A[j].second] = true;
				}
			}
		}
		rep(j, m) {
			int num = 0;
			Rep1(k, z[j].l, z[j].r) {
				if (use[k])num ^= 1;
			}
			if (z[j].t^num)B[i - r].first[j] = 1;
		}
	}
	int rn=Gaussian_elimination(B);
	fill(use, use + m, false);
	//bが存在してるか判定
	bitset<300> b_;
	rep(i, rn) {
		rep(j, m) {
			if (b[j] && B[i].first[j])use[i] = true;
		}
		if (use[i]) {
			rep(j, m) {
				b_[j] = b_[j]^B[i].first[j];
			}
		}
	}

	/*rep(i, n-r) {
		rep(j,m){
			cout << B[i].first[j];
		}
		cout << endl;
	}
*/
	//rep(i, 300) {
	//	cout << b[i];
	//}
	//cout << endl;
	//rep(i, 300) {
	//	cout << b_[i];
	//}
	//cout << endl;
	if (b != b_) {
		cout << 0 << endl; return;
	}
	ll ans = 1;
	int d = n - r - rn;
	rep(i, d) {
		ans = ans * 2 % mod;
	}
	cout << ans << endl;
}
void check() {
	int n, m, x; cin >> n >> m >> x;
	vector<int> a(n);
	rep(i, n) {
		cin >> a[i];
	}
	rep(i, m) {
		int t, l, r; cin >> t >> l >> r;
	}
	int cnt = 0;
	rep(i, (1 << n)) {
		int s = 0;
		rep(j, n) {
			if (i&(1 << j)) {
				s ^= a[j];
			}
		}
		if (s == x)cnt++;
	}
	cout << cnt << endl;
}
int main() {
	solve();
	//check();
	//stop
	return 0;
}
0