题目
题目描述
Hanks 博士是BT (Bio-Tech,生物技术) 领域的知名专家。现在,他正在为一个细胞实 验做准备工作:培养细胞样本。
Hanks 博士手里现在有N 种细胞,编号从1~N,一个第i 种细胞经过1 秒钟可以分裂为 Si 个同种细胞(Si 为正整数)。现在他需要选取某种细胞的一个放进培养皿,让其自由分裂, 进行培养。一段时间以后,再把培养皿中的所有细胞平均分入M 个试管,形成M 份样本, 用于实验。Hanks 博士的试管数M 很大,普通的计算机的基本数据类型无法存储这样大的 M 值,但万幸的是,M 总可以表示为m1 的m2 次方,即M =m1^m2 ,其中m1,m2 均为基本 数据类型可以存储的正整数。
注意,整个实验过程中不允许分割单个细胞,比如某个时刻若培养皿中有4 个细胞, Hanks 博士可以把它们分入2 个试管,每试管内2 个,然后开始实验。但如果培养皿中有5 个细胞,博士就无法将它们均分入2 个试管。此时,博士就只能等待一段时间,让细胞们继 续分裂,使得其个数可以均分,或是干脆改换另一种细胞培养。
为了能让实验尽早开始,Hanks 博士在选定一种细胞开始培养后,总是在得到的细胞“刚 好可以平均分入M 个试管”时停止细胞培养并开始实验。现在博士希望知道,选择哪种细 胞培养,可以使得实验的开始时间最早。
输入
第一行有一个正整数 N,代表细胞种数。
第二行有两个正整数 m1,m2,以一个空格隔开,m1^ m2 即表示试管的总数M。
第三行有 N 个正整数,第i 个数Si 表示第i 种细胞经过1 秒钟可以分裂成同种细胞的个 数。
输出共一行,为一个整数,表示从开始培养细胞到实验能够开始所经过的 最少时间(单位为秒)。 如果无论 Hanks 博士选择哪种细胞都不能满足要求,则输出整数-1。
输出
共一行,为一个整数,表示从开始培养细胞到实验能够开始所经过的 最少时间(单位为秒)。 如果无论 Hanks 博士选择哪种细胞都不能满足要求,则输出整数-1。
样例
样例输入 | 样例输出 |
| |
| |
数据范围
对于 50%的数据,有m1^m2 ≤ 30000。
对于所有的数据,有1 ≤N≤ 10000,1 ≤m1 ≤ 30000,1 ≤m2 ≤ 10000,1 ≤ Si ≤ 2,000,000,000。
提示
【输入输出样例1 说明】 经过 1 秒钟,细胞分裂成3 个,经过2 秒钟,细胞分裂成9 个,……,可以看出无论怎么分 裂,细胞的个数都是奇数,因此永远不能分入2 个试管。
【输入输出样例2 说明】 第 1 种细胞最早在3 秒后才能均分入24 个试管,而第2 种最早在2 秒后就可以均分(每 试管144/(241)=6 个)。故实验最早可以在2 秒后开始。
题目大意
问时ans的最小值。
算法
这道题是数论。
我们把m1分解质因数:×××……×,即
×××……×
于是,我们发现,如果si含有每一个,那么i细胞才能分裂成的倍数个。
如果可以,那么就要求要分裂多少次了。我们设x表示si包含多少个,t表示至少分裂多少次,则 ,即t=max{t,}。
最后我们用t更新ans即可。
此算法时间复杂度为。
注意事项
1.不开long long会WA。
2.质因数没有1,不然死循环。
代码
#include <cstdio>
#include <cmath>
#define ll long long
using namespace std;
int max(int a,int b);
int num;
ll m1,m2,p[30],a[30];
void split(ll x)
{
num=0;
for (int i=2; i<=sqrt(x); i++)
{
if (x%i!=0) continue;
num++;
p[num]=i;
a[num]=0;
while (x%i==0)
{
a[num]++;
x/=i;
}
a[num]*=m2;
}
if (x>1)
{
num++;
p[num]=x;
a[num]=m2;
}
}
int main()
{
int n,t;
ll s,ans=-1;
freopen("cell.in","r",stdin);
freopen("cell.out","w",stdout);
scanf("%d",&n);
scanf("%lld %lld",&m1,&m2);
split(m1);
for (int i=1; i<=n; i++)
{
t=0;
bool adj=true;
scanf("%lld",&s);
for (int j=1; j<=num; j++)
{
if (s%p[j]!=0)
{
adj=false;
break;
}
int x=0;
while (s%p[j]==0)
{
x++;
s/=p[j];
}
t=max(t,ceil((a[j]-1)/x)+1);
}
if (adj)
if (ans==-1 || t<ans)
ans=t;
}
printf("%lld",ans);
}
int max(int a,int b)
{
if (a>b) return a;
else return b;
}