結果

問題 No.375 立方体のN等分 (1)
ユーザー kongarishisyamo
提出日時 2016-06-07 02:15:59
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 2,573 bytes
コンパイル時間 573 ms
コンパイル使用メモリ 60,084 KB
実行使用メモリ 6,816 KB
最終ジャッジ日時 2024-10-08 17:03:55
合計ジャッジ時間 3,674 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 31 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<cmath>

using namespace std;

#define MAX 316228

long long N;
bool arr[MAX+1];
long long insuu[MAX];
long long insuun=0;
long long pn[MAX];
long long pc=0;
long long ansmin;

long long min(long long a,long long b){
	return (a<b)?a:b;
}

long long max(long long a,long long b){
	return (a>b)?a:b;
}

void search1(long long a,long long b,long long c,int sn,int pl){
    if(a+b+c>=ansmin) return;
    if(sn==-1||sn==insuun){
        ansmin=min(ansmin,a+b+c);
        return;
    }
    if(max(a,max(b,c))==a){
        search1(a,b*insuu[sn],c,sn+pl,pl);
        search1(a,b,c*insuu[sn],sn+pl,pl);
    }
    else if(max(a,max(b,c))==b){
        search1(a*insuu[sn],b,c,sn+pl,pl);
        search1(a,b,c*insuu[sn],sn+pl,pl);
    }
    else{
        search1(a*insuu[sn],b,c,sn+pl,pl);
        search1(a,b*insuu[sn],c,sn+pl,pl);
    }
}

void search2(long long a,long long b,long long c,int sn,int pl){
    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);
    if(a+b+c>=ansmin) return;
    search2(a,b*insuu[sn],c,sn+pl,pl);
    if(a+b+c>=ansmin) return;
    search2(a,b,c*insuu[sn],sn+pl,pl);
}

int main(){

	cin>>N;

	long long ansmax=N-1;

	for(long long i=0;i<=MAX;i++){
		arr[i]=true;
	}

	arr[0]=false;
	arr[1]=false;

	for(long long i=2;i<=MAX;i++){
		if(arr[i]){
			pn[pc]=i;
			pc++;
			for(long long j=i*2;j<=MAX;j+=i){
				arr[j]=false;
			}
		}
	}

	for(long long i=0;i<pc;i++) insuu[i]=0;

	for(long long i=0;i<pc;i++){
		while(N%pn[i]==0){
			N/=pn[i];
			insuu[insuun]=pn[i];
			insuun++;
		}
		if(N==1) break;
	}

	if(N!=1) insuu[insuun]=N,insuun++;
	
	long long ni[3]={1,1,1};
    long long mini=0,insuui=pc-1;
	long long insuun2=insuun-1;

	ansmin=ansmax+3;

    while(insuun2!=-1){
        mini=0;
        for(int i=1;i<3;i++){
            if(ni[i]<ni[mini]) mini=i;
        }
        ni[mini]*=insuu[insuun2];
        insuun2--;
		if(insuun2<=25){
			search1(ni[0],ni[1],ni[2],insuun2,-1);
		}
		if(insuun2<=16){
			search2(ni[0],ni[1],ni[2],insuun2,-1);
		}
    }
	
	ansmin=min(ansmin,ni[0]+ni[1]+ni[2]);

	ni[0]=1,ni[1]=1,ni[2]=1;
    mini=0,insuui=pc-1;
	insuun2=0;

    while(insuun2!=insuun){
        mini=0;
        for(int i=1;i<3;i++){
            if(ni[i]<ni[mini]) mini=i;
        }
        ni[mini]*=insuu[insuun2];
        insuun2++;
		if(insuun-insuun2<=25){
			search1(ni[0],ni[1],ni[2],insuun2,1);
		}
		if(insuun-insuun2<=16){
			search2(ni[0],ni[1],ni[2],insuun2,1);
		}
    }
	cout<<ansmin-3<<" "<<ansmax<<endl;
}
0