#include #include #include #include #include #include using namespace std; using i64 = int64_t; vector factorize(int n) { vector v(n+1); for(int d=2; d*d<=n; ++d) { while(n % d == 0) { ++v[d]; n /= d; } } if(n > 1) { ++v[n]; } return v; } template vector sieve(Int n) { vector is_prime(n+1, true); is_prime[0] = is_prime[1] = false; Int sq = Int(sqrt(n)); for(Int i=2; i<=sq; ++i) { if(is_prime[i]) { for(Int j=i*i; j<=n; j+=i) { is_prime[j] = false; } } } vector res; for(Int i=2; i<=n; ++i) { if(is_prime[i]) { res.push_back(i); } } return res; } void debug(const vector &vec_) { for(size_t i=0, n=vec_.size(); i primes = sieve(x); int n = static_cast(primes.size()); vector fac = factorize(x); // 約数の個数(dcnt), 共通の素因数の個数(ccnt), i番目の素数, k個使う, 値 queue> que; que.emplace(1, 0, 0, 0, 1); i64 res = inf; int max_dcnt = 0; while(!que.empty()) { int dcnt, ccnt, i, k; i64 val; tie(dcnt, ccnt, i, k, val) = que.front(); que.pop(); if(val >= x) { continue; } if(ccnt >= K) { if(max_dcnt < dcnt) { max_dcnt = dcnt; res = val; } else if(max_dcnt == dcnt) { res = min(res, val); } } // i番目をk個からk+1個にする int p = primes[i]; bool have = fac[p] >= k+1; que.emplace(dcnt*(k+2)/(k+1), ccnt+have, i, k+1, val*p); // j番目をはじめて使う for(int j=i+1; j= 1; que.emplace(dcnt*2, ccnt+have, j, 1, val*p); } } return res; } int main(void) { int N, K; scanf("%d%d", &N, &K); i64 res = f(N, K); printf("%ld\n", res); return 0; }