#include using namespace std; typedef long long ll; template class Vec { protected: using iterator = typename std::vector::iterator; using const_iterator = typename std::vector::const_iterator; using reference = T &; using const_reference = const T &; std::vector v; template Vec unop_new(Unop op) const { Vec res(v.size()); transform(begin(v), end(v), res.begin(), op); return res; } template Vec &binop(const Vec &r, Binop op) { transform(r.begin(), r.end(), v.begin(), v.begin(), op); return *this; } template Vec binop_new(const Vec &r, Binop op) const { Vec res(v.size()); transform(r.begin(), r.end(), v.begin(), res.begin(), op); return res; } public: Vec(int n) : v(n) {} Vec(int n, const T &val) : v(n, val) {} Vec(const std::vector &w) : v(w) {} int size() const noexcept { return v.size(); } const_iterator begin() const noexcept { return v.begin(); } const_iterator end() const noexcept { return v.end(); } iterator begin() noexcept { return v.begin(); } iterator end() noexcept { return v.end(); } reference operator[](int i) { return v[i]; } const_reference operator[](int i) const { return v[i]; } Vec operator-() const { return unop_new([](T val) { return -val; }); }; Vec &operator+=(const Vec &r) { return binop(r, [](T x, T y) { return x + y; }); } Vec &operator-=(const Vec &r) { return binop(r, [](T x, T y) { return x - y; }); } Vec operator+(const Vec &r) const { return binop_new(r, [](T x, T y) { return x + y; }); } Vec operator-(const Vec &r) const { return binop_new(r, [](T x, T y) { return x - y; }); } T dot(const Vec &r) const { return inner_product(v.begin(), v.end(), r.begin(), T(0)); } T norm() const { return this->dot(v); } void push_back(const T &r) { v.push_back(r); } void concat(const Vec &r) { v.insert(v.end(), r.begin(), r.end()); } }; template class Matrix : public Vec> { public: using Vec>::Vec; Matrix(int n, int m, const T &val) : Vec>::Vec(n, Vec(m, val)) {} int Y() const { return this->size(); } int X() const { return (*this)[0].size(); } Matrix transpose() const { const int row = Y(), col = X(); Matrix res(col, row); for (int j = 0; j < col; ++j) { for (int i = 0; i < row; ++i) { res[j][i] = (*this)[i][j]; } } return res; } Matrix operator*(const Matrix &r) const { Matrix tr = r.transpose(); const int row = Y(), col = tr.Y(); assert(X() == tr.X()); Matrix res(row, col); for (int i = 0; i < row; ++i) { for (int j = 0; j < col; ++j) { res[i][j] = (*this)[i].dot(tr[j]); } } return res; } Vec operator*(const Vec &r) const { const int row = Y(), col = r.Y(); assert(r.size() == col); Vec res(row); for (int i = 0; i < row; ++i) { res[i] = (*this)[i].dot(r); } return res; } Matrix &operator*=(const Matrix &r) { return *this = *this * r; } Matrix operator^(ll n) const { const int m = Y(); assert(m == X()); Matrix A = *this, res(m, m, 0); for (int i = 0; i < m; ++i) res[i][i] = 1; while (n > 0) { if (n % 2) res *= A; A = A * A; n /= 2; } return res; } void concat_right(const Vec &r) { const int n = Y(); assert(n == r.size()); for (int i = 0; i < n; ++i) { (*this)[i].push_back(r[i]); } } void concat_right(const Matrix &r) { const int n = Y(); assert(n == r.Y()); for (int i = 0; i < n; ++i) { (*this)[i].concat(r[i]); } } void concat_below(const Vec &r) { assert(Y() == 0 || X() == r.size()); this->push_back(r); } void concat_below(const Matrix &r) { assert(Y() == 0 || X() == r.X()); for (Vec i : r) (*this).push_back(i); } int rank() const { Matrix A = *this; if (Y() == 0) return 0; const int n = Y(), m = X(); int r = 0; for (int i = 0; r < n && i < m; ++i) { int pivot = r; for (int j = r + 1; j < n; ++j) { if (abs(A[j][i]) > abs(A[pivot][i])) pivot = j; } std::swap(A[pivot], A[r]); if (A[r][i] == 0) continue; for (int k = m - 1; k >= i; --k) A[r][k] = A[r][k] / A[r][i]; for (int j = r + 1; j < n; ++j) { for (int k = m - 1; k >= i; --k) { A[j][k] ^= A[r][k] * A[j][i]; } } ++r; } return r; } }; Vec create(int x) { Vec res(30); for (int i = 0; i < 30; i++) { res[i] = x >> i & 1; } return res; } int A[310]; template class Modulo { int n; static typename std::enable_if::type inv(ll a, ll p) { return (a == 1 ? 1 : (1 - p * inv(p % a, a)) / a + p); } public: Modulo() : n(0) { ; } Modulo(int m) : n(m) { if (n >= M) n %= M; else if (n < 0) n = (n % M + M) % M; } Modulo(ll m) { if (m >= M) m %= M; else if (m < 0) m = (m % M + M) % M; n = m; } explicit operator int() const { return n; } explicit operator ll() const { return n; } bool operator==(const Modulo &a) const { return n == a.n; } Modulo &operator+=(const Modulo &a) { n += a.n; if (n >= M) n -= M; return *this; } Modulo &operator-=(const Modulo &a) { n -= a.n; if (n < 0) n += M; return *this; } Modulo &operator*=(const Modulo &a) { n = (ll(n) * a.n) % M; return *this; } Modulo operator+(const Modulo &a) const { Modulo res = *this; return res += a; } Modulo operator-(const Modulo &a) const { Modulo res = *this; return res -= a; } Modulo operator-() const { return Modulo(0) - *this; } Modulo operator*(const Modulo &a) const { Modulo res = *this; return res *= a; } Modulo operator^(ll m) const { if (m == 0) return Modulo(1); const Modulo a = *this; Modulo res = (a * a) ^(m / 2); return m % 2 ? res * a : res; } typename std::enable_if::type operator/(const Modulo &a) const { return *this * inv(ll(a), M); } typename std::enable_if::type operator/=(const Modulo &a) { return *this *= inv(ll(a), M); } friend bool is_zero(const Modulo &x) { return int(x) == 0; } friend int abs(const Modulo &x) { return int(x); } static Modulo fact(int n, bool sw = true) { static std::vector v1 = {1}, v2 = {1}; if (n >= (int) v1.size()) { const int from = v1.size(), to = n + 1024; v1.reserve(to); v2.reserve(to); for (int i = from; i < to; ++i) { v1.push_back(v1.back() * Modulo(i)); v2.push_back(v2.back() / Modulo(i)); } } return sw ? v1[n] : v2[n]; } static Modulo comb(int a, int b) { if (b < 0 || b > a) return 0; return Modulo::fact(a, true) * Modulo::fact(b, false) * Modulo::fact(a - b, false); } }; typedef Modulo<1000000007, true> mInt; int main() { int N, M, X; cin >> N >> M >> X; auto xv = create(X); for (int i = 0; i < N; i++) cin >> A[i]; Matrix mat1(30 + M, N); for (int j = 0; j < 30; j++) { mat1[j][N] = xv[j]; } for (int i = 0; i < N; i++) { auto r = create(A[i]); for (int j = 0; j < 30; j++) { mat1[j][i] = r[j]; } } for (int i = 30; i < 30 + M; i++) { int t, a, b; cin >> t >> a >> b; --a; for (int j = a; j < b; j++) { mat1[i][j] = 1; } mat1[i][N] = t; } mInt ans = mInt(2) ^ (N - mat1.rank()); cout << (int) ans << endl; }