2019HDU多校 6223 （最小质因子指数幂）

传送门

题解：考虑先筛出$1e4$以内的素数，对于每个输入的 $n$ 先除掉 $1e4$ 以内的素因子（如果存在的话），并且更新答案。

如果除完后$n$ 还大于1，说明存在大于$1e4$的素因子，而且这些素因子的最多只会有$4$次方，因为更大的话就会超过
$1e18$,因此剩下的$n$只会是 $p$² , $p$³, $p$⁴。

$p$²和$p$⁴可以一起判断， $p$³则可以通过二分判断

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<stdio.h>
#include<string.h>
#include<queue>
#include<cmath>
#include<map>
#include<set>
#include<vector> 
using namespace std;
#define inf 0x3f3f3f3f
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define mem(a,b) memset(a,b,sizeof(a));
#define lowbit(x)  x&-x;  
#define debugint(name,x) printf("%s: %d\n",name,x);
#define debugstring(name,x) printf("%s: %s\n",name,x);
typedef long long ll;
typedef unsigned long long ull;
const double eps = 1e-6;
const int maxn = 1e4+5;
const int mod = 1e9+7;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int prime[maxn],tot = 0;
bool pri[maxn];
void init(){
    pri[1] = 0;
    for(int i = 2; i < maxn; i++){
        if(!pri[i]){
            prime[tot++] = i;
            for(int j = i+i; j < maxn; j+=i){
                if(!pri[j]){
                    pri[j] = 1;
                }
            }
        }
    }
}
ll n;
int main() {
    init();
    int t;
    cin>>t;
    while(t--){
        scanf("%lld",&n);
        int ans = 100;
        for(int i = 0; i < tot && n >= prime[i]; i++){
            if(n%prime[i]==0){
                int cnt = 0;
                while(n%prime[i]==0){
                    cnt++;
                    n/=prime[i];
                }
                ans = min(ans,cnt);
            }
        }
        if(n > 1){
			ll c1 = (ll)sqrt(1.0*n);
			if(c1*c1 == n){
				//ans = min(ans,2);
				ll c2 = (ll)sqrt(1.0*c1);
				if(c2*c2==c1)
					ans = min(ans,4);
				else
					ans = min(ans,2);
			}else{
				ll l = 1, r = (ll)pow(1.0*n,1.0/3.0)+1, mid;
				ll res = 0; 
				while(l <= r){
					mid = l+r>>1;
					if(mid*mid*mid > n){
						r = mid-1;	
					}
					else if(mid*mid*mid < n){
						l = mid+1;
					}else
						break;
				}
				if(mid*mid*mid == n) ans = min(ans,3);
				else ans = 1; 
			}
        }
        printf("%d\n",ans);
    }
}