Pagini recente » Cod sursa (job #1958269) | Cod sursa (job #393902) | Cod sursa (job #2429245) | Cod sursa (job #1395113) | Cod sursa (job #2067535)
#include <iostream>
#include<cstdio>
#include <utility>
using namespace std;
int p,q,fina;
pair <int,int> divizori[100];
int legendre(long long n,int d){
int rasp=0;
while(d<=n){
rasp+=(n/=d);
}
return rasp;
}
bool ver(long long n){
long long i;
for(i=1;i<=fina;i++){
if(legendre(n,divizori[i].first)<divizori[i].second)
return false;
}
return true;
}
long long cautbin(){
int p2=1<<30;
long long pas=0;
while(p2>0){
if(ver(pas+p2)==false)
pas+=1LL*p2;
p2/=2;
}
return pas+1;
}
int main()
{
FILE*fin,*fout;
fin=fopen("gfact.in","r");
fout=fopen("gfact.out","w");
fscanf(fin,"%d%d",&p,&q);
long long d=2,i=0;
while(d*d<=p){
if(p%d==0){
divizori[++i].first=d;
}
while(p%d==0){
p/=d;
divizori[i].second++;
}
divizori[i].second*=q;
d++;
}
if(p>1){
divizori[++i].first=p;
divizori[i].second=q;
}
fina=i;
fprintf(fout,"%lld",cautbin());
return 0;
}