#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; const ll mod = 1000000007; const ll INF = (ll)1000000007 * 1000000007; typedef pair P; #define rep(i,n) for(int i=0;i=0;i--) #define Rep(i,sta,n) for(int i=sta;i=sta;i--) typedef long double ld; const ld eps = 1e-8; const ld pi = acos(-1.0); typedef pair LP; int dx[4]={1,-1,0,0}; int dy[4]={0,0,1,-1}; templatebool chmax(T &a, const T &b) {if(abool chmin(T &a, const T &b) {if(b struct SegmentTree { public: SegmentTree() : SegmentTree(0) {} SegmentTree(int n) : SegmentTree(std::vector(n, e())) {} SegmentTree(const std::vector& v) : _n(int(v.size())) { log = ceil_pow2(_n); size = 1 << log; d = std::vector(2 * size, e()); for (int i = 0; i < _n; i++) d[size + i] = v[i]; for (int i = size - 1; i >= 1; i--) { update(i); } } void set_val(int p, S x) { assert(0 <= p && p < _n); p += size; d[p] = x; for (int i = 1; i <= log; i++) update(p >> i); } S get(int p) { assert(0 <= p && p < _n); return d[p + size]; } S query(int l, int r) { assert(0 <= l && l <= r && r <= _n); S sml = e(), smr = e(); l += size; r += size; while (l < r) { if (l & 1) sml = op(sml, d[l++]); if (r & 1) smr = op(d[--r], smr); l >>= 1; r >>= 1; } return op(sml, smr); } S all_query() { return d[1]; } template int max_right(int l) {//f(op([l,r)))==trueを満たす最大のr return max_right(l, [](S x) { return f(x); }); } template int max_right(int l, F f) { assert(0 <= l && l <= _n); assert(f(e())); if (l == _n) return _n; l += size; S sm = e(); do { while (l % 2 == 0) l >>= 1; if (!f(op(sm, d[l]))) { while (l < size) { l = (2 * l); if (f(op(sm, d[l]))) { sm = op(sm, d[l]); l++; } } return l - size; } sm = op(sm, d[l]); l++; } while ((l & -l) != l); return _n; } template int min_left(int r) {//f(op([l,r)))==trueを満たす最小のl return min_left(r, [](S x) { return f(x); }); } template int min_left(int r, F f) { assert(0 <= r && r <= _n); assert(f(e())); if (r == 0) return 0; r += size; S sm = e(); do { r--; while (r > 1 && (r % 2)) r >>= 1; if (!f(op(d[r], sm))) { while (r < size) { r = (2 * r + 1); if (f(op(d[r], sm))) { sm = op(d[r], sm); r--; } } return r + 1 - size; } sm = op(d[r], sm); } while ((r & -r) != r); return 0; } private: int _n, size, log; std::vector d; void update(int k) { d[k] = op(d[2 * k], d[2 * k + 1]); } int ceil_pow2(int n) { int x = 0; while ((1U << x) < (unsigned int)(n)) x++; return x; } }; int n; vector A,inv; int F(int a,int b){return max(a,b);} int e(){return -mod;} void solve(){ cin >> n;A.resize(n);inv.resize(n); rep(i,n) { cin >> A[i];A[i]--; inv[A[i]]=i; } SegmentTree seg1(n+1); seg1.set_val(n,0); per(i,n){ seg1.set_val(A[i],seg1.query(A[i]+1,n+1)+1); } // rep(i,n){ // cout << seg.get(i) << " "; // } // cout << "" << endl; int K=seg1.all_query(); vector L1; rep(i,n){ // rep(i,n){ // cout << seg1.get(i) << " "; // } // cout << "" << endl; int p=seg1.get(A[i]); if(p+L1.size()==K){ if(!L1.empty()) if(A[L1.back()]>A[i]){ seg1.set_val(A[i],e()); continue; }; L1.push_back(i); } seg1.set_val(A[i],e()); } // rep(i,K){ // cout << L1[i] << " "; // } // cout << "" << endl; reverse(A.begin(),A.end()); SegmentTree seg2(n+1); seg2.set_val(0,0); per(i,n){ seg2.set_val(A[i]+1,seg2.query(0,A[i]+1)+1); } vector L2; rep(i,n){ // rep(i,n){ // cout << seg2.get(i) << " "; // } // cout << "" << endl; int p=seg2.get(A[i]+1); if(p+L2.size()==K){ if(!L2.empty()) if(A[L2.back()]