OJ P1200 【noip2009 普及】细胞问题

#include<iostream>
#include<iomanip>
#include<cmath>
#include<cstdio>
#include<queue>
#include<algorithm>
#include<string>
#include<cstring>
using namespace std;
int m1[10000],m2[10000],a1[10000],a2[10000];
int n;
int mm1,mm2;
int ans=-1;
int k=0,kk=0;
int main()
{
 
 cin>>n;
 cin>>mm1>>mm2;
 
 int i=2;
 while (mm1>1)
 {
  if(mm1%i==0) 
  {
   m1[++k]=i;
   m2[k]=0;
   while(mm1%i==0) {m2[k]++;mm1/=i;}
  }
  else i++;
 }
 for(int i=1;i<=k;i++)
  m2[i]*=mm2;
 //将m1^m2进行质因数分解
 
 
 for(int i=1;i<=n;i++)//读一个细胞分裂速度就运算一个
 {
  bool q=1;//判断是不是所有m1^m2的质因数都在aaa中
  long long  aaa;
  cin>>aaa;
  memset(a1,0,sizeof(a1));
  memset(a2,0,sizeof(a2));
  kk=0;
  for(int i=1;i<=k;)
  {
   
   if(aaa%m1[i]==0) 
   {
    a1[++kk]=m1[i];
    a2[kk]=0;
    while(aaa%m1[i]==0) {a2[kk]++;aaa/=m1[i];}
   }
   else {q=0;break;}
   i++;
  }//用m1^m2的所有质因数检验aaa
  if(!q) continue;
  
  int aans=0;
  for(int i=1;i<=k;i++)
  {
   if(m2[i]%a2[i]!=0 && (aans<(m2[i]/a2[i]+1) || aans==0)) aans=m2[i]/a2[i]+1;
   if(m2[i]%a2[i]==0 && (aans<(m2[i]/a2[i]) || aans==0)) aans=m2[i]/a2[i];
  }
  if(ans==-1 || ans>aans)
   ans=aans;//记录单个最大但总体最小值
 }
 cout<<ans<<endl;
return 0;
}