2019HDU多校集训——Minimal Power of Prime

题目大意,输出一个数n的质因数的最小次幂。想到了先筛出小部分的质因数,后面又想到了开方之类的操作,但是脑子不正常,把开方这里想歪了。。。。。。感觉最近状态太差了,有时候想题的时候像个憨憨一样。好像以前也有这种情况出现过,然后多睡觉就好了??

解题思路:根据唯一分解定理,一个合数必然是某些素数次方的乘积。所以先线性筛 筛出1e4范围内的素数,再去分解n,同时得到一个当前的最小次幂ans,如果最后n的值为1,那么很显然也就不用继续了,只需要判断最小次幂就行了。但是重点很显然是分解完以后n不为1,也就是1e4内的素数不能完全的分解n,换句话说就是n有那种大于1e4的质因子。

分类讨论:设v1,v2,v3,v4,分别是大于1e4的质数。

那么n可能是v1^4,这种情况结果为min(ans,4)

或者(v1^2)*(v2*v3),这种情况结果为1

或者v1*v2*v3*v4,这种情况结果为1

或者(v1^2)*(v2^2),这种情况结果为min(ans,2)

或者v1,这种情况为1

或者v1^2,这种情况结果为min(ans,2)

或者v1^3,这种情况结果为min(ans,3)。

其实就是判断能不能开2,3,4次方,否则结果为1。我用的是二分判断的,可能用pow再写一点处理精度误差的代码应该也可以。

最后,代码:

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e4;
int prime[maxn+10];
bool vis[maxn+10];
int tot;
void eular()
{
    for(int i=2;i<=maxn;++i)
    {
        if(!vis[i])
            prime[tot++]=i;
        for(int j=0;j<tot&&i*prime[j]<=maxn;++j)
        {
            vis[prime[j]*i]=true;
            if(i%prime[j]==0)
                break;
        }
    }
}
bool check(int k,long long int val)
{
    long long int l=1,r=pow(val,1.0/k)+1;
    long long int mid;
    long long int temp;
    bool flag;
    while(l<=r)
    {
        mid=l+r>>1;
        temp=1;
        for(int i=0;i<k;++i)
        {
            temp*=mid;
            if(temp>val)
            {
                break;
            }
        }
        if(temp>val)
        {
            r=mid-1;
        }
        else if(temp<val)
        {
            l=mid+1;
        }
        else
            return true;
    }
    return false;
}
int main()
{
    eular();
    ios::sync_with_stdio(false);
    int T;
    scanf("%d",&T);
    while(T--)
    {
        long long int n;
        scanf("%lld",&n);
        int ans=100;
        bool flag;
        for(int i=0;i<tot;++i)
        {
            int cnt=0;
            flag=false;
            while(n%prime[i]==0)
            {
                ++cnt;
                n/=prime[i];
                flag=true;
            }
            if(flag)
                ans=min(ans,cnt);
        }
        if(n==1||ans==1)
        {
            printf("%d\n",ans);
        }
        else
        {
            if(check(4,n))
                ans=min(ans,4);
            else if(check(3,n))
                ans=min(ans,3);
            else if(check(2,n))
                ans=min(ans,2);
            else
                ans=1;
            printf("%d\n",ans);
        }
    }
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值