結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
![]() |
提出日時 | 2020-04-26 01:14:18 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,482 ms / 2,000 ms |
コード長 | 7,186 bytes |
コンパイル時間 | 515 ms |
コンパイル使用メモリ | 82,100 KB |
実行使用メモリ | 172,160 KB |
最終ジャッジ日時 | 2024-09-16 14:05:53 |
合計ジャッジ時間 | 25,197 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
ソースコード
import sys,bisect,string,math,time,functools,randomfrom heapq import heappush,heappop,heapifyfrom collections import deque,defaultdict,Counterfrom itertools import permutations,combinations,groupbydef Golf():*a,=map(int,open(0))def I():return int(input())def S_():return input()def IS():return input().split()def LS():return [i for i in input().split()]def LI():return [int(i) for i in input().split()]def LI_():return [int(i)-1 for i in input().split()]def NI(n):return [int(input()) for i in range(n)]def NI_(n):return [int(input())-1 for i in range(n)]def StoLI():return [ord(i)-97 for i in input()]def ItoS(n):return chr(n+97)def LtoS(ls):return ''.join([chr(i+97) for i in ls])def GI(V,E,ls=None,Directed=False,index=1):org_inp=[];g=[[] for i in range(V)]FromStdin=True if ls==None else Falsefor i in range(E):if FromStdin:inp=LI()a,b,c=(inp+[1])[:3]org_inp.append(inp)else:index=0a,b,c=(ls[i]+[1])[:3]if index==1:a-=1;b-=1aa=(a,c);bb=(b,c);g[a].append(bb)if not Directed:g[b].append(aa)return g,org_inpdef GGI(h,w,search=None,replacement_of_found='.',mp_def={'#':1,'.':0}):#h,w,g,sg=GGI(h,w,search=['S','G'],replacement_of_found='.',mp_def={'#':1,'.':0}) # sample usagemp=[1]*(w+2);found={}for i in range(h):s=input()for char in search:if char in s:found[char]=((i+1)*(w+2)+s.index(char)+1)mp_def[char]=mp_def[replacement_of_found]mp+=[1]+[mp_def[j] for j in s]+[1]mp+=[1]*(w+2)return h+2,w+2,mp,founddef TI(n):return GI(n,n-1)def bit_combination(k,n=2):rt=[]for tb in range(n**k):s=[tb//(n**bt)%n for bt in range(k)];rt+=[s]return rtdef show(*inp,end='\n'):if show_flg:print(*inp,end=end)YN=['YES','NO'];Yn=['Yes','No']mo=10**9+7inf=float('inf')l_alp=string.ascii_lowercase#sys.setrecursionlimit(10**7)input=lambda: sys.stdin.readline().rstrip()class Comb:def __init__(self,n,mo=10**9+7):self.fac=[0]*(n+1)self.inv=[1]*(n+1)self.fac[0]=1self.fact(n)for i in range(1,n+1):self.fac[i]=i*self.fac[i-1]%moself.inv[n]*=iself.inv[n]%=moself.inv[n]=pow(self.inv[n],mo-2,mo)for i in range(1,n):self.inv[n-i]=self.inv[n-i+1]*(n-i+1)%moreturndef fact(self,n):return self.fac[n]def invf(self,n):return self.inv[n]def comb(self,x,y):if y<0 or y>x:return 0return self.fac[x]*self.inv[x-y]*self.inv[y]%moshow_flg=Falseshow_flg=True## Segment Tree #### Initializer Template ### Range Sum: sg=SegTree(n)# Range Minimum: sg=SegTree(n,inf,min,inf)class SegTree:def __init__(self,n,init_val=0,function=lambda a,b:a+b,ide=0):self.n=nself.ide_ele=ide_ele=ideself.num=num=2**(n-1).bit_length()self.seg=seg=[self.ide_ele]*2*self.numself.lazy=lazy=[self.ide_ele]*2*self.numself.segfun=segfun=function#set_valfor i in range(n):self.seg[i+self.num-1]=init_val#builtfor i in range(self.num-2,-1,-1):self.seg[i]=self.segfun(self.seg[2*i+1],self.seg[2*i+2])def update(self,k,x):k += self.num-1self.seg[k] = xwhile k:k = (k-1)//2self.seg[k] = self.segfun(self.seg[k*2+1],self.seg[k*2+2])def evaluate(k,l,r): #遅延評価処理if lazy[k]!=0:node[k]+=lazy[k]if(r-l>1):lazy[2*k+1]+=lazy[k]//2lazy[2*k+2]+=lazy[k]//2lazy[k]=0def query(self,p,q):if q<=p:return self.ide_elep += self.num-1q += self.num-2res=self.ide_elewhile q-p>1:if p&1 == 0:res = self.segfun(res,self.seg[p])if q&1 == 1:res = self.segfun(res,self.seg[q])q -= 1p = p//2q = (q-1)//2if p == q:res = self.segfun(res,self.seg[p])else:res = self.segfun(self.segfun(res,self.seg[p]),self.seg[q])return resdef find_min_index(self,p,q,m):if q<=p:return self.ide_elep += self.num-1q += self.num-2res=self.ide_elewhile q-p>1:if p&1 == 0:res = self.segfun(res,self.seg[p])if q&1 == 1:res = self.segfun(res,self.seg[q])q -= 1p = p//2q = (q-1)//2if p == q:res = self.segfun(res,self.seg[p])else:res = self.segfun(self.segfun(res,self.seg[p]),self.seg[q])return resdef __str__(self):# 生配列を表示rt=self.seg[self.num-1:self.num-1+self.n]return str(rt)class Tree:def __init__(self,inp_size=None,init=True):if init:self.stdin(inp_size)returndef stdin(self,inp_size=None):if inp_size==None:self.size=int(input())else:self.size=inp_sizeself.edges,_=GI(self.size,self.size-1)returndef listin(self,ls):self.size=len(ls)+1self.edges,_=GI(self.size,self.size-1,ls)returndef __str__(self):return str(self.edges)def dfs(self,x,func=lambda prv,nx,dist:prv+dist,root_v=0):q=deque()q.append(x)v=[-1]*self.sizev[x]=root_vwhile q:c=q.pop()for nb,d in self.edges[c]:if v[nb]==-1:q.append(nb)v[nb]=func(v[c],nb,d)return vdef EulerTour(self,x,func=lambda prv,nx,dist:prv+dist,root_v=0):q=deque()q.append((-1,x))v=[None]*self.sizev[x]=root_vet=[]while q:cb,ce=q.pop()et.append(ce)for nb,d in self.edges[ce]:if v[nb]==None:q.append((nb,ce))q.append((ce,nb))v[nb]=func(v[ce],nb,d)vid=[[-1,-1]for i in range(self.size)]for i,j in enumerate(et):if vid[j][0]==-1:vid[j][0]=ielse:vid[j][1]=ireturn v,et,viddef LCA_init(self,depth,et):self.st=SegTree(self.size*2-1,func=min,ide=inf)for i,j in enumerate(et):self.st.update(i,j)self.LCA_init_stat==Truereturndef LCA(self,root,x,y):if self.LCA_init_stat==False:depth,et,vid=self.EulerTour(root)self.LCA_init(depth,et)return self.st.query(x,y+1)n=I()a=LI()sg=SegTree(n,function=math.gcd)for i in range(n):sg.update(i,a[i])ans=0r=0for l in range(n):r=max(l,r)while sg.query(l,r+1)>1 and r<n:r=min(r+1,n)ans+=n-rprint(ans)