結果

問題 No.811 約数の個数の最大化
ユーザー mine691mine691
提出日時 2019-04-12 23:11:02
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,230 bytes
コンパイル時間 1,964 ms
コンパイル使用メモリ 180,240 KB
実行使用メモリ 4,356 KB
最終ジャッジ日時 2023-10-13 09:08:25
合計ジャッジ時間 2,899 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge15
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,348 KB
testcase_01 AC 1 ms
4,352 KB
testcase_02 AC 3 ms
4,348 KB
testcase_03 AC 2 ms
4,348 KB
testcase_04 AC 1 ms
4,352 KB
testcase_05 WA -
testcase_06 WA -
testcase_07 AC 2 ms
4,348 KB
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 AC 4 ms
4,352 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;

vector<int> v;

int ans = 1;
int d = 1;
map<int, int> mp;
queue<int> que;

//map型
#include <map>
template <typename T>
map<T, int> prime_factor(T n)
{
    map<T, int> res;
    for (T i = 2; i * i <= n; i++)
    {
        bool f = false;
        while (n % i == 0)
        {
            f = true;
            res[i]++;
            n /= i;
        }
        if (f == true)
            v.push_back(i);
    }
    if (n != 1)
    {
        res[n]++; //nは素数
        v.push_back(n);
    }
    return res;
}

// 𝑂(𝑁*log(log 𝑁))
vector<int> eratosthenes(int N)
{
    vector<int> res;
    vector<bool> prime(N + 1, 1);
    prime[0] = prime[1] = false;
    for (int i = 2; i * i <= N; ++i)
    {
        if (prime[i] == true)
        {
            for (int j = 2 * i; j <= N; j += i)
            {
                prime[j] = false;
            }
        }
    }
    for (int i = 0; i <= N; ++i)
    {
        if (prime[i] == true)
            res.push_back(i);
    }
    return res;
}

int n, k;

int main()
{
    cin >> n >> k;
    auto p = prime_factor(n);
    int vs = v.size();
    for (int i = 0; i < vs; i++)
    {
        que.push(v[i]);
    }
    for (int i = 0; i < k; i++)
    {
        int j = que.front();
        que.pop();
        p[j]--, mp[j]++;
        if (p[j] > 0)
            que.push(j);
    }
    for (auto e : v)
    {
        for (int i = 0; i < mp[e]; i++)
        {
            ans *= e;
        }
    }
    //k個選んだ
    //cout << ans << endl;
    auto prime = eratosthenes(n);

    for (int i = 0; i < prime.size(); i++)
        d *= (mp[prime[i]] + 1);

    while (ans * 2 < n)
    {
        int m = 0, t = 1;
        for (int i = 0; i < prime.size() && prime[i] * ans < n && prime[i] < 30; i++)
        {
            //cout << d / (mp[prime[i]] + 1) * (mp[prime[i]] + 2) << endl;
            if (d / (mp[prime[i]] + 1) * (mp[prime[i]] + 2) > m)
                m = d / (mp[prime[i]] + 1) * (mp[prime[i]] + 2), t = prime[i];
        }
        d = m, ans *= t, mp[t]++;
        //cout << ans << endl;
    }
    cout << ans << endl;
}
0