码迷,mamicode.com
首页 > 其他好文 > 详细

noip 2009 细胞分裂

时间:2016-08-06 23:24:15      阅读:369      评论:0      收藏:0      [点我收藏+]

标签:

/*数论题 考察唯一分解定理 当然用到一些技巧*/
#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;
}

 

noip 2009 细胞分裂

标签:

原文地址:http://www.cnblogs.com/yanlifneg/p/5745037.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!