Pagini recente » Cod sursa (job #1631673) | Borderou de evaluare (job #867939) | Cod sursa (job #2805070) | Cod sursa (job #56832) | Cod sursa (job #796608)
Cod sursa(job #796608)
#include <iostream>
#include <fstream>
#include <vector>
#include <cmath>
#include <cstdio>
#include <set>
using namespace std;
ifstream inp("frac.in");
ofstream outp("frac.out");
int main()
{
unsigned long long n,nn,nsr2,p;
inp >> n >> p;
nn=n;
nsr2 = sqrt(n)+1;
unsigned long long di=2;
vector<unsigned long long> divisors;
while (nn!=1)
{
if (nn%di==0)
divisors.push_back(di);
while (nn%di==0)
nn/=di;
di++;
if (di>nsr2)
{
if (nn!=1)
divisors.push_back(nn);
break;
}
}
vector<unsigned long long> combFact;
vector<unsigned long long> nr1;
for (unsigned long long i=0;i<1<<divisors.size();i++)
{
unsigned long long pr=1;
int nr0=0;
for (int j = 0;j<divisors.size();j++)
{
if (i &(1<<j))
{
pr*=divisors[j];
nr0++;
}
}
if (nr0)
{
combFact.push_back(pr);
nr1.push_back(nr0);
}
}
unsigned long long res=0;
{
unsigned long long ll = 0;
unsigned long long rr=0x2000000000000001;
unsigned long long final=0;
while (ll<rr)
{
unsigned long long mm = (ll+rr)/2;
res=0;
for (int i=0;i<combFact.size();i++)
if (nr1[i]%2==1)
{
res+=mm/combFact[i];
}
else
{
res-=mm/combFact[i];
}
res = mm - res;
if (res == p)
{
final=mm;
rr = mm;
}
else
if (res < p)
{
ll = mm + 1;
}
else
{
rr = mm;
}
if (ll==rr) break;
}
outp << final<<endl;
}
}
/*
10
+5
+3
-1
7
*/
//1/13
//2/13