結果
問題 | No.2697 Range LIS Query |
ユーザー |
![]() |
提出日時 | 2024-03-22 22:10:06 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 7,030 ms / 10,000 ms |
コード長 | 4,975 bytes |
コンパイル時間 | 495 ms |
コンパイル使用メモリ | 82,176 KB |
実行使用メモリ | 274,228 KB |
最終ジャッジ日時 | 2024-09-30 11:39:59 |
合計ジャッジ時間 | 67,053 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 15 |
ソースコード
class lazy_segtree():def update(self,k):self.d[k]=self.op(self.d[2*k],self.d[2*k+1])def all_apply(self,k,f):self.d[k]=self.mapping(f,self.d[k])if (k<self.size):self.lz[k]=self.composition(f,self.lz[k])def push(self,k):self.all_apply(2*k,self.lz[k])self.all_apply(2*k+1,self.lz[k])self.lz[k]=self.identitydef __init__(self,V,OP,E,MAPPING,COMPOSITION,ID):self.n=len(V)self.log=(self.n-1).bit_length()self.size=1<<self.logself.d=[E for i in range(2*self.size)]self.lz=[ID for i in range(self.size)]self.e=Eself.op=OPself.mapping=MAPPINGself.composition=COMPOSITIONself.identity=IDfor i in range(self.n):self.d[self.size+i]=V[i]for i in range(self.size-1,0,-1):self.update(i)def set(self,p,x):assert 0<=p and p<self.np+=self.sizefor i in range(self.log,0,-1):self.push(p>>i)self.d[p]=xfor i in range(1,self.log+1):self.update(p>>i)def get(self,p):assert 0<=p and p<self.np+=self.sizefor i in range(self.log,0,-1):self.push(p>>i)return self.d[p]def prod(self,l,r):assert 0<=l and l<=r and r<=self.nif l==r:return self.el+=self.sizer+=self.sizefor i in range(self.log,0,-1):if (((l>>i)<<i)!=l):self.push(l>>i)if (((r>>i)<<i)!=r):self.push(r>>i)sml,smr=self.e,self.ewhile(l<r):if l&1:sml=self.op(sml,self.d[l])l+=1if r&1:r-=1smr=self.op(self.d[r],smr)l>>=1r>>=1return self.op(sml,smr)def all_prod(self):return self.d[1]def apply_point(self,p,f):assert 0<=p and p<self.np+=self.sizefor i in range(self.log,0,-1):self.push(p>>i)self.d[p]=self.mapping(f,self.d[p])for i in range(1,self.log+1):self.update(p>>i)def apply(self,l,r,f):assert 0<=l and l<=r and r<=self.nif l==r:returnl+=self.sizer+=self.sizefor i in range(self.log,0,-1):if (((l>>i)<<i)!=l):self.push(l>>i)if (((r>>i)<<i)!=r):self.push((r-1)>>i)l2,r2=l,rwhile(l<r):if (l&1):self.all_apply(l,f)l+=1if (r&1):r-=1self.all_apply(r,f)l>>=1r>>=1l,r=l2,r2for i in range(1,self.log+1):if (((l>>i)<<i)!=l):self.update(l>>i)if (((r>>i)<<i)!=r):self.update((r-1)>>i)def max_right(self,l,g):assert 0<=l and l<=self.nassert g(self.e)if l==self.n:return self.nl+=self.sizefor i in range(self.log,0,-1):self.push(l>>i)sm=self.ewhile(1):while(l%2==0):l>>=1if not(g(self.op(sm,self.d[l]))):while(l<self.size):self.push(l)l=(2*l)if (g(self.op(sm,self.d[l]))):sm=self.op(sm,self.d[l])l+=1return l-self.sizesm=self.op(sm,self.d[l])l+=1if (l&-l)==l:breakreturn self.ndef min_left(self,r,g):assert (0<=r and r<=self.n)assert g(self.e)if r==0:return 0r+=self.sizefor i in range(self.log,0,-1):self.push((r-1)>>i)sm=self.ewhile(1):r-=1while(r>1 and (r%2)):r>>=1if not(g(self.op(self.d[r],sm))):while(r<self.size):self.push(r)r=(2*r+1)if g(self.op(self.d[r],sm)):sm=self.op(self.d[r],sm)r-=1return r+1-self.sizesm=self.op(self.d[r],sm)if (r&-r)==r:breakreturn 0e=[0]*17def f(x,y):return 4*x+ydef op(x,y):z=[0]*17z[16]=x[16]+y[16]for l1 in range(4):for r1 in range(l1,4):for l2 in range(r1,4):for r2 in range(l2,4):z[f(l1,r2)]=max(z[f(l1,r2)],x[f(l1,r1)]+y[f(l2,r2)])return zID=-1def mapping(val,x):if val==-1:return xres=[0]*17res[16]=x[16]res[f(val,val)]=x[16]return resdef composition(F,G):if F==ID:return Greturn Ffrom sys import stdininput=lambda :stdin.readline()[:-1]N=int(input())A=list(map(lambda x:int(x)-1,input().split()))init=[[0]*17 for i in range(N)]for i in range(N):init[i][16]=1init[i][f(A[i],A[i])]=1seg=lazy_segtree(init,op,e,mapping,composition,ID)for _ in range(int(input())):query=list(map(int,input().split()))if query[0]==1:l,r=query[1:]print(max(seg.prod(l-1,r)[:16]))else:l,r,x=query[1:]seg.apply(l-1,r,x-1)