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;

            else

                temp = 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;

}