結果
問題 | No.811 約数の個数の最大化 |
ユーザー | どらら |
提出日時 | 2019-04-12 22:07:05 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 144 ms / 2,000 ms |
コード長 | 1,135 bytes |
コンパイル時間 | 1,969 ms |
コンパイル使用メモリ | 178,788 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-09-14 18:59:15 |
合計ジャッジ時間 | 3,422 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge6 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 12 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define REP(i,a,n) for(int i=(a); i<(int)(n); i++) #define rep(i,n) REP(i,0,n) #define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it) #define ALLOF(c) (c).begin(), (c).end() typedef long long ll; typedef unsigned long long ull; map<int,int> factorize(int x){ map<int,int> ret; int d, q; while(x>=4 && x%2==0){ret[2]++; x/=2;} d = 3; q = x/d; while(q>=d){ if(x%d==0){ ret[d]++; x=q; }else{ d+=2; } q=x/d; } ret[x]++; return ret; } vector<int> divisor(int n){ vector<int> ret; for(int i=1; i*i<=n; i++){ if(n%i==0){ ret.push_back(i); if(i!=n/i) ret.push_back(n/i); } } //sort(ret.begin(), ret.end()); return ret; } int main(){ int N, K; cin >> N >> K; map<int,int> f = factorize(N); int cnt = 0; int M = -1; REP(i,1,N){ map<int,int> g = factorize(i); vector<int> d = divisor(i); int cmn = 0; FOR(it,g){ cmn += min(f[it->first], it->second); } if(cmn >= K && d.size() > cnt){ cnt = d.size(); M = i; } } cout << M << endl; return 0; }