Problem: You must have heard all about the Foolland on your Geography lessons. Specifically, you must know that federal structure of this country has been the same for many centuries. The country consists of n cities, some pairs of cities are connected by bidirectional roads, each road is described by its length li.

The fools lived in their land joyfully, but a recent revolution changed the king. Now the king is Vasily the Bear. Vasily divided the country cities into regions, so that any two cities of the same region have a path along the roads between them and any two cities of different regions don't have such path. Then Vasily decided to upgrade the road network and construct exactly p new roads in the country. Constructing a road goes like this:

  1. We choose a pair of distinct cities u, v that will be connected by a new road (at that, it is possible that there already is a road between these cities). 
  2. We define the length of the new road: if cities u, v belong to distinct regions, then the length is calculated as min(109, S + 1) (S — the total length of all roads that exist in the linked regions), otherwise we assume that the length equals 1000. 
  3. We build a road of the specified length between the chosen cities. If the new road connects two distinct regions, after construction of the road these regions are combined into one new region. 



Vasily wants the road constructing process to result in the country that consists exactly of q regions. Your task is to come up with such road constructing plan for Vasily that it meets the requirement and minimizes the total length of the built roads.

Input

The first line contains four integers n (1 ≤ n ≤ 105), m (0 ≤ m ≤ 105), p (0 ≤ p ≤ 105), q (1 ≤ q ≤ n) — the number of cities in the Foolland, the number of existing roads, the number of roads that are planned to construct and the required number of regions.

Next m lines describe the roads that exist by the moment upgrading of the roads begun. Each of these lines contains three integers xi, yi, li: xi, yi — the numbers of the cities connected by this road (1 ≤ xi, yi ≤ n, xi ≠ yi), li — length of the road (1 ≤ li ≤ 109). Note that one pair of cities can be connected with multiple roads.

Output

If constructing the roads in the required way is impossible, print a single string "NO" (without the quotes). Otherwise, in the first line print word "YES" (without the quotes), and in the next p lines print the road construction plan. Each line of the plan must consist of two distinct integers, giving the numbers of the cities connected by a road. The road must occur in the plan in the order they need to be constructed. If there are multiple optimal solutions, you can print any of them.

Examples

Input

9 6 2 2
1 2 2
3 2 1
4 6 20
1 3 8
7 8 3
5 7 2


Output

YES
9 5
1 9


Input

2 0 1 2


Output

NO


Input

2 0 0 2


Output

YES

Note

