結果

問題 No.811 約数の個数の最大化
ユーザー chocorusk
提出日時 2019-04-12 21:33:05
言語 C++11(廃止可能性あり)
(gcc 13.3.0)
結果
AC  
実行時間 88 ms / 2,000 ms
コード長 1,603 bytes
コンパイル時間 1,453 ms
コンパイル使用メモリ 114,392 KB
実行使用メモリ 20,992 KB
最終ジャッジ日時 2024-09-14 17:23:02
合計ジャッジ時間 2,355 ms
ジャッジサーバーID
(参考情報)
judge1 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#define popcount __builtin_popcount
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
const int MAX=100001;
vector<ll> prime;
bool isprime[MAX];
void sieve(){
	for(ll i=3; i<MAX; i+=2){
		isprime[i]=1;
	}
	isprime[2]=1;
	prime.push_back(2);
	for(ll i=3; i<MAX; i++){
		if(isprime[i]){
			prime.push_back(i);
			for(ll j=2*i; j<MAX; j+=i){
				isprime[j]=0;
			}
		}
	}
	return;
}
int main()
{
    int n, k; cin>>n>>k;
    map<int, int> mp[100001];
    sieve();
    for(int i=0; i<prime.size(); i++){
        int p=prime[i];
        if(p>n) break;
        for(int i=p; i<=n; i+=p){
            int x=i, e=0;
            while(x%p==0){
                x/=p; e++;
            }
            mp[i][p]=e;
        }
    }
    vector<P> ans;
    for(int i=2; i<n; i++){
        int ct=0;
        for(auto pe:mp[n]){
            int p=pe.first, e=pe.second;
            auto itr=mp[i].find(p);
            if(itr==mp[i].end()) continue;
            else{
                ct+=min(itr->second, e);
            }
        }
        if(ct<k) continue;
        int d=1;
        for(auto pe:mp[i]){
            d*=(pe.second+1);
        }
        ans.push_back(P(-d, i));
    }
    sort(ans.begin(), ans.end());
    cout<<ans[0].second<<endl;
    return 0;
}
0