結果
| 問題 |
No.375 立方体のN等分 (1)
|
| コンテスト | |
| ユーザー |
kongarishisyamo
|
| 提出日時 | 2016-06-07 16:35:25 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
AC
|
| 実行時間 | 416 ms / 5,000 ms |
| コード長 | 1,619 bytes |
| コンパイル時間 | 1,327 ms |
| コンパイル使用メモリ | 105,672 KB |
| 実行使用メモリ | 79,896 KB |
| 最終ジャッジ日時 | 2024-10-11 21:05:34 |
| 合計ジャッジ時間 | 7,981 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 32 |
ソースコード
#include<iostream>
#include<cmath>
#include<algorithm>
#include<map>
#include<vector>
using namespace std;
#define MAX 10000001
long long N;
vector<long long> insuu;
long long insuun=0;
long long pc=0;
long long ansmin;
long long cou=0;
map<long long,map<long long,map<long long,map<long long,long long>>>> memo;
void search2(long long a,long long b,long long c,int sn,int pl){
if(memo.find(a)!=memo.end()&&memo[a].find(b)!=memo[a].end()&&memo[a][b].find(c)!=memo[a][b].end()&&memo[a][b][c].find(sn)!=memo[a][b][c].end()){
return;
}
else{
memo[a][b][c][sn]=a+b+c;
memo[a][c][b][sn]=a+b+c;
memo[b][a][c][sn]=a+b+c;
memo[b][c][a][sn]=a+b+c;
memo[c][b][a][sn]=a+b+c;
memo[c][a][b][sn]=a+b+c;
}
if(a+b+c>=ansmin) return;
if(sn==-1||sn==insuun){
ansmin=min(ansmin,a+b+c);
return;
}
search2(a*insuu[sn],b,c,sn+pl,pl);
search2(a,b*insuu[sn],c,sn+pl,pl);
search2(a,b,c*insuu[sn],sn+pl,pl);
}
int main(){
vector<bool> arr;
vector<long long> pn;
cin>>N;
long long ansmax=N-1;
for(long long i=0;i<=MAX;i++){
arr.push_back(true);
}
arr[0]=false;
arr[1]=false;
for(long long i=2;i<=MAX;i++){
if(arr[i]){
pn.push_back(i);
pc++;
for(long long j=i*2;j<=MAX;j+=i){
arr[j]=false;
}
}
}
while(!arr.empty()) arr.pop_back();
for(long long i=0;i<pc;i++){
while(N%pn[i]==0){
N/=pn[i];
insuu.push_back(pn[i]);
insuun++;
}
if(N==1) break;
}
while(!pn.empty()) pn.pop_back();
if(N!=1) insuu.push_back(N),insuun++;
ansmin=ansmax+3;
search2(1,1,1,insuun-1,-1);
ansmin-=3;
cout<<ansmin<<" "<<ansmax<<endl;
}
kongarishisyamo