結果

問題 No.803 Very Limited Xor Subset
ユーザー heno239heno239
提出日時 2019-03-14 20:57:56
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 4,562 bytes
コンパイル時間 1,237 ms
コンパイル使用メモリ 122,720 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-09-09 13:45:00
合計ジャッジ時間 3,171 ms
ジャッジサーバーID
(参考情報)
judge13 / judge11
このコードへのチャレンジ(β)

テストケース

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

ソースコード

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;
		rep(j, 30) {
			if (a&(1 << j)) {
				A[i].first[j] = 1;
			}
		}
		A[i].second = i;
		/*cout << i << " ";
		rep(j, 30) {
			cout << A[i].first[j];
		}
		cout << endl;*/
	}
	int r=Gaussian_elimination(A);
	/*cout << r << "hello?" << endl;
	rep(i, n) {
		cout << A[i].second << " ";
		rep(j, 30) {
			cout << A[i].first[j];
		}
		cout << endl;
	}*/
	//debug
	/*rep(i, n) {
		rep(j, 30) {
			cout << A[i].first[j];
		}
		cout << endl;
	}
*/
	int s = 0;
	rep(i, r) {
		rep(j, 30) {
			if (A[i].first[j] && (x&(1 << j))) {
				use[A[i].second] = true;
			}
		}
		if (use[A[i].second]) {
			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;
	}
	B.resize(n-r+1);
	bitset<300> b = 0;
	rep(i, m) {
		int num = 0;
		Rep1(j, z[i].l, z[i].r) {
			if (use[j])num ^= 1;
		}
		//cout << num << endl;
		if (num)b[i] = 1;
	}
	B[n - r].first.resize(m);
	rep(j, m) {
		B[n - r].first[j] = b[j];
	}
	rep(j, m)b[j] = 0;
	rep(j, m) {
		if (z[j].t)b[j] = 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 (num)B[i - r].first[j] = 1;
		}
		/*cout << "hello!!! I am "<<A[i].second<< endl;
		rep(j, n) {
			cout << use[j];
		}
		cout << endl;*/
	}
	//cout << n - r << endl;
	/*rep(i, n - r) {
		rep(j, m) {
			cout << B[i].first[j];
		}
		cout << endl;
	}*/
	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];
			}
		}
	}
	//cout << n - r << endl;
	/*rep(i, n-r) {
		rep(j,m){
			cout << B[i].first[j];
		}
		cout << endl;
	}
	rep(i, m) {
		cout << b[i];
	}
	cout << endl;
	rep(i, m) {
		cout << b_[i];
	}
	cout << endl;*/
	if (b != b_) {
		cout << 0 << endl; return;
	}
	ll ans = 1;
	int d = n - r + 1 - rn; d--;
	rep(j, m) {
		if (b_[j]) {
			d++; break;
		}
	}
	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) {
		cin >> z[i].t >> z[i].l >> z[i].r;
		z[i].l--; z[i].r--;
	}
	int cnt = 0;
	rep(i, (1 << n)) {
		int u[15] = {};
		int s = 0;
		rep(j, n) {
			if (i&(1 << j)) {
				s ^= a[j];
				u[j] = 1;
			}
		}
		if (s == x) {
			bool f = true;
			rep(j, m) {
				int num = 0;
				Rep1(k, z[j].l, z[j].r) {
					if (i&(1 << k))num ^= 1;
				}
				if (num^z[j].t)f = false;
			}
			if(f)cnt++;
		}
	}
	cout << cnt << endl;
}
int main() {
	solve();
	//check();
	//stop
	return 0;
}
0