NOIP2009普及组细胞分裂(数论)——yhx
2016-05-02 16:45
381 查看
题目描述
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,以一个空格隔开,
即表示试管的总数 M = m1^m2。
第三行有 N 个正整数,第 i 个数 Si表示第 i 种细胞经过 1 秒钟可以分裂成同种细胞的个
数。
输出格式:
输出文件 cell.out 共一行,为一个整数,表示从开始培养细胞到实验能够开始所经过的
最少时间(单位为秒)。
如果无论 Hanks 博士选择哪种细胞都不能满足要求,则输出整数-1。
题目可以转化成求方程ax=km1m2(k为整数)的最小整数解。【我也来试一试打公式】
先把m1质因数分解,再把每个次数都乘上m2就是题中的M质因数分解的结果。
方程如果有解,那么必须m1的每个因数都是a的因数,这样a经过幂之后才可能成为m1的倍数。在此基础上,取各个因数中次数相差最多(指的是倍数)的一个的倍数,就是x的最小值。
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,以一个空格隔开,
即表示试管的总数 M = m1^m2。
第三行有 N 个正整数,第 i 个数 Si表示第 i 种细胞经过 1 秒钟可以分裂成同种细胞的个
数。
输出格式:
输出文件 cell.out 共一行,为一个整数,表示从开始培养细胞到实验能够开始所经过的
最少时间(单位为秒)。
如果无论 Hanks 博士选择哪种细胞都不能满足要求,则输出整数-1。
#include<cstdio> #include<cmath> struct prm { int x,t; }m[10010]; int a[10010]; int main() { int i,j,k,n,p,q,m1,m2,x,y,z,cntm,ans,t; bool ok; scanf("%d",&n); scanf("%d%d",&m1,&m2); cntm=0; for (i=2;m1!=1;i++) if (m1%i==0) { cntm++; m[cntm].x=i; while (m1%i==0) { m[cntm].t++; m1/=i; } m[cntm].t*=m2; } ans=-1; for (i=1;i<=n;i++) scanf("%d",&a[i]); for (i=1;i<=n;i++) { p=0; ok=1; for (j=1;j<=cntm;j++) if (a[i]%m[j].x==0) { t=0; while (a[i]%m[j].x==0) { a[i]/=m[j].x; t++; } if (m[j].t%t==0) t=m[j].t/t; else t=m[j].t/t+1; if (t>p) p=t; } else { ok=0; break; } if (ok&&(ans==-1||p<ans)) ans=p; } printf("%d\n",ans); }
题目可以转化成求方程ax=km1m2(k为整数)的最小整数解。【我也来试一试打公式】
先把m1质因数分解,再把每个次数都乘上m2就是题中的M质因数分解的结果。
方程如果有解,那么必须m1的每个因数都是a的因数,这样a经过幂之后才可能成为m1的倍数。在此基础上,取各个因数中次数相差最多(指的是倍数)的一个的倍数,就是x的最小值。
相关文章推荐
- uva 11572 unique snowflakes——yhx
- uva 1471 defence lines——yhx
- NOIP2010提高组 关押罪犯 -SilverN
- NOIP2011提高组 聪明的质监员 -SilverN
- NOIP2010普及组T3 接水问题 ——S.B.S.
- NOIP2010普及组T4 三国游戏——S.B.S.
- 【NOIP提高组2015D2T1】uva 714 copying books【二分答案】——yhx
- 【NOIP合并果子】uva 10954 add all【贪心】——yhx
- #include <NOIP2009 Junior> 细胞分裂 ——using namespace wxl;
- #include <NOIP2008 Junior> 双栈排序 ——using namespace wxl;
- NOIP2010普及组 三国游戏 -SilverN
- NOIP2009 提高组T3 机器翻译 解题报告-S.B.S
- NOIP2010提高组乌龟棋 -SilverN
- NOIP2010提高组 机器翻译 -SilverN
- UVa 297 Quadtrees -SilverN
- UVa 673 Parentheses Balance -SilverN
- uva 1606 amphiphilic carbon molecules【把缩写写出来,有惊喜】(滑动窗口)——yhx
- uva 11054 wine trading in gergovia (归纳【好吧这是我自己起的名字】)——yhx
- NOIP2008提高组(前三题) -SilverN
- uva 11134 fabled rooks (贪心)——yhx