結果
問題 | No.2697 Range LIS Query |
ユーザー |
|
提出日時 | 2024-03-23 00:27:42 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 5,340 ms / 10,000 ms |
コード長 | 4,804 bytes |
コンパイル時間 | 2,440 ms |
コンパイル使用メモリ | 221,544 KB |
実行使用メモリ | 134,580 KB |
最終ジャッジ日時 | 2024-09-30 12:51:07 |
合計ジャッジ時間 | 49,519 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 15 |
ソースコード
#include<bits/stdc++.h>using namespace std;//#pragma GCC optimize("Ofast")#define rep(i,n) for(ll i=0;i<n;i++)#define repl(i,l,r) for(ll i=(l);i<(r);i++)#define per(i,n) for(ll i=(n)-1;i>=0;i--)#define perl(i,r,l) for(ll i=r-1;i>=l;i--)#define fi first#define se second#define pb push_back#define ins insert#define pqueue(x) priority_queue<x,vector<x>,greater<x>>#define all(x) (x).begin(),(x).end()#define CST(x) cout<<fixed<<setprecision(x)#define vtpl(x,y,z) vector<tuple<x,y,z>>#define rev(x) reverse(x);using ll=long long;using vl=vector<ll>;using vvl=vector<vector<ll>>;using pl=pair<ll,ll>;using vpl=vector<pl>;using vvpl=vector<vpl>;const ll MOD=1000000007;const ll MOD9=998244353;const int inf=1e9+10;const ll INF=4e18;const ll dy[9]={1,0,-1,0,1,1,-1,-1,0};const ll dx[9]={0,1,0,-1,1,-1,1,-1,0};template<class T> inline bool chmin(T& a, T b) {if (a >= b) {a = b;return true;}return false;}template<class T> inline bool chmax(T& a, T b) {if (a < b) {a = b;return true;}return false;}template <typename X, typename M>//faがint lenになってるstruct SegTreeLazy {//遅延セグ木 単位元に注意(updateなら選ばれない数、affineなら(1,0))using FX = function<X(X, X)>;using FA = function<X(X, M, int)>;using FM = function<M(M, M)>;int n;FX fx;FA fa;FM fm;const X ex;const M em;vector<X> dat;vector<M> lazy;SegTreeLazy(int n_, FX fx_, FA fa_, FM fm_, X ex_, M em_): n(), fx(fx_), fa(fa_), fm(fm_), ex(ex_), em(em_), dat(n_ * 4, ex), lazy(n_ * 4, em) {int x = 1;while (n_ > x) x *= 2;n = x;}void set(int i, X x) { dat[i + n - 1] = x; }void build() {for (int k = n - 2; k >= 0; k--) dat[k] = fx(dat[2 * k + 1], dat[2 * k + 2]);}/* lazy eval */void eval(int k, int len) {if (lazy[k] == em) return; // 更新するものが無ければ終了if (k < n - 1) { // 葉でなければ子に伝搬lazy[k * 2 + 1] = fm(lazy[k * 2 + 1], lazy[k]);lazy[k * 2 + 2] = fm(lazy[k * 2 + 2], lazy[k]);}// 自身を更新dat[k] = fa(dat[k],lazy[k],len);//fa(dat[k], fp(lazy[k], len));lazy[k] = em;}void update(int a, int b, M x, int k, int l, int r) {eval(k, r - l);if (a <= l && r <= b) { // 完全に内側の時lazy[k] = fm(lazy[k], x);eval(k, r - l);} else if (a < r && l < b) { // 一部区間が被る時update(a, b, x, k * 2 + 1, l, (l + r) / 2); // 左の子update(a, b, x, k * 2 + 2, (l + r) / 2, r); // 右の子dat[k] = fx(dat[k * 2 + 1], dat[k * 2 + 2]);}}void update(int a, int b, M x) { update(a, b, x, 0, 0, n); }X query_sub(int a, int b, int k, int l, int r) {eval(k, r - l);if (r <= a || b <= l) { // 完全に外側の時return ex;} else if (a <= l && r <= b) { // 完全に内側の時return dat[k];} else { // 一部区間が被る時X vl = query_sub(a, b, k * 2 + 1, l, (l + r) / 2);X vr = query_sub(a, b, k * 2 + 2, (l + r) / 2, r);return fx(vl, vr);}}X query(int a, int b) { return query_sub(a, b, 0, 0, n); }X operator[](int i){return query(i,i+1);}};int main(){ll n;cin >> n;vvl emp(4,vl(4));auto fx=[](vvl a,vvl b){vvl c(4,vl(4));rep(i,4){repl(j,i,4){repl(k,j,4){repl(l,k,4){chmax(c[i][l],a[i][j]+b[k][l]);}}}}return c;};auto fa=[](vvl a,ll nb,int len){vvl b(4,vl(4));b[nb][nb]=len;/*vvl c(4,vl(4));rep(i,4){repl(j,i,4){repl(k,j,4){repl(l,k,4){chmax(c[i][l],a[i][j]+b[k][l]);}}}}*/return b;};auto fm=[](ll a,ll b){return b;};SegTreeLazy<vvl,ll> st(n,fx,fa,fm,emp,-1);vector<vvl> sg(4,vvl(4,vl(4)));rep(i,4)sg[i][i][i]=1;rep(i,n){ll a;cin >> a;a--;st.set(i,sg[a]);}st.build();ll q;cin >> q;while(q--){ll t;cin >> t;if(t==1){ll l,r;cin >> l >> r;l--;auto f=st.query(l,r);ll ans=0;rep(i,4)rep(j,4)chmax(ans,f[i][j]);cout << ans << endl;}else{ll l,r,x;cin >> l >> r >> x;l--;x--;st.update(l,r,x);}}}