J - Mysterious Bacteria(求最大的可行幂,质因子分解+负数判断)

题目

题意：给出一个数 x，求满足 x=p^k 的 k 值，其中 p、k 为整数。
思路: x=p1^e1 * p2^e2 * … * pn^en;
k=gcd(p1,p2…pn);
假如x为负数时,求出的k为偶数，这样是不可能成立的，因为任何数的偶次幂非负。此时将求出的k一直/2,直到k为奇数，就是x为负数时的答案。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#define en '\n'
using namespace std;
typedef long long ll;
const int N=65536+5;
int is[N],prime[N],tot;
void getprime()//素数打表
{
    tot=0;
    for(int i=2;i<=N-5;++i)
    {
        if(!is[i]) prime[++tot]=i;
        for(int j=1;j<=tot;++j)
        {
            if((ll)i*prime[j]>N-5) break;
            is[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
int gcd(int a,int b){
    return b?gcd(b,a%b):a;
}
int solve(ll n)//用打表后的素数，质因子分解
{
    if(!n) return 0;
    int res=1,flag=n<0;
    if(flag) n*=-1;
    for(int i=1;(ll)prime[i]*prime[i]<=n&&i<=tot;++i)
    {
        int num=0;
        while(n%prime[i]==0) ++num,n/=prime[i];
        if(num)
        {
            if(res==1) res=num;
            else
            {
                res=gcd(res,num);
                if(res==1) break;
            }
        }
    }
    if(flag)//负数时。
        while(!(res&1)) res>>=1;
    return res;
}
int main()
{
    getprime();
    int T,cas=0;scanf("%d",&T);
    while(T--)
    {
        ll n;scanf("%lld",&n);
        printf("Case %d: %d\n",++cas,solve(n));
    }
}