Language: C++

#include <bits/stdc++.h>
using namespace std;
#ifndef ONLINE_JUDGE
	#include <boost/multiprecision/cpp_int.hpp>
	#include <boost/multiprecision/cpp_dec_float.hpp>
	using bll = boost::multiprecision::cpp_int;
	using bdouble = boost::multiprecision::cpp_dec_float_100;
#endif
#ifdef LOCAL_DEV
	void debug_impl() { std::cerr << std::endl; }
	template<typename Head, typename... Tail> void debug_impl(Head head, Tail... tail) { std::cerr << " " << head << (sizeof...(tail) ? "," : ""); debug_impl(tail...); }
	#define debug(...) { std::cerr << std::boolalpha << "(" << #__VA_ARGS__ << ") ="; debug_impl(__VA_ARGS__); std::cerr << std::noboolalpha; }
#else
	#define debug(...) {}
#endif
#ifdef LOCAL_TEST
	#define BOOST_STACKTRACE_USE_ADDR2LINE
	#define BOOST_STACKTRACE_ADDR2LINE_LOCATION /usr/local/opt/binutils/bin/addr2line
	#define _GNU_SOURCE
	#include <boost/stacktrace.hpp>
	namespace std {
		template<typename T> class dvector : public std::vector<T> {
		public:
			dvector() : std::vector<T>() {}
			explicit dvector(size_t n, const T& value = T()) : std::vector<T>(n, value) {}
			dvector(const std::vector<T>& v) : std::vector<T>(v) {}
			dvector(const std::initializer_list<T> il) : std::vector<T>(il) {}
			dvector(const typename std::vector<T>::iterator first, const typename std::vector<T>::iterator last) : std::vector<T>(first, last) {}
			dvector(const std::string::iterator first, const std::string::iterator last) : std::vector<T>(first, last) {}
			T& operator[](size_t n) {
				try { return this->at(n); } catch (const std::exception& e) {
					std::cerr << boost::stacktrace::stacktrace() << std::endl; return this->at(n);
				}
			}
			const T& operator[](size_t n) const {
				try { return this->at(n); } catch (const std::exception& e) {
					std::cerr << boost::stacktrace::stacktrace() << std::endl; return this->at(n);
				}
			}
		};
	}
	class dbool {
	private:
		bool boolvalue;
	public:
		dbool() : boolvalue(false) {}
		dbool(bool b) : boolvalue(b) {}
		dbool(const dbool &b) : boolvalue(b.boolvalue) {}
		operator bool&() { return boolvalue; }
		operator const bool&() const { return boolvalue; }
	};
	template<typename T> std::ostream& operator<<(std::ostream& s, const dvector<T>& v) {
		for (int i = 0, len = v.size(); i < len; ++i){ s << v[i]; if (i < len - 1) s << "\t"; } return s; }
	template<typename T> std::ostream& operator<<(std::ostream& s, const dvector< dvector<T> >& vv) {
		for (int i = 0, len = vv.size(); i < len; ++i){ s << vv[i] << std::endl; } return s; }
	template<typename T> std::ostream& operator<<(std::ostream& s, const std::set<T>& se) {
		s << "{ "; for (auto itr = se.begin(); itr != se.end(); ++itr){ s << (*itr) << "\t"; } s << "}"; return s; }
	template<typename T> std::ostream& operator<<(std::ostream& s, const std::multiset<T>& se) {
		s << "{ "; for (auto itr = se.begin(); itr != se.end(); ++itr){ s << (*itr) << "\t"; } s << "}"; return s; }
	template<typename T1, typename T2> std::ostream& operator<<(std::ostream& s, const std::map<T1, T2>& m) {
		s << "{" << std::endl; for (auto itr = m.begin(); itr != m.end(); ++itr){ s << "\t" << (*itr).first << " : " << (*itr).second << std::endl; } s << "}"; return s; }
	template<typename T1, typename T2> std::ostream& operator<<(std::ostream& s, const std::pair<T1, T2>& p) { 
		return s << "(" << p.first << ", " << p.second << ")"; }
	#define vector dvector
	#define bool dbool
	class SIGFPE_exception : std::exception {};
	class SIGSEGV_exception : std::exception {};
	void catch_SIGFPE(int e) { std::cerr << boost::stacktrace::stacktrace() << std::endl; throw SIGFPE_exception();	}
	void catch_SIGSEGV(int e) { std::cerr << boost::stacktrace::stacktrace() << std::endl; throw SIGSEGV_exception(); }
	signed convertedmain();
	signed main() { signal(SIGFPE, catch_SIGFPE); signal(SIGSEGV, catch_SIGSEGV); return convertedmain(); }
	#define main() convertedmain()
