杭电2019多校第四场 Minimal Power of Prime (数论+思维)

                                    Minimal Power of Prime

                      Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
                                              Total Submission(s): 5331    Accepted Submission(s): 328

                                                                 Statistic | Submit | Clarifications | Back

题意：把一个数质因数分解，求所有质因数中的最小指数

题解：打印1-10000之间的素数表然后质因数分解，分解完剩下的那个数，最多是某个质数的4次方（因为这个质数必然大于1e4），所以讨论一下就好了，当然要注意：先看是4次方再看2次方（因为如果满足4次方一定满足2次方，但是满足4次方也满足2次方，2次方的话就不是质因数了），3次方无所谓，因为开3次方会损失精度，所以就二分一下~~最后如果2,3,4次方都不满足，那么它就是一个质数

为什么是打印1-10000不是别的区间，因为1-20000就超时

 

比赛的时候想的2，3次方没有清楚为啥最多到几次方有点坑

/*
Niubility_n
*/
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int maxn=1e4+10;
ll prime[maxn];
int vis[maxn];
int cnt;
void isprime()
{
    vis[0]=vis[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!vis[i])prime[cnt++]=i;
        for(int j=0;j<cnt && i*prime[j]<maxn;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)break;
        }
    }
}
int judge(ll n)
{
    ll l=10010,r=1000000;
    while(l<=r)
    {
        ll mid=(l+r)>>1;
        if((1ll*mid*mid*mid)>n)r=mid-1;
        else if((1ll*mid*mid*mid)<n)l=mid+1;
        else return 1;
    }
    return 0;
}
int main()
{
    int t;
    isprime();
    /*
    for(int i=0;i<cnt;i++)printf("%lld ",prime[i]);
    printf("\n");
    */
    scanf("%d",&t);
    while(t--)
    {
        ll n;
        scanf("%lld",&n);
        int ans=1000;
        for(int i=0;i<cnt;i++)
        {
            if(n%prime[i]==0)
            {
                int sum=0;
                while(n%prime[i]==0)
                {
                    sum++;
                    n/=prime[i];
                }
                ans=min(sum,ans);
            }
            if(n==1)break;
        }
        if(n>10009)
        {
            ll ans1=(ll)sqrt(n);
            ll ans2=(ll)sqrt((ll)sqrt(n));
            if(ans2*ans2*ans2*ans2==n)ans=min(ans,4);
            else if(ans1*ans1==n)ans=min(ans,2);
            else if(judge(n))ans=min(ans,3);
            else ans=min(ans,1);
        }
        printf("%d\n",ans);
    }
    return 0;
}