Pagini recente » Cod sursa (job #2806369) | Cod sursa (job #469388) | Cod sursa (job #321408) | Cod sursa (job #550370) | Cod sursa (job #443222)
Cod sursa(job #443222)
#include <cstdio>
using namespace std;
FILE* fin=fopen("gfact.in","r");
FILE* fout=fopen("gfact.out","w");
typedef unsigned long long uint64;
int p,q,decomp[30][2],nf=0;
#define max(a,b) (((a)>(b))?(a):(b))
inline int num(uint64 a,uint64 k){
int nr=0;
for(int i=a;i<=k;i*=a){
nr+=k/i;
}
return nr;
}
inline uint64 search(uint64 a,uint64 b){
uint64 beg=0LL,end=10009676333LL,last=0LL,mdl;
while(beg<=end){
mdl=beg+((end-beg)>>1LL);
if(num(a,mdl)>=b){
last=mdl,end=mdl-1LL;
}else{
beg=mdl+1LL;
}
}
return last;
}
int main(){
fscanf(fin,"%d %d\n",&p,&q);
for(int i=2,k;i*i<=p;++i){
if(p%i==0){
k=0;
while(p%i==0){
k++,p/=i;
}
decomp[nf][0]=i;
decomp[nf][1]=k*q;
nf++;
}
}
if(p>1){
decomp[nf][0]=p;
decomp[nf][1]=q;
nf++;
}
uint64 maxv=0,c=0;
for(int i=0;i<nf;i++){
c=search(decomp[i][0],decomp[i][1]);
maxv=max(c,maxv);
}
fprintf(fout,"%llu ",maxv);
fclose(fin);
fclose(fout);
return 0;
}