#endif
//#define int long long
using ll = long long;
//constexpr int INF = 1e9;//INT_MAX=(1<<31)-1=2147483647
constexpr ll INF = (ll)1e18;//(1LL<<63)-1=9223372036854775807
constexpr ll MOD = (ll)1e9 + 7;
constexpr double EPS = 1e-9;
constexpr ll dx[4] = {1LL, 0LL, -1LL, 0LL};
constexpr ll dy[4] = {0LL, 1LL, 0LL, -1LL};
constexpr ll dx8[8] = {1LL, 0LL, -1LL, 0LL, 1LL, 1LL, -1LL, -1LL};
constexpr ll dy8[8] = {0LL, 1LL, 0LL, -1LL, 1LL, -1LL, 1LL, -1LL};
#define rep(i, n)   for(ll i=0, i##_length=(n); i< i##_length; ++i)
#define repeq(i, n) for(ll i=1, i##_length=(n); i<=i##_length; ++i)
#define rrep(i, n)   for(ll i=(n)-1; i>=0; --i)
#define rrepeq(i, n) for(ll i=(n)  ; i>=1; --i)
#define all(v) (v).begin(), (v).end()
#define rall(v) (v).rbegin(), (v).rend()
void p() { std::cout << '\n'; }
template<typename Head, typename... Tail> void p(Head head, Tail... tail) { std::cout << head << (sizeof...(tail) ? " " : ""); p(tail...); }
template<typename T> inline void pv(std::vector<T> &v) { for(ll i=0, N=v.size(); i<N; i++) std::cout << v[i] << " \n"[i==N-1]; }
template<typename T> inline T gcd(T a, T b) { return b ? gcd(b,a%b) : a; }
template<typename T> inline T lcm(T a, T b) { return a / gcd(a,  b) * b; }
template<typename T> inline bool chmax(T &a, T b) { return a < b && (a = b, true); }
template<typename T> inline bool chmin(T &a, T b) { return a > b && (a = b, true); }
template<typename T> inline void uniq(std::vector<T> &v) { v.erase(unique(v.begin(), v.end()), v.end()); }

/*-----8<-----template-----8<-----*/

/*
実装では木を 1-indexed の配列で表現している。
ノード k について, 親ノードは k/2, 子ノードは 2k, 2k+1 である。

LazySegmentTree(n, f, g, h, M1, OM0):= サイズ n の初期化。
	f は2つの区間の要素をマージする二項演算, 
	g は要素と作用素をマージする二項演算(第三引数は対応する区間の長さ), 
	h は作用素同士をマージする二項演算, M1 はモノイドの単位元, 
	OM0 は作用素の単位元である。
set(k,x):= k 番目の要素に x を代入する。
build():= セグメント木を構築する。
query(a,b):= 区間 [a,b) に対して二項演算した結果を返す。
update(a, b, x) := 区間 [a,b) に作用素 x を適用する。
operator[k] := k 番目の要素を返す。
find_first(a, check) := [a,x) が check を満たす最初の要素位置 x を返す。見つからない場合は -1 を返す。
find_last(b, check) := [x,b) が check を満たす最後の要素位置 x を返す。見つからない場合は -1 を返す。
*/
template< typename Monoid, typename OperatorMonoid, typename F , typename G , typename H >
class LazySegmentTree {
public:
	ll sz, height;
	vector< Monoid > data;
	vector< OperatorMonoid > lazy;
	const F f;
	const G g;
	const H h;
	const Monoid M1;
	const OperatorMonoid OM0;

