Cod sursa(job #126465)

Utilizator ciprianfFarcasanu Alexandru Ciprian ciprianf Data 22 ianuarie 2008 11:50:17
Problema GFact Scor 75
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.63 kb
#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;
}