結果
問題 | No.1250 汝は倍数なりや? |
ユーザー | kohei2019 |
提出日時 | 2020-12-17 01:58:33 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 382 ms / 1,000 ms |
コード長 | 992 bytes |
コンパイル時間 | 136 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 31,364 KB |
最終ジャッジ日時 | 2024-09-20 05:58:41 |
合計ジャッジ時間 | 5,583 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 49 |
ソースコード
import sys import copy N,H = map(int,input().split()) lsA = list(map(int,input().split())) def factorization(a):#素因数分解√N arr = [] temp = a for i in range(2, int(-(-a**0.5//1))+1): if temp%i==0: cnt=0 while temp%i==0: cnt+=1 temp //= i arr.append([i, cnt]) if temp!=1: arr.append([temp, 1]) if arr==[]: arr.append([a, 1]) return arr #[素因数、個数] if 0 in lsA: print('YES') sys.exit() if H == 1: print('YES') sys.exit() lsp = factorization(H) for i in range(N): if not lsp: print('YES') sys.exit() k = lsA[i] for j in range(len(lsp)): p = lsp[j][0] while k % p == 0: lsp[j][1] -= 1 k //= p lsp2 = [] for j in range(len(lsp)): if lsp[j][1] > 0: lsp2.append(lsp[j]) lsp = copy.copy(lsp2) if not lsp: print('YES') sys.exit() print('NO')