Pagini recente » Cod sursa (job #1581571) | Cod sursa (job #2407693) | Cod sursa (job #2423083) | Cod sursa (job #1914503) | Cod sursa (job #66932)
Cod sursa(job #66932)
//100 de puncte
#include<stdio.h>
#include<values.h>
int main()
{
long long unsigned p,q,a,sol,j,nr,r,k=2,x,max=0,pr,i,m;
FILE*f=fopen("gfact.in","r");
FILE*g=fopen("gfact.out","w");
fscanf(f,"%llu %llu",&p,&q);
for(k=2;k*k<=p;++k)
{
r=0;
while(p%k==0) {++r; p/=k;}
if(r!=0)
{
r=r*q;
i=1; j=2*r*k;
while(i<=j)
{
m=(i+j)/2;
nr=0; pr=k;
while((m/pr)!=0) { nr+=(m/pr); pr*=k; }
if(nr==r)
{
while(m%k!=0) m--;
break;
}
else if(nr>r) j=m-1;
else i=m+1;
}
if(m>max) max=m;
}
}
if(p>1)
{k=p; r=q;
i=1; j=2*r*k;
while(i<=j)
{
m=(i+j)/2;
nr=0; pr=k;
while((m/pr)!=0) { nr+=(m/pr); pr*=k; }
if(nr==r)
{
while(m%k!=0) m--;
break;
}
else if(nr>r) j=m-1;
else i=m+1;
}
if(m>max) max=m;
}
fprintf(g,"%llu",max);
return 0;
}