結果
問題 |
No.2972 確率的素数判定
|
ユーザー |
|
提出日時 | 2025-02-04 00:52:01 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 232 ms / 2,000 ms |
コード長 | 2,193 bytes |
コンパイル時間 | 1,008 ms |
コンパイル使用メモリ | 82,476 KB |
実行使用メモリ | 91,384 KB |
最終ジャッジ日時 | 2025-02-04 00:52:09 |
合計ジャッジ時間 | 6,618 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
import sys import math import bisect from heapq import heapify, heappop, heappush from collections import deque, defaultdict, Counter from functools import lru_cache from itertools import accumulate, combinations, permutations, product sys.set_int_max_str_digits(10 ** 6) sys.setrecursionlimit(1000000) MOD = 10 ** 9 + 7 MOD99 = 998244353 input = lambda: sys.stdin.readline().strip() NI = lambda: int(input()) NMI = lambda: map(int, input().split()) NLI = lambda: list(NMI()) SI = lambda: input() SMI = lambda: input().split() SLI = lambda: list(SMI()) EI = lambda m: [NLI() for _ in range(m)] # 高速エラストテネス sieve[n]はnの最小の素因数 def make_prime_table(n): sieve = list(range(n + 1)) sieve[0] = -1 sieve[1] = -1 for i in range(4, n + 1, 2): sieve[i] = 2 for i in range(3, int(n ** 0.5) + 1, 2): if sieve[i] != i: continue for j in range(i * i, n + 1, i * 2): if sieve[j] == j: sieve[j] = i return sieve prime_table = make_prime_table(100001) # 素数列挙 is_prime = [int(i==p) for i, p in enumerate(prime_table)] # 素因数分解 上のprime_tableと組み合わせて使う def prime_factorize(n): result = [] while n != 1: p = prime_table[n] e = 0 while n % p == 0: n //= p e += 1 result.append((p, e)) return result # Nの素因数分解を辞書で返す(単体) def prime_fact(n): root = int(n**0.5) + 1 prime_dict = {} for i in range(2, root): cnt = 0 while n % i == 0: cnt += 1 n = n // i if cnt: prime_dict[i] = cnt if n != 1: prime_dict[n] = 1 return prime_dict # 約数列挙(単体) def divisors(x): res = set() for i in range(1, int(x**0.5) + 2): if x % i == 0: res.add(i) res.add(x//i) return res def main(): prime_count = list(accumulate(is_prime)) T = NI() for _ in range(T): N, P, Q = NMI() K = prime_count[N] ans = K * P / (K*P + (N-K)*(100-Q)) print(ans) if __name__ == "__main__": main()