#include <bits/stdc++.h>
using namespace std;

#define rep(i, n) for (int i = 0; i < (int)n; i++)

struct LR {
  int l;
  int r;
  LR() {}
  LR(int l, int r) : l(l), r(r) {}
  bool operator==(const LR& rhs) const {
    return (l == rhs.l && r == rhs.r);
  }
};

int main() {
  int N, Q;
  cin >> N >> Q;
  vector<int> A(N);
  rep(i, N) {
    cin >> A[i];
  }
  vector<int> ki(Q);
  vector<LR> lri(Q);
  rep(i, Q) {
    int L, R, K;
    cin >> L >> R >> K;
    L--; R--;
    lri[i] = LR(L, R);
    ki[i] = K;
  }
  map<int, vector<int>> mp;
  vector<vector<int>> div(N);
  rep(i, N) {
    for (int j = 1; j*j <= A[i]; j++) {
      if (A[i] % j == 0) {
	div[i].push_back(j);
	if (j*j != A[i])
	  div[i].push_back(A[i]/j);
      }
    }
  }
  rep(i, N) {
    for (int j : div[i]) {
      mp[j].push_back(i);
    }
  }
  rep(i, Q) {
    int n = lower_bound(mp[ki[i]].begin(),
			mp[ki[i]].end(),
			lri[i].l) - mp[ki[i]].begin();
    int m = upper_bound(mp[ki[i]].begin(),
			mp[ki[i]].end(),
			lri[i].r) - mp[ki[i]].begin();;
    cout << m - n << endl;
  }
  return 0;
}