結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #

//愚直に素因数分解 & 最大ケース埋め込み

#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;
}
0