結果
| 問題 |
No.811 約数の個数の最大化
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2022-10-02 20:15:48 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 11 ms / 2,000 ms |
| コード長 | 852 bytes |
| コンパイル時間 | 2,182 ms |
| コンパイル使用メモリ | 195,976 KB |
| 最終ジャッジ日時 | 2025-02-07 20:38:40 |
|
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 12 |
ソースコード
#include <bits/stdc++.h>
#define rep(i,n) for(int i = 0; i < (n); i++)
using namespace std;
typedef long long ll;
int main(){
cin.tie(0);
ios::sync_with_stdio(0);
// 共通の素因数の数がK -> 最大公約数の素因数の数がK
// O(NlogN)
int N,K; cin >> N >> K;
vector<int> div_cnt(N + 1, 2);
div_cnt[1] = 1;
for(int i = 2; i <= N; i++)
for(int j = i * 2; j <= N; j += i)
div_cnt[j]++;
vector<int> pf_cnt(N + 1, 1);
pf_cnt[1] = 0;
for(int i = 2; i * i <= N; i++) {
if(pf_cnt[i] == 1) {
for(int j = 2; i * j <= N; j++) {
pf_cnt[i * j] = pf_cnt[i] + pf_cnt[j];
}
}
}
int M = 0;
for(int i = 2; i < N; i++)
if(pf_cnt[gcd(i, N)] >= K && div_cnt[M] < div_cnt[i]) M = i;
cout << M << endl;
}