Pagini recente » Cod sursa (job #790784) | Cod sursa (job #188959) | Cod sursa (job #2134464) | Cod sursa (job #1693846) | Cod sursa (job #2149888)
#include<cstdio>
using namespace std;
int cnt=0,ciu[1000007],prim[100007];
void ciur()
{
for(int i=2;i<=1000007;i++)
if(ciu[i]==0)
{
cnt++;
prim[cnt]=i;
for(int j=i;j<=1000000;j+=i) ciu[j]=1;
}
}
long long pinex(long long a,long long b)
{
int cnt2=0,val[100007];
long long bb=b,aa=a;
for(int i=1;i<=cnt;i++)
if(bb%prim[i]==0)
{
cnt2++;
val[cnt2]=prim[i];
while(bb%prim[i]==0) bb/=prim[i];
}
if(bb!=1) {cnt2++;val[cnt2]=bb;}
//printf("%d\n",cnt2);
for(long i=1;i<=((1<<cnt2)-1);i++)
{
long long x=i,p=1;
int cnt3=0;
for(int j=1;j<=cnt2;j++)
{
if(x&1) {cnt3++;p*=val[j];}
x>>=1;
}
if(cnt3&1) aa-=(a/p);
else aa+=(a/p);
}
return aa;
}
long long cautbin(long long p, long long n)
{
int step=(1<<30), start = 1;
for(;step;step>>=1)
{
int index=step+start;
if(pinex(index,n)<p)
start=index;
}
return start;
}
int main()
{
freopen("frac.in","r",stdin);
freopen("frac.out","w",stdout);
ciur();
long long n,p;
scanf("%lld %lld",&n,&p);
long long x=cautbin(p,n);
printf("%lld",x+1);
}