2020牛客寒假算法基础集训营6 E-立方数 - 质因子分解

题目链接：https://ac.nowcoder.com/acm/contest/3007/E
题目大意：


我们要求最大的A。
我们根据唯一分解定理：
$$N=e_1^{p_1}\ast e_2^{p_2}\ast e_2^{p_2}\ast e_2^{p_2}\ast ...e_k^{p_k}$$
那么
$$如果pi\%3==0那么A\ast =e_i\ast (pi/3)$$
现在我们就是要$$枚举ei来得到A$$
我们发现$$e_i的最大值为N^{1/3}$$
我们枚举$$N^{1/3}里的素数那么复杂度是O(T{N}^{(1/3)}/ln(N^{(1/3)}))$$
我们可以进一步优化：$$如果枚举到N^{1/4},那么我们把N在N^{1/4}里的素数全部筛出后得到X$$
$$那么X的形式为e^p\ast B,并且e,B>N^{1/4}，那么可以知道：p<4。如果p>=4那么e^p\ast B>N$$
$$\begin{array}{l}1:p=3那么B一定等于1证明N^{3/4}\ast B,而且B>N^{1/4}那么N^{3/4}\ast B>N不可能,所以B=1\\ \\ 2:p<3那么ei对结果没有贡献\end{array}$$
$$所以我们枚举N^{1/4}的素数,再判断X是否为立方数，复杂度:O(T{N}^{(1/4)}/ln(N^{(1/4)}))$$

#include<bits/stdc++.h>
#define LL long long
using namespace std;
 
#define MAXN 1000005
#define MAXL 5299710
LL prime[MAXN];
LL Lf[MAXN];//存储立方
LL check[MAXL];
LL tot = 0;
void get_prime(LL n){
    memset(check, 0, sizeof(check));
    for(LL i = 2; i <= n; ++i){
      if (!check[i])
        prime[tot++] = i, Lf[tot-1]=i*i*i;
 
      for (LL j = 0; j < tot; ++j){
        if (i * prime[j] > n)
          break;
 
        check[i*prime[j]] = 1;
        if (i % prime[j] == 0)
          break;
      }
    }
}
 
int main(){
 
    get_prime(31622);
    int t;scanf("%d", &t);
    while(t--){
        LL n, A=1;scanf("%lld", &n);
        for(int i=0; i<tot&&prime[i]<=n; i++){
            while(n%Lf[i]==0){//立方筛完
                A*=prime[i];
                n/=Lf[i];
            }
            while(n%prime[i]==0){
                n/=prime[i];
            }
        }
 
        int lb = 1, rb = 1000000;//二分判断X是不是立方数
        while(lb <= rb){
            int md = lb + rb >> 1;
            if((LL)md * md * md < n)
                lb = md + 1;
            else
                rb = md - 1;
        }
        if((LL)lb * lb * lb == n)
            A *= lb;
        printf("%lld\n", A);
    }
 
    return 0;
}