题目描述
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
说明/提示
【输入输出说明】
经过 11 1秒钟,细胞分裂成3 33 个,经过2 2 2秒钟,细胞分裂成9 9 9个,……,可以看出无论怎么分裂,细胞的个数都是奇数,因此永远不能分入 22 2个试管。
【输入输出样例22 2说明】
第 111 种细胞最早在333 秒后才能均分入242424 个试管,而第222 种最早在222 秒后就可以均分(每试管144/(241)=6144/(241)=6144/(241)=6 个)。故实验最早可以在222 秒后开始。
【数据范围】
对于 50%的数据,有m1m2≤30000m_1^{m_2} ≤ 30000m1m2≤30000。
对于所有的数据,有1≤N≤10000,1≤m1≤30000,1≤m2≤10000,1≤Si≤2,000,000,0001 ≤N≤ 10000,1 ≤m_1 ≤ 30000,1 ≤m_2 ≤ 10000,1 ≤ S_i ≤ 2,000,000,0001≤N≤10000,1≤m1≤30000,1≤m2≤10000,1≤Si≤2,000,000,000。
NOIP 2009 普及组 第三题
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <set>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long ll;
const ll inf = 1e9;
const int maxn = 1e5 + 10;
int prim[maxn];
bool flag[maxn];
ll a[maxn];
ll b[maxn];
int n, m1, m2;
int tot;
ll getValue(ll x) {
memset(b, 0, tot * sizeof(int));
ll sqr = sqrt(x);
for(int i = 0; i < tot; i++) {
while(x % prim[i] == 0) {
x /= prim[i];
b[i]++;
}
}
// for(int i = 0; i < tot; i++)
// if(b[i]) cout << prim[i] << " " << b[i] << endl;
ll tem = 0;
for(int i = 0; i < tot; i++) {
if(a[i] == 0) continue;
if(b[i] == 0) return inf;
else tem = max(tem, max(0ll, a[i] / b[i] + bool(a[i] % b[i]) ));//, cout << prim[i] << " " << a[i] - b[i] << endl;
}
return tem;
}
void getPrim() {
int max_num = 30010;
for(int i = 2; i <= max_num; i++) {
if(flag[i] == 0) prim[tot++] = i;
for(int j = 0; j < tot; j++) {
if(prim[j] * i > max_num) break;
flag[i * prim[j]] = true;
if(i % prim[j] == 0) break;
}
}
}
int main() {
getPrim();
cin >> n >> m1 >> m2;
for(int i = 0; prim[i] <= m1; i++) {
while(m1 % prim[i] == 0) {
m1 /= prim[i];
a[i] += m2;
}
}
// for(int i = 0; i < tot; i++)
// if(a[i]) cout << prim[i] << " " << a[i] << endl;
ll ans = inf, data;
for(int i = 0; i < n; i++) {
cin >> data;
ans = min(ans, getValue(data));
}
if(ans == inf) cout << -1 << endl;
else cout << ans << endl;
return 0;
}