牛客数学 最小的指数

博客围绕最小的指数题目展开,牛妹有 T 个超级大的数,牛牛想将数分解质因数并记录质因子指数最小值,数多处理繁琐便请人帮忙。给出了输入输出描述及样例,还有数据范围 1≤T≤100000,1≤n≤1018。

最小的指数题目链接
题目描述
牛妹手里捧着 T T T 个数,但每个数都超级超级大。
牛牛有强迫症,一看到数就想把它分解质因数,并记录下它的质因子指数的最小值。
形式化地,对于数 x = p 1 a 1 p 2 a 2 . . . p k a k x=p_1^{a_1}p_2^{a_2}...p_k^{a_k} x=p1a1p2a2...pkak ( a 1 , a 2 , . . . , a k > 0 ) (a_1,a_2,...,a_k > 0) (a1,a2,...,ak>0) ,牛牛会记录 m i n { a 1 , a 2 , . . . , a k } min\{a_1,a_2,...,a_k\} min{a1,a2,...,ak}
但是数一多,牛牛感觉处理起来过于繁琐,所以想请你来帮他快速计算。
特别地,如果 x = 1 x = 1 x=1 ,则最小值为 0 。
输入描述
第一行,输入 T T T
接下来 T T T 行,每行输入一个数 x x x
输出描述
输出 T T T 行,每行输出答案。
输入样例

8
6
9
12
108
97
100
2333
2147483648

输出样例

1
2
1
2
1
2
1
31

1 ≤ T ≤ 100000 , 1 ≤ n ≤ 1 0 18 1≤T≤100000,1≤n≤10^{18} 1T100000,1n1018
在这里插入图片描述

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
typedef long long ll;
const int maxn = 1e4 + 7;
int primes[maxn];
int cnt = 0;
bool st[maxn];
void get_primes()
{
    for(int i = 2; i <= maxn; i++)
    {
        if(!st[i]) primes[cnt++] = i;
        for(int j = 0; primes[j] <= maxn/i; j++)
        {
            st[primes[j]*i] = 1;
            if(i%primes[j] == 0) break;
        }
    }
}
int main()
{
    get_primes();
    int t;
    scanf("%d",&t);
    while(t--)
    {
        ll x;
        scanf(" %lld",&x);
        if(x == 1)
        {
            puts("0");
            continue;
        }
        ll ans = 1e9, ci = 0;
        for(int i = 0; i < cnt; i++)
        {
            if(x%primes[i] == 0)
            {
                ci = 0;
                while(x%primes[i] == 0)
                {
                    ci++;
                    x /= primes[i];
                }
                ans = min(ans,ci);
            }
        }
        ll line4, line3, line2;
        if(x > 1)
        {
            line2 = (ll) sqrt(x*1.0);
            if(line2*line2 == x)
            {
                line4 = (ll) sqrt(1.0*line2);
                if(line4*line4 == line2) ci = 4;
                else ci = 2;
            }
            else
            {
                ll l = 1, r = (ll) pow(1.0*x, 1.0/3.0) + 1, mid;
                while( l < r)
                {
                    mid = l + r >> 1;
                    if(mid * mid * mid > x) r = mid - 1;
                    else if(mid * mid * mid < x) l = mid + 1;
                    else break;
                }
                if(l*l*l == x) ci = 3;
                else ci = 1;
            }
            ans = min(ans,ci);
        }
        printf("%lld\n",ans);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

幸愉信奥

谢谢亲的支持,我会继续努力啦~

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值