import sys readline=sys.stdin.readline def SA_IS(lst,compressed=False): N=len(lst) if not compressed: decomp=sorted(list(set(lst))) comp={x:i for i,x in enumerate(decomp)} lst=[comp[s]+1 for s in lst]+[0] else: lst=[s+1 for s in lst]+[0] def induced_sort(N,lst,type,LMS,bucket_count,sorted_LMS_index=None): buckets_left=[[] for b in range(bucket_count)] buckets_right=[[] for b in range(bucket_count)] if sorted_LMS_index==None: for i in LMS: buckets_right[lst[i]].append(i) else: for i in sorted_LMS_index[::-1]: buckets_right[lst[LMS[i]]].append(LMS[i]) for b in range(bucket_count): for i in buckets_left[b]: if i and type[i-1]=="L": buckets_left[lst[i-1]].append(i-1) for i in buckets_right[b][::-1]: if i and type[i-1]=="L": buckets_left[lst[i-1]].append(i-1) buckets_right[b]=[] for b in range(bucket_count-1,-1,-1): for i in buckets_right[b]: if i and type[i-1]=="S": buckets_right[lst[i-1]].append(i-1) for i in buckets_left[b][::-1]: if i and type[i-1]=="S": buckets_right[lst[i-1]].append(i-1) suffix_array=[N] for b in range(1,bucket_count): suffix_array+=buckets_left[b]+buckets_right[b][::-1] return suffix_array stack=[] while N>=1: type=[None]*(N+1) type[N]="S" for i in range(N-1,-1,-1): if lst[i]lst[i+1]: type[i]="L" else: type[i]=type[i+1] LMS=[i for i in range(1,N) if type[i-1]=="L" and type[i]=="S"]+[N] bucket_count=max(lst)+1 stack.append((N,lst,type,LMS,bucket_count)) suffix_array=induced_sort(N,lst,type,LMS,bucket_count) LMS_substring=[None]*(N+1) for i in range(len(LMS)-1): LMS_substring[LMS[i]]=lst[LMS[i]:LMS[i+1]] LMS_substring[N]=lst[N:N+1] num=0 prev=[0] for i in suffix_array: if LMS_substring[i]!=None: if prev!=LMS_substring[i]: prev=LMS_substring[i] num+=1 LMS_substring[i]=num lst=[LMS_substring[i] for i in LMS] N=len(lst)-1 sorted_LMS_index=[0] for N,lst,type,LMS,bucket_count in stack[::-1]: suffix_array=induced_sort(N,lst,type,LMS,bucket_count,sorted_LMS_index) sorted_LMS_index=suffix_array suffix_array=suffix_array[1:] return suffix_array def LCP(lst,suffix_array): N=len(lst) idx=[None]*N for i in range(N): idx[suffix_array[i]]=i lcp=[None]*(N-1) r=0 for l in range(N): if idx[l]==N-1: continue ll=suffix_array[idx[l]+1] while r1: i>>= 1 self.segment_tree[i]=self.f(self.segment_tree[i<<1],self.segment_tree[i<<1|1]) def Build(self,lst): for i,x in enumerate(lst,self.N): self.segment_tree[i]=x for i in range(self.N-1,0,-1): self.segment_tree[i]=self.f(self.segment_tree[i<<1],self.segment_tree[i<<1|1]) def Fold(self,L=None,R=None): if L==None: L=self.N else: L+=self.N if R==None: R=self.N*2 else: R+=self.N vL=self.e vR=self.e while L>=1 R>>=1 return self.f(vL,vR) def Fold_Index(self,L=None,R=None): if L==None: L=self.N else: L+=self.N if R==None: R=self.N*2 else: R+=self.N if L==R: return None x=self.Fold(L-self.N,R-self.N) while L>=1 R>>=1 while i>=1 r>>=1 if f(self.f(vl,vr)): return self.N v=self.e while True: while L%2==0: L>>=1 vv=self.f(v,self.segment_tree[L]) if f(vv): v=vv L+=1 else: while L>=1 r>>=1 if f(self.f(vl,vr)): return 0 v=self.e while True: R-=1 while R>1 and R%2: R>>=1 vv=self.f(self.segment_tree[R],v) if f(vv): v=vv else: while R=R-L) r=ST_min.Bisect_Right(i,lambda le:le>=R-L)+1 ans=ST_sum.Fold(0,l)+(l-i)*(R-L-1)+R-L-1+(r-i-1)*(R-L-1) print(ans)