#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include 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 P; #define stop char nyaa;cin>>nyaa; #define rep(i,n) for(int i=0;i=0;i--) #define Rep(i,sta,n) for(int i=sta;i=1;i--) #define Rep1(i,sta,n) for(int i=sta;i<=n;i++) typedef complex Point; const ld eps = 1e-8; const ld pi = acos(-1.0); typedef pair LDP; typedef pair LP; typedef vector vec; typedef pair vect; typedef vector 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 "< 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 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; }