Cod sursa(job #2635095)

Utilizator Sho10Andrei Alexandru Sho10 Data 13 iulie 2020 12:26:59
Problema Frac Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.21 kb
#include <bits/stdc++.h> //Andrei Alexandru a.k.a Sho10
#define ll long long int
#pragma GCC optimize("O3")
#pragma GCC optimize("Ofast")
#define all(a) (a).begin(), (a).end()
#define sz size
#define f first
#define s second
#define pb push_back
#define er erase
#define in insert
#define mp make_pair
#define pi pair
#define rc(s) return cout<<s,0
#define endl '\n'
#define mod 1000000007
#define PI 3.14159265359
#define MAXN 100005
#define INF 1000000005
#define LINF 1000000000000000005
#define CODE_START  ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);
using namespace std;
ll n,k,nr,rs;
vector<ll>v;
void ans(ll pos,ll sign,ll y){
if(pos==v.size()){
    rs+=sign*(nr/y);
}else {
ans(pos+1,sign,y);
ans(pos+1,-sign,y*v[pos]);
}
}
ll calc(ll val){
    rs=0;
    nr=val;
    ans(0,1,1);
    return rs;
}
int32_t main(){
CODE_START;
ifstream cin("frac.in");
ofstream cout("frac.out");
cin>>n>>k;
for(ll i=2;i*i<=n;i++)
{
    ll s1=0;
    while(n%i==0){
        s1=1;
        n/=i;
    }
    if(s1==1){
        v.pb(i);
    }
}
if(n!=1){
    v.pb(n);
}
ll r=0,x=(1ll<<60);
while(x){
    if(calc(r+x)<k){
        r+=x;
    }
    x=x/2;
}
r++;
cout<<r<<endl;
}