Pagini recente » Cod sursa (job #740894) | Cod sursa (job #2026673) | Cod sursa (job #406736) | Cod sursa (job #2750749) | Cod sursa (job #2144263)
#include <fstream>
using namespace std;
ifstream in("gfact.in");
ofstream out("gfact.out");
const int N = 120001;
long long p, q, b, prim[N];
long long putere(int poz, long long x) {
long long nr = 0;
while (x != 0) {
x /= poz;
nr += x;
}
return nr;
}
void cautbin(int poz) {
long long r = 0, pas = 1LL << 52;
while (pas != 0) {
if (r + pas <= q * poz * prim[poz] && putere(poz, r + pas) < q * prim[poz]) {
r += pas;
}
pas >>= 1;
}
b = max(b, r + 1);
}
void precalc_prime(int nr) {
bool isprim = true;
long long cp = p;
for (int d = 2; d * d <= p; d++) {
while (cp % d == 0) {
prim[d]++;
cp /= d;
isprim = false;
}
if (prim[d] != 0) cautbin(d);
}
if (isprim) {
prim[p]++;
cautbin(p);
}
}
int main()
{
in >> p >> q;
precalc_prime(p);
out << b;
out.close();
}