Nu exista pagina, dar poti sa o creezi ...
Cod sursa(job #2353345)
Utilizator | Data | 24 februarie 2019 11:11:55 | |
---|---|---|---|
Problema | GFact | Scor | 65 |
Compilator | cpp-64 | Status | done |
Runda | Arhiva de probleme | Marime | 1.1 kb |
#include <fstream>
#include <climits>
#include <cmath>
using namespace std;
ifstream fin("gfact.in");
ofstream fout("gfact.out");
int Div[50005], p, q;
bool test(unsigned long long val)
{
int r = sqrt(p) + 1;
for(int i = 2; i <= r; ++i)
{
if(Div[i])
{
unsigned long long rez = 0;
int p = i;
while(val / p != 0){
rez += val / p;
p *= i;
}
if(rez < Div[i])
return 0;
}
}
return 1;
}
int main()
{
fin >> p >> q;
int d = 2;
int cp = p;
while(p > 1)
{
int e = 0;
while(p % d == 0)
{
++e;
p /= d;
}
Div[d] = e * q;
++d;
}
p = cp;
unsigned long long rez = 0;
unsigned long long st = 1;
unsigned long long dr = ULLONG_MAX - 3;
while(st <= dr)
{
unsigned long long mij = (st + dr) / 2;
if(test(mij))
{
rez = mij;
dr = mij - 1;
}
else st = mij + 1;
}
fout << rez << '\n';
return 0;
}