Pagini recente » Cod sursa (job #1332326) | Cod sursa (job #3277548) | Cod sursa (job #1473221) | Cod sursa (job #1361357) | Cod sursa (job #354152)
Cod sursa(job #354152)
#include<cstdio>
const long long MAXIM=(long long)1<<60;
int p,q,n,a[1001],b[1001];
void dfp(int m)
{
for(int i=2;i*i<=m;i++)
{
if(m%i==0)
{
++n;
a[n]=i;
}
while(m%i==0)
{
++b[n];
m/=i;
}
}
if(m!=1)
{
++n;
a[n]=m;
b[n]=1;
}
}
bool ver(long long x, int y, int z)
{
int nrfac=0;
while(x)
{
nrfac+=x/y;
x/=y;
}
if(nrfac>=z)
return true;
return false;
}
bool ok(long long x)
{
for(int i=1;i<=n;i++)
if(ver(x,a[i],b[i]*q)==false)
return false;
return true;
}
int main()
{
freopen("gfact.in","r",stdin);
freopen("gfact.out","w",stdout);
long long i,pas;
scanf("%d%d",&p,&q);
dfp(p);
pas=MAXIM;
for(i=1;pas;pas>>=1)
if(!ok(i+pas))
i+=pas;
printf("%lld ",i+1);
return 0;
}