結果
問題 | No.1809 Divide NCK |
ユーザー |
|
提出日時 | 2022-01-14 22:41:55 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 93 ms / 2,000 ms |
コード長 | 994 bytes |
コンパイル時間 | 372 ms |
コンパイル使用メモリ | 82,232 KB |
実行使用メモリ | 77,724 KB |
最終ジャッジ日時 | 2024-11-20 13:01:03 |
合計ジャッジ時間 | 4,771 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 39 |
ソースコード
# 素数列挙from collections import Counterdef prime_numbers(x: int) -> list:if x < 2:return []prime_numbers = [i for i in range(x)]prime_numbers[1] = 0for prime_number in prime_numbers:if prime_number ** 2 > x:breakif prime_number == 0:continuefor composite_number in range(2 * prime_number, x, prime_number):prime_numbers[composite_number] = 0return [prime_number for prime_number in prime_numbers if prime_number != 0]# 素因数分解def prime_factorize(n: int) -> list:res = []while n % 2 == 0:res.append(2)n //= 2f = 3while f ** 2 <= n:if n % f == 0:res.append(f)n //= felse:f += 2if n != 1:res.append(n)return resn, k, m = map(int, input().split())primes = prime_numbers(10 ** 6 + 10)p = prime_factorize(m)c = Counter(p)f = set(p)ans = 10 ** 18for i in f:cnt = 0cur = iwhile cur <= n:cnt += n // cur - (n - k) // cur - k // curcur *= ians = min(ans, cnt // c[i])print(ans)