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 GCC optimize("O0")
#pragma GCC optimize("Ofast")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
#pragma GCC optimize("-ffloat-store")
using namespace std;
long long n, m, a, b, k, temp, x, y;
const int MAXN = 100000 + 11;
inline void read(vector<long long> &v);
inline long long max(long long a, long long b);
inline long long min(long long a, long long b);
inline long long gcd(long long a, long long b);
class DisjSet {
  int *rank, *parent, n, *sze, comp;
  long long int *roadsum;

 public:
  DisjSet(int n) {
    rank = new int[n];
    parent = new int[n];
    sze = new int[n];
    comp = n;
    roadsum = new long long int[n];
    this->n = n;
    makeSet();
  }
  void makeSet() {
    for (int i = 0; i < n; i++) {
      parent[i] = i;
      sze[i] = 1;
      roadsum[i] = 0;
    }
  }
  int find(int x) {
    if (parent[x] != x) {
      parent[x] = find(parent[x]);
    }
    return parent[x];
  }
  int getsize(int x) { return sze[find(x)]; }
  long long getsum(int x) { return roadsum[find(x)]; }
  void Union(int x, int y, long long int len) {
    int xset = find(x);
    int yset = find(y);
    roadsum[yset] += len;
    if (xset == yset) return;
    comp--;
    if (rank[xset] < rank[yset]) {
      parent[xset] = yset;
      sze[yset] += sze[xset];
      roadsum[yset] += roadsum[xset];
    } else if (rank[xset] > rank[yset]) {
      parent[yset] = xset;
      sze[xset] += sze[yset];
      roadsum[xset] += roadsum[yset];
    } else {
      parent[yset] = xset;
      sze[xset] += sze[yset];
      roadsum[xset] += roadsum[yset];
      rank[xset] = rank[xset] + 1;
    }
  }
  int getcomps() { return comp; }
};
void solveforthiscase(const int &test) {
  int p, q;
  cin >> n >> m >> p >> q;
  class DisjSet c(n);
  while (m--) {
    long long l;
    cin >> a >> b >> l;
    c.Union(a - 1, b - 1, l);
  }
  if (c.getcomps() < q || (p && q == n)) {
    cout << "NO\n";
    return;
  }
  set<pair<long long, long long> > s;
  for (int i = 0; i < n; i++) s.insert({c.getsum(i), c.find(i)});
  vector<pair<long long, long long> > ans;
  while (s.size() > q && p) {
    p--;
    if (false) cerr << "p" << ' ' << '=' << ' ' << p << '\n';
    int u = (*s.begin()).second;
    long long c1 = (*s.begin()).first;
    s.erase(s.begin());
    int v = (*s.begin()).second;
    long long c2 = (*s.begin()).first;
    s.erase(s.begin());
    if (false) cerr << "u + 1" << ' ' << '=' << ' ' << u + 1 << '\n';
    if (false) cerr << "v + 1" << ' ' << '=' << ' ' << v + 1 << '\n';
    c.Union(u, v, min(c1 + c2 + 1, 1e9));
    s.insert({c.getsum(u), c.find(u)});
    ans.push_back({u + 1, v + 1});
    if (false) cerr << "p" << ' ' << '=' << ' ' << p << '\n';
    if (false) cerr << "s.size()" << ' ' << '=' << ' ' << s.size() << '\n';
  }
  while (p--) {
    if (ans.size())
      ans.push_back(ans.back());
    else {
      vector<int> vis(MAXN, -1);
      for (int i = 0; i < n; i++)
        if (vis[c.find(i)] == -1)
          vis[c.find(i)] = i + 1;
        else {
          ans.push_back({i + 1, vis[c.find(i)]});
          break;
        }
    }
  }
  if (s.size() != q) {
    cout << "NO\n";
    return;
  }
  cout << "YES\n";
  for (auto x : ans) cout << x.first << ' ' << x.second << '\n';
  return;
}
int main() {
  ios_base::sync_with_stdio(false);
  cout << fixed << setprecision(10);
  cin.tie(NULL);
  cout.tie(NULL);
  clock_t start, end;
  start = clock();
  int test = 1;
  for (int i = 1; i <= test; i++) {
    solveforthiscase(i);
  }
  end = clock();
  long double time_taken =
      (long double)(end - start) / (long double)(CLOCKS_PER_SEC);
  if (false)
    cerr << "Time taken by program is : " << (long double)time_taken * 1000
         << " milisec ";
}
inline long long max(long long a, long long b) { return ((a > b) ? a : b); }
inline long long min(long long a, long long b) { return ((a > b) ? b : a); }
inline long long gcd(long long a, long long b) {
  if (b == 0) return a;
  a %= b;
  return gcd(b, a);
}
inline void read(vector<long long> &v) {
  for (int i = 0; i < v.size(); i++) cin >> v[i];
}