結果
| 問題 |
No.811 約数の個数の最大化
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2020-07-03 20:23:57 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 16 ms / 2,000 ms |
| コード長 | 1,053 bytes |
| コンパイル時間 | 1,890 ms |
| コンパイル使用メモリ | 200,344 KB |
| 最終ジャッジ日時 | 2025-01-11 14:10:54 |
|
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 12 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
#define REP(i,n) for(int i=0; i<(int)(n); i++)
#define FOR(i,b,e) for (int i=(int)(b); i<(int)(e); i++)
#define ALL(x) (x).begin(), (x).end()
const double PI = acos(-1);
map<long long, int> factorize(long long n) {
map<long long, int> ret;
for (long long i = 2; i * i <= n; i++) {
while (n % i == 0) {
++ret[i];
n /= i;
}
}
if (n != 1)
ret[n] = 1;
return ret;
}
int divs[100000+1];
int main() {
ios_base::sync_with_stdio(0);
cin.tie(0);
for (int i = 1; i <= 100000; i++) {
for (int j = i; j <= 100000; j += i)
divs[j]++;
}
int n, k;
cin >> n >> k;
auto fac = factorize(n);
int ret = 0;
for (int i = 1; i < n; i++) {
int x = i;
int cnt = 0;
for (auto &pr: fac) {
REP (j, pr.second) {
if (x % pr.first == 0) {
++cnt;
x /= pr.first;
}
}
}
if (cnt < k) continue;
if (ret == 0 || divs[i] > divs[ret])
ret = i;
}
cout << ret << endl;
return 0;
}