题目
NOIP2009-普及组T3-细胞分裂 |
---|
难度级别:A; 运行时间限制:1000ms; 运行空间限制:51200KB; 代码长度限制:2000000B |
试题描述 |
H
a
n
k
s
Hanks
Hanks 博士是
B
T
BT
BT (
B
i
o
−
T
e
c
h
Bio-Tech
Bio−Tech,生物技术) 领域的知名专家。现在,他正在为一个细胞实验做准备工作:培养细胞样本。 H a n k s Hanks Hanks 博士手里现在有 N N N 种细胞,编号从 1 − N 1-N 1−N,一个第 i i i 种细胞经过 1 1 1 秒钟可以分裂为 S i S_i Si 个同种细胞( S i S_i Si 为正整数)。现在他需要选取某种细胞的一个放进培养皿,让其自由分裂,进行培养。一段时间以后,再把培养皿中的所有细胞平均分入 M M M 个试管,形成 M M M 份样本,用于实验。 H a n k s Hanks Hanks 博士的试管数 M M M 很大,普通的计算机的基本数据类型无法存储这样大的 M M M 值,但万幸的是, M M M 总可以表示为 m 1 m_1 m1 的 m 2 m_2 m2 次方,即 M = m 1 m 2 M = m_1^{m_2} M=m1m2 ,其中 m 1 , m 2 m_1,m_2 m1,m2 均为基本数据类型可以存储的正整数。 注意,整个实验过程中不允许分割单个细胞,比如某个时刻若培养皿中有 4 4 4 个细胞, H a n k s Hanks Hanks 博士可以把它们分入 2 2 2 个试管,每试管内 2 2 2 个,然后开始实验。但如果培养皿中有 5 5 5 个细胞,博士就无法将它们均分入 2 2 2 个试管。此时,博士就只能等待一段时间,让细胞们继续分裂,使得其个数可以均分,或是干脆改换另一种细胞培养。 为了能让实验尽早开始, H a n k s Hanks Hanks 博士在选定一种细胞开始培养后,总是在得到的细胞“刚好可以平均分入 M M M 个试管”时停止细胞培养并开始实验。现在博士希望知道,选择哪种细胞培养,可以使得实验的开始时间最早。 |
输入 |
第一行,有一个正整数
N
N
N ,代表细胞种数。 第二行,有两个正整数 m 1 , m 2 m_1,m_2 m1,m2 ,以一个空格隔开,即表示试管的总数 M = m 1 m 2 M = m_1^{m_2} M=m1m2 。 第三行有 N N N 个正整数,第 i i i 个数 S i S_i Si表示第 i i i 种细胞经过 1 1 1 秒钟可以分裂成同种细胞的个数。 |
输出 |
一个整数,表示从开始培养细胞到实验能够开始所经过的最少时间(单位为秒)。 如果无论 H a n k s Hanks Hanks 博士选择哪种细胞都不能满足要求,则输出整数 − 1 -1 −1 。 |
输入示例 |
1 2 1 3 |
输出示例 |
-1 |
输入示例2 |
2 24 1 30 12 |
输出示例2 |
2 |
数据说明 |
#1:经过
1
1
1 秒钟,细胞分裂成
3
3
3 个,经过
2
2
2 秒钟,细胞分裂成
9
9
9 个,……,可以看出无论怎么分裂,细胞的个数都是奇数,因此永远不能分入
2
2
2 个试管。 #2:第 1 1 1 种细胞最早在 3 3 3 秒后才能均分入 24 24 24 个试管,而第 2 2 2 种最早在 2 2 2 秒后就可以均分(每试管 144 / 2 4 1 = 6 144/24^1=6 144/241=6 个)。故实验最早可以在 2 2 2 秒后开始。 |
数据范围 |
对于 50%的数据,有
m
1
m
2
≤
30000
m_1^{m_2}≤30000
m1m2≤30000 对于所有的数据,有 1 ≤ N ≤ 10000 , 1 ≤ m 1 ≤ 30000 , 1 ≤ m 2 ≤ 10000 , 1 ≤ S i ≤ 2 , 000 , 000 , 000 1≤N≤10000,1≤m_1≤30000,1≤m_2≤10000,1≤S_i≤2,000,000,000 1≤N≤10000,1≤m1≤30000,1≤m2≤10000,1≤Si≤2,000,000,000 |
题解
思路
一道数学题。(可以欣赏美丽的 LaTeX \LaTeX LATEX 公式了)
先来看题目要做的是什么:枚举 S i n S_i^n Sin ( n n n为这里的秒数),使得保证 m 1 m 2 ∣ S i n m_1^{m_2}|S_i^n m1m2∣Sin (即 S i n S_i^n Sin 被 m 1 m 2 m_1^{m_2} m1m2 整除)时, n n n 最小。
显然对于 , 1 ≤ m 1 ≤ 30000 , 1 ≤ m 2 ≤ 10000 ,1≤m_1≤30000,1≤m_2≤10000 ,1≤m1≤30000,1≤m2≤10000 这样的数据,我们是存不下 M M M 的数值的,更不用说对 M M M 进行一系列的操作和计算了。但是题目中给出的 M = m 1 m 2 M=m_1^{m_2} M=m1m2 的条件很有意思,结合着我们需要满足 m 1 m 2 ∣ S i n m_1^{m_2}|S_i^n m1m2∣Sin 的目的,不难看出这个条件是暗示让我们对 M M M 和 S i S_i Si 进行质因数分解,然后再通过质因数和其指数判断是否能满足 m 1 m 2 ∣ S i n m_1^{m_2}|S_i^n m1m2∣Sin 。
而这样,我们就成功地将一个乘方的计算,转换成了指数相乘的计算。( 若 m 1 = p 1 i p 2 j ⋯ p n k , 则 m 1 m 2 = p 1 m 2 i p 2 m 2 j ⋯ p n m 2 k 若m_1=p_1^ip_2^j \cdots p_n^k,则m_1^{m_2}=p_1^{m_2i}p_2^{m_2j} \cdots p_n^{m_2k} 若m1=p1ip2j⋯pnk,则m1m2=p1m2ip2m2j⋯pnm2k)
怎么判断是否满足 m 1 m 2 ∣ S i n m_1^{m_2}|S_i^n m1m2∣Sin 呢,举个例子:
p 1 p_1 p1 | p 2 p_2 p2 | p 3 p_3 p3 | p 4 p_4 p4 | p 5 p_5 p5 | p 6 p_6 p6 | |
---|---|---|---|---|---|---|
M M M | 2 2 2 | 3 3 3 | 3 3 3 | 3 3 3 | 4 4 4 | 0 0 0 |
S i n S_i^n Sin | 1 ∗ n 1*n 1∗n | 2 ∗ n 2*n 2∗n | 1 ∗ n 1*n 1∗n | 1 ∗ n 1*n 1∗n | 1 ∗ n 1*n 1∗n | 0 ∗ n 0*n 0∗n |
只要满足对于 M M M 的每一个质因数 p p p ,它的指数都要比 S i n S_i^n Sin 来的小即可。在这里, 1 ∗ n ≥ 5 1*n\ge5 1∗n≥5 即可,取最小的 n n n 为 5 5 5。我们可以循环每一个 M M M 的质因数,然后让 n n n 取能够满足规定指数大小关系的最小值,最后对于每一个 S i S_i Si 打擂台求出最后 n n n 的最小值。
什么时候永远也不能满足要求呢?那就是当 M M M 中有 S i S_i Si 没有的质因数的时候。因为无论给 S i S_i Si 乘多少次方,都不可能创造新的质因数。就比如若 M M M 的质因子 p 6 p_6 p6 指数为 1 1 1 时,就不可能满足要求。
最后,怎么把质因数和它的指数存下来呢,我们注意到题目规定 1 ≤ m 1 ≤ 30000 1≤m_1≤30000 1≤m1≤30000 ,如果开 30000 30000 30000 大小的数组就造成了浪费,小于 30000 30000 30000 的质数只有 3300 3300 3300 左右个,我们可以在初始化时先跑一遍质数存下来,以便后期使用。
对了,还有一个特判,也就是当 m 1 = 1 m_1=1 m1=1 时,怎样都可以均分,应该直接输出 0 0 0 ,否则进入后面由于 1 1 1 不是质数,会输出 − 1 -1 −1
完整代码
感觉代码写的有点臃肿,见谅,读者读懂其意思即可
#include<bits/stdc++.h>
using namespace std;
int prime[3310],mx[3310],sx[3310],s[10010];//存下质因数,M的指数,Si的指数,还有所有的s
bool isprime(int n)//判断质数
{
for(int i=2;i*i<=n;i++)//用i*i<=n免去了计算sqrt(n)的过程
if(n%i==0)return false;
return true;
}
int main()
{
int n,m1,m2,tot=0,minn=1<<30,flag=0,np,ff=0;//tot记录质数总数,minn打擂台,flag记录某一个Si是否有结果,ff记录所有Si是否有结果
prime[++tot]=2;
for(int i=3;i<=30000;i+=2)//跳着判断,心里觉得会快点,实际上差别不大
if(isprime(i))
prime[++tot]=i;
scanf("%d %d %d",&n,&m1,&m2);
for(int i=1;i<=n;i++)scanf("%d",&s[i]);
if(m1==1)
{
printf("0");
return 0;
}
for(int i=1;prime[i]<=m1;i++)//对m1进行质因数分解
{
while(m1%prime[i]==0)
{
m1/=prime[i];
mx[i]++;
}
mx[i]*=m2;//在这里就相当于m1^m2的操作
}
for(int i=1;i<=n;i++)
{
memset(sx,0,sizeof(sx));//每次清空si的指数的数组备用
for(int j=1;j<=tot;j++)//对Si进行质因数分解,由于m<30000,只需要判断Si是否有30000以下的质因数就行了
{ //而且我们只存了30000以下的质数,用for循环判断prime[i]<=s[i]会RE
while(s[i]%prime[j]==0)
{
s[i]/=prime[j];
sx[j]++;
}
}
np=1;//每一次的n
flag=0;
for(int j=1;j<=tot;j++)
{
if(mx[j]&&sx[j])//如果这个质因数两者都有
{
flag=1;
if(mx[j]>sx[j]*np)np=mx[j]/sx[j];
if(mx[j]>sx[j]*np)np++;
}
else if(mx[j]&&!sx[j]){flag=false;break;}//没有则不可能满足要求,退出循环
}
if(flag)minn=min(minn,np),ff=1;//打擂台
}
if(ff)printf("%d",minn);
else printf("-1");
return 0;
}
写在后面
2009 2009 2009 年的题目好喜欢考数学题,普及和提高都有一道,还都是 H a n k s Hanks Hanks 家族的,有点意思
总结:注意特判,否则又要扣分啊,这里是 m 1 = 1 m_1=1 m1=1 时需要特判