Cod sursa(job #2661245)

Utilizator AlexNicuNicu Alexandru AlexNicu Data 21 octombrie 2020 17:17:21
Problema GFact Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.46 kb
#include <fstream>
#include <cstdio>

using namespace std;

ifstream cin ("gfact.in");
ofstream cout ("gfact.out");

#define PMAX 2000000000
#define QMAX 30000
#define MAXDIV 15
long long divprim[MAXDIV], putere_prim[MAXDIV];

long long legendre( long long fact, long long nrprim ) {
    long long ans = 0;
    long long putere = nrprim;
    while ( putere <= fact ) {
      ans += (fact / putere);
      putere *= nrprim;
    }
    return ans;
}

bool verif( long long fact, long long nrdivprimi ) {
  int i;
  for ( i = 0; i < nrdivprimi; i++ )
    if ( legendre( fact, divprim[i] ) < putere_prim[i] )
      return 0;
  return 1;
}

int main()
{
    long long p, q, div, nrdivprimi, putere;
    long long st, dr, mij;
    cin >> p >> q;
    div = 2;
    nrdivprimi = 0;
    while ( div * div <= p ) {
        putere = 0;
        while ( p % div == 0 ) {
          putere++;
          p /= div;
        }
        if ( putere > 0 ) {
          divprim[nrdivprimi] = div;
          putere_prim[nrdivprimi] = putere * q;
          nrdivprimi++;
        }
        div++;
    }
    if ( p > 1 ) {
        divprim[nrdivprimi] = p;
        putere_prim[nrdivprimi] = q;
        nrdivprimi++;
    }
    st = 0;
    dr = (long long)PMAX * QMAX + 1000000;
    while ( dr - st > 1 ) {
      mij = ( st + dr ) / 2;
      if ( verif( mij, nrdivprimi ) == 0 )
        st = mij;
      else
        dr = mij;
    }
    cout << dr;
    return 0;
}