結果
| 問題 |
No.811 約数の個数の最大化
|
| コンテスト | |
| ユーザー |
mine691
|
| 提出日時 | 2019-04-12 23:36:20 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 107 ms / 2,000 ms |
| コード長 | 1,078 bytes |
| コンパイル時間 | 1,552 ms |
| コンパイル使用メモリ | 168,792 KB |
| 実行使用メモリ | 5,376 KB |
| 最終ジャッジ日時 | 2024-09-15 06:22:23 |
| 合計ジャッジ時間 | 2,732 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 12 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const ll mod = 1e9 + 7;
const ll inf = (1 << 30) - 1;
const ll infll = (1LL << 61) - 1;
int n, k, ans;
template <typename T>
T gcd(T a, T b)
{
return (b == 0 ? a : gcd(b, a % b));
}
template <typename T>
int divisor(T n)
{
int d = 0;
for (T i = 1; i * i <= n; i++)
{
if (n % i == 0)
{
d++;
if (i * i != n)
d++;
}
}
return d;
}
template <typename T>
vector<T> prime_factor(T n)
{
vector<T> res;
for (T i = 2; i * i <= n; i++)
{
while (n % i == 0)
{
res.push_back(i);
n /= i;
}
}
if (n != 1)
res.push_back(n); //nは素数
return res;
}
int main()
{
cin >> n >> k;
int dm = 0;
for (int i = 1; i < n; i++)
{
int d = divisor(i);
//cout << d << endl;
int g = gcd(i, n);
auto p = prime_factor(g);
if (p.size() >= k && d > dm)
dm = d, ans = i;
}
cout << ans << endl;
}
mine691