#line 1 "A.cpp" // #pragma GCC target("avx2") // #pragma GCC optimize("O3") // #pragma GCC optimize("unroll-loops") #include using namespace std; using ll = long long; using ull = unsigned long long; template using pq = priority_queue; template using qp = priority_queue, greater>; #define vec(T, A, ...) vector A(__VA_ARGS__); #define vvec(T, A, h, ...) vector> A(h, vector(__VA_ARGS__)); #define vvvec(T, A, h1, h2, ...) vector>> A(h1, vector>(h2, vector(__VA_ARGS__))); #ifndef RIN__LOCAL #define endl "\n" #endif #define spa ' ' #define len(A) A.size() #define all(A) begin(A), end(A) #define fori1(a) for(ll _ = 0; _ < (a); _++) #define fori2(i, a) for(ll i = 0; i < (a); i++) #define fori3(i, a, b) for(ll i = (a); i < (b); i++) #define fori4(i, a, b, c) for(ll i = (a); ((c) > 0 || i > (b)) && ((c) < 0 || i < (b)); i += (c)) #define overload4(a, b, c, d, e, ...) e #define fori(...) overload4(__VA_ARGS__, fori4, fori3, fori2, fori1)(__VA_ARGS__) template vector> ENUMERATE(vector &A, ll s = 0){ vector> ret(A.size()); for(int i = 0; i < A.size(); i++) ret[i] = {i + s, A[i]}; return ret; } vector> ENUMERATE(string &A, ll s = 0){ vector> ret(A.size()); for(int i = 0; i < A.size(); i++) ret[i] = {i + s, A[i]}; return ret; } #define enum1(A) fori(A.size()) #define enum2(a, A) for(auto a:A) #define enum3(i, a, A) for(auto&& [i, a]: ENUMERATE(A)) #define enum4(i, a, A, s) for(auto&& [i, a]: ENUMERATE(A, s)) #define enum(...) overload4(__VA_ARGS__, enum4, enum3, enum2, enum1)(__VA_ARGS__) template vector> ZIP(vector &A, vector &B){ int n = min(A.size(), B.size()); vector> ret(n); for(int i = 0; i < n; i++) ret[i] = {A[i], B[i]}; return ret; } template vector> ENUMZIP(vector &A, vector &B, ll s = 0){ int n = min(A.size(), B.size()); vector> ret(n); for(int i = 0; i < n; i++) ret[i] = {i + s, A[i], B[i]}; return ret; } #define zip4(a, b, A, B) for(auto&& [a, b]: ZIP(A, B)) #define enumzip5(i, a, b, A, B) for(auto&& [i, a, b]: ENUMZIP(A, B)) #define enumzip6(i, a, b, A, B, s) for(auto&& [i, a, b]: ENUMZIP(A, B, s)) #define overload6(a, b, c, d, e, f, g, ...) g #define zip(...) overload6(__VA_ARGS__, enumzip6, enumzip5, zip4, _, _, _)(__VA_ARGS__) vector stoc(string &S){ int n = S.size(); vector ret(n); for(int i = 0; i < n; i++) ret[i] = S[i]; return ret; } #define INT(...) int __VA_ARGS__; inp(__VA_ARGS__); #define LL(...) ll __VA_ARGS__; inp(__VA_ARGS__); #define STRING(...) string __VA_ARGS__; inp(__VA_ARGS__); #define CHAR(...) char __VA_ARGS__; inp(__VA_ARGS__); #define VEC(T, A, n) vector A(n); inp(A); #define VVEC(T, A, n, m) vector> A(n, vector(m)); inp(A); const ll MOD1 = 1000000007; const ll MOD9 = 998244353; template auto min(const T& a){ return *min_element(all(a)); } template auto max(const T& a){ return *max_element(all(a)); } template auto clamp(T &a, const S &l, const S &r) { return (a > r ? r : a < l ? l : a); } template inline bool chmax(T &a, const S &b) { return (a < b ? a = b, 1 : 0); } template inline bool chmin(T &a, const S &b) { return (a > b ? a = b, 1 : 0); } template inline bool chclamp(T &a, const S &l, const S &r) { auto b = clamp(a, l, r); return (a != b ? a = b, 1 : 0); } void FLUSH(){cout << flush;} void print(){cout << endl;} template void print(Head &&head, Tail &&... tail) { cout << head; if (sizeof...(Tail)) cout << spa; print(forward(tail)...); } template void print(vector &A){ int n = A.size(); for(int i = 0; i < n; i++){ cout << A[i]; if(i != n - 1) cout << ' '; } cout << endl; } template void print(vector> &A){ for(auto &row: A) print(row); } template void print(pair &A){ cout << A.first << spa << A.second << endl; } template void print(vector> &A){ for(auto &row: A) print(row); } template void prisep(vector &A, S sep){ int n = A.size(); for(int i = 0; i < n; i++){ cout << A[i]; if(i == n - 1) cout << endl; else cout << sep; } } template void priend(T A, S end){ cout << A << end; } template void priend(T A){ priend(A, spa); } template void inp(T&... a){ (cin >> ... >> a); } template void inp(vector &A){ for(auto &a:A) cin >> a; } template void inp(vector> &A){ for(auto &row:A) inp(row); } template void inp(pair &A){ inp(A.first, A.second); } template void inp(vector> &A){ for(auto &row: A) inp(row.first, row.second); } template T sum(vector &A){ T tot = 0; for(auto a:A) tot += a; return tot; } template vector compression(vector X){ sort(all(X)); X.erase(unique(all(X)), X.end()); return X; } vector> read_edges(int n, int m, bool direct=false, int indexed=1){ vector> edges(n, vector()); for(int i = 0; i < m; i++){ INT(u, v); u -= indexed; v -= indexed; edges[u].push_back(v); if(!direct) edges[v].push_back(u); } return edges; } vector> read_tree(int n, int indexed=1){ return read_edges(n, n - 1, false, indexed); } template vector>> read_wedges(int n, int m, bool direct=false, int indexed=1){ vector>> edges(n, vector>()); for(int i = 0; i < m; i++){ INT(u, v); T w; inp(w); u -= indexed; v -= indexed; edges[u].push_back({v, w}); if(!direct) edges[v].push_back({u, w}); } return edges; } template vector>> read_wtree(int n, int indexed=1){ return read_wedges(n, n - 1, false, indexed); } struct UnionFind{ int n; vector par; int group; UnionFind() = default; UnionFind(int n) : n(n){ par.assign(n, -1); group = n; } int find(int x){ if(par[x] < 0) return x; par[x] = find(par[x]); return par[x]; } bool unite(int x, int y){ x = find(x); y = find(y); if(x == y) return false; if(par[x] > par[y]) swap(x, y); group--; par[x] += par[y]; par[y] = x; return true; } bool same(int x, int y){ return find(x) == find(y); } int size(int x){ return -par[find(x)]; } vector roots(){ vector ret; for(int i = 0; i < n; i++){ if(i == find(i)) ret.push_back(i); } return ret; } }; #line 2 "Library/C++/other/Modint.hpp" template struct Modint{ int x; Modint() : x(0){} Modint(int64_t y){ if(y >= 0) x = y % MOD; else x = (y % MOD + MOD) % MOD; } Modint &operator+=(const Modint &p){ x += p.x; if(x >= MOD) x -= MOD; return *this; } Modint &operator-=(const Modint &p){ x -= p.x; if(x < 0) x += MOD; return *this; } Modint &operator*=(const Modint &p){ x = int(1LL * x * p.x % MOD); return *this; } Modint &operator/=(const Modint &p){ *this *= p.inverse(); return *this; } Modint &operator%=(const Modint &p){ assert(p.x == 0); return *this; } Modint operator-() const{ return Modint(-x); } Modint& operator++() { x++; if (x == MOD) x = 0; return *this; } Modint& operator--() { if (x == 0) x = MOD; x--; return *this; } Modint operator++(int) { Modint result = *this; ++*this; return result; } Modint operator--(int) { Modint result = *this; --*this; return result; } friend Modint operator+(const Modint &lhs, const Modint &rhs){ return Modint(lhs) += rhs; } friend Modint operator-(const Modint &lhs, const Modint &rhs){ return Modint(lhs) -= rhs; } friend Modint operator*(const Modint &lhs, const Modint &rhs){ return Modint(lhs) *= rhs; } friend Modint operator/(const Modint &lhs, const Modint &rhs){ return Modint(lhs) /= rhs; } friend Modint operator%(const Modint &lhs, const Modint &rhs){ assert(rhs.x == 0); return Modint(lhs); } bool operator==(const Modint &p) const{ return x == p.x; } bool operator!=(const Modint &p) const{ return x != p.x; } bool operator<(const Modint &rhs) const{ return x < rhs.x; } bool operator<=(const Modint &rhs) const{ return x <= rhs.x; } bool operator>(const Modint &rhs) const{ return x > rhs.x; } bool operator>=(const Modint &rhs) const{ return x >= rhs.x; } Modint inverse() const{ int a = x, b = MOD, u = 1, v = 0, t; while(b > 0){ t = a / b; a -= t * b; u -= t * v; swap(a, b); swap(u, v); } return Modint(u); } Modint pow(int64_t k) const{ Modint ret(1); Modint y(x); while(k > 0){ if(k & 1) ret *= y; y *= y; k >>= 1; } return ret; } friend ostream &operator<<(ostream &os, const Modint &p){ return os << p.x; } friend istream &operator>>(istream &is, Modint &p){ int64_t y; is >> y; p = Modint(y); return (is); } static int get_mod(){ return MOD; } }; struct Arbitrary_Modint{ int x; static int MOD; static void set_mod(int mod){ MOD = mod; } Arbitrary_Modint() : x(0){} Arbitrary_Modint(int64_t y){ if(y >= 0) x = y % MOD; else x = (y % MOD + MOD) % MOD; } Arbitrary_Modint &operator+=(const Arbitrary_Modint &p){ x += p.x; if(x >= MOD) x -= MOD; return *this; } Arbitrary_Modint &operator-=(const Arbitrary_Modint &p){ x -= p.x; if(x < 0) x += MOD; return *this; } Arbitrary_Modint &operator*=(const Arbitrary_Modint &p){ x = int(1LL * x * p.x % MOD); return *this; } Arbitrary_Modint &operator/=(const Arbitrary_Modint &p){ *this *= p.inverse(); return *this; } Arbitrary_Modint &operator%=(const Arbitrary_Modint &p){ assert(p.x == 0); return *this; } Arbitrary_Modint operator-() const{ return Arbitrary_Modint(-x); } Arbitrary_Modint& operator++() { x++; if (x == MOD) x = 0; return *this; } Arbitrary_Modint& operator--() { if (x == 0) x = MOD; x--; return *this; } Arbitrary_Modint operator++(int) { Arbitrary_Modint result = *this; ++*this; return result; } Arbitrary_Modint operator--(int) { Arbitrary_Modint result = *this; --*this; return result; } friend Arbitrary_Modint operator+(const Arbitrary_Modint &lhs, const Arbitrary_Modint &rhs){ return Arbitrary_Modint(lhs) += rhs; } friend Arbitrary_Modint operator-(const Arbitrary_Modint &lhs, const Arbitrary_Modint &rhs){ return Arbitrary_Modint(lhs) -= rhs; } friend Arbitrary_Modint operator*(const Arbitrary_Modint &lhs, const Arbitrary_Modint &rhs){ return Arbitrary_Modint(lhs) *= rhs; } friend Arbitrary_Modint operator/(const Arbitrary_Modint &lhs, const Arbitrary_Modint &rhs){ return Arbitrary_Modint(lhs) /= rhs; } friend Arbitrary_Modint operator%(const Arbitrary_Modint &lhs, const Arbitrary_Modint &rhs){ assert(rhs.x == 0); return Arbitrary_Modint(lhs); } bool operator==(const Arbitrary_Modint &p) const{ return x == p.x; } bool operator!=(const Arbitrary_Modint &p) const{ return x != p.x; } bool operator<(const Arbitrary_Modint &rhs) { return x < rhs.x; } bool operator<=(const Arbitrary_Modint &rhs) { return x <= rhs.x; } bool operator>(const Arbitrary_Modint &rhs) { return x > rhs.x; } bool operator>=(const Arbitrary_Modint &rhs) { return x >= rhs.x; } Arbitrary_Modint inverse() const{ int a = x, b = MOD, u = 1, v = 0, t; while(b > 0){ t = a / b; a -= t * b; u -= t * v; swap(a, b); swap(u, v); } return Arbitrary_Modint(u); } Arbitrary_Modint pow(int64_t k) const{ Arbitrary_Modint ret(1); Arbitrary_Modint y(x); while(k > 0){ if(k & 1) ret *= y; y *= y; k >>= 1; } return ret; } friend ostream &operator<<(ostream &os, const Arbitrary_Modint &p){ return os << p.x; } friend istream &operator>>(istream &is, Arbitrary_Modint &p){ int64_t y; is >> y; p = Arbitrary_Modint(y); return (is); } static int get_mod(){ return MOD; } }; int Arbitrary_Modint::MOD = 998244353; using modint9 = Modint<998244353>; using modint1 = Modint<1000000007>; using modint = Arbitrary_Modint; #line 368 "A.cpp" using mint = modint9; template struct WeightedUnionFind{ int n; vector par; vector D; vector W; vector sum; int group; WeightedUnionFind(int n, vector W) : n(n), W(W){ par.assign(n, -1); group = n; D.assign(n, T(0)); sum.assign(n, T(0)); } WeightedUnionFind(int n): WeightedUnionFind(n, vector(n, T(0))) {} int find(int x){ if(par[x] < 0) return x; int p = find(par[x]); D[x] += D[par[x]]; par[x] = find(p); return p; } bool unite(int x, int y, T d){ // x = y + d T dd = d; int xp = find(x); int yp = find(y); d -= D[x]; x = xp; d += D[y]; y = yp; if(x == y){ assert(d == 0); return false; } if(par[x] > par[y]){ swap(x, y); d *= -1; } sum[x] += dd * mint(par[x]) * mint(par[y]); group--; par[x] += par[y]; D[y] = -d; par[y] = x; return true; } bool same(int x, int y){ return find(x) == find(y); } int size(int x){ return -par[find(x)]; } vector roots(){ vector ret; for(int i = 0; i < n; i++){ if(i == find(i)) ret.push_back(i); } return ret; } T diff(int x, int y){ assert(same(x, y)); return D[x] - D[y]; } void all_add(int a, T b){ a = find(a); W[a] += b; } T get(int a){ int p = find(a); return W[p] + diff(a, p); } }; ////////// https://tiramister.net/blog/posts/persistent-unionfind/ const int INF = 1 << 25; // 頂点数 const int V_NUM = 200100; class persistentUF { public: int now; // 現在時刻 int par[V_NUM]; // 各頂点の親 int rank[V_NUM]; // その頂点を根とする木の深さ int time[V_NUM]; // 親がいつ更新されたか mint sum[V_NUM]; explicit persistentUF() { now = 0; for (int i = 0; i < V_NUM; ++i) { par[i] = i; } fill(rank, rank + V_NUM, 1); fill(sum, sum + V_NUM, 0); fill(time, time + V_NUM, INF); // 自身が根の間は便宜上INFとする } // 時刻tにおけるxの親を返す int find(int x, int t) { if (t < time[x]) { return x; } else { return find(par[x], t); } } // 時刻tにてxとyが同じ木に属するか判定 bool same(int x, int y, int t) { return find(x, t) == find(y, t); } // 頂点xとyを繋げる // 繋げた直後の時刻を返す int unite(int x, int y, mint d) { ++now; // 時間を進める x = find(x, now); y = find(y, now); if (x == y) return now; // rank[x] >= rank[y]にする if (rank[x] < rank[y]) swap(x, y); // rankの大きい方、つまりxにyをくっつける par[y] = x; time[y] = now; // timeに時刻を記録 sum[x] += d * rank[x] * rank[y]; rank[x] += rank[y]; return now; } }; void solve(){ INT(n, x, Q); persistentUF UF; vec(int, times, 1, 0); vec(int, D, 1, 0); fori(Q){ INT(t); // print(t, x); if(t == 1){ INT(v, w); if(UF.same(x, v, times.back())) continue; int res = UF.unite(x, v, w); times.push_back(res); D.push_back(w); } else if(t == 2){ INT(u, v); if(!UF.same(u, v, times.back())) print(-1); else if(u == v){ print(0); } else{ int l = 0; int r = times.size(); while(r - l > 1){ int mid = (l + r) / 2; if(UF.same(u, v, mid)) r = mid; else l = mid; } ll d = D[r]; print(d); x = (x + d) % n; } } else if(t == 3){ INT(v); print(UF.sum[UF.find(v, times.back())]); } else{ INT(v); x = (x + v) % n; } } } int main(){ cin.tie(0)->sync_with_stdio(0); // cout << fixed << setprecision(12); int t; t = 1; // cin >> t; while(t--) solve(); return 0; }