import sys,collections,bisect,operator def MillerRabin_64bit(x): if x <= 1: return False if x == 2: return True if x % 2 == 0: return False if x < 1 << 30: test = [2, 7, 61] else: test = [2,325,9375,28178,450775,9780504,1795265022] s = 0 d = x - 1 while d & 1 == 0: s += 1 d >>= 1 for i in test: if x <= i: return True y = pow(i,d,x) if y == 1 or y == x - 1: continue c = 0 for j in range(s-1): y = y * y % x if y == x - 1: c = 1 break if not c: return False return True def eratosthenes(n): is_prime = [False, False] + [True] * (n-1) for p in range(2, n+1): if not(is_prime[p]): continue for k in range(p*2, n+1, p): is_prime[k] = False return is_prime #https://github.com/shakayami/ACL-for-python/blob/master/segtree.py class segtree(): n=1 size=1 log=2 d=[0] op=None e=10**15 def __init__(self,V,OP,E): self.n=len(V) self.op=OP self.e=E self.log=(self.n-1).bit_length() self.size=1<>i) def get(self,p): assert 0<=p and p>=1 r>>=1 return self.op(sml,smr) def all_prod(self): return self.d[1] def max_right(self,l,f): assert 0<=l and l<=self.n assert f(self.e) if l==self.n: return self.n l+=self.size sm=self.e while(1): while(l%2==0): l>>=1 if not(f(self.op(sm,self.d[l]))): while(l1 and (r%2)): r>>=1 if not(f(self.op(self.d[r],sm))): while(r