#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
ll N, Q;
// w, c, d
vector<vector<tuple<ll, ll, ll>>> Adj;
const ll LEVEL = 20;
vector<vector<ll>> Anc;
vector<ll> Dist, Depth;
// query id, coefficient, x, y
vector<vector<tuple<ll, ll, ll, ll>>> Calc;
// sum of lengths, number
vector<pair<ll, ll>> ColorInfo;
vector<ll> R;

void init_dfs(ll v, ll par, ll depth, ll dist) {
  Anc[v][0] = par;
  Depth[v] = depth;
  Dist[v] = dist;
  for (auto edge : Adj[v]) {
    ll w = get<0>(edge);
    if (w == par) continue;
    ll d = get<2>(edge);
    init_dfs(w, v, depth + 1, dist + d);
  }
}
ll lca(ll u, ll v) {
  if (Depth[u] > Depth[v]) swap(u, v);
  for (ll k = LEVEL - 1; k >= 0; --k) {
    if ((Depth[v] - Depth[u]) & (1 << k)) v = Anc[v][k];
  }
  if (u == v) return u;
  for (ll k = LEVEL - 1; k >= 0; --k) {
    if (Anc[u][k] != Anc[v][k]) {
      u = Anc[u][k];
      v = Anc[v][k];
    }
  }
  return Anc[u][0];
}
void dfs(ll v, ll par) {
  for (auto calc : Calc[v]) {
    ll i = get<0>(calc), coeff = get<1>(calc);
    ll x = get<2>(calc), y = get<3>(calc);
    R[i] += coeff * (Dist[v] - ColorInfo[x].first + y * ColorInfo[x].second);
  }
  for (auto edge : Adj[v]) {
    ll w = get<0>(edge);
    if (w == par) continue;
    ll c = get<1>(edge), d = get<2>(edge);
    ColorInfo[c].first += d;
    ++ColorInfo[c].second;
    dfs(w, v);
    ColorInfo[c].first -= d;
    --ColorInfo[c].second;
  }
}
int main() {
  cin >> N >> Q;
  Adj.resize(N);
  for (ll i = 0; i < N - 1; ++i) {
    ll a, b, c, d;
    cin >> a >> b >> c >> d;
    --a;
    --b;
    --c;
    Adj[a].push_back(make_tuple(b, c, d));
    Adj[b].push_back(make_tuple(a, c, d));
  }
  Anc.resize(N);
  for (auto &anc : Anc) anc.resize(LEVEL, -1);
  Depth.resize(N);
  Dist.resize(N);
  init_dfs(0, -1, 0, 0);
  for (ll k = 1; k < LEVEL; ++k) {
    for (ll v = 0; v < N; ++v) {
      if (Anc[v][k - 1] != -1) Anc[v][k] = Anc[Anc[v][k - 1]][k - 1];
    }
  }
  Calc.resize(N);
  for (ll i = 0; i < Q; ++i) {
    ll x, y, u, v;
    cin >> x >> y >> u >> v;
    --x;
    --u;
    --v;
    ll w = lca(u, v);
    Calc[u].push_back(make_tuple(i, 1, x, y));
    Calc[v].push_back(make_tuple(i, 1, x, y));
    Calc[w].push_back(make_tuple(i, -2, x, y));
  }
  ColorInfo.resize(N - 1);
  R.resize(Q);
  dfs(0, -1);
  for (ll r : R) cout << r << endl;
  return 0;
}

### Prompt

Your challenge is to write a cpp solution to the following problem:

There is a tree with N vertices numbered 1 to N. The i-th edge in this tree connects Vertex a_i and Vertex b_i, and the color and length of that edge are c_i and d_i, respectively. Here the color of each edge is represented by an integer between 1 and N-1 (inclusive). The same integer corresponds to the same color, and different integers correspond to different colors.

Answer the following Q queries:

* Query j (1 \leq j \leq Q): assuming that the length of every edge whose color is x_j is changed to y_j, find the distance between Vertex u_j and Vertex v_j. (The changes of the lengths of edges do not affect the subsequent queries.)

Constraints

