Pagini recente » Cod sursa (job #1184949) | Cod sursa (job #2950178) | Cod sursa (job #1740074) | Cod sursa (job #409238) | Cod sursa (job #126465)
Cod sursa(job #126465)
#include <stdio.h>
int main()
{ int p,q,st,dr,m,nr,cp,d,c=0,factmax=0,qmax=0,fact=1,dmax;
FILE*f=fopen("gfact.in","r");
FILE*g=fopen("gfact.out","w");
fscanf(f,"%d%d",&p,&q);
for (d=2;d*d<=p;d++)
{
for(c=0,fact=1;p%d==0;c++,p/=d,fact*=d);
if(fact>factmax) { factmax=fact; qmax=c; dmax=d;}
}
if(p>1 && factmax<p) { qmax=1; dmax=p;}
q=q*qmax;
st=1;
dr=2000000000;
d=dmax;
while(st<=dr)
{
nr=0;
cp=d;
m=(st+dr)/2;
while(m/cp)
{
nr=nr+m/cp;
cp=cp*d;
}
if (nr==q)
{
while(m%d) m--;
fprintf(g,"%d",m);
return 0;
}
if (nr>q) dr=m-1;
if (nr<q) st=m+1;
}
fprintf(g,"%d",st);
return 0;
}