結果
| 問題 |
No.811 約数の個数の最大化
|
| コンテスト | |
| ユーザー |
utouto97
|
| 提出日時 | 2019-08-18 12:59:53 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,355 bytes |
| コンパイル時間 | 1,681 ms |
| コンパイル使用メモリ | 165,656 KB |
| 実行使用メモリ | 6,820 KB |
| 最終ジャッジ日時 | 2024-10-01 13:55:53 |
| 合計ジャッジ時間 | 2,446 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 4 WA * 8 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define rep(i,l,r) for(int i=(int)(l);i<(int)(r);i++)
#define all(x) (x).begin(),(x).end()
template<class T>bool chmax(T &a,T b){if(a<b){a=b;return 1;}return 0;}
template<class T>bool chmin(T &a,T b){if(a>b){a=b;return 1;}return 0;}
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef vector<vi> vvi;
const int inf = 1LL<<60;
const int mod = 1e9 + 7;
const double eps = 1e-9;
/*{
}*/
int Pow(int b, int p){
if(p == 0) return 1;
return Pow(b*b, p/2)*(p%2 ? b : 1);
}
signed main(){
int n, k;
cin >> n >> k;
int t = n;
vector<pii> pf;
rep(i, 2, t/i+1){
int c = 0;
while(t%i == 0){ c++; t /= i; }
if(c) pf.emplace_back(i, c);
}
if(t != 1) pf.emplace_back(t, 1);
int m = pf.size();
vi v(m), z(n+1);
int cur = 0;
while(k){
if(v[cur] < pf[cur].second){
v[cur]++;
z[pf[cur].first]++;
k--;
}
cur = (cur+1)%m;
}
int ans = 1;
rep(i, 0, m){
ans *= Pow(pf[i].first, v[i]);
}
vi p(n+1, 1);
rep(i, 2, n+1){
rep(j, 2, n/i+1){
p[i*j] = 0;
}
}
rep(i, 0, 30){
for(int j = 2; ans * j < n; j++){
if(p[j] and z[j] <= i){
ans *= j;
z[j]++;
}
}
}
cout << ans << endl;
// rep(i, 0, n+1) if(z[i]) cout << i << " " << z[i] << endl;
return 0;
}
utouto97