Time Limit: 1 second
Memory Limit: 128 MB

【问题描述】

Hanks博士是BT(Bio-Tech，生物技术)领域的知名专家。现在，他正在为一个细胞实验做准备工作：培养细胞样本。
Hanks博士手里现在有N种细胞，编号从1~N，一个第i种细胞经过1秒钟可以分裂为Si个同种细胞（Si为正整数）。现在他需要选取某种细胞的一个放进培养皿，让其自由分裂，进行培养。一段时间以后，再把培养皿中的所有细胞平均分入M个试管，形成M份样本，用于实验。Hanks博士的试管数M很大，普通的计算机的基本数据类型无法存储这样大的M值，但万幸的是，M总可以表示为m1的m2次方，即M=m1m2，其中m1，m2均为基本数据类型可以存储的正整数。
注意，整个实验过程中不允许分割单个细胞，比如某个时刻若培养皿中有4个细胞，Hanks博士可以把它们分入2个试管，每个试管内2个，然后开始实验。但如果培养皿中有5个细胞，博士就无法将它们均分入2个试管。此时，博士就只能等待一段时间，让细胞们继续分裂，使得其个数可以均分，或是干脆改换另一种细胞培养。
为了能让实验尽早开始，Hanks博士在选定一种细胞开始培养后，总是在得到的细胞“刚好可以平均分入M个试管”时停止细胞培养并开始实验。现在博士希望知道，选择哪种细胞培养，可以使得实验的开始时间最早。

【输入格式】

共有3行。
第一行有1个整数N，代表细胞种数。
第二行有两个正整数m1，m2，以一个空格隔开，m1m2即表示试管的总数M。
第三行有N个正整数，第i个数表示第i种细胞经过1秒钟可以分裂为同种细胞的个数。

【输出格式】

共1行，为一个整数，表示从开始培养细胞到实验能够开始所经过的最小时间（单位为秒）。
如果无论Hanks博士选择哪种细胞都不能满足要求，则输出整数-1。

【输入样例1】

1
2 1
3

【输出样例1】

-1

【输入输出样例1说明】

经过1秒钟，细胞分裂成3个，经过2秒钟，细胞分裂成9个，……，可以看出无论怎么分裂，细胞的个数都是奇数，因此永远不能分入2个试管。

【输入样例2】

2
24 1
30 12

【输出样例2】

2

【输入输出样例2说明】

第1种细胞最早在3秒后才能均分入24个试管，而第2种最早在2秒后就可以均分(每试管144/(241)=6个)。故实验最早可以在2秒后开始。

【数据范围】

对于50%的数据，有m1m2≤30000。
对于所有的数据，有1≤N≤10000，1≤m1≤30000，1≤m2≤10000，1≤Si≤2,000,000,000。

【题目链接】:http://noi.qz5z.com/viewtask.asp?id=p093

【题意】

【题解】

/*qi^(ti)%(m1^m2)==0;ti要求最小;然后求min(ti);对m1质因数分解;分成a1*a2..ak的形式把a1..ak它们的指数处理出来->b[1..k];然后b[1..k]都乘上m2;对于n个数字;q[i]如果q[i]不能被a1..ak中的某个数字整除;则这个细胞无解;如果能,则处理出这个数字q[i]的质因数a1..ak的各个指数w[1..k]对于w[i]c[i] = (b[i]-1)/w[i] + 1;则time[第i个细胞]=max(c[1..k]);然后找time最小的细胞就好;或者记录无解;m1==1的情况特判一下。
*/

【完整代码】

#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define LL long long
#define rep1(i,a,b) for (int i = a;i <= b;i++)
#define rep2(i,a,b) for (int i = a;i >= b;i--)
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define rei(x) scanf("%d",&x)
#define rel(x) scanf("%lld",&x)typedef pair<int, int> pii;
typedef pair<LL, LL> pll;const int dx[9] = { 0,1,-1,0,0,-1,-1,1,1 };
const int dy[9] = { 0,0,0,-1,1,-1,1,-1,1 };
const double pi = acos(-1.0);
const int N = 1e4+100;
const int M = 3e4 + 100;int n,m1,m2,k;
int q[N],w[M],a[M],b[M],mi[N];
bool can[N];int main()
{//freopen("F:\\rush.txt", "r", stdin);rei(n);rei(m1), rei(m2);if (m1 == 1){puts("0");return 0;}rep1(i, 1, n)rei(q[i]);for (int i = 2;i<=m1;i++)if (m1%i == 0){a[++k] = i;while (m1%i == 0){m1 /= i;b[k]++;}}rep1(i, 1, k)b[i] *= m2;memset(can, true, sizeof can);rep1(i, 1, n){rep1(j, 1, k)w[j] = 0;int ju = -1;rep1(j, 1, k){if (q[i] % a[j] != 0){ju = j;break;}int temp = q[i];while (temp % a[j] == 0){temp /= a[j];w[j]++;}}if (ju != -1){can[i] = false;continue;}int d = -1;rep1(j, 1, k){int temp = (b[j] - 1) / w[j] + 1;if (d == -1)d = temp;elsetemp = max(temp, d);}mi[i] = d;}int ans = -1;rep1(i, 1, n)if (can[i] && (ans == -1 || mi[i] < ans))ans = mi[i];printf("%d\n", ans);return 0;
}