Problem: Note that this is the second problem of the two similar problems. You can hack this problem if you solve it. But you can hack the previous problem only if you solve both problems.

You are given a tree with n nodes. In the beginning, 0 is written on all edges. In one operation, you can choose any 2 distinct leaves u, v and any integer number x and add x to values written on all edges on the simple path between u and v. Note that in previous subtask x was allowed to be any real, here it has to be integer.

For example, on the picture below you can see the result of applying two operations to the graph: adding 2 on the path from 7 to 6, and then adding -1 on the path from 4 to 5. 

<image>

You are given some configuration of nonnegative integer pairwise different even numbers, written on the edges. For a given configuration determine if it is possible to achieve it with these operations, and, if it is possible, output the sequence of operations that leads to the given configuration. Constraints on the operations are listed in the output format section.

Leave is a node of a tree of degree 1. Simple path is a path that doesn't contain any node twice.

Input

The first line contains a single integer n (2 ≤ n ≤ 1000) — the number of nodes in a tree.

Each of the next n-1 lines contains three integers u, v, val (1 ≤ u, v ≤ n, u ≠ v, 0 ≤ val ≤ 10 000), meaning that there is an edge between nodes u and v with val written on it. It is guaranteed that these edges form a tree. It is guaranteed that all val numbers are pairwise different and even. 

Output

If there aren't any sequences of operations which lead to the given configuration, output "NO".

If it exists, output "YES" in the first line. In the second line output m — number of operations you are going to apply (0 ≤ m ≤ 10^5). Note that you don't have to minimize the number of the operations!

In the next m lines output the operations in the following format:

u, v, x (1 ≤ u, v ≤ n, u not = v, x — integer, -10^9 ≤ x ≤ 10^9), where u, v — leaves, x — number we are adding. 

It is guaranteed that if there exists a sequence of operations producing given configuration, then there exists a sequence of operations producing given configuration, satisfying all the conditions above.

Examples

Input


5
1 2 2
2 3 4
3 4 10
3 5 18


Output


NO

Input


6
1 2 6
1 3 8
1 4 12
2 5 2
2 6 4


Output


YES
4
3 6 1
4 6 3
3 4 7
4 5 2

Note

The configuration from the first sample is drawn below, and it is impossible to achieve.

<image>

The sequence of operations from the second sample is illustrated below.

<image>

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int MOD = (int)1e9 + 7;
const int FFTMOD = 119 << 23 | 1;
const int INF = (int)1e9 + 23111992;
const long long LINF = (long long)1e18 + 23111992;
const long double PI = acos((long double)-1);
const long double EPS = 1e-9;
inline long long gcd(long long a, long long b) {
  long long r;
  while (b) {
    r = a % b;
    a = b;
    b = r;
  }
  return a;
}
inline long long lcm(long long a, long long b) { return a / gcd(a, b) * b; }
inline long long fpow(long long n, long long k, int p = MOD) {
  long long r = 1;
  for (; k; k >>= 1) {
    if (k & 1) r = r * n % p;
    n = n * n % p;
  }
  return r;
}
template <class T>
inline int chkmin(T& a, const T& val) {
  return val < a ? a = val, 1 : 0;
}
template <class T>
inline int chkmax(T& a, const T& val) {
  return a < val ? a = val, 1 : 0;
}
inline unsigned long long isqrt(unsigned long long k) {
  unsigned long long r = sqrt(k) + 1;
  while (r * r > k) r--;
  return r;
}
inline long long icbrt(long long k) {
  long long r = cbrt(k) + 1;
  while (r * r * r > k) r--;
  return r;
}
inline void addmod(int& a, int val, int p = MOD) {
  if ((a = (a + val)) >= p) a -= p;
}
inline void submod(int& a, int val, int p = MOD) {
  if ((a = (a - val)) < 0) a += p;
}
inline int mult(int a, int b, int p = MOD) { return (long long)a * b % p; }
inline int inv(int a, int p = MOD) { return fpow(a, p - 2, p); }
inline int sign(long double x) { return x < -EPS ? -1 : x > +EPS; }
inline int sign(long double x, long double y) { return sign(x - y); }
mt19937 mt(chrono::high_resolution_clock::now().time_since_epoch().count());
inline int mrand() { return abs((int)mt()); }
void solve() {
  int n;
  cin >> n;
  vector<vector<pair<int, int>>> adj(n);
  vector<int> ws(n);
  for (int i = (0); i < (n - 1); ++i) {
    int u, v;
    cin >> u >> v >> ws[i];
    u--, v--;
    adj[u].push_back(make_pair(v, i)), adj[v].push_back(make_pair(u, i));
  }
  for (int u = (0); u < (n); ++u)
    if (int((adj[u]).size()) == 2) {
      cout << "NO\n";
      return;
    }
  function<int(int, int)> find = [&](int u, int p) {
    if (int((adj[u]).size()) == 1) return u;
    for (auto e : adj[u]) {
      int v = e.first;
      if (v ^ p) {
        return find(v, u);
      }
    }
    assert(0);
    return -1;
  };
  function<int(int, int, int, int)> upd = [&](int u, int p, int v, int w) {
    if (u == v) {
      return 1;
    }
    for (int i = (0); i < (int((adj[u]).size())); ++i) {
      int nu = adj[u][i].first;
      int ix = adj[u][i].second;
      if (nu != p && upd(nu, u, v, w)) {
        ws[ix] -= w;
        return 1;
      }
    }
    return 0;
  };
  vector<tuple<int, int, int>> res;
  function<void(int, int)> dfs = [&](int u, int p) {
    for (int i = (0); i < (int((adj[u]).size())); ++i) {
      int v = adj[u][i].first;
      int ix = adj[u][i].second;
      if (v ^ p) {
        vector<int> ver;
        for (auto e2 : adj[u]) {
          int v2 = e2.first;
          if (v2 != v) {
            ver.push_back(v2);
          }
        }
        int x = ver[0];
        int y = ver[1];
        int t1 = find(v, u);
        int t2 = find(x, u);
        int t3 = find(y, u);
        int t = ws[ix];
        res.push_back(make_tuple(t1, t2, t / 2));
        res.push_back(make_tuple(t1, t3, t / 2));
        res.push_back(make_tuple(t2, t3, -t / 2));
        assert(upd(t1, -1, t2, t / 2));
        assert(upd(t1, -1, t3, t / 2));
        assert(upd(t2, -1, t3, -t / 2));
        dfs(v, u);
      }
    }
  };
  if (n == 2) {
    res.push_back(make_tuple(0, 1, ws[adj[0][0].second]));
  } else {
    for (int rt = (0); rt < (n); ++rt)
      if (3 <= int((adj[rt]).size())) {
        dfs(rt, -1);
        break;
      }
  }
  cout << "YES\n";
  cout << int((res).size()) << "\n";
  for (auto e : res) {
    int u, v, w;
    tie(u, v, w) = e;
    cout << u + 1 << " " << v + 1 << " " << w << "\n";
  }
}
int main(int argc, char* argv[]) {
  ios_base::sync_with_stdio(0), cin.tie(0);
  if (argc > 1) {
    assert(freopen(argv[1], "r", stdin));
  }
  if (argc > 2) {
    assert(freopen(argv[2], "wb", stdout));
  }
  solve();
  cerr << "\nTime elapsed: " << 1000 * clock() / CLOCKS_PER_SEC << "ms\n";
  return 0;
}