結果
問題 | No.1809 Divide NCK |
ユーザー | ktr216 |
提出日時 | 2022-01-14 23:12:56 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 13 ms / 2,000 ms |
コード長 | 1,288 bytes |
コンパイル時間 | 1,684 ms |
コンパイル使用メモリ | 170,944 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-11-20 14:05:50 |
合計ジャッジ時間 | 2,969 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 39 |
ソースコード
#include <bits/stdc++.h>using namespace std;using ll = long long;using VI = vector<int>;using P = pair<int, int>;#define REP(i, n) for (int i = 0; i < (int)(n); i++)#define FOR(i, a, b) for (ll i = a; i < (ll)(b); i++)#define ALL(a) (a).begin(),(a).end()constexpr int INF = 1001001001;constexpr ll LINF = 1001001001001001001ll;constexpr int DX[] = {1, 0, -1, 0};constexpr int DY[] = {0, 1, 0, -1};template< typename T1, typename T2>inline bool chmax(T1 &a, T2 b) {return a < b && (a = b, true); }template< typename T1, typename T2>inline bool chmin(T1 &a, T2 b) {return a > b && (a = b, true); }ll f(ll x, ll y) {ll ret = 0;while (x > 0) {ret += x / y;x /= y;}return ret;}int main() {ll N, K, M;cin >> N >> K >> M;vector<pair<ll, int>> A;ll P = M;for (ll i = 2; i * i <= P; i++) {int count = 0;while (P % i == 0) {count++;P /= i;}if (count) A.push_back(make_pair(i, count));}if (P > 1) A.push_back(make_pair(P, 1));ll ans = LINF;for (auto i : A) {ll x = i.first, y = i.second;//cout << x << " " << y << endl;chmin(ans, (f(N, x) - f(N - K, x) - f(K, x)) / y);}cout << ans << endl;}