信息学奥赛一本通 1947：【09NOIP普及组】细胞分裂 | 洛谷 P1069 [NOIP2009 普及组] 细胞分裂

【题目链接】

ybt 1947：【09NOIP普及组】细胞分裂
洛谷 P1069 [NOIP2009 普及组] 细胞分裂

【题目考点】

1. 质因数分解

每个整数x都可以进行质因数分解，写成公式为$x=p_1^{a_1}{p}_2^{a_2}...p_n^{a_n}$，其中$p_1\sim p_n$都为s的质因数。

【解题思路】

本题可以抽象为：有数字$s_1,s_2,...,s_n$，s指其中的某个数字。看在这个数字序列中，哪个个数字$s_j$的幂$s_j^i$能整除$m_1^{m_2}$，且幂指数最小。
根据质因数分解公式$x=p_1^{a_1}{p}_2^{a_2}...p_n^{a_n}$，我们可以设计一个表示幂指数的结构，里面包含底数与指数。
将$m_1$分解质因数，得到$m_1=p_1^{a_1}{p}_2^{a_2}...p_u^{a_u}$，
那么$m_1^{m_2}=(p_1^{a_1}{p}_2^{a_2}...p_u^{a_u}{)}^{m_2}=p_1^{a_1{m}_2}{p}_2^{a_2{m}_2}...p_u^{a_u{m}_2}$
将每个s都分解质因数，得到$s_j=q_1^{b_1}{q}_2^{b_2}...q_t^{b_t}$

如果$s_j^i$能整除$m_1^{m_2}$，那么$s_j$的质因数必须包含$m_1$的所有质因数。
在分解$m_1$的质因数的同时，将$m_1$的所有质因数的1次方乘积，得到数字$m_0$。如果$s_j$包含$m_1$的所有质因数，那么$s_j$一定可以整除$m_0$
如果每个$s_j$的质因数都不能包含$m_1$的所有质因数，则输出-1。

接下来求$s_j$的指数i：
对于$m_1$中的每个质因数$p_k$，它在m1中的指数为$a_k$，那么在$m_1^{m_2}$中$p_k$的指数为$a_k{m}_2$，在$s_j$中$p_k$的指数为$b_k$，那么$s_j^i$中$p_k$的指数为$i\cdot b_k$。要想让$s_j^i$能整除$m_1^{m_2}$，必须有$i\cdot b_k\ge a_k{m}_2$，所以$i\ge \frac{a_k{m}_2}{b_k}$，i最小为$\lceil \frac{a_k{m}_2}{b_k}\rceil$。
遍历所有可能的k，求$\lceil \frac{a_k{m}_2}{b_k}\rceil$的最大值，即为能使$s_j^i$整除$m_1^{m_2}$的i的最小值。
对每个$s_j$都求出其对应的能使$s_j^i$整除$m_1^{m_2}$的i的最小值，比较得到其中i的最小值，即为问题的结果。

特殊地，如果$m_1$是1，那么任何1个细胞在0秒后都有1个，是符合要求的，输出0。

【题解代码】

解法1：使用结构体及vector

#include <bits/stdc++.h>
using namespace std;
#define N 10005
#define INF 0x3f3f3f3f
struct Pow
{
    int b, i;//b底数，i指数  
    Pow(){}
    Pow(int x, int y):b(x),i(y){}
};
int n, m1, m2, m0 = 1, s[N], mini = INF;
vector<Pow> v_m;
void initVM()
{
    int a = 2, x = m1;
    while(a <= x)
    {
        if(x % a == 0)
        {
            if(v_m.empty() == false && v_m.back().b == a)
                v_m.back().i++;
            else
                v_m.push_back(Pow(a,1));
            x /= a;
        }
        else
            a++;
    }
    for(i = 0; i < v_m.size(); ++i)
    {
        m0 *= v_m[i].b;//所有m1的质因子的一次方的乘积 
        v_m[i].i *= m2;//m的所有质数乘m2 
    }
}
int main()
{
    cin >> n >> m1 >> m2;
    for(int i = 1; i <= n; ++i)
        cin >> s[i];
    if(m1 == 1)//任何细胞0秒时都有1个 
    {
        cout << 0;
        return 0;
    }
    initVM();
    for(int j = 1; j <= n; ++j)
    {
        if(s[j] % m0 != 0)//s[i]不包含m1中的某些质因子 
            continue;
        int mxi = 0;//m的质因子个数除s[j]的质因子个数（向上取整）得到的最大倍数 
        for(int i = 0; i < v_m.size(); ++i)
        {
            int a = s[j], k = 0;//s[j]中包含k个v_m[i].b 
            while(a % v_m[i].b == 0)
            {
                a /= v_m[i].b;
                k++;
            }
            mxi = max(mxi, (int)ceil((double)v_m[i].i/k));
        }
        mini = min(mini, mxi);
    }
    if(mini == INF)
        cout << -1;
    else
        cout << mini;
    return 0;
}

解法2：只使用数组

#include <bits/stdc++.h>
using namespace std;
#define N 10005
#define INF 0x3f3f3f3f
int s[N], n, m1, m2, m0 = 1, mini = INF, m_b[N], m_i[N], mn;//m_b[i]:m的第i个质因数，m_i[i]:m的第i个质数因数的指数  mn:个数 
int main()
{
    int a, x, k, mxi;
    cin >> n >> m1 >> m2;
    for(int i = 1; i <= n; ++i)
        cin >> s[i];
    if(m1 == 1)//任何细胞0秒时都有1个 
    {
        cout << 0;
        return 0;
    }
    a = 2, x = m1;//对m做质因数分解 
    while(a <= x)
    {
        if(x % a == 0)
        {
            if(m_b[mn] == a)
                m_i[mn]++;
            else
            {
                m_i[mn] *= m2;//m的所有质数乘m2         
                m_b[++mn] = a;
                m_i[mn] = 1;
                m0 *= a;//所有m1的质因子的一次方的乘积 
            }
            x /= a;
        }
        else
            a++;
    }
    m_i[mn] *= m2;
    for(int j = 1; j <= n; ++j)
    {
        if(s[j] % m0 != 0)//s[i]不包含m1中的某些质因子 
            continue;
        mxi = 0;//m的质因子个数除s[j]的质因子个数（向上取整）得到的最大倍数 
        for(int i = 1; i <= mn; ++i)
        {
            a = s[j], k = 0;//s[j]中包含k个v_m[i].b 
            while(a % m_b[i] == 0)
            {
                a /= m_b[i];
                k++;
            }
            mxi = max(mxi, (int)ceil((double)m_i[i]/k));
        }
        mini = min(mini, mxi);
    }
    cout << (mini == INF ? -1 : mini);
    return 0;
}