noip 2009 细胞分裂

/*数论题 考察唯一分解定理 当然用到一些技巧*/
#include<iostream>
#include<cstdio>
#include<cstring>
#define maxn 30010
using namespace std;
int n,m1,m2,prime[maxn],num,mi[maxn],S,s[maxn],ans=0x7fffffff,cnt;
bool f[maxn];
void Get_prime(int x)
{
    for(int i=2;i<=x;i++)
      {
          if(f[i]==0)prime[++num]=i;
          for(int j=1;j<=num;j++)
            {
                if(i*prime[j]>x)break;
                f[i*prime[j]]=1;
                if(i%prime[j]==0)break;
          }
      }
}
void Get_mi(int x,int y)
{
    if(x<2)return;//m1==1 的时候会re 
    for(int i=1;prime[i]<=x;i++)
      {
          while(x%prime[i]==0)mi[i]++,x/=prime[i];
          if(x==1)break;
      }
    for(int i=1;i<=m2;i++)
      mi[i]*=y;
}
int cla(int a,int b)
{
    if(a<b)return 0;
    int t=a/b;
    if(t*b==a)return t;
    else return t+1;
}
int main()
{
    scanf("%d%d%d",&n,&m1,&m2);
    Get_prime(m1);
    Get_mi(m1,m2);
    for(int k=1;k<=n;k++)
      {
          scanf("%d",&S);int falg=0,A=0;
          for(int i=1;prime[i]<=m1;i++)//相当于只分解m1 那这m1的因子来分解Si 
            {
                cnt=0;if(!prime[i])break;
                while(S%prime[i]==0)cnt++,S/=prime[i];
                if(mi[i]&&!cnt)
                  {
                      falg=1;break;
              }
            if(mi[i]==0&&cnt==0)continue;
            A=max(A,cla(mi[i],cnt));
          }
        if(falg==0)ans=min(ans,A);
      }
    if(ans>=0x7fffff)printf("-1\n");
    else printf("%d\n",ans);
    return 0;
}

 

posted @ 2016-08-06 22:23  一入OI深似海  阅读(343)  评论(0编辑  收藏  举报