题目描述
HanksHanksHanks 博士是 BTBTBT (Bio−TechBio-TechBio−Tech,生物技术) 领域的知名专家。现在,他正在为一个细胞实验做准备工作:培养细胞样本。
HanksHanksHanks 博士手里现在有 NN N种细胞,编号从 1−N1-N1−N,一个第 ii i种细胞经过 111 秒钟可以分裂为SiS_iSi个同种细胞(SiS_iSi为正整数)。现在他需要选取某种细胞的一个放进培养皿,让其自由分裂,进行培养。一段时间以后,再把培养皿中的所有细胞平均分入M M M个试管,形成M M M份样本,用于实验。HanksHanksHanks 博士的试管数M M M很大,普通的计算机的基本数据类型无法存储这样大的MM M值,但万幸的是,MMM 总可以表示为m1 m_1m1的m2 m_2m2次方,即M=m1m2M = m_1^{m_2}M=m1m2,其中 m1,m2m_1,m_2m1,m2均为基本数据类型可以存储的正整数。
注意,整个实验过程中不允许分割单个细胞,比如某个时刻若培养皿中有 44 4个细胞,
HanksHanks Hanks博士可以把它们分入 222 个试管,每试管内2 22 个,然后开始实验。但如果培养皿中有5 55个细胞,博士就无法将它们均分入2 22 个试管。此时,博士就只能等待一段时间,让细胞们继续分裂,使得其个数可以均分,或是干脆改换另一种细胞培养。
为了能让实验尽早开始,HanksHanks Hanks博士在选定一种细胞开始培养后,总是在得到的细胞“刚好可以平均分入 MM M个试管”时停止细胞培养并开始实验。现在博士希望知道,选择哪种细胞培养,可以使得实验的开始时间最早。
输入格式
第一行,有一个正整数 NNN,代表细胞种数。
第二行,有两个正整数 m1,m2m_1,m_2m1,m2,以一个空格隔开,即表示试管的总数 M=m1m2M = m_1^{m_2}M=m1m2.
第三行有 N 个正整数,第 i 个数 Si表示第 i 种细胞经过 1 秒钟可以分裂成同种细胞的个数。
输出格式
一个整数,表示从开始培养细胞到实验能够开始所经过的最少时间(单位为秒)。
如果无论HanksHanksHanks博士选择哪种细胞都不能满足要求,则输出整数−1-1−1。
输入输出样例
输入 #1
1 2 1 3
输出 #1
-1
输入 #2
2 24 1 30 12
输出 #2
2
说明/提示
【输入输出说明】
经过 1秒钟,细胞分裂成3 个,经过2秒钟,细胞分裂成9个,……,可以看出无论怎么分裂,细胞的个数都是奇数,因此永远不能分入 2个试管。
【输入输出样例22 2说明】
第 1种细胞最早在3秒后才能均分入24个试管,而第2种最早在2 秒后就可以均分(每试管144/(241)=6144/(241)=6144/(241)=6 个)。故实验最早可以在2 秒后开始。
【数据范围】
对于 50%的数据,有m1m2≤30000m_1^{m_2} ≤ 30000m1m2≤30000。
对于所有的数据,有1≤N≤10000,1≤m1≤30000,1≤m2≤10000,1≤Si≤2,000,000,000。
NOIP 2009 普及组 第三题
题解
这是一道数论题,质因数分解。。。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
const int M=1<<30;
using namespace std;
const int N=10005;
int n,m1,m2,s,a[N],sta[N],pop=0;
void init()
{
int f=m1,i;
for(i=2;i<=f;i++)
{
if(m1%i==0&&!a[i])
{
pop++;
sta[pop]=i;
}
while(m1%i==0)
{
a[i]+=m2;
m1/=i;
}
}
}
int work(int x)
{
int sum=0;
int f=s;
while(f%x==0)
{
sum++;
f/=x;
}
return sum;
}
inline int max(int x,int y)
{
return x>y?x:y;
}
int main()
{
scanf("%d",&n);
scanf("%d%d",&m1,&m2);
if(m1==1)
{
printf("0\n");
return 0;
}
init();
int ans=M,i,j;
for(i=1;i<=n;i++)
{
scanf("%d",&s);
bool flag=true;
int Max=-1;
for(j=1;j<=pop;j++)
{
int k=work(sta[j]);
if(!k)
{
flag=false;
break;
}
if(k>=a[sta[j]])
Max=max(Max,1);
else
Max=max(Max,a[sta[j]]/k+(a[sta[j]]%k==0?0:1));
}
if(flag)
ans=min(ans,Max);
}
if(ans==M)
printf("-1\n");
else
printf("%d\n",ans);
return 0;
}