#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; constexpr int INF = 1001001001; // constexpr int mod = 1000000007; constexpr int mod = 998244353; template inline bool chmax(T& x, T y){ if(x < y){ x = y; return true; } return false; } template inline bool chmin(T& x, T y){ if(x > y){ x = y; return true; } return false; } struct mint { int x; mint() : x(0) {} mint(int64_t y) : x(y >= 0 ? y % mod : (mod - (-y) % mod) % mod) {} mint& operator+=(const mint& p){ if((x += p.x) >= mod) x -= mod; return *this; } mint& operator-=(const mint& p){ if((x -= p.x) < 0) x += mod; return *this; } mint& operator*=(const mint& p){ x = (int)(1LL * x * p.x % mod); return *this; } mint& operator/=(const mint& p){ *this *= p.inverse(); return *this; } mint operator-() const { return mint(-x); } mint operator+(const mint& p) const { return mint(*this) += p; } mint operator-(const mint& p) const { return mint(*this) -= p; } mint operator*(const mint& p) const { return mint(*this) *= p; } mint operator/(const mint& p) const { return mint(*this) /= p; } bool operator==(const mint& p) const { return x == p.x; } bool operator!=(const mint& p) const { return x != p.x; } mint pow(int64_t n) const { mint res = 1, mul = x; while(n > 0){ if(n & 1) res *= mul; mul *= mul; n >>= 1; } return res; } mint inverse() const { return pow(mod - 2); } friend ostream& operator<<(ostream& os, const mint& p){ return os << p.x; } friend istream& operator>>(istream& is, mint& p){ int64_t val; is >> val; p = mint(val); return is; } }; int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); int N, Q; cin >> N >> Q; vector A(N); for(int i = 0; i < N; ++i) cin >> A[i]; sort(begin(A), end(A)); vector prod(N + 1, 1); for(int i = 0; i < N; ++i) prod[i + 1] = prod[i] * A[i]; vector> dp(N + 1, vector(N + 1)); dp[N][0] = 1; for(int i = N - 1; i >= 0; --i){ for(int j = 0; j <= N - 1 - i; ++j){ dp[i][j] += dp[i + 1][j] * (A[i] - 1); dp[i][j + 1] += dp[i + 1][j]; } } for(int q = 0; q < Q; ++q){ int l, r, p; cin >> l >> r >> p; int ans = 0; for(int c = l; c <= r; ++c){ int i = lower_bound(begin(A), end(A), c) - begin(A); (ans ^= (prod[i] * dp[i][p]).x) %= mod; } cout << ans << '\n'; } return 0; }