Cod sursa(job #2285614)

Utilizator victor1306Victor Mihaila victor1306 Data 18 noiembrie 2018 20:21:25
Problema GFact Scor 0
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.15 kb
#include <fstream>
#include<cmath>
using namespace std;
ifstream cin("gfact.in");
ofstream cout("gfact.out");
const int L=31;
int divizor[100],putere[100],nd,p,q;
void descompunere(int n)
{
    int dv = 2;
    while (dv * dv <= n)
    {
        if (n % dv == 0)
        {
            divizor[nd] = dv;
            while (n % dv == 0)
            {
                putere[nd]++;
                n /= dv;
            }
            nd++;
        }
        dv++;
    }
    if (n > 1)
    {
        divizor[nd] = n;
        putere[nd++] = 1;
    }
}

long long nrp(long long n, int p)
{

    long long nr=0;
    while (n>=p)
    {
        nr+=n/p;
        n/=p;
    }
    return nr;
}
bool  df(long long n)
{
    for (int i=0; i<nd; i++)
    {
        if (nrp(n,divizor[i])<putere[i] * q)
        {
            return 0;
        }
    }
    return 1;
}
long long caut()
{
    long long r = 0, pas = 1<< L;
    while (pas != 0)
    {
        if (!df(r + pas))
        {
            r += pas;
        }
        pas /= 2;
    }
    return r+1;
}
int main()
{
    cin>>p>>q;
    descompunere(p);
    cout<<caut();
}