Consider the first sample. Before the reform the Foolland consists of four regions. The first region includes cities 1, 2, 3, the second region has cities 4 and 6, the third region has cities 5, 7, 8, the fourth region has city 9. The total length of the roads in these cities is 11, 20, 5 and 0, correspondingly. According to the plan, we first build the road of length 6 between cities 5 and 9, then the road of length 23 between cities 1 and 9. Thus, the total length of the built roads equals 29.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
#pragma comment(linker, "/STACK:66777216")
#pragma hdrstop
using namespace std;
inline void sIO() {}
inline void iIO() {
  freopen("input.txt", "r", stdin);
  freopen("output.txt", "w", stdout);
}
inline void fIO(string fn) {
  freopen((fn + ".in").c_str(), "r", stdin);
  freopen((fn + ".out").c_str(), "w", stdout);
}
inline void TM() {}
inline void swap(short& a, short& b) { b ^= a ^= b ^= a; }
inline void swap(int& a, int& b) { b ^= a ^= b ^= a; }
inline void swap(char& a, char& b) { b ^= a ^= b ^= a; }
inline void swap(long long& a, long long& b) { b ^= a ^= b ^= a; }
template <class T>
inline T abs(T x) {
  return x < 0 ? -x : x;
}
template <class T>
inline T sqr(T x) {
  return x * x;
}
template <class T>
inline T min(T& a, T& b) {
  return a < b ? a : b;
}
template <class T>
inline T max(T& a, T& b) {
  return a > b ? a : b;
}
template <class T>
inline T gcd(T a, T b) {
  if (a < b) swap(a, b);
  while (b) {
    a %= b;
    swap(a, b);
  }
  return a;
}
template <class T>
inline T lcm(T a, T b) {
  return a / gcd(a, b) * b;
}
template <class T>
inline bool isPrime(T n) {
  if (n < 2) return false;
  T kk = (T)sqrt(n + 0.);
  for (T i = 2; i <= kk; ++i)
    if (!(n % i)) return false;
  return true;
}
template <class T>
inline string toa(T x) {
  stringstream ss;
  ss << x;
  string ret;
  ss >> ret;
  return ret;
}
template <class T>
inline T ppow(T a, long long b) {
  T ret = 1;
  while (b) {
    if (b & 1) ret *= a;
    a *= a;
    b >>= 1;
  }
  return ret;
}
template <class T>
inline T ppow(T a, long long b, long long md) {
  T ret = 1;
  a %= md;
  while (b) {
    if (b & 1) ret = ret * a % md;
    a = a * a % md;
    b >>= 1;
  }
  return ret % md;
}
inline int toi(string s) {
  stringstream ss;
  ss << s;
  int ret;
  ss >> ret;
  return ret;
}
inline long long tol(string s) {
  stringstream ss;
  ss << s;
  long long ret;
  ss >> ret;
  return ret;
}
inline int Random() { return ((rand() << 16) | rand()); }
inline int Random(int l, int r) {
  if (l == r) return l;
  return Random() % (r - l) + l;
}
inline char upperCase(char ch) {
  return (ch >= 'a' && ch <= 'z') ? ch ^ 32 : ch;
}
inline char lowerCase(char ch) {
  return (ch >= 'A' && ch <= 'Z') ? ch ^ 32 : ch;
}
inline string upperCase(string s) {
  int ls = s.length();
  for (int i = 0; i < ls; ++i)
    if (s[i] >= 'a' && s[i] <= 'z') s[i] ^= 32;
  return s;
}
inline string lowerCase(string s) {
  int ls = s.length();
  for (int i = 0; i < ls; ++i)
    if (s[i] >= 'A' && s[i] <= 'Z') s[i] ^= 32;
  return s;
}
inline int dig(char ch) { return ch - 48; }
inline bool isAlpha(char ch) {
  return (ch >= 'A' && ch <= 'Z' || ch >= 'a' && ch <= 'z');
}
inline bool isDigit(char ch) { return (ch >= '0' && ch <= '9'); }
inline bool isLowerCase(char ch) { return (ch >= 'a' && ch <= 'z'); }
inline bool isUpperCase(char ch) { return (ch >= 'A' && ch <= 'Z'); }
int __;
const int INF = 0x3f3f3f3f;
const long long LINF = 0x3f3f3f3f3f3f3f3fLL;
const double EPS = 1e-9;
const int MD = 1000000009;
const int MAXN = 1000000000;
int n, m, p, q, r[111111], cityLen[111111], x, y, z, w[111111], need;
set<pair<int, int> > regions;
pair<int, int> cur1, cur2;
bool was[111111];
long long sum = 0;
int findSet(int v) {
  if (v != r[v]) return r[v] = findSet(r[v]);
  return v;
}
void unite(int a, int b) {
  if ((a = findSet(a)) != (b = findSet(b))) {
    if (rand() & 1) swap(a, b);
    r[a] = b;
  }
}
bool debug;
int main() {
  sIO();
  srand(time(NULL));
  scanf("%d%d%d%d", &n, &m, &p, &q), debug = (n == 1000);
  for (int i = 0; i < n; ++i) r[i] = i;
  memset(cityLen, 0, sizeof(cityLen));
  for (int i = 0; i < m; ++i)
    scanf("%d%d%d", &x, &y, &z), unite(--x, --y),
        cityLen[x] = min(MAXN, cityLen[x] + z),
        cityLen[y] = min(MAXN, cityLen[y] + z);
  memset(w, 0, sizeof(w)), memset(was, 0, sizeof(was));
  for (int i = 0; i < n; ++i)
    r[i] = findSet(r[i]),
    w[r[i]] = min((long long)(MAXN << 1), (long long)w[r[i]] + cityLen[i]),
    was[r[i]] = true;
  for (int i = 0; i < n; ++i)
    if (was[i]) regions.insert(make_pair(w[i] >> 1, i));
  if ((int)regions.size() < q || q == n && p || (int)regions.size() - p > q) {
    printf("NO");
    return 0;
  }
  printf("YES\n"), need = p - (int)regions.size() + q;
  while ((int)regions.size() > q) {
    cur1 = *regions.begin(), regions.erase(cur1), cur2 = *regions.begin(),
    regions.erase(cur2);
    unite(cur1.second, cur2.second);
    printf("%d %d\n", cur1.second + 1, cur2.second + 1);
    sum += (cur1.first + cur2.first + min(MAXN, 1 + cur1.first + cur2.first));
    regions.insert(make_pair(min(MAXN, min(MAXN, (cur1.first << 1)) +
                                           min(MAXN, (cur2.first << 1)) + 1),
                             cur1.second));
  }
  memset(w, -1, sizeof(w));
  for (int i = 0; i < n; ++i) {
    r[i] = findSet(r[i]);
    if (w[r[i]] != -1) {
      for (int j = 0; j < need; ++j) printf("%d %d\n", w[r[i]] + 1, i + 1);
      return 0;
    }
    w[r[i]] = i;
  }
  return 0;
}