#include #include using namespace std; #define SZ(x) (int) (x).size() #define REP(i, n) for(int i = 0; i < (n); i++) #define FOR(i, a, b) for(auto i = (a); i < (b); i++) #define For(i, a, b, c) \ for(auto i = (a); i != (b); i += (c)) #define REPR(i, n) for(auto i = (n) -1; i >= 0; i--) #define ALL(s) (s).begin(), (s).end() #define so(V) sort(ALL(V)) #define rev(V) reverse(ALL(V)) #define uni(v) v.erase(unique(ALL(v)), (v).end()) #define eb emplace_back typedef unsigned long long ull; typedef long long ll; typedef vector vi; typedef vector vll; typedef vector vb; typedef vector vvi; typedef vector vvll; typedef pair PI; typedef pair PL; const double EPS = 1e-6; const int MOD = 1e9 + 7; const int INF = (1 << 30); const ll LINF = 1e18; const double math_PI = acos(-1); template vector make_v(size_t a) { return vector(a); } template auto make_v(size_t a, Ts... ts) { return vector(ts...))>( a, make_v(ts...)); } template typename enable_if::value == 0>::type fill_v( T& t, const V& v) { t = v; } template typename enable_if::value != 0>::type fill_v( T& t, const V& v) { for(auto& e: t) fill_v(e, v); } template bool chmax(T& a, const T& b) { if(a < b) { a = b; return true; } return false; } template bool chmin(T& a, const T& b) { if(a > b) { a = b; return true; } return false; } template istream& operator>>(istream& is, pair& p) { cin >> p.first >> p.second; return is; } template istream& operator>>(istream& is, vector& vec) { for(T& x: vec) is >> x; return is; } template string join(vector& vec, string splitter) { stringstream ss; REP(i, SZ(vec)) { if(i != 0) ss << splitter; ss << vec[i]; } return ss.str(); } template ostream& operator<<(ostream& os, vector& vec) { os << join(vec, " "); return os; } class HLdecomposition { private: int V; vector> G; vector stsize, parent, pathtop, in, out; int root; void BuildStsize(int u, int p) { stsize[u] = 1, parent[u] = p; for(int& v: G[u]) { if(v == p) { if(v == G[u].back()) break; else swap(v, G[u].back()); } BuildStsize(v, u); stsize[u] += stsize[v]; if(stsize[v] > stsize[G[u][0]]) { swap(v, G[u][0]); } } } void BuildPath(int u, int p, int& tm) { in[u] = tm++; for(int v: G[u]) { if(v == p) continue; pathtop[v] = (v == G[u][0] ? pathtop[u] : v); BuildPath(v, u, tm); } out[u] = tm; } public: void add_edge(int u, int v) { G[u].push_back(v), G[v].push_back(u); } void build(int _root = 0) { root = _root; int tm = 0; BuildStsize(root, -1); pathtop[root] = root; BuildPath(root, -1, tm); } //蜈・・鬆らせ縺ョ繧、繝ウ繝・ャ繧ッ繧ケ縺ョ驟榊・荳翫〒縺ョid繧定ソ斐☆ inline int get(int a) { return in[a]; } int lca(int a, int b) { int pa = pathtop[a], pb = pathtop[b]; while(pathtop[a] != pathtop[b]) { if(in[pa] > in[pb]) { a = parent[pa], pa = pathtop[a]; } else { b = parent[pb], pb = pathtop[b]; } } if(in[a] > in[b]) swap(a, b); return a; } void subtree_query( int a, const function& func) { func(in[a], out[a]); } // 萓・ hl.path_query(q, r, [&](int l, int r){ seg.range(l, r, s); }) void path_query(int a, int b, const function& func) { int pa = pathtop[a], pb = pathtop[b]; while(pathtop[a] != pathtop[b]) { if(in[pa] > in[pb]) { func(in[pa], in[a] + 1); a = parent[pa], pa = pathtop[a]; } else { func(in[pb], in[b] + 1); b = parent[pb], pb = pathtop[b]; } } if(in[a] > in[b]) swap(a, b); func(in[a], in[b] + 1); } HLdecomposition(int node_size) : V(node_size), G(V), stsize(V, 0), parent(V, -1), pathtop(V, -1), in(V, -1), out(V, -1) { } }; using namespace atcoder; using mint = modint1000000007; using S = vector>; S op(S a, S b) { S c = S(2, vector(2)); REP(i, 2) { REP(j, 2) { REP(k, 2) { c[i][j] += a[i][k] * b[k][j]; } } } return c; } S e() { S c = S(2, vector(2)); c[0][0] = 1; c[1][1] = 1; return c; } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); int N; cin >> N; HLdecomposition hld(N + N - 1); REP(i, N - 1) { int A, B; cin >> A >> B; hld.add_edge(A, N + i); hld.add_edge(N + i, B); } hld.build(); segtree seg(2 * N - 1); int Q; cin >> Q; REP(q, Q) { char c; cin >> c; if(c == 'x') { int i; cin >> i; int x1, x2, x3, x4; cin >> x1 >> x2 >> x3 >> x4; seg.set(hld.get(N + i), vector> { { mint(x1), mint(x2) }, { mint(x3), mint(x4) } }); } else { int i, j; cin >> i >> j; auto c = e(); vector V; hld.path_query(i, j, [&](int l, int r) { V.eb(l, r); }); rev(V); for(auto [l, r]: V) { c = op(c, seg.prod(l, r)); } cout << c[0][0].val() << " " << c[0][1].val() << " " << c[1][0].val() << " " << c[1][1].val() << endl; } } return 0; }