結果

問題 No.375 立方体のN等分 (1)
コンテスト
ユーザー h_noson
提出日時 2016-06-04 23:23:06
言語 C++11(old_compat)
(gcc 12.4.0 + boost 1.89.0)
コンパイル:
g++-12 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -include bits/stdc++.h -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 1,496 ms / 5,000 ms
コード長 1,476 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,516 ms
コンパイル使用メモリ 177,052 KB
実行使用メモリ 134,528 KB
最終ジャッジ日時 2026-03-08 16:06:33
合計ジャッジ時間 5,932 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <vector>
#include <map>
#include <algorithm>
using namespace std;

#define REP(i,s,e) for (i = s; i <= e; i++)
#define rep(i,n) REP (i,0,(int)(n)-1)
#define RREP(i,s,e) for (i = s; i >= e; i--)
#define rrep(i,n) RREP (i,(int)(n)-1,0)
#define INF (int)1e8
#define MOD (int)(1e9+7)

typedef long long ll;

ll a[3];
bool prime[1000000];
map<tuple<ll,ll,ll>,ll> mp;

ll dfs(ll n) {
    tuple<ll,ll,ll> tp = make_tuple(a[0],a[1],a[2]);
    if (mp.count(tp))
        return mp[tp];
    if (n == 1) 
        return a[0] + a[1] + a[2] - 3;
    ll ret = 0;
    int i;
    for (i = sqrt(n); i >= 2; i--) {
        if (prime[i] && n % i == 0) {
            a[0] *= i;
            ret = dfs(n/i);
            a[0] /= i;
            a[1] *= i;
            ret = min(ret,dfs(n/i));
            a[1] /= i;
            a[2] *= i;
            ret = min(ret,dfs(n/i));
            a[2] /= i;
            return mp[tp] = ret;
        }
    }
    a[0] *= n;
    ret = dfs(1);
    a[0] /= n;
    a[1] *= n;
    ret = min(ret,dfs(1));
    a[1] /= n;
    a[2] *= n;
    ret = min(ret,dfs(1));
    a[2] /= n;
    return mp[tp] = ret;
}

int main(void) {
    ll i, j, n;
    cin >> n;

    rep (i,1000000) prime[i] = true;
    for (i = 2; i * i <= 1000000; i++) {
        if (prime[i]) {
            for (j = i*2; j <= 1000000; j+=i)
                prime[j] = false;
        }
    }
    a[0] = a[1] = a[2] = 1;
    cout << dfs(n) << " " << n - 1 << endl;
    return 0;
}
0