* 2 \leq N \leq 10^5
* 1 \leq Q \leq 10^5
* 1 \leq a_i, b_i \leq N
* 1 \leq c_i \leq N-1
* 1 \leq d_i \leq 10^4
* 1 \leq x_j \leq N-1
* 1 \leq y_j \leq 10^4
* 1 \leq u_j < v_j \leq N
* The given graph is a tree.
* All values in input are integers.

Input

Input is given from Standard Input in the following format:


N Q
a_1 b_1 c_1 d_1
:
a_{N-1} b_{N-1} c_{N-1} d_{N-1}
x_1 y_1 u_1 v_1
:
x_Q y_Q u_Q v_Q


Output

Print Q lines. The j-th line (1 \leq j \leq Q) should contain the answer to Query j.

Example

Input

5 3
1 2 1 10
1 3 2 20
2 4 4 30
5 2 1 40
1 100 1 4
1 100 1 5
3 1000 3 4


Output

130
200
60

### Response

```cpp
#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
ll N, Q;
// w, c, d
vector<vector<tuple<ll, ll, ll>>> Adj;
const ll LEVEL = 20;
vector<vector<ll>> Anc;
vector<ll> Dist, Depth;
// query id, coefficient, x, y
vector<vector<tuple<ll, ll, ll, ll>>> Calc;
// sum of lengths, number
vector<pair<ll, ll>> ColorInfo;
vector<ll> R;

void init_dfs(ll v, ll par, ll depth, ll dist) {
  Anc[v][0] = par;
  Depth[v] = depth;
  Dist[v] = dist;
  for (auto edge : Adj[v]) {
    ll w = get<0>(edge);
    if (w == par) continue;
    ll d = get<2>(edge);
    init_dfs(w, v, depth + 1, dist + d);
  }
}
ll lca(ll u, ll v) {
  if (Depth[u] > Depth[v]) swap(u, v);
  for (ll k = LEVEL - 1; k >= 0; --k) {
    if ((Depth[v] - Depth[u]) & (1 << k)) v = Anc[v][k];
  }
  if (u == v) return u;
  for (ll k = LEVEL - 1; k >= 0; --k) {
    if (Anc[u][k] != Anc[v][k]) {
      u = Anc[u][k];
      v = Anc[v][k];
    }
  }
  return Anc[u][0];
}
void dfs(ll v, ll par) {
  for (auto calc : Calc[v]) {
    ll i = get<0>(calc), coeff = get<1>(calc);
    ll x = get<2>(calc), y = get<3>(calc);
    R[i] += coeff * (Dist[v] - ColorInfo[x].first + y * ColorInfo[x].second);
  }
  for (auto edge : Adj[v]) {
    ll w = get<0>(edge);
    if (w == par) continue;
    ll c = get<1>(edge), d = get<2>(edge);
    ColorInfo[c].first += d;
    ++ColorInfo[c].second;
    dfs(w, v);
    ColorInfo[c].first -= d;
    --ColorInfo[c].second;
  }
}
int main() {
  cin >> N >> Q;
  Adj.resize(N);
  for (ll i = 0; i < N - 1; ++i) {
    ll a, b, c, d;
    cin >> a >> b >> c >> d;
    --a;
    --b;
    --c;
    Adj[a].push_back(make_tuple(b, c, d));
    Adj[b].push_back(make_tuple(a, c, d));
  }
  Anc.resize(N);
  for (auto &anc : Anc) anc.resize(LEVEL, -1);
  Depth.resize(N);
  Dist.resize(N);
  init_dfs(0, -1, 0, 0);
  for (ll k = 1; k < LEVEL; ++k) {
    for (ll v = 0; v < N; ++v) {
      if (Anc[v][k - 1] != -1) Anc[v][k] = Anc[Anc[v][k - 1]][k - 1];
    }
  }
  Calc.resize(N);
  for (ll i = 0; i < Q; ++i) {
    ll x, y, u, v;
    cin >> x >> y >> u >> v;
    --x;
    --u;
    --v;
    ll w = lca(u, v);
    Calc[u].push_back(make_tuple(i, 1, x, y));
    Calc[v].push_back(make_tuple(i, 1, x, y));
    Calc[w].push_back(make_tuple(i, -2, x, y));
  }
  ColorInfo.resize(N - 1);
  R.resize(Q);
  dfs(0, -1);
  for (ll r : R) cout << r << endl;
  return 0;
}

```