	LazySegmentTree(ll n, const F f, const G g, const H h, const Monoid &M1, const OperatorMonoid OM0)
			: f(f), g(g), h(h), M1(M1), OM0(OM0) {
		sz = 1;
		height = 0;
		while(sz < n) sz <<= 1, height++;
		data.assign(2 * sz, M1);
		lazy.assign(2 * sz, OM0);
	}

	void set(ll k, const Monoid &x) {
		data[k + sz] = x;
	}

	void build() {
		for(ll k = sz - 1; k > 0; k--) {
			data[k] = f(data[2 * k + 0], data[2 * k + 1]);
		}
	}

	inline void propagate(ll k) {
		if(lazy[k] != OM0) {
			lazy[2 * k + 0] = h(lazy[2 * k + 0], lazy[k]);
			lazy[2 * k + 1] = h(lazy[2 * k + 1], lazy[k]);
			data[k] = reflect(k);
			lazy[k] = OM0;
		}
	}

	inline Monoid reflect(ll k) {
		return lazy[k] == OM0 ? data[k] : g(data[k], lazy[k]);
	}

	inline void recalc(ll k) {
		while(k >>= 1) data[k] = f(reflect(2 * k + 0), reflect(2 * k + 1));
	}

	inline void thrust(ll k) {
		for(ll i = height; i > 0; i--) propagate(k >> i);
	}

	void update(ll a, ll b, const OperatorMonoid &x) {
		thrust(a += sz);
		thrust(b += sz - 1);
		for(ll l = a, r = b + 1; l < r; l >>= 1, r >>= 1) {
			if(l & 1) lazy[l] = h(lazy[l], x), ++l;
			if(r & 1) --r, lazy[r] = h(lazy[r], x);
		}
		recalc(a);
		recalc(b);
	}

	Monoid query(ll a, ll b) {
		thrust(a += sz);
		thrust(b += sz - 1);
		Monoid L = M1, R = M1;
		for(ll l = a, r = b + 1; l < r; l >>= 1, r >>= 1) {
			if(l & 1) L = f(L, reflect(l++));
			if(r & 1) R = f(reflect(--r), R);
		}
		return f(L, R);
	}

	Monoid operator[](const ll &k) {
		return query(k, k + 1);
	}

	template< typename C >
	ll find_subtree(ll a, const C &check, Monoid &M, bool type) {
		while(a < sz) {
			propagate(a);
			Monoid nxt = type ? f(reflect(2 * a + type), M) : f(M, reflect(2 * a + type));
			if(check(nxt)) a = 2 * a + type;
			else M = nxt, a = 2 * a + 1 - type;
		}
		return a - sz;
	}

	template< typename C >
	ll find_first(ll a, const C &check) {
		Monoid L = M1;
		if(a <= 0) {
			if(check(f(L, reflect(1)))) return find_subtree(1, check, L, false);
			return -1;
		}
		thrust(a + sz);
		ll b = sz;
		for(a += sz, b += sz; a < b; a >>= 1, b >>= 1) {
			if(a & 1) {
				Monoid nxt = f(L, reflect(a));
				if(check(nxt)) return find_subtree(a, check, L, false);
				L = nxt;
				++a;
			}
		}
		return -1;
	}

	template< typename C >
	ll find_last(ll b, const C &check) {
		Monoid R = M1;
		if(b >= sz) {
			if(check(f(reflect(1), R))) return find_subtree(1, check, R, true);
			return -1;
		}
		thrust(b + sz - 1);
		ll a = sz;
		for(b += sz; a < b; a >>= 1, b >>= 1) {
			if(b & 1) {
				Monoid nxt = f(reflect(--b), R);
				if(check(nxt)) return find_subtree(b, check, R, true);
				R = nxt;
			}
		}
		return -1;
	}
};

/*
最小共通祖先(Doubling-Lowest-Common-Ancestor)
DoublingLowestCommonAncestor(g):= 木 g で初期化する。
build(v):= vを根としたLCAを構築する。
query(u, v):= 頂点 u と v の最小共通祖先を返す。
*/
template< typename G >
class DoublingLowestCommonAncestor {
public:
	const G &g;
	vector< int > depth;//depth[v]=根からvまでの距離
	const int LOG;
	vector< vector< int > > table;

