HDU 6623 Minimal Power of Prime

米勒羅賓超時，別試着優化了
比較容易想到先把小區間的素數篩出來，然後除掉這部分，然後剩下的因子肯定不會很多，直接特判就行了，我看的題解是用的[1-1e4]的區間去篩，然後剩下的質因子不會超過四個。
則：

• $n= x^{4}$，ans = 4
• $n= x^{3}$，ans = 3
• $n= x^{2}$，ans = 2

否則 ans =1 ，可以嘗試用手寫一遍有哪幾種情況.

然後問題變爲了如何判斷　n 爲一個數的幾次方，很容易想到二分（這個很穩定，也好寫），或者 pow 函數，但是　pow 會有誤差，畢竟　double 和 long double　精度都沒有long long 高，並且　pow，真香

注：我起初用的是（　判斷 n 是否爲　x　的3次方　）

LL s=pow(n+0.5,1.0/3);

然後ＷＡ到自閉，然後用

LL s=pow(n,1.0/4)+0.3;

就過了，經過漫長的debug,後面發現 當 上面這個式子誤差第一個式子誤差挺大的，比如$n＝100007^{3}$，兩個s不一樣!!!，就這樣有了較大誤差(玄學誤差)

#include<vector>
#include<algorithm>
#include<iostream>
#include<cstdio>
#include<map>
#include<cmath>
#include<queue>
#include<cstring>
#define LL long long
using namespace std;
const int N=1e4+5;
const int M=3e6+5;
const int INF=1e5;
int pri[N],len;
void getpri(){
    for(int i=2;i<N;i++)
    if(!pri[i]){
        pri[++len]=i;
        for(int j=i+i;j<N;j+=i)
            pri[j]=1;
    }
}
int main(){

    getpri();
    int t;cin>>t;
    while(t--){
        LL n;scanf("%lld",&n);
        int ans=INF;
        for(int i=1;i<=len;i++)
        {
            if(n==1)break;
            if(n%pri[i])continue;
            int x=0;
            while(n%pri[i]==0)n/=pri[i],x++;
            ans=min(x,ans);
        }
        if(n==1){printf("%d\n",ans);continue;}

        LL s=pow(n,1.0/4)+0.5;
        if(s*s*s*s==n){printf("%d\n",min(ans,4));continue;}
        s=pow(n,1.0/3)+0.5;
        if(s*s*s==n){printf("%d\n",min(ans,3));continue;}
        s=pow(n,1.0/2)+0.5;
        if(s*s==n){printf("%d\n",min(ans,2));continue;}
        printf("1\n");
    }

}