Cod sursa(job #3156639)

Utilizator Alex_BerbescuBerbescu Alexandru Alex_Berbescu Data 11 octombrie 2023 23:11:50
Problema GFact Scor 0
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.3 kb
#pragma GCC optimize("O3")
#pragma GCC optimize("fast-math")
#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
#define INF 10000000
#define ll long long
using namespace std;
ll p, q, st, dr, mid, sol;
vector<pair<ll, ll> >primes;
ll legendre(int val, int baza)
{
    ll ans = 0;
    while(val > 1)
    {
        ans += 1ll * val / baza;
        val /= baza;
    }
    return ans;
}
ifstream fin("gfact.in");
ofstream fout("gfact.out");
int32_t main(int argc, char * argv[])
{
    fin >> p >> q;
    int d = 2, fm = 0;
    do
    {
        fm = 0;
        while(p % d == 0)
        {
            fm++;
            p /= d;
        }
        if(fm > 0)
        {
            primes.push_back({d, fm * q});
        }
        d = d + 1 + (d & 1);
        if(d * d > p)
        {
            d = p;
        }
    }
    while(p > 1);
    st = 1, dr = INF, mid = 0;
    while(st <= dr)
    {
        mid = (st + dr) / 2;
        bool verif = true;
        for(int j = 0; j < primes.size(); ++j)
        {
            if(legendre(mid, primes[j].first) < primes[j].second)
            {
                verif = false;
            }
        }
        if(verif)
        {
            sol = mid;
            dr = mid - 1;
        }
        else
        {
            st = mid + 1;
        }
    }
    fout << sol;
    return 0;
}