Problem: Petya sometimes has to water his field. To water the field, Petya needs a tank with exactly V ml of water.

Petya has got N tanks, i-th of them initially containing ai ml of water. The tanks are really large, any of them can contain any amount of water (no matter how large this amount is).

Also Petya has got a scoop that can contain up to K ml of water (initially the scoop is empty). This scoop can be used to get some water from some tank, and after that pour it all into some tank (it is impossible to get water from multiple tanks without pouring it, or leave some water in the scoop when pouring it). When Petya tries to get some water from a tank, he gets min(v, K) water, where v is the current volume of water in the tank.

Is it possible to obtain a tank with exactly V ml of water using these operations? If it is possible, print a sequence of operations that allows to do it. If there are multiple ways to obtain needed amount of water in some tank, print any of them.

Input

The first line contains 3 integers: N (2 ≤ N ≤ 5000), K (1 ≤ K ≤ 5000), and V (0 ≤ V ≤ 109) — the number of tanks, the maximum volume of water the scoop can contain, and the required amount of water in some tank, respectively.

The second line contains N integers ai (0 ≤ ai ≤ 105), where ai is initial volume of water in i-th tank.

Output

If it is impossible to obtain a tank with exactly V ml of water, print NO. 

Otherwise print YES in the first line, and beginning from the second line, print the sequence of operations in the following format: 

Each line has to contain 3 numbers denoting a compressed operation: "cnt x y" (1 ≤ cnt ≤ 109, 1 ≤ x, y ≤ N), where x is the index of the tank where we get water, y is the index of the tank where we pour water, and cnt is the number of times we transfer water from tank x to tank y. 

The number of these lines must not exceed N + 5.

Examples

Input

2 3 5
2 3


Output

YES
1 2 1


Input

2 3 4
2 3


Output

NO


Input

5 2 0
1 3 5 7 9


Output

YES
2 2 1
3 3 1
4 4 1
5 5 1

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
int dp[5005][5005];
int a[5005];
void solve() {
  int n, k, v;
  scanf("%d %d %d ", &n, &k, &v);
  long long sum = 0;
  for (int i = 1; i <= n; ++i) {
    scanf(" %d ", &a[i]);
    sum += a[i];
  }
  if (sum < v) {
    printf("NO\n");
    return;
  }
  if (v % k == 0) {
    printf("YES\n");
    for (int i = 2; i <= n; ++i) {
      int cnt = a[i] / k;
      int mod = a[i] % k;
      if (mod > 0) {
        ++cnt;
      }
      if (cnt > 0) printf("%d %d %d\n", cnt, i, 1);
    }
    if (v / k > 0) printf("%d %d %d ", v / k, 1, 2);
    return;
  }
  for (int j = 0; j < k; ++j) {
    dp[0][j] = -1;
  }
  dp[0][0] = 0;
  for (int i = 1; i <= n; ++i) {
    for (int j = 0; j < k; ++j) {
      dp[i][j] = -1;
    }
    int rem = a[i] / k;
    int mod = a[i] % k;
    for (int j = 0; j < k; ++j) {
      if (dp[i - 1][j] != -1) {
        dp[i][j] = dp[i - 1][j] + rem;
      }
    }
    for (int j = 0; j < k; ++j) {
      if (dp[i - 1][j] != -1) {
        dp[i][(j + mod) % k] =
            max(dp[i][(j + mod) % k],
                dp[i - 1][j] + rem + ((j + mod) >= k ? 1 : 0));
      }
    }
  }
  int g = (sum - v) % k;
  if (dp[n][v % k] >= v / k) {
    int curMod = v % k;
    int curRem = v / k;
    int lastRem = v / k;
    curRem = dp[n][curMod];
    printf("YES\n");
    vector<int> vertex;
    vector<int> nv;
    for (int i = n; i >= 1; --i) {
      int rem = a[i] / k;
      int mod = a[i] % k;
      if (dp[i - 1][curMod] + rem == curRem) {
        curRem -= rem;
        nv.push_back(i);
      } else {
        int cnt = rem;
        if (curMod - mod < 0) {
          --curRem;
          ++cnt;
        }
        curMod = (curMod - mod + k) % k;
        vertex.push_back(i);
        curRem -= rem;
      }
    }
    assert(vertex.size() > 0);
    for (int i = 1; i < vertex.size(); ++i) {
      int node = vertex[i];
      int rem = a[node] / k;
      int mod = a[node] % k;
      int cnt = rem;
      if (mod > 0) {
        ++cnt;
      }
      printf("%d %d %d\n", cnt, node, vertex[0]);
    }
    for (int node : nv) {
      if (a[node] / k > 0) printf("%d %d %d\n", a[node] / k, node, vertex[0]);
    }
    int un = 1;
    if (un == vertex[1]) {
      un = 2;
    }
    int cnt = dp[n][v % k] - lastRem;
    assert(cnt >= 0);
    if (cnt > 0) printf("%d %d %d\n", cnt, vertex[0], un);
  } else if (dp[n][g] != -1) {
    printf("YES\n");
    vector<int> vertex;
    vector<int> nv;
    int curMod = g;
    int curRem = dp[n][g];
    long long tsum = 0;
    for (int i = n; i >= 1; --i) {
      int rem = a[i] / k;
      int mod = a[i] % k;
      if (dp[i - 1][curMod] + rem == curRem) {
        curRem -= rem;
        nv.push_back(i);
      } else {
        int cnt = rem;
        if (curMod - mod < 0) {
          --curRem;
          ++cnt;
        }
        tsum += a[i];
        curMod = (curMod - mod + k) % k;
        vertex.push_back(i);
        curRem -= rem;
      }
    }
    int un = 1;
    if (un == vertex[0]) {
      un = 2;
    }
    for (int i = 1; i < vertex.size(); ++i) {
      int node = vertex[i];
      printf("%d %d %d\n", (a[node] + k - 1) / k, node, vertex[0]);
    }
    for (int node : nv) {
      if ((a[node] + k - 1) / k > 0 && un != node) {
        printf("%d %d %d\n", (a[node] + k - 1) / k, node, un);
      }
    }
    if (tsum / k > 0) {
      printf("%lld %d %d\n", tsum / k, vertex[0], un);
    }
    tsum = sum - (sum - v) % k;
    assert((tsum - v) % k == 0);
    if ((tsum - v) / k > 0)
      printf("%lld %d %d\n", (tsum - v) / k, un, vertex[0]);
  } else {
    printf("NO\n");
  }
}
int main() {
  solve();
  return 0;
}