結果
問題 | No.2570 最大最大公約数 |
ユーザー |
![]() |
提出日時 | 2023-12-02 16:38:34 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,019 ms / 2,000 ms |
コード長 | 1,513 bytes |
コンパイル時間 | 275 ms |
コンパイル使用メモリ | 82,468 KB |
実行使用メモリ | 73,856 KB |
最終ジャッジ日時 | 2024-11-13 17:14:33 |
合計ジャッジ時間 | 6,191 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 28 |
ソースコード
#緑以下Ex-B#素因数分解 ライブラリclass PrimeFact:def fact(self,N): #Nを素因数分解し、(素因数, 次数) のリストを返すL=[]for i in range(2,N):if i**2>N: breakif N%i==0:c=0while N%i==0: c+=1; N//=iL.append((i,c))return L if N==1 else L+[(N,1)]def div(self,N): #Nの約数を全列挙するLow,High=[],[]for i in range(1,N+1):if i**2>N: breakif N%i==0:Low.append(i)if i**2!=N: High.append(N//i)return Low+High[::-1]def euler(self,N): #オイラーのファイ関数: X**(euler(M))≡1 mod M if gcd(X,M)=1L=self.fact(N); ans=Nfor p,_ in L: ans=ans*(p-1)//preturn ansimport math#入力受取N,K = map(int,input().split())A = list(map(int,input().split()))PF = PrimeFact()ans = 1for first_ope in range(-K,K+1):G = A[0]+first_opeif G<=0: continueD = PF.div(G)nokori = K - abs(first_ope)#判定問題: 残りleft回の操作で、gcd(A) = X にできるか?for X in D:left = nokorifor next in A[1:]:diff = next % Xcnt = 10**18#減らす操作if next - diff > 0:cnt = min(cnt, diff)#増やす操作cnt = min(cnt, X-diff)left -= cntif left>=0: ans = max(ans,X)print(ans)