LightOJ1220 Mysterious Bacteria(思维+唯一分解)

题目大意:

给你一个整数n(可能为负数),让你求满足a^p=n的最大的p

思路:

当n是正数时,直接对n进行素因子分解,在对它的素因子的个数进行gcd,比如12=2^2*3,gcd(2,1)就是最大的p;

当n是负数时,则p的值一定是奇数,因为一个数的偶数次方一定为整数,因此需要将它的素因子个数全都化为奇数


这题主要在于n为负数时的思考!


#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<climits>
#include<cmath>
using namespace std;
typedef long long ll;
const int maxn = 1000100;
bool vis[maxn], flag;
ll prim[maxn/10], n;
int m;
void initPrime()
{
    m = 0;
    for(ll i = 2; i < maxn; i++)
    {
        if(!vis[i])
            prim[m++] = i;
        for(ll j = i*i; j < maxn; j += i)
            vis[j] = true;
    }
}
int solve(ll x)
{
    int ans = 0;
    for(int i = 0; i < m && prim[i]*prim[i] <= x; i++)
    {
        int res = 0;
        while(x % prim[i] == 0&&x)
            {x /= prim[i], res++;}
        if(flag)
            while(res%2==0&&res)
                res /= 2;
        if(res && !ans)
            ans = res;
        else if(res)
            ans = __gcd(ans, res);
    }
    if(x > 1)
        ans = 1;
    return ans;
}
int main()
{
    int T;
    initPrime();
    scanf("%d", &T);
    for(int kase = 1; kase <= T; kase++)
    {
        scanf("%lld", &n);
        flag = false;
        if(n < 0)
            n = -n, flag = true;
        printf("Case %d: %d\n", kase, solve(n));
    }
    return 0;
}



Sigma函数是指一个数字的所有因子之和。给定一个数字n,需要求出有多少个数字的Sigma函数是偶数。\[2\] 为了解决这个问题,可以先筛选出n范围内的素数(范围在10^6即可),然后对n进行素因子分解。对于每个因子,如果它的Sigma函数中连乘的每一项都是偶数,那么整个Sigma函数就是偶数。具体实现中,可以判断每个因子的平方根是否为偶数,如果是偶数,则减去(平方根+1)/2。\[1\] 另外,还可以使用O(1)的做法来解决这个问题。根据观察,所有的完全平方数及其两倍的值都会导致Sigma函数为偶数。因此,可以直接计算n的平方根,然后减去(平方根+1)/2即可得到结果。\[3\] #### 引用[.reference_title] - *1* [Sigma Function](https://blog.csdn.net/PNAN222/article/details/50938232)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* *3* [【LightOJ1336】Sigma Function(数论)](https://blog.csdn.net/qq_30974369/article/details/79009498)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值