Pagini recente » Cod sursa (job #618214) | Cod sursa (job #2671039) | Cod sursa (job #3176667) | Cod sursa (job #2248143) | Cod sursa (job #353475)
Cod sursa(job #353475)
#include <cstdio>
int p,q,f[1000],c[1000],nr;
void desc(int y) // descompunere in fact primi
{
nr=0;
for (int i=2;i*i<=y;i++)
if (y%i==0)
{
nr++;
f[nr]=i;
c[i]=0;
while (y%i==0)
{
y=y/i;
c[i]++;
}
}
if (y!=1)
{
nr++;
f[nr]=y;
c[nr]=1;
}
}
bool nrz(int n,int a,int b)
{
int b2=0;
while (n)
b2+=(n/=a);
return b2>=b;
}
bool div(int n)
{
int i;
for (i=1;i<=nr;i++)
if (!nrz(n,f[i],c[i]))
return 0;
return 1;
}
int cautbin()
{
int i,pas=1<<30;
for (i=0;pas;pas>>=1) // pas=pas/2;
if (!div(i+pas))
i=i+pas;
return i+1;
}
int main()
{
int i;
freopen("gfact.in","r",stdin);
freopen("gfact.out","w",stdout);
scanf("%d%d",&p,&q);
desc(p);
for (i=1;i<=nr;i++) c[i]=c[i]*q;
int x=cautbin();
printf("%d",x);
return 0;
}