結果
問題 | No.811 約数の個数の最大化 |
ユーザー |
|
提出日時 | 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; | ^~~
ソースコード
#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;}