結果

問題 No.2751 429-like Number
ユーザー miho-4miho-4
提出日時 2024-05-10 22:03:46
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 1,039 bytes
コンパイル時間 178 ms
コンパイル使用メモリ 82,368 KB
実行使用メモリ 156,996 KB
最終ジャッジ日時 2024-05-10 22:04:08
合計ジャッジ時間 16,760 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 440 ms
141,764 KB
testcase_01 AC 435 ms
142,300 KB
testcase_02 RE -
testcase_03 RE -
testcase_04 AC 442 ms
141,692 KB
testcase_05 AC 469 ms
142,284 KB
testcase_06 RE -
testcase_07 RE -
testcase_08 AC 738 ms
156,268 KB
testcase_09 AC 572 ms
153,220 KB
testcase_10 AC 545 ms
152,556 KB
testcase_11 WA -
testcase_12 RE -
testcase_13 RE -
testcase_14 AC 576 ms
152,360 KB
testcase_15 AC 725 ms
155,248 KB
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import random

def mr_test_normal(n: int):
    assert n > 1
    if n == 2:
        return True
    elif n%2 == 0:
        return False
    d = n - 1
    s = 0
    while d & 1 == 0:
        s += 1; d >>= 1
    
    k = 20
    for _ in range(k):
        a = random.randint(2, n-1)
        x = pow(a, d, n) # a^d mod n
        if x == 1 or x == n-1:
            continue
        for _ in range(s):
            x = pow(x, 2, n)
            if x == n-1: # a^(d*2^r) mod n
                break
        else: return False
    return True

P=[]
prime = [True for i in range(10**7+1)]
p = 2

while (p * p <= 10**7+1):
  if prime[p] ==1:
    P.append(p)
    for i in range(p * 2, 10**7+1, p):
      prime[i] = 0
  p += 1

S=set(P)
l=len(P)
t=int(input())
for _ in range(t):
  n=int(input())
  flag=0
  for i in range(l):
    if n%P[i]==0:
      n//=P[i]
      for j in range(i,l):
        if n%P[j]==0:
          n//=P[j]
          if mr_test_normal(n):
            flag=1
            break
    if flag:
      break
  print("Yes" if flag else "No")
0