#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
constexpr int Inf = 1000000030;
constexpr ll INF= 2000000000000000000;
constexpr ll MOD = 1000000007;
const double PI = 3.1415926535897;
typedef pair<int,int> P;
typedef pair<int,P> PP;

template<class T> inline bool chmax(T &a, T b) {
    if (a < b) {
        a = b;
        return 1;
    }
    return 0;
}

template<class T> inline bool chmin(T &a, T b) {
    if (a > b) {
        a = b;
        return 1;
    }
    return 0;
}

ll mod(ll val, ll M) {
    val = val % M;
    if(val < 0) {
        val += M;
    }
    return val;
}

template<typename T>
T RS(T N, T P, T M){
    if(P == 0) {
        return 1;
    }
    if(P < 0) {
        return 0;
    }
    if(P % 2 == 0){
        ll t = RS(N, P/2, M);
        if(M == -1) return t * t;
        return t * t % M;
    }
    if(M == -1) {
        return N * RS(N,P - 1,M);
    }
    return N * RS(N, P-1, M) % M;
}

int dp[5050][5050];

int main() {
    int N,Q,P;
    cin >> N >> Q >> P;
    vector<int> vec(N);
    for(int i = 0;i < N;i++) {
        cin >> vec.at(i);
        assert(0 <= vec.at(i) && vec.at(i) < P);
    }
    for(int i = 0;i < 5050;i++) {
        for(int j = 0;j < 5050;j++) {
            dp[i][j] = -Inf;
        }
    }
    for(int i = 0;i < N;i++) {
        dp[i + 1][vec.at(i)] = i + 1;
        for(int j = 0;j < 5050;j++) {
            if(j < P) chmax(dp[i + 1][(j * vec.at(i)) % P],dp[i][j]);
        }
        for(int j = 0;j < P;j++) {
            chmax(dp[i + 1][j],dp[i][j]);
        }
    }
    for(int i = 0;i < Q;i++) {
        int l,r,k;
        cin >> l >> r >> k;
        assert(1 <= l && l <= r && r <= N);
        assert(0 <= k && k < P);
        if(dp[r][k] >= l) cout << "Yes" << endl;
        else cout << "No" << endl;
    }
}