結果
| 問題 |
No.1809 Divide NCK
|
| コンテスト | |
| ユーザー |
MasKoaTS
|
| 提出日時 | 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;
}
MasKoaTS