Cod sursa(job #468089)

Utilizator unknownliviuMaria Liviu Valentin unknownliviu Data 2 iulie 2010 11:29:36
Problema GFact Scor 95
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.9 kb
#include<fstream>
using namespace std;
ifstream in("gfact.in");
ofstream out("gfact.out");
long long p,q,k;
struct diviz{long long div,put;};
diviz v[10];

void read()
{
	in>>p>>q;
}

long long putere(long long a,long long b)//puterea la care apare b in a!
{
	long long rez=0;
	while(a>=b)
		rez+=(a/=b);
	return rez;
}
void descompun(long long p)
{
	int m;
	for(int i=2;i*i<=p;i++)
	{
		m=0;
		while(p%i==0)
		{
			m++;
			p/=i;
		}
		if(m)
		{
			k++;
			v[k].div=i;
			v[k].put=m;
		}
	}
	if(p!=1)
	{
		++k;
		v[k].div=p;
		v[k].put=1;
	}
}
bool sol(long long n)
{
	for(int i=1;i<=k;i++)
		if(putere(n,v[i].div)<q*v[i].put)
			return false;
	return true;
}
long long caut()
{
	long long i,pas=(long long)1<<40;
	for(i=0;pas!=0;pas>>=1)
		if(sol(i+pas)==false)
			i+=pas;
	return i+1;
}
int main()
{
	read();
	descompun(p);
	out<<caut();
	return 0;
}