Cod sursa(job #373024)

Utilizator Teodor94Teodor Plop Teodor94 Data 12 decembrie 2009 14:20:36
Problema GFact Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.93 kb
#include<cstdio>

long long nr,prim[5000],put[5000];

void desc(long long x)
{
	if (x%2==0)
	{
		nr++;
		prim[nr]=2;
		while (x%2==0)
		{
			put[nr]++;
			x/=2;
		}
	}
	for (long long i=3;i*i<=x;i+=2)
		if (x%i==0)
		{
			nr++;
			prim[nr]=i;
			while (x%i==0)
			{
				put[nr]++;
				x/=i;
			}
		}
	if (x!=1)
	{
		nr++;
		prim[nr]=x;
		put[nr]=1;
	}
}

bool ok(long long x)
{
	bool ok=true;
	for (int i=1;i<=nr;i++)
	{
		long long s=0,x1=x;
		while (x1)
			s += (x1/=prim[i]);
		if (s<put[i])
		{
			ok=false;
			break;
		}
	}
	return ok;
}

long long cautbin()
{
	long long i,pas=(long long)1<<60;
	for (i=0;pas;pas>>=1)
		if (!ok(i+pas))
			i+=pas;
	return i+1;
}

int main()
{
	freopen("gfact.in","r",stdin);
	freopen("gfact.out","w",stdout);
	long long p=0,q=0;
	scanf("%lld%lld",&p,&q);
 	desc(p);
	for (int i=1;i<=nr;i++)
		put[i]*=q;
	printf("%lld",cautbin());
	return 0;
}