結果
問題 | No.376 立方体のN等分 (2) |
ユーザー |
|
提出日時 | 2016-07-04 12:03:05 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2,826 ms / 5,000 ms |
コード長 | 2,106 bytes |
コンパイル時間 | 1,236 ms |
コンパイル使用メモリ | 103,096 KB |
実行使用メモリ | 8,848 KB |
最終ジャッジ日時 | 2024-11-06 23:31:24 |
合計ジャッジ時間 | 17,449 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 38 |
ソースコード
#include <iostream> #include <vector> #include <algorithm> #include <cmath> #include <map> #include <functional> #define repeat(i,n) for (int i = 0; (i) < (n); ++(i)) typedef long long ll; template <class T> bool setmin(T & l, T const & r) { if (not (r < l)) return false; l = r; return true; } using namespace std; vector<int> sieve_of_eratosthenes(int n) { // enumerate primes in [2,n] with O(n log log n) vector<bool> is_prime(n+1, true); is_prime[0] = is_prime[1] = false; for (int i = 2; i*i <= n; ++i) if (is_prime[i]) for (int k = i+i; k <= n; k += i) is_prime[k] = false; vector<int> primes; for (int i = 2; i <= n; ++i) if (is_prime[i]) primes.push_back(i); return primes; } map<ll,int> factors(ll n, vector<int> const & primes) { map<ll,int> result; for (int p : primes) { if (n < p *(ll) p) break; while (n % p == 0) { result[p] += 1; n /= p; } } if (n != 1) result[n] += 1; return result; } const ll inf = 1e18+9; template <typename F> ll go(ll n, map<ll,int> & ps, map<ll,int>::iterator it, ll acc, F cont) { if (it == ps.end()) { return cont(n, ps, acc); } else { ll p; int cnt; tie(p, cnt) = *it; ++ it; ll ans = inf; repeat (i,cnt+1) { setmin(ans, go(n, ps, it, acc, cont)); ps[p] -= 1; acc *= p; n /= p; } ps[p] += cnt+1; return ans; } } template <typename F> ll go(ll n, map<ll,int> ps, F cont) { return go(n, ps, ps.begin(), 1, cont); } int main() { ll n; cin >> n; vector<int> primes = sieve_of_eratosthenes(ceil(sqrt(n))); map<ll,ll> memo; auto ps = factors(n, primes); ll ans = go(n, ps, [&](ll n, map<ll,int> ps, ll a) { if (not memo.count(n)) { memo[n] = go(n, ps, [&](ll c, map<ll,int> ps, ll b) { return (b - 1) + (c - 1); }); } return (a - 1) + memo[n]; }); cout << ans << ' ' << n-1 << endl; return 0; }