Cod sursa(job #2067535)

Utilizator ivan.tudorIvan Tudor ivan.tudor Data 16 noiembrie 2017 16:24:31
Problema GFact Scor 90
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.99 kb
#include <iostream>
#include<cstdio>
#include <utility>
using namespace std;
int p,q,fina;
pair <int,int> divizori[100];
int legendre(long long n,int d){
  int rasp=0;
  while(d<=n){
    rasp+=(n/=d);
  }
  return rasp;
}
bool ver(long long n){
  long long i;
  for(i=1;i<=fina;i++){
    if(legendre(n,divizori[i].first)<divizori[i].second)
      return false;
  }
  return true;
}

long long cautbin(){
  int p2=1<<30;
  long long pas=0;
  while(p2>0){
    if(ver(pas+p2)==false)
      pas+=1LL*p2;
    p2/=2;
  }
  return pas+1;
}
int main()
{
    FILE*fin,*fout;
    fin=fopen("gfact.in","r");
    fout=fopen("gfact.out","w");
    fscanf(fin,"%d%d",&p,&q);
    long long d=2,i=0;
    while(d*d<=p){
      if(p%d==0){
        divizori[++i].first=d;
      }
      while(p%d==0){
        p/=d;
        divizori[i].second++;
      }
      divizori[i].second*=q;
      d++;
    }
    if(p>1){
      divizori[++i].first=p;
      divizori[i].second=q;
    }
    fina=i;
    fprintf(fout,"%lld",cautbin());
    return 0;
}