結果
問題 | No.2809 Sort Query |
ユーザー |
|
提出日時 | 2024-07-12 23:20:32 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 516 ms / 2,000 ms |
コード長 | 5,611 bytes |
コンパイル時間 | 2,641 ms |
コンパイル使用メモリ | 216,756 KB |
最終ジャッジ日時 | 2025-02-22 05:30:14 |
ジャッジサーバーID (参考情報) |
judge5 / judge6 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 71 |
ソースコード
#include <bits/stdc++.h>using namespace std;vector<long long> zaatull(vector<long long> &A){vector<long long> B = A;sort(B.begin(),B.end());B.erase(unique(B.begin(),B.end()),B.end());return B;}using SS = long long;class SegmentTree{public:long long siz = 1,n;vector<SS> dat;SS op(SS a, SS b){return a+b;}SS e(){return 0;}void renew (SS &a,SS x){a = op(a,x);//a = x;//その他.}void make(long long N){while(siz < N) siz *= 2;n = N; dat.resize(siz*2,e());}void make2(long long N,vector<SS> &A){make(N);for(long long i=0; i<N; i++){long long pos = i+siz;dat.at(pos) = A.at(i);}for(long long i=siz-1; i>0; i--) dat.at(i) = op(dat.at(i*2),dat.at(i*2+1));}void update(long long pos,SS x){pos = pos+siz;renew(dat.at(pos),x);while(pos != 1){pos = pos/2;dat.at(pos) = op(dat.at(pos*2),dat.at(pos*2+1));}}SS findans(long long l, long long r){SS retl = e(),retr = e();l += siz,r += siz;while(l < r){if(l&1) retl = op(retl,dat.at(l++));if(r&1) retr = op(dat.at(--r),retr);l >>= 1; r >>= 1;}return op(retl,retr);}SS get(long long pos){return dat.at(pos+siz);}SS rangeans(long long l, long long r){return findans(l,r);}SS allrange(){return dat.at(1);}//rightは) leftは[で 渡す&返す.long long maxright(const function<bool(SS)> f,long long l = 0){l += siz; long long r = n+siz;vector<long long> ls,rs;while(l < r){if(l&1) ls.push_back(l++);if(r&1) rs.push_back(--r);l >>= 1; r >>= 1;}SS okl = e();for(long long i=0; i<ls.size(); i++){l = ls.at(i);SS now = op(okl,dat.at(l));if(!f(now)){while(l < siz){l <<= 1;now = op(okl,dat.at(l));if(f(now)){okl = now; l++;}}return l-siz;}okl = now;}for(long long i=rs.size()-1; i>=0; i--){l = rs.at(i);SS now = op(okl,dat.at(l));if(!f(now)){while(l < siz){l <<= 1;now = op(okl,dat.at(l));if(f(now)){okl = now; l++;}}return l-siz;}okl = now;}return n;}long long minleft(const function<bool(SS)> f,long long r = -1){if(r == -1) r = n;long long l = siz; r += siz;vector<long long> ls,rs;while(l < r){if(l&1) ls.push_back(l++);if(r&1) rs.push_back(--r);l >>= 1; r >>= 1;}SS okr = e();for(long long i=0; i<rs.size(); i++){r = rs.at(i);SS now = op(dat.at(r),okr);if(!f(now)){while(r < siz){r <<= 1; r++;now = op(dat.at(r),okr);if(f(now)){okr = now; r--;}}return r+1-siz;}}for(long long i=ls.size()-1; i>=0; i--){r = ls.at(i);SS now = op(dat.at(r),okr);if(!f(now)){while(r < siz){r <<= 1; r++;now = op(dat.at(r),okr);if(f(now)){okr = now; r--;}}return r+1-siz;}}return 0;}//ノーマルセグ木.一年の時を経て漸く二分探索実装した.};int main() {ios_base::sync_with_stdio(false);cin.tie(nullptr);long long N,Q; cin >> N >> Q;vector<long long> B(N,-1);vector<long long> St;vector<long long> A(N),all;for(auto &a : A) cin >> a;all = A;vector<tuple<long long,long long,long long>> Query(Q,{-1,-1,-1});for(auto &[t,k,x] : Query){cin >> t;if(t == 1) cin >> k >> x,k--,all.push_back(x);if(t == 3) cin >> k,k--;}all = zaatull(all);long long n = all.size();for(auto &a : A) a = lower_bound(all.begin(),all.end(),a)-all.begin();for(auto &[t,k,x] : Query) if(t == 1) x = lower_bound(all.begin(),all.end(),x)-all.begin();SegmentTree Z; Z.make(n);for(auto &a : A) Z.update(a,1);long long border = -1;auto f = [&](SS x){return x<border;};bool sorted = false;for(auto &[t,k,x] : Query){if(t == 1){St.push_back(k);B.at(k) = x;}if(t == 2){vector<pair<int,int>> change;for(auto pos : St){if(B.at(pos) == -1) continue;border = pos+1;long long a1 = Z.maxright(f),a2 = B.at(pos);if(sorted == false) a1 = A.at(pos);change.push_back({a1,a2});B.at(pos) = -1;}for(auto &[a1,a2] : change) Z.update(a1,-1),Z.update(a2,1);sorted = true;St.clear();}if(t == 3){if(B.at(k) != -1) cout << all.at(B.at(k)) << "\n";else if(sorted == false) cout << all.at(A.at(k)) << "\n";else{border = k+1;cout << all.at(Z.maxright(f)) << "\n";}}}}