結果

問題 No.1809 Divide NCK
ユーザー MasKoaTS
提出日時 2021-12-08 22:06:34
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 11 ms / 2,000 ms
コード長 919 bytes
コンパイル時間 785 ms
コンパイル使用メモリ 74,592 KB
最終ジャッジ日時 2025-01-26 06:58:39
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

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) {
	VV<ll> prime = {};
	ll f = 2;
	while (f * 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 legendre(ll N, ll K, V<ll>& p) {
	ll R = N - K;
	ll cnt = 0;
	while (N) {
		N /= p[0];
		cnt += N;
	}
	while (K) {
		K /= p[0];
		cnt -= K;
	}
	while (R) {
		R /= p[0];
		cnt -= R;
	}
	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, legendre(N, K, p));
	}

	cout << ans << endl;
	return 0;
}
0