結果

問題 No.811 約数の個数の最大化
ユーザー betrue12betrue12
提出日時 2019-04-12 21:41:04
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 53 ms / 2,000 ms
コード長 1,308 bytes
コンパイル時間 1,902 ms
コンパイル使用メモリ 181,744 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-14 17:46:06
合計ジャッジ時間 2,890 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 12
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:59:13: warning: 'ans' may be used uninitialized [-Wmaybe-uninitialized]
   59 |     cout << ans << endl;
      |             ^~~
main.cpp:34:9: note: 'ans' was declared here
   34 |     int ans;
      |         ^~~

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

vector<int64_t> sieve(int64_t n){
    vector<int64_t> ret;
    vector<bool> is_prime(n+1);
    for(int64_t i=2; i<=n; i++) is_prime[i] = true;
    for(int64_t i=2; i<=n; i++){
        if(is_prime[i]){
            ret.push_back(i);
            for(int64_t j=2*i; j<=n; j+=i) is_prime[j] = false;
        }
    }
    return ret;
}

int main(){
    int N, K;
    cin >> N >> K;

    auto primes = sieve(100000);
    map<int, int> pn;
    int num = N;
    for(int64_t p : primes){
        if(p*p > num) break;
        while(num % p == 0){
            pn[p]++;
            num /= p;
        }
    }
    if(num > 1) pn[num]++;

    int mx = 0;
    int ans;
    for(int i=1; i<N; i++){
        map<int, int> pi;
        int num = i;
        for(int64_t p : primes){
            if(p*p > num) break;
            while(num % p == 0){
                pi[p]++;
                num /= p;
            }
        }
        if(num > 1) pi[num]++;

        int k = 0;
        for(auto& p : pn){
            k += min(p.second, pi[p.first]);
        }
        if(k < K) continue;
        int result = 1;
        for(auto& p : pi) result *= 1+p.second;
        if(mx < result){
            mx = result;
            ans = i;
        }
    }
    cout << ans << endl;
    return 0;
}
0