結果
問題 |
No.2972 確率的素数判定
|
ユーザー |
![]() |
提出日時 | 2024-11-29 22:15:36 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 437 ms / 2,000 ms |
コード長 | 885 bytes |
コンパイル時間 | 389 ms |
コンパイル使用メモリ | 82,664 KB |
実行使用メモリ | 78,600 KB |
最終ジャッジ日時 | 2024-11-29 22:15:42 |
合計ジャッジ時間 | 5,859 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
class Prime(): def __init__(self, N): self.N = N self.Is_prime = [1 for i in range(N + 1)] self.Is_prime[0] = self.Is_prime[1] = 0 self.P = [] for p in range(2, N + 1): if self.Is_prime[p] == 0: continue self.P.append(p) for d in range(2, N//p + 1): self.Is_prime[p * d] = 0 q = 1 r = p while (p * d) % (p * r) == 0: q += 1 r *= p def is_prime(self, p): return self.Is_prime[p] def primelist(self): return self.P P = Prime(100000) C = [0 for i in range(100001)] for i in range(2, 100001): if P.is_prime(i): C[i] = 1 C[i] += C[i - 1] for i in range(int(input())): n, p, q = map(int, input().split()) if n == 1: print(0) continue if p == 0 and q == 100: print(0) continue x = p * C[n] y = (100 - q) * (n - C[n]) print(x / (x + y))