結果
問題 |
No.811 約数の個数の最大化
|
ユーザー |
![]() |
提出日時 | 2020-07-20 07:25:19 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 93 ms / 2,000 ms |
コード長 | 957 bytes |
コンパイル時間 | 2,306 ms |
コンパイル使用メモリ | 200,028 KB |
最終ジャッジ日時 | 2025-01-12 01:12:16 |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 12 |
ソースコード
#include <bits/stdc++.h> using namespace std; constexpr long lim = 100001; vector<vector<int>> pr(lim), id(lim); vector<bool> is_pr(lim, true); void all_fac() { is_pr.at(0) = 0; is_pr.at(1) = 0; for (int i = 2; i < lim; i++) { if (is_pr.at(i)) { pr.at(i).push_back(i); id.at(i).push_back(1); for (int j = 2 * i; j < lim; j += i) { is_pr.at(j) = false; int cnt = 0; int nw = j; while(nw % i == 0) { nw /= i; cnt++; } pr.at(j).push_back(i); id.at(j).push_back(cnt); } } } } int main() { all_fac(); int N, K; cin >> N >> K; int mx = 0, ans = 0; for (int i = 1; i < N; i++) { int tmp = gcd(i, N); int tmp2 = accumulate(id.at(tmp).begin(), id.at(tmp).end(), 0); if (tmp2 >= K) { int tmp3 = 1; for (auto a : id.at(i)) tmp3 *= (a + 1); if (tmp3 > mx) mx = tmp3, ans = i; } } cout << ans << "\n"; }