个人训练 6 H

H. Ancient Wisdom

题目链接
题目大意：给出一个式子CD^3 = A^2,1<=c<1e18,求出D的最小值。A和D均为整数。
思路：
将式子转化为C=A^2/ D^3 = CD^3/ D^3.要令C * D^3的开方是一个整数，即只需算sqrt(CD)是一个整数即可。
令X是一个非0正整数，有CD=X^2,已知D的最大值定为C，令D=C，
然后对D进行因数分解，能够开方出整数的因数就移到到左边，变为C * D/T=A^2/T,(sqrt(T)为整数），直至最后分解出这类因数，D/T就是最终结果。也就是要分解C的质因数来解题。
然后就是处理因数了，我们可以算出3e6以内的所有质因数，然后开始分解因数，结果用ans来存。质因数从小开始枚举，每当出现奇数个相同的质因数peime[i]，就令ans=ansprime[i]。接下来解释一下为什么只需算3e6以内的质因数，因为（3e6）^3>=1e18的，因此肯定不会存在三个大于3e6的质因数，所以最后剩余的是两个大质数之积或者一个大质数，最后只需判断剩余的是否可以开方出整数即可，如果不可以，那么ans就是C。
代码：

#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <queue>
#include <vector>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
ll pri[3000005],vis[3000005];
int tot;
void get_prime(){
    tot = 0;
    for (int i=2;i<=3000000;i++) {
        if (!vis[i]) pri[tot++]=i;
        for (int j=0;j<tot&&i*pri[j]<=3000000;j++) {
            vis[i*pri[j]]=1;
            if (i%pri[j]==0) break;
        }
    }
}
int main(){
	get_prime();
	ll res=1,c;
	cin>>c;
	for(int i=0;i<tot;i++){
		ll x=pri[i];
		if(c%x==0){
			int cnt=0;
			while(c%x==0){
				cnt++;
				c/=x;
			}
			if(cnt%2==1) res*=x;
		}
	}
	ll d=sqrt(c);
	if(d*d!=c&&d!=0){
		res*=c;
	}
	cout<<res<<endl;
	return 0;
}