#include using namespace std; template struct BIT { int n; vector dat; BIT(int n=0){ initialize(n); } void initialize(int nin){ n = nin; dat.resize(n, 0); } T sum(int i){ T s = 0; while(i >= 0){ s += dat[i]; i = (i & (i+1)) - 1; } return s; } T sum_between(int i, int j){ return sum(j) - sum(i-1); } void plus(int i, T x){ while(i < n){ dat[i] += x; i |= i+1; } } // a[0]+...+a[ret] >= x int lower_bound(T x){ if(x < 0) return -1; int ret = -1; int k = 1; while(2*k <= n) k <<= 1; for( ;k>0; k>>=1){ if(ret+k < n && dat[ret+k] < x){ x -= dat[ret+k]; ret += k; } } return ret + 1; } }; template struct Segtree { int n, n_org; T e; vector dat; typedef function Func; Func f; Segtree(){} Segtree(int n_input, Func f_input, T e_input){ initialize(n_input, f_input, e_input); } void initialize(int n_input, Func f_input, T e_input){ n_org = n_input; f = f_input; e = e_input; n = 1; while(n < n_input) n <<= 1; dat.resize(2*n-1, e); } void build(vector& A){ for(int k=0; k=0; k--) dat[k] = f(dat[2*k+1], dat[2*k+2]); } void update(int k, T a){ k += n - 1; dat[k] = a; while(k > 0){ k = (k - 1)/2; dat[k] = f(dat[2*k+1], dat[2*k+2]); } } T get(int k){ return dat[k+n-1]; } T between(int a, int b){ return query(a, b+1, 0, 0, n); } T query(int a, int b, int k, int l, int r){ if(r<=a || b<=l) return e; if(a<=l && r<=b) return dat[k]; T vl = query(a, b, 2*k+1, l, (l+r)/2); T vr = query(a, b, 2*k+2, (l+r)/2, r); return f(vl, vr); } // [S, t] が条件checkを満たす最大のtを求める int bisect(int S, function check){ T val = get(S); int k = S+n-1, l = S, r = S+1; while(true){ while(k%2) k = (k-1)/2, r += r-l; T val2 = f(val, dat[k]); if(check(val2)){ if(r == n) return n_org-1; val = val2; k++; int d = r-l; l += d, r += d; }else{ break; } } while(k> N; vector P(N); for(int i=0; i> P[i], P[i]--; int64_t ans = 0; for(int t=0; t<2; t++){ priority_queue> que; vector> out(N); for(int i=0; i P[i]){ out[i].push_back(que.top().second); que.pop(); } que.emplace(P[i], i); } BIT bit(N); Segtree st(N, [](int a, int b){ return max(a, b); }, -1); for(int i=0; i