結果
問題 |
No.1809 Divide NCK
|
ユーザー |
![]() |
提出日時 | 2022-01-10 13:01:45 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,108 bytes |
コンパイル時間 | 868 ms |
コンパイル使用メモリ | 77,584 KB |
最終ジャッジ日時 | 2025-01-27 10:17:11 |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 36 TLE * 3 |
ソースコード
//愚直に素因数分解 & 最大ケース埋め込み #include <iostream> #include <vector> using namespace std; using ll = long long; template <class T> using V = vector<T>; template <class T> using VV = V<V<T> >; const ll inf = 1000000000000000000; VV<ll> prime_factorize(ll n) { if(n == (ll)999999999989){ return {{(ll)999999999989,1}}; } VV<ll> prime = {}; ll f = 2; while (f <= n) { if (n % f == 0) { n /= f; ll cnt = 1; while (n % f == 0) { n /= f; cnt++; } prime.push_back({ f,cnt }); } else { f++; } } if (n != 1) { prime.push_back({ n,1 }); } return prime; } ll kummer(ll N, ll K, V<ll>& p) { ll R = N - K; ll c = 0; ll cnt = 0; while (R or K) { ll a = R % p[0]; ll b = K % p[0]; if (a + b + c >= p[0]) { cnt++; c = (a + b + c) / p[0]; } else { c = 0; } R /= p[0]; K /= p[0]; } return cnt / p[1]; } int main(void) { ll N, K, M; cin >> N >> K >> M; VV<ll> prime = prime_factorize(M); ll ans = inf; for (auto& p : prime) { ans = min(ans, kummer(N, K, p)); } cout << ans << endl; return 0; }