Cod sursa(job #468087)

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

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

int putere(int a,int b)//puterea la care apare b in a!
{
	int rez=0;
	while(a>=b)
		rez+=(a/=b);
	return rez;
}
void descompun(int 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(int n)
{
	for(int i=1;i<=k;i++)
		if(putere(n,v[i].div)<q*v[i].put)
			return false;
	return true;
}
int caut()
{
	int i,pas=1<<25;
	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;
}