結果
問題 |
No.1809 Divide NCK
|
ユーザー |
![]() |
提出日時 | 2022-01-12 21:48:11 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 56 ms / 2,000 ms |
コード長 | 667 bytes |
コンパイル時間 | 428 ms |
コンパイル使用メモリ | 82,080 KB |
実行使用メモリ | 57,728 KB |
最終ジャッジ日時 | 2024-11-15 14:02:49 |
合計ジャッジ時間 | 3,400 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 39 |
ソースコード
N,K,M = map(int,input().split()) if not(1 <= K <= N <= 10 ** 18 and 2 <= M <= 10 ** 12): exit(1) def prime_factorize(n): prime = [] f = 2 while(f * f <= n): if(n % f == 0): n //= f cnt = 1 while(n % f == 0): n //= f cnt += 1 prime.append((f,cnt)) else: f += 1 if(n != 1): prime.append((n,1)) return prime def kummer(N,K,p,q): R = N - K cnt = c = 0 while(R or K): a = R % p b = K % p if(a + b + c >= p): cnt += 1 c = (a + b + c) // p else: c = 0 R //= p K //= p return cnt // q prime = prime_factorize(M) ans = 10 ** 18 for p,c in prime: cnt = kummer(N,K,p,c) if(ans > cnt): ans = cnt print(ans)