#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define endl "\n" #define ll long long #define ch char #define vec vector #define vll vector #define sll set #define pll pair #define mkp make_pair #define mll map #define puf push_front #define pub push_back #define pof pop_front() #define pob pop_back() #define em empty() #define fi first #define se second #define fr front() #define ba back() #define be begin() #define rbe rbegin() #define en end() #define ren rend() #define all(x) x.begin(),x.end() #define rall(x) x.rbegin(),x.rend() #define fo(i,x,y) for(ll i=x;i<=y;++i) #define fa(i,v) for(auto &i:v) #define re return #define sz size() #define so(v) sort(all(v)) #define pop_count(x) __builtin_popcount(x) #define rso(v) sort(rall(v)) #define rev(v) reverse(all(v)) #define i(x) for(ll i=0;i void ff(T a,A... b){ cout< void cc(T a,A... b){ cout< void cl(T a,A... b){ cout< void cn(T a,A... b){ cout< void ci(A&... a){ (cin>>...>>a); }; templatevoid ou(T v){fa(i,v)cout<void oun(T v){fa(i,v)cout<void ouu(T v){fa(i,v){fa(j,i)cout< void oul(T v){fa(i,v)cout<void in(T &v){fa(i,v)cin>>i;} templatevoid inn(T &v){fa(i,v)fa(j,i)cin>>j;} templatevoid oump(T &v){fa(i,v)ff(i.fi,i.se);} templatevoid pi(pair &p){ci(p.fi,p.se);} templatevoid po(pair &p){ff(p.fi,p.se);} template void fl(T a,A... b){ cout< dist1(1.0, 100000); i(10000){ // 各分布法に基いて乱数を生成 ll n = dist1(engine); } return 0; } */ mll to_prime(ll x){ mll mp; for(ll i=2;i*i<=x;++i){ while(x%i==0){ ++mp[i]; x/=i; } } if(x!=1) ++mp[x]; re mp; } #define acc(v) accumulate(v.begin(),v.end(),0LL) #define acci(v,i) accumulate(v.begin(),v.begin()+i,0LL) #define dll deque int main(void){ init(); solve(); return 0; } template class pnt{ public: T x,y; pnt(T x=0,T y=0):x(x),y(y){} pnt operator + (const pnt r)const { return pnt(x+r.x,y+r.y);} pnt operator - (const pnt r)const { return pnt(x-r.x,y-r.y);} pnt operator * (const pnt r)const { return pnt(x*r.x,y*r.y);} pnt operator / (const pnt r)const { return pnt(x/r.x,y/r.y);} pnt &operator += (const pnt r){ x+=r.x;y+=r.y;return *this;} pnt &operator -= (const pnt r){ x-=r.x;y-=r.y;return *this;} pnt &operator *= (const pnt r){ x*=r.x;y*=r.y;return *this;} pnt &operator /= (const pnt r){ x/=r.x;y/=r.y;return *this;} ll dist(const pnt r){ re (x-r.x)*(x-r.x)+(y-r.y)*(y-r.y); } ll man(const pnt r){ re abs(x-r.x)+abs(y-r.y); } pnt rot(const dou theta){ T xx,yy; xx=cos(theta)*x-sin(theta)*y; yy=sin(theta)*x+cos(theta)*y; return pnt(xx,yy); } }; istream &operator >> (istream &is,pnt &r){is>>r.x>>r.y;return is;} ostream &operator << (ostream &os,pnt &r){os<bool chmaxeq(T& a, const T& b) { if (a <= b) { a = b; return 1; } return 0; } templatebool chmineq(T& a, const T& b) { if (b <= a) { a = b; return 1; } return 0; } templatebool chmax(T& a, const T& b) { if (a < b) { a = b; return 1; } return 0; } templatebool chmin(T& a, const T& b) { if (b < a) { a = b; return 1; } return 0; } struct Trie{ struct Node{ vll nxt; vec done; ll dep,cnt=0; Node(ll c_):nxt(30),dep(c_){} }; ll root=0; vectree={Node(root)}; void ins(st s){ ll c=0; for(ll i=0;i1)++ans; else break; c=to; } re ans; } }; #define fo(i,x,y) for(ll i=x;i<=y;++i) #define rfo(_ii,_xx,_yy) for(ll _ii=_xx;_ii>=_yy;--_ii) #define qll queue template using pq= priority_queue; template using pqg= priority_queue,greater>; vec>rle(st s){//run_length_encoding ll n=s.sz; vec>ans; for(ll i=0;i> mat_mul(vector> a, vector> b, ll mod) { // 行列乗算 int n = a.size(); vector> res(n, vector(n)); for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { for (int k = 0; k < n; k++) { res[i][j] += a[i][k] * b[k][j]; res[i][j] %= mod; } } } return res; } vector> mat_pow(vector> a, ll b, ll mod) { // 行列累乗 int n = a.size(); vector> res(n, vector(n)); for (int i = 0; i < n; i++) res[i][i] = 1; while (b) { if (b & 1) res = mat_mul(res, a, mod); a = mat_mul(a, a, mod); b >>= 1; } return res; } void Yes(bool f){ ff(f?"Yes":"No");re; } void yes(bool f){ ff(f?"yes":"no");re; } void YES(bool f){ ff(f?"YES":"NO");re; } void sub(); void solve(){ //ge(ll,t); ll t=1; xx(t){ sub(); } } template struct Edge { int rv, from, to; // rev:逆向きの辺の番号 T cap, original_cap; Edge(int f, int t, T c,int r ) : rv(r), from(f), to(t), cap(c), original_cap(c) {} }; template struct Graph { vector>> G; Graph(int n = 0) : G(n) {} vector>& operator[](int i) { return G[i]; }//G[i]でi番目の辺を返す const size_t size() const { return G.size(); }// Edge& redge(Edge e) { // 逆向きの辺を返す return G[e.to][e.rv]; // 自己ループはないと仮定(あれば G[e.to][e.rev + 1] とする必要がある) } void add_edg(int from, int to, T cap) { // 有向辺を加える G[from].push_back( Edge( from, to, cap,(int)G[to].size() )); G[to].push_back(Edge( to, from, 0 ,(int)G[from].size() - 1)); } }; /* FordFulkerson: Ford-Fulkersonのアルゴリズムで最大流を求める構造体 max_flow(G,s,t):sからtへのグラフGの最大流を求める 副作用:G は最大流の残余ネットワークになる 計算量: O(|f*||E|) (f*:最大流) (この最悪ケースになることはほぼ無い) */ template struct Ford { const T INF = 1e9; vector used; Ford(){}; T dfs(Graph& G, int v, int t, T f) { // 増加可能経路を見つけて増加分のフローを返す if (v == t) return f; used[v] = true; for (auto& e : G[v]) { if (!used[e.to] && e.cap > 0) { T d = dfs(G, e.to, t, min(f, e.cap)); if (d > 0) { e.cap -= d; G.redge(e).cap += d; return d; } } } return 0; } T max_flow(Graph& G, int s, int t) { T flow = 0; while (true) { used.assign(G.size(), 0); T f = dfs(G, s, t, INF); // f が増加分のフロー if (f == 0) { return flow; } else { flow += f; } } return 0; } /* Fordf; ff(f.max_flow(g,s,sink)); */ }; struct UF{ vll par,rk,siz; UF(ll n):par(n+5,-1),rk(n+5,0){ } ll root(ll x){ if(par[x]<0)return x; else return par[x]=root(par[x]); } bool same(ll x,ll y){ return root(x)==root(y); } bool unite(ll x,ll y){ ll rx=root(x),ry=root(y); if(rx==ry) return 0; if(rk[rx]>e(n); vll dis(n+50,MAX); pqgq;q.push({0ll,start}); dis[start]=0; while(q.em^1){ auto [d,now]=q.top();q.pop(); if(dis[now]>cost){//O((n+m)log(n)) vll dis(n+5,MAX); pqgq;q.push({0ll,start}); dis[start]=0; while(q.em^1){ auto [d,now]=q.top();q.pop(); if(dis[now]> g; vector vid, head, heavy, parent, depth, inv, in, out; HLDecomposition() {} HLDecomposition(int n) { init(n); } void init(int n) { g.resize(n); vid.resize(n, -1); head.resize(n); heavy.resize(n, -1); parent.resize(n); depth.resize(n); inv.resize(n); in.resize(n); out.resize(n); } void add(int u, int v) { g[u].insert(v); g[v].insert(u); } void build(int root) { dfs(root, -1); t = 0; dfs_hld(root); } int dfs(int curr, int prev) { parent[curr] = prev; int sub = 1, max_sub = 0; heavy[curr] = -1; for (int next : g[curr]) if (next != prev) { depth[next] = depth[curr] + 1; int sub_next = dfs(next, curr); sub += sub_next; if (max_sub < sub_next) max_sub = sub_next, heavy[curr] = next; }return sub; } int t = 0; void dfs_hld(int v = 0) { vid[v] = in[v] = t; t++; inv[in[v]] = v; if (0 <= heavy[v]) { head[heavy[v]] = head[v]; dfs_hld(heavy[v]); } for (auto u : g[v]) if (depth[v] < depth[u]) if (u != heavy[v]) { head[u] = u; dfs_hld(u); } out[v] = t; } void foreach(int u, int v, function f) { // [x,y] if (vid[u] > vid[v]) swap(u, v); f(max(vid[head[v]], vid[u]), vid[v]); if (head[u] != head[v]) foreach(u, parent[head[v]], f); } int ancestor(int from, int times) { while (true) { if (depth[head[from]] > depth[from] - times) { times -= depth[from] - depth[head[from]] + 1; if (head[from] == 0)return -1; from = parent[head[from]]; } else return inv[vid[from] - times]; } } int lca(int u, int v) {//lowest common ancestor 最近共通祖先 if (vid[u] > vid[v]) swap(u, v); if (head[u] == head[v]) return u; return lca(u, parent[head[v]]); } int distance(int u, int v) { return depth[u] + depth[v] - 2 * depth[lca(u, v)]; } int child(int parent, int child, int times) { assert(depth[parent] < depth[child]); int d = distance(parent, child); assert(times - 1 <= d); return ancestor(child, d - times); } int go(int from, int to, int times) { int d = distance(from, to); assert(0 <= times and times <= d); int lc = lca(from, to); if (lc == to)return ancestor(from, times); if (lc == from)return child(from, to, times); int dd = distance(from, lc); if (dd <= times)return go(lc, to, times - dd); return ancestor(from, times); } }; vll z_algo(string pattern, string text) {//z algorithm 最長共通接頭辞 string str = pattern + "$" + text; ll n = str.length(),np=pattern.sz; vll Z(n+5),ret(n+5); for (ll i = 1,L=0,R=0,j; i < n; ++i) { if (Rnp) ret[i-np-1]=Z[i]; R--; }else{ j= i-L; if (Z[j]np) ret[i-np-1]=Z[i]; } } //ff("z"); ou(Z); re ret; } // 2^10 = 1024 //vll dy={-1,-1,-1,0,0,1,1,1},dx={-1,0,1,-1,1,-1,0,1}; /* O(2*10^8) 9*10^18 1LL<<62 4*10^18 ~~(v.be,v.be+n,x); not include v.be+n set.lower_bound(x); ->. *++~ ! /%* +- << < == & && +=?: */ // 12345678901234567890 //vll dy={-1,0,0,1},dx={0,-1,1,0}; ll mod = 1000000007; const ll INF = mod * mod; typedef pairP; #define per(i,n) for(int i=n-1;i>=0;i--) #define Rep(i,sta,n) for(int i=sta;i=1;i--) #define Rep1(i,sta,n) for(int i=sta;i<=n;i++) ll mod_pow(ll x, ll n, ll m = mod) { if (n < 0) { ll res = mod_pow(x, -n, m); return mod_pow(res, m - 2, m); } if (abs(x) >= m)x %= m; if (x < 0)x += m; //if (x == 0)return 0; ll res = 1; while (n) { if (n & 1)res = res * x % m; x = x * x % m; n >>= 1; } return res; } //mod should be <2^31 struct modint { int n; modint() :n(0) { ; } modint(ll m) { if (m < 0 || mod <= m) { m %= mod; if (m < 0)m += mod; } n = m; } operator int() { return n; } }; bool operator==(modint a, modint b) { return a.n == b.n; } bool operator<(modint a, modint b) { return a.n < b.n; } modint operator+=(modint& a, modint b) { a.n += b.n; if (a.n >= mod)a.n -= (int)mod; return a; } modint operator-=(modint& a, modint b) { a.n -= b.n; if (a.n < 0)a.n += (int)mod; return a; } modint operator*=(modint& a, modint b) { a.n = ((ll)a.n * b.n) % mod; return a; } modint operator+(modint a, modint b) { return a += b; } modint operator-(modint a, modint b) { return a -= b; } modint operator*(modint a, modint b) { return a *= b; } modint operator^(modint a, ll n) { if (n == 0)return modint(1); modint res = (a * a) ^ (n / 2); if (n % 2)res = res * a; return res; } ll inv(ll a, ll p) { return (a == 1 ? 1 : (1 - p * inv(p % a, a)) / a + p); } modint operator/(modint a, modint b) { return a * modint(inv(b, mod)); } modint operator/=(modint& a, modint b) { a = a / b; return a; } const int max_n = 1 << 20; modint fact[max_n], factinv[max_n]; void init_f() { fact[0] = modint(1); for (int i = 0; i < max_n - 1; i++) { fact[i + 1] = fact[i] * modint(i + 1); } factinv[max_n - 1] = modint(1) / fact[max_n - 1]; for (int i = max_n - 2; i >= 0; i--) { factinv[i] = factinv[i + 1] * modint(i + 1); } } modint comb(int a, int b) { if (a < 0 || b < 0 || a < b)return 0; return fact[a] * factinv[b] * factinv[a - b]; } modint combP(int a, int b) { if (a < 0 || b < 0 || a < b)return 0; return fact[a] * factinv[a - b]; } ll gcd(ll a, ll b) { a = abs(a); b = abs(b); if (a < b)swap(a, b); while (b) { ll r = a % b; a = b; b = r; } return a; } //(1)区間に一様加算 (2)区間の合計値 template class segtree{ public: vector val,to; segtree(){val.resize(NV*2); to.resize(NV*2);}; ll comp(int l,int r){ return l+r;}; ll query(int x,int y,int l=0,int r=NV,int k=1) { if(r<=x || y<=l) return 0; if(x<=l && r<=y) return to[k]; x=max(x,l); y=min(y,r); if(val[k]>=0) return val[k]?(y-x):0; return comp(query(x,y,l,(l+r)/2,k*2),query(x,y,(l+r)/2,r,k*2+1)); } void update(int x,int y,int v,int l=0,int r=NV,int k=1) { if(l>=r) return; if(x<=l && r<=y) { val[k]=v; to[k]=val[k]?(r-l):0; } else if(l < y && x < r) { if(val[k]!=-1) { val[k*2]=val[k*2+1]=val[k]; to[k*2]=to[k*2+1]=val[k]?(r-l)/2:0; val[k]=-1; } update(x,y,v,l,(l+r)/2,k*2); update(x,y,v,(l+r)/2,r,k*2+1); to[k]=comp(to[k*2],to[k*2+1]); } } }; void sub() { ge(ll,n,q); vec> s(2); vec v={0,0}; xx(q) { ge(ll,x,l,r); if(x==0) { vll c(2); i(2) c[i]=s[i].query(l,r+1); i(2) if(c[i^1]