	DoublingLowestCommonAncestor(const G &g) : g(g), depth(g.size()), LOG(32 - __builtin_clz(g.size())) {
		table.assign(LOG, vector< int >(g.size(), -1));
	}

	void dfs(int idx, int par, int d) {
		table[0][idx] = par;
		depth[idx] = d;
		for(auto &to : g[idx]) {
			if(to != par) dfs(to, idx, d + 1);
		}
	}

	void build(int v) {
		dfs(v, -1, 0);
		for(int k = 0; k + 1 < LOG; k++) {
			for(int i = 0,tablesize = table[k].size(); i < tablesize; i++) {
				if(table[k][i] == -1) table[k + 1][i] = -1;
				else table[k + 1][i] = table[k][table[k][i]];
			}
		}
	}

	int query(int u, int v) {
		if(depth[u] > depth[v]) swap(u, v);
		for(int i = LOG - 1; i >= 0; i--) {
			if(((depth[v] - depth[u]) >> i) & 1) v = table[i][v];
		}
		if(u == v) return u;
		for(int i = LOG - 1; i >= 0; i--) {
			if(table[i][u] != table[i][v]) {
				u = table[i][u];
				v = table[i][v];
			}
		}
		return table[0][u];
	}
};

////

class EulerTourForEdge{
private:
	vector<ll> downindex, upindex, depth, bottomindex;

	void dfs(ll v,ll p,ll d){
		depth[v]=d;
		for(ll u:g[v]){
			if(u==p) continue;
			downindex[u]=bottomindex.size();
			bottomindex.emplace_back(u);
			dfs(u,v,d+1);
			upindex[u]=bottomindex.size();
			bottomindex.emplace_back(u);
		}
	}
public:
	vector<vector<ll>> g;

	EulerTourForEdge() {}
	EulerTourForEdge(vector<vector<ll>> &g): downindex(g.size()), upindex(g.size()), depth(g.size()), g(g) {}
	EulerTourForEdge(ll n): downindex(n), upindex(n), depth(n), g(n) {}
	
	//u,v間に無向辺を追加
	void add_edge(ll u, ll v){
		g[u].emplace_back(v);
		g[v].emplace_back(u);
	}

	//オイラーツアー構築
	void build(ll v=0){
		bottomindex.clear();
		downindex[v]=bottomindex.size();
		bottomindex.emplace_back(v);
		dfs(v,-1,0);
		upindex[v]=bottomindex.size();
		bottomindex.emplace_back(v);
	}

	//u,vのうち、より葉に近い頂点を返す
	ll child(ll u, ll v){
		return depth[u]<depth[v] ? v : u;
	}

	//bottomindex=オイラーツアーの各頂点の行きがけ帰りがけ順を記録した配列
	ll bottom(ll e){
		return bottomindex[e];
	}

	//頂点u-v間に対して関数fを適用する
	template<typename F>
	void query(ll u, ll v, F f){
		auto querylambda = [](ll u, ll v, F f){
			if(depth[u]>depth[v]) swap(u,v);
			f(downindex[u]+1, downindex[v]+1);
		};
		ll lcav=lca.query(u, v);
		if(lcav == u || lcav == v){
			querylambda(u, v, f);
		} else {
			querylambda(lcav, u, f);
			querylambda(lcav, v, f);
		}
	}

	//頂点vの部分木全体へ値valueのupdateを行う関数gを適用する
	template<typename T, typename G>
	void update(ll v, T value, G g){
		g(downindex[v], value);
		g(upindex[v], -value);
	}
};

/*-----8<-----library-----8<-----*/

