結果
問題 | No.1809 Divide NCK |
ユーザー |
![]() |
提出日時 | 2022-01-14 21:47:11 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,461 bytes |
コンパイル時間 | 1,979 ms |
コンパイル使用メモリ | 198,304 KB |
最終ジャッジ日時 | 2025-01-27 11:14:35 |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 38 WA * 1 |
ソースコード
#include <bits/stdc++.h>using namespace std;#include <atcoder/modint>using namespace atcoder;using mint = modint998244353;//using mint = modint1000000007;using ll = long long;using ld = long double;using pll = pair<ll, ll>;using tlll = tuple<ll, ll, ll>;constexpr ll INF = 1LL << 60;template<class T> bool chmin(T& a, T b) {if (a > b) {a = b; return true;} return false;}template<class T> bool chmax(T& a, T b) {if (a < b) {a = b; return true;} return false;}ll safemod(ll A, ll M) {return (A % M + M) % M;}ll divfloor(ll A, ll B) {if (B < 0) {return divfloor(-A, -B);} return (A - safemod(A, B)) / B;}ll divceil(ll A, ll B) {if (B < 0) {return divceil(-A, -B);} return divfloor(A + B - 1, B);}#define FINALANS(A) do {cout << (A) << '\n'; exit(0);} while (false)vector<pll> factorize(ll n){vector<pll> ret;for (ll p = 2; p * p <= n; p++){ll e = 0;while (n % p == 0){n /= p;e++;}if (e > 0)ret.push_back(make_pair(p, e));}if (n > 1)ret.push_back(make_pair(n, 1));return ret;}ll f(ll n, ll p) // n! は p で何回割り切れるか{ll ret = 0;for (ll k = p; k <= n; k *= p){ret += n / k;}return ret;}int main(){ll N, K, M;cin >> N >> K >> M;auto PE = factorize(M);ll ans = INF;for (auto &[p, e] : PE){ll e2 = f(N, p) - f(N - K, p) - f(K, p);ll tmp = max(0LL, e2 / e);chmin(ans, tmp);}cout << ans << endl;}