結果

問題 No.1250 汝は倍数なりや?
コンテスト
ユーザー shotoyoo
提出日時 2020-10-09 22:05:22
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 1,526 bytes
コンパイル時間 194 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 104,064 KB
最終ジャッジ日時 2024-07-20 11:31:12
合計ジャッジ時間 5,539 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 32 RE * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = lambda : sys.stdin.readline().rstrip()
sys.setrecursionlimit(max(1000, 10**9))
write = lambda x: sys.stdout.write(x+"\n")


n,h = list(map(int, input().split()))
a = list(map(int, input().split()))
def hurui(n):
    pl = []
    mpf = [None]*(n+1)
    for d in range(2,n):
        if mpf[d] is None:
            mpf[d] = d
            pl.append(d)
        for p in pl:
            if p*d>n or p>mpf[d]:
                break
            mpf[p*d] = p
    return pl, mpf
from collections import defaultdict
def factor2(n, m=None):
    # mを与えると、高々その素因数まで見て、残りは分解せずにそのまま出力する
    f = {}
    tmp = n
    M = int(-(-n**0.5//1))+1
    if m is not None:
        M = min(m+1, M)
    for i in range(2, M+1):
        if tmp<i:
            break
        if tmp%i==0:
            cnt=0
            while tmp%i==0:
                cnt+=1
                tmp //= i
            f[i] = cnt
    if tmp!=1:
        f[tmp] = 1
    if not f:
        f[n] = 1
    return f

if 0 in a or h==1:
    ans = "YES"
else:
    def factor(num):
        d = defaultdict(int)
        while num>1:
            d[mpf[num]] += 1
            num //= mpf[num]
        return d
    pl, mpf = hurui(max(abs(item) for item in a)+10)
    d = defaultdict(int)
    for num in a:
        dd = factor(abs(num))
        for k,v in dd.items():
            d[k] += v
    d2 = factor2(h)
    if all(d2[k]<=d[k] for k in d2.keys()):
        ans = "YES"
    else:
        ans = "NO"
print(ans)
0