//#define _GLIBCXX_DEBUG //#pragma GCC target("avx2") //#pragma GCC optimize("O3") //#pragma GCC optimize("unroll-loops") #include using namespace std; #ifdef LOCAL #include #define OUT(...) debug_print::multi_print(#__VA_ARGS__, __VA_ARGS__) #else #define OUT(...) (static_cast(0)) #endif #define endl '\n' #define lfs cout<= (ll)(n); i--) namespace template_tute{ using ll = long long; using ld = long double; const ll MOD1 = 1e9+7; const ll MOD9 = 998244353; const ll INF = 4.1e18; using P = pair; template using PQ = priority_queue; template using QP = priority_queue,greater>; templatebool chmin(T1 &a,T2 b){if(a>b){a=b;return true;}else return false;} templatebool chmax(T1 &a,T2 b){if(a({a,b,c})-min({a,b,c});} void ans1(bool x){if(x) cout<<"Yes"<void ans(bool x,T1 y,T2 z){if(x)cout<void anss(T1 x,T2 y,T3 z){ans(x!=y,x,z);}; templatevoid debug(const T &v,ll h,ll w,string sv=" "){for(ll i=0;ivoid debug(const T &v,ll n,string sv=" "){if(n!=0)cout<void debug(const vector&v){debug(v,v.size());} templatevoid debug(const vector>&v){for(auto &vv:v)debug(vv,vv.size());} templatevoid debug(stack st){while(!st.empty()){cout<void debug(queue st){while(!st.empty()){cout<void debug(deque st){while(!st.empty()){cout<void debug(PQ st){while(!st.empty()){cout<void debug(QP st){while(!st.empty()){cout<void debug(const set&v){for(auto z:v)cout<void debug(const multiset&v){for(auto z:v)cout<void debug(const array &a){for(auto z:a)cout<void debug(const map&v){for(auto z:v)cout<<"["<vector>vec(ll x, ll y, T w){vector>v(x,vector(y,w));return v;} vectordx={1,-1,0,0,1,1,-1,-1};vectordy={0,0,1,-1,1,-1,1,-1}; templatevector make_v(size_t a,T b){return vector(a,b);} templateauto make_v(size_t a,Ts... ts){return vector(a,make_v(ts...));} templateostream &operator<<(ostream &os, const pair&p){return os << "(" << p.first << "," << p.second << ")";} templateostream &operator<<(ostream &os, const vector &v){os<<"[";for(auto &z:v)os << z << ",";os<<"]"; return os;} templatevoid rearrange(vector&ord, vector&v){ auto tmp = v; for(int i=0;ivoid rearrange(vector&ord,Head&& head, Tail&&... tail){ rearrange(ord, head); rearrange(ord, tail...); } template vector ascend(const vector&v){ vectorord(v.size());iota(ord.begin(),ord.end(),0); sort(ord.begin(),ord.end(),[&](int i,int j){return make_pair(v[i],i) vector descend(const vector&v){ vectorord(v.size());iota(ord.begin(),ord.end(),0); sort(ord.begin(),ord.end(),[&](int i,int j){return make_pair(v[i],-i)>make_pair(v[j],-j);}); return ord; } template vector inv_perm(const vector&ord){ vectorinv(ord.size()); for(int i=0;i0);return n>=0?n/div:(n-div+1)/div;} ll CEIL(ll n,ll div){assert(div>0);return n>=0?(n+div-1)/div:n/div;} ll digitsum(ll n){ll ret=0;while(n){ret+=n%10;n/=10;}return ret;} ll modulo(ll n,ll d){return (n%d+d)%d;}; templateT min(const vector&v){return *min_element(v.begin(),v.end());} templateT max(const vector&v){return *max_element(v.begin(),v.end());} templateT acc(const vector&v){return accumulate(v.begin(),v.end(),T(0));}; templateT reverse(const T &v){return T(v.rbegin(),v.rend());}; //mt19937 mt(chrono::steady_clock::now().time_since_epoch().count()); int popcount(ll x){return __builtin_popcountll(x);}; int poplow(ll x){return __builtin_ctzll(x);}; int pophigh(ll x){return 63 - __builtin_clzll(x);}; templateT poll(queue &q){auto ret=q.front();q.pop();return ret;}; templateT poll(priority_queue &q){auto ret=q.top();q.pop();return ret;}; templateT poll(QP &q){auto ret=q.top();q.pop();return ret;}; templateT poll(stack &s){auto ret=s.top();s.pop();return ret;}; ll MULT(ll x,ll y){if(LLONG_MAX/x<=y)return LLONG_MAX;return x*y;} ll POW2(ll x, ll k){ll ret=1,mul=x;while(k){if(mul==LLONG_MAX)return LLONG_MAX;if(k&1)ret=MULT(ret,mul);mul=MULT(mul,mul);k>>=1;}return ret;} ll POW(ll x, ll k){ll ret=1;for(int i=0;isputn(d, len) != len) { dest.setstate(std::ios_base::badbit); } } return dest; } namespace converter{ int dict[500]; const string lower="abcdefghijklmnopqrstuvwxyz"; const string upper="ABCDEFGHIJKLMNOPQRSTUVWXYZ"; const string digit="0123456789"; const string digit1="123456789"; void regi_str(const string &t){ for(int i=0;ito_int(const string &s,const string &t){ regi_str(t); vectorret(s.size()); for(int i=0;iinput(const string &t){ string s;cin>>s; return to_int(s,t); } vectorto_int(const string &s){ auto t=s; sort(t.begin(),t.end()); t.erase(unique(t.begin(),t.end()),t.end()); return to_int(s,t); } vectorinput(){ string s;cin>>s; return to_int(s); } vector>to_int(const vector&s,const string &t){ regi_str(t); vector>ret(s.size(),vector(s[0].size())); for(int i=0;i>input(int h,const string &t){ vectors(h);for(int i=0;i>s[i]; return to_int(s,t); } vector>to_int(const vector&s){ string t; for(int i=0;i>input(int h){ vectors(h);for(int i=0;i>s[i]; return to_int(s); } string to_str(const vector&s,const string &t){ regi_int(t); string ret; for(auto z:s)ret+=dict[z]; return ret; } vector to_str(const vector>&s,const string &t){ regi_int(t); vectorret(s.size()); for(int i=0;i struct edge { int to; T cost; int id; edge():to(-1),id(-1){}; edge(int to, T cost = 1, int id = -1):to(to), cost(cost), id(id){} operator int() const { return to; } }; template using Graph = vector>>; template Graphrevgraph(const Graph &g){ Graphret(g.size()); for(int i=0;i Graph readGraph(int n,int m,int indexed=1,bool directed=false,bool weighted=false){ Graph ret(n); for(int es = 0; es < m; es++){ int u,v; T w=1; cin>>u>>v;u-=indexed,v-=indexed; if(weighted)cin>>w; ret[u].emplace_back(v,w,es); if(!directed)ret[v].emplace_back(u,w,es); } return ret; } template Graph readParent(int n,int indexed=1,bool directed=true){ Graphret(n); for(int i=1;i>p; p-=indexed; ret[p].emplace_back(i); if(!directed)ret[i].emplace_back(p); } return ret; } } using namespace template_tute; struct SegmentTreeBeats { static const long long INF = 1LL << 60; struct Node { Node *left, *right; long long max_v, smax_v, max_c; long long min_v, smin_v, min_c; long long sum; long long len, ladd, lval; Node() : left(0), right(0), ladd(0), lval(INF) {} void init(long long x) { max_v = min_v = sum = x; smax_v = -INF; smin_v = INF; max_c = min_c = 1; } void init_empty() { max_v = smax_v = -INF; min_v = smin_v = INF; max_c = min_c = 0; } void update_max(long long x) { sum += (x - max_v) * max_c; /* */if (max_v == min_v) max_v = min_v = x; else if (max_v == smin_v) max_v = smin_v = x; else max_v = x; if (lval != INF and x < lval) lval = x; } void update_min(long long x) { sum += (x - min_v) * min_c; /* */if (max_v == min_v) max_v = min_v = x; else if (max_v == smin_v) min_v = smax_v = x; else min_v = x; if (lval != INF and lval < x) lval = x; } void addall(long long x) { max_v += x; if (smax_v != -INF) smax_v += x; min_v += x; if (smin_v != INF) smin_v += x; sum += len * x; if (lval != INF) lval += x; else ladd += x; } void updateall(long long x) { max_v = min_v = x; smax_v = -INF; smin_v = INF; max_c = min_c = len; sum = len * x; lval = x; ladd = 0; } void push() { if (lval != INF) { left -> updateall(lval); right -> updateall(lval); lval = INF; return; } if (ladd != 0) { left -> addall(ladd); right -> addall(ladd); ladd = 0; } if (max_v < left -> max_v) left -> update_max(max_v); if (left -> min_v < min_v) left -> update_min(min_v); if (max_v < right -> max_v) right -> update_max(max_v); if (right -> min_v < min_v) right -> update_min(min_v); } void update() { sum = left -> sum + right -> sum; if (left -> max_v < right -> max_v) { max_v = right -> max_v; max_c = right -> max_c; smax_v = max(left -> max_v, right -> smax_v); } else if (left -> max_v > right -> max_v) { max_v = left -> max_v; max_c = left -> max_c; smax_v = max(left -> smax_v, right -> max_v); } else { max_v = left -> max_v; max_c = left -> max_c + right -> max_c; smax_v = max(left -> smax_v, right -> smax_v); } if (left -> min_v < right -> min_v) { min_v = left -> min_v; min_c = left -> min_c; smin_v = min(left -> smin_v, right -> min_v); } else if (left -> min_v > right -> min_v) { min_v = right -> min_v; min_c = right -> min_c; smin_v = min(left -> min_v, right -> smin_v); } else { min_v = left -> min_v; min_c = left -> min_c + right -> min_c; smin_v = min(left -> smin_v, right -> smin_v); } } }; int n, n0; Node *root; void _update_min(long long x, int a, int b, Node *nd, int l, int r) { if (b <= l || r <= a || nd -> max_v <= x) return; if (a <= l && r <= b && nd -> smax_v < x) { nd -> update_max(x); return; } nd -> push(); _update_min(x, a, b, nd -> left, l, (l + r) >> 1); _update_min(x, a, b, nd -> right, (l + r) >> 1, r); nd -> update(); } void _update_max(long long x, int a, int b, Node *nd, int l, int r) { if (b <= l || r <= a || x <= nd -> min_v) return; if (a <= l && r <= b && x < nd -> smin_v) { nd -> update_min(x); return; } nd -> push(); _update_max(x, a, b, nd -> left, l, (l + r) >> 1); _update_max(x, a, b, nd -> right, (l + r) >> 1, r); nd -> update(); } void _add_val(long long x, int a, int b, Node *nd, int l, int r) { if (b <= l || r <= a) return; if (a <= l && r <= b) { nd -> addall(x); return; } nd -> push(); _add_val(x, a, b, nd -> left, l, (l + r) >> 1); _add_val(x, a, b, nd -> right, (l + r) >> 1, r); nd -> update(); } void _update_val(long long x, int a, int b, Node *nd, int l, int r) { if (b <= l || r <= a) return; if (a <= l && r <= b) { nd -> updateall(x); return; } nd -> push(); _update_val(x, a, b, nd -> left, l, (l + r) >> 1); _update_val(x, a, b, nd -> right, (l + r) >> 1, r); nd -> update(); } long long _query_max(int a, int b, Node *nd, int l, int r) { if (b <= l || r <= a) return -INF; if (a <= l && r <= b) return nd -> max_v; nd -> push(); long long lv = _query_max(a, b, nd -> left, l, (l + r) >> 1); long long rv = _query_max(a, b, nd -> right, (l + r) >> 1, r); return max(lv, rv); } long long _query_min(int a, int b, Node *nd, int l, int r) { if (b <= l || r <= a) return INF; if (a <= l && r <= b) return nd -> min_v; nd -> push(); long long lv = _query_min(a, b, nd -> left, l, (l + r) >> 1); long long rv = _query_min(a, b, nd -> right, (l + r) >> 1, r); return min(lv, rv); } long long _query_sum(int a, int b, Node *nd, int l, int r) { if (b <= l || r <= a) return 0; if (a <= l && r <= b) return nd -> sum; nd -> push(); long long lv = _query_sum(a, b, nd -> left, l, (l + r) >> 1); long long rv = _query_sum(a, b, nd -> right, (l + r) >> 1, r); return lv + rv; } SegmentTreeBeats(int n, long long val = 0) : n(n) { n0 = 1; while (n0 < n) n0 <<= 1; Node *nds = new Node[2 * n0]; root = nds; nds[0].len = n0; for (int i = 0; i < n0 - 1; ++i) { nds[i].left = &nds[2 * i + 1]; nds[i].right = &nds[2 * i + 2]; nds[2 * i + 1].len = nds[2 * i + 2].len = nds[i].len >> 1; } for (int i = 0; i < n; ++i) nds[n0 - 1 + i].init(val); for (int i = n; i < n0; ++i) nds[n0 - 1 + i].init_empty(); for (int i = n0 - 2; i >= 0; i--) nds[i].update(); } SegmentTreeBeats(vector a) : n(a.size()) { n0 = 1; while (n0 < n) n0 <<= 1; Node *nds = new Node[2 * n0]; root = nds; nds[0].len = n0; for (int i = 0; i < n0 - 1; ++i) { nds[i].left = &nds[2 * i + 1]; nds[i].right = &nds[2 * i + 2]; nds[2 * i + 1].len = nds[2 * i + 2].len = nds[i].len >> 1; } for (int i = 0; i < n; ++i) nds[n0 - 1 + i].init(a[i]); for (int i = n; i < n0; ++i) nds[n0 - 1 + i].init_empty(); for (int i = n0 - 2; i >= 0; i--) nds[i].update(); } void update_min(int a, int b, long long x) { _update_min(x, a, b, root, 0, n0); } void update_max(int a, int b, long long x) { _update_max(x, a, b, root, 0, n0); } void add_val(int a, int b, long long x) { _add_val(x, a, b, root, 0, n0); } void update_val(int a, int b, long long x) { _update_val(x, a, b, root, 0, n0); } long long query_max(int a, int b) { return _query_max(a, b, root, 0, n0); } long long query_min(int a, int b) { return _query_min(a, b, root, 0, n0); } long long query_sum(int a, int b) { return _query_sum(a, b, root, 0, n0); } void print(){ for(int i=0;i>n; vectora(n); rep(i,0,n)cin>>a[i]; vectorb(n); rep(i,0,n)cin>>b[i]; SegmentTreeBeats sa(a); SegmentTreeBeats sb(b); vectorret; rep(i,0,n){ sa.update_min(0,i+1,a[i]); sb.update_min(0,i+1,b[i]); ll ok=-1,ng=i+1; while(ng-ok>=2){ ll mid=(ok+ng)/2; if(sa.query_max(mid,mid+1)>T; while(T--){ solve(); } return 0; }