void solve() {
	ll N, Q;
	vector<ll> a,b,c,d;
	vector<ll> x,y,u,v;
	scanf("%lld%lld", &N, &Q);
	a.assign(N, 0);
	b.assign(N, 0);
	c.assign(N, 0);
	d.assign(N, 0);
	rep(i, N-1){
		scanf("%lld%lld%lld%lld", &a[i], &b[i], &c[i], &d[i]);
		a[i]--;b[i]--;
	}
	x.assign(Q, 0);
	y.assign(Q, 0);
	u.assign(Q, 0);
	v.assign(Q, 0);
	rep(i, Q) {
		scanf("%lld%lld%lld%lld", &x[i], &y[i], &u[i], &v[i]);
		u[i]--;v[i]--;
	}

	//colorp[x]=[i,...] : 色xになっている辺の一覧
	//colorq[x]=[i,...] : 色xが必要なクエリ番号の一覧
	vector<vector<ll>> colorp(N),colorq(N);
	rep(i,N)colorp[c[i]].emplace_back(i);
	rep(i,Q)colorq[x[i]].emplace_back(i);

	//グラフ作る
	vector<vector<ll>> g(N);
	rep(i,N-1){
		ll from=a[i];
		ll to=b[i];
		g[from].push_back(to);
		g[to].push_back(from);
	}

	DoublingLowestCommonAncestor< vector<vector<ll>> > lca(g);
	lca.build(0);

	EulerTourForEdge euler(g);
	euler.build();

	//distanceseg=根からの距離を管理するセグ木
	//upとdownの分でsizeは2倍必要
	auto F=[](ll x,ll y){return x+y;};
	LazySegmentTree<ll,ll,decltype(F),decltype(F),decltype(F)> distanceseg(2*N, F, F, F, 0LL, 0LL);

	auto updatedistanceseg=[&](ll v,ll value){
		distanceseg.update(v, v+1, value);
	};
	
	vector<ll> ans(Q,0);

	//辺の両端の頂点のうち、降りた先の頂点に降りた分の距離d[i]を加算
	//euler.query(u, v, [&](ll l, ll r){ ll dis=distanceseg.query(l,r); });
	//　で頂点u-v間の距離disが得られる
	rep(i,N-1){
		ll ch=euler.child(a[i], b[i]);
		euler.update(ch, d[i], updatedistanceseg);
	}

	//ans[i] = u-lca-v 間の距離
	rep(i,Q){
		auto f=[&](ll l, ll r){
			ans[i]+=distanceseg.query(l, r);
		};
		euler.query(u[i], v[i], f);
	}


	//countseg:色xの辺の数、sumseg:色xの距離の和を管理するセグ木
	//auto F=[](ll x,ll y){return x+y;};
	LazySegmentTree<ll,ll,decltype(F),decltype(F),decltype(F)> countseg(2*N, F, F, F, 0LL, 0LL);
	LazySegmentTree<ll,ll,decltype(F),decltype(F),decltype(F)> sumseg(2*N, F, F, F, 0LL, 0LL);

	auto updatecountseg=[&](ll v,ll value){countseg.update(v, v+1, value);};
	auto updatesumseg=[&](ll v,ll value){sumseg.update(v, v+1, value);};

	//色ごとにループ
	rep(t,N){
		//i:クエリ番号
		//色tによる距離変更を適用
		for(ll i : colorp[t]){
			ll ch=euler.child(a[i], b[i]);
			euler.update(ch, +1, updatecountseg);
			euler.update(ch, +d[i], updatesumseg);
		}

		//ans[i]に距離増減を反映
		for(ll i : colorq[t]){
			ll res=0,num=0;
			auto f=[&](ll l, ll r){
				assert(l<=r);
				num+=countseg.query(l, r);
				res+=sumseg.query(l, r);
			};
			euler.query(u[i], v[i], f);
			ans[i]+=num*y[i];
			ans[i]-=res;
		}

		//色tによる距離変更を元に戻す
		for(ll i : colorp[t]){
			ll ch=euler.child(a[i], b[i]);
			euler.update(ch, -1, updatecountseg);
			euler.update(ch, -d[i], updatesumseg);
		}
	}

	rep(i,Q)printf("%lld\n", ans[i]);
}

// https://atcoder.jp/contests/abc133/tasks/abc133_f
signed main() {
	solve();
	return 0;
}


Compiler output:
a.cc:4:18: fatal error: boost/multiprecision/cpp_int.hpp: No such file or directory
    4 |         #include <boost/multiprecision/cpp_int.hpp>
      |                  ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
compilation terminated.