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 dx[4] = {1, 0, -1, 0}, dy[4] = {0, 1, 0, -1};
const long long linf = 4000000000000000000LL;
const long long inf = 1000000007;
const long double pi = 3.1415926535;
void pv(vector<int> a) {
  for (auto& x : a) cout << x << " ";
  cout << '\n';
}
void pv(vector<long long> a) {
  for (auto& x : a) cout << x << " ";
  cout << '\n';
}
void pv(vector<vector<int>> a) {
  for (int i = (0); i < (int(a.size())); ++i) {
    cout << i << '\n';
    pv(a[i]);
    cout << '\n';
  }
}
void pv(vector<vector<long long>> a) {
  for (int i = (0); i < (int(a.size())); ++i) {
    cout << i << '\n';
    pv(a[i]);
  }
  cout << '\n';
}
void pv(vector<string> a) {
  for (auto& x : a) cout << x << '\n';
  cout << '\n';
}
void build_primes(vector<int>& primes, int size) {
  vector<int> visited;
  visited.resize(size, 0);
  for (int i = (2); i < (size); ++i) {
    if (visited[i] == 0) {
      primes.push_back(i);
      int a = i;
      while (a < size) {
        visited[a] = 1;
        a += i;
      }
    }
  }
}
vector<vector<long long>> matrix_mult(vector<vector<long long>>& a,
                                      vector<vector<long long>>& b) {
  int n = a.size();
  vector<vector<long long>> answer;
  answer.resize(n);
  for (int i = 0; i < n; i++) answer[i].resize(n, 0);
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < n; j++) {
      for (int k = 0; k < n; k++)
        answer[i][j] = (answer[i][j] + a[i][k] * b[k][j]) % inf;
    }
  }
  return answer;
}
int modInverse(int a, int m) {
  int m0 = m;
  int y = 0, x = 1;
  if (m == 1) return 0;
  while (a > 1) {
    int q = a / m;
    int t = m;
    m = a % m, a = t;
    t = y;
    y = x - q * y;
    x = t;
  }
  if (x < 0) x += m0;
  return x;
}
long long power(long long x, long long y) {
  long long k = 1LL << 60;
  long long z = 1;
  while (k != 0) {
    z *= z;
    z %= inf;
    if (y >= k) {
      z *= x;
      z %= inf;
      y -= k;
    }
    k >>= 1;
  }
  return z;
}
struct point {
  long double x, y;
  bool operator<(const point& rhs) const {
    if (x == rhs.x) return y < rhs.y;
    return x < rhs.x;
  }
};
struct pt {
  long long x, y;
  bool operator<(const point& rhs) const {
    if (x == rhs.x) return y < rhs.y;
    return x < rhs.x;
  }
};
long double area(point x, point y, point z) {
  return (x.y * y.x + y.y * z.x + z.y * x.x - x.x * y.y - y.x * z.y -
          z.x * x.y) /
         2.0;
}
bool clockwise(point x, point y, point z) { return area(x, y, z) > 0; }
long double area(pt x, pt y, pt z) {
  return (x.y * y.x + y.y * z.x + z.y * x.x - x.x * y.y - y.x * z.y -
          z.x * x.y) /
         2.0;
}
bool clockwise(pt x, pt y, pt z) { return area(x, y, z) > 0; }
long long gcd(long long a, long long b) {
  if (a > b) swap(a, b);
  if (a == 0) return b;
  return gcd(a, b % a);
}
int popcount(long long a) {
  int count = 0;
  while (a) {
    count += (a & 1);
    a >>= 1;
  }
  return count;
}
long long choose(long long n, long long r) {
  long long p = 1, k = 1;
  if (n - r < r) r = n - r;
  if (r != 0) {
    while (r) {
      p *= n;
      k *= r;
      long long m = gcd(p, k);
      p /= m;
      k /= m;
      n--;
      r--;
    }
  } else
    p = 1;
  return p;
}
vector<long long> prefix_hash(string& a, vector<long long>& powers) {
  int n = int(a.size());
  vector<long long> prefix(n + 1);
  for (int i = (0); i < (n); ++i)
    prefix[i + 1] = (prefix[i] + powers[i] * (a[i] - 'a' + 1)) % inf;
  return prefix;
}
struct custom_hash {
  static uint64_t splitmix64(uint64_t x) {
    x += 0x9e3779b97f4a7c15;
    x = (x ^ (x >> 30)) * 0xbf58476d1ce4e5b9;
    x = (x ^ (x >> 27)) * 0x94d049bb133111eb;
    return x ^ (x >> 31);
  }
  uint64_t operator()(uint64_t x) const {
    static const uint64_t FIXED_RANDOM =
        chrono::steady_clock::now().time_since_epoch().count();
    return splitmix64(x + FIXED_RANDOM);
  }
};
struct custom_hash_fast {
  uint64_t operator()(uint64_t x) const {
    static const uint64_t FIXED_RANDOM =
        chrono::steady_clock::now().time_since_epoch().count();
    x ^= FIXED_RANDOM;
    return x ^ (x >> 16);
  }
};
void setIO(string second) {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
}
int n;
vector<vector<pair<int, int>>> adj;
vector<int> t;
vector<int> tin, tout;
int timenow;
void dfs(int i, int p, int val) {
  if (int(t.size()) == val) return;
  if (int(adj[i].size()) == 1) {
    t.push_back(i);
  }
  for (auto& x : adj[i]) {
    if (x.first != p) {
      dfs(x.first, i, val);
    }
  }
}
void preorder(int i, int p) {
  tin[i] = timenow++;
  for (auto& x : adj[i]) {
    if (x.first != p) {
      preorder(x.first, i);
    }
  }
  tout[i] = timenow++;
}
struct type {
  int a, b, c;
};
vector<type> ans;
void add(int v, int leaf, int val) {
  if (int(adj[v].size()) == 1) {
    ans.push_back({1 + v, 1 + leaf, val});
    return;
  }
  t.resize(0);
  preorder(v, -1);
  for (auto& x : adj[v]) {
    if (int(t.size()) == 2) break;
    if (tin[x.first] <= tin[leaf] && tout[x.first] >= tout[leaf]) continue;
    dfs(x.first, v, int(t.size()) + 1);
  }
  ans.push_back({t[0] + 1, leaf + 1, val / 2});
  ans.push_back({t[1] + 1, leaf + 1, val / 2});
  ans.push_back({t[0] + 1, t[1] + 1, -val / 2});
}
int main() {
  setIO("");
  cin >> n;
  adj.resize(n);
  for (int i = (0); i < (n - 1); ++i) {
    int a, b, c;
    cin >> a >> b >> c;
    adj[--a].push_back({--b, c});
    adj[b].push_back({a, c});
  }
  for (int i = (0); i < (n); ++i) {
    if (int(adj[i].size()) == 2) {
      cout << "NO" << '\n';
      return 0;
    }
  }
  cout << "YES" << '\n';
  tin.resize(n);
  tout.resize(n);
  for (int i = (0); i < (n); ++i) {
    for (auto& x : adj[i]) {
      if (x.first > i) continue;
      if (int(adj[x.first].size()) == 1) {
        add(i, x.first, x.second);
      } else {
        t.resize(0);
        dfs(x.first, i, 1);
        int node = t[0];
        add(i, node, x.second);
        add(x.first, node, -x.second);
      }
    }
  }
  cout << int(ans.size()) << '\n';
  for (auto& x : ans) cout << x.a << " " << x.b << " " << x.c << '\n';
}