結果
問題 | No.875 Range Mindex Query |
ユーザー |
👑 ![]() |
提出日時 | 2021-02-10 16:49:43 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 4,912 bytes |
コンパイル時間 | 163 ms |
コンパイル使用メモリ | 82,076 KB |
実行使用メモリ | 104,448 KB |
最終ジャッジ日時 | 2024-07-08 03:12:44 |
合計ジャッジ時間 | 3,866 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 4 WA * 14 |
ソースコード
class Segment_Tree():"""このプログラム内は1-index"""def __init__(self,L,calc,unit,index):"""calcを演算とするリストLのSegment Treeを作成calc:演算(2変数関数,モノイド)unit:モノイドcalcの単位元 (xe=ex=xを満たすe)index:数列の第1要素のindex"""self.calc=calcself.unit=unitself.index=indexN=len(L)d=max(1,(N-1).bit_length())k=1<<dself.data=[unit]*k+L+[unit]*(k-len(L))self.N=kself.depth=dfor i in range(k-1,0,-1):self.data[i]=self.calc(self.data[i<<1],self.data[i<<1|1])def get(self,k,index=1):"""第k要素を取得"""assert 0<=k-index<self.N,"添字が範囲外"return self.data[k-index+self.N]def update(self,k,x,index=1):"""第k要素をxに変え,更新を行う.k:数列の要素x:更新後の値"""assert 0<=k-index<self.N,"添字が範囲外"m=(k-index)+self.Nself.data[m]=xwhile m>1:m>>=1self.data[m]=self.calc(self.data[m<<1],self.data[m<<1|1])def product(self,From,To,index=1,left_closed=True,right_closed=True):L=(From-index)+self.N+(not left_closed)R=(To-index)+self.N+(right_closed)vL=self.unitvR=self.unitwhile L<R:if L&1:vL=self.calc(vL,self.data[L])L+=1if R&1:R-=1vR=self.calc(self.data[R],vR)L>>=1R>>=1return self.calc(vL,vR)def all_product(self):return self.data[1]def max_right(self,left,cond,index=1):"""以下の2つをともに満たすxの1つを返す.\n(1) r=left or cond(data[left]*data[left+1]*...*d[right-1]):True(2) r=N+index or cond(data[left]*data[left+1]*...*data[right]):False※condが単調減少の時,cond(data[left]*...*data[right-1])を満たす最大のrightとなる.cond:関数(引数が同じならば結果も同じ)cond(unit):Trueindex<=left<=r<n+index"""left-=indexassert 0<=left<=self.N,"添字が範囲外"assert cond(self.unit),"単位元が条件を満たさない."if left==self.N:return self.N+indexleft+=self.N-(index-1)sm=self.unitcalc=self.calcfirst=Truewhile first or (left & (-left))!=left:first=Falsewhile left%2==0:left>>=1if not cond(calc(sm,self.data[left])):while left<self.N:left<<=1if cond(self.calc(sm,self.data[left])):sm=self.calc(sm,self.data[left])left+=1return left-self.N+indexsm=self.calc(sm,self.data[left])left+=1return self.N+indexdef min_left(self,right,cond,index=1):"""以下の2つをともに満たすxの1つを返す.\n(1) l=right or cond(data[left]*data[left+1]*...*d[right]):True(2) l=index or cond(data[left-1]*data[left+1]*...*data[right]):False※condが単調減少の時,cond(data[left]*...*data[right-1])を満たす最大のrightとなる.cond:関数(引数が同じならば結果も同じ)cond(unit):Trueindex<=left<=r<n+index"""right-=indexassert 0<=right<=self.N,"添字が範囲外"assert cond(self.unit),"単位元が条件を満たさない."if right==0:return indexright+=self.Nsm=self.unitcalc=self.calcfirst=1while first or (right & (-right))!=right:first=0right-=1while right>1 and right&1:right>>=1if not cond(calc(self.data[right],sm)):while right<self.N:right=2*right+1if cond(calc(self.data[right],sm)):sm=calc(self.data[right],sm)right-=1return right+1-self.N+indexsm=calc(self.data[right],sm)return index#================================================import sysinput=sys.stdin.readlinewrite=sys.stdout.writeN,Q=map(int,input().split())A=["*"]+list(map(int,input().split()))A_inv=[0]*(N+1)for i,a in enumerate(A[1:],1):A_inv[a]=iX=[]S=Segment_Tree(A[1:],min,N+1,1)for _ in range(Q):t,l,r=map(int,input().split())if t==1:p,q=A[l],A[r]A_inv[p],A_inv[q]=r,lS.update(l,q,1)S.update(r,p,1)else:alpha=S.product(l,r)X.append(A_inv[alpha])write("\n".join(map(str,X)))