Mysterious Bacteria LightOJ - 1220 (惟一分解定理）

题目链接在这里
x=b^p,给出x，求p的最大值。

x只有b一项组成，b需要包括所有的p1,p2,…,pn。b=p1^c1p2^c2…pn^cn。x=b^p=(p1^c1p2^c2…pn^cn)^p。那么p是p1,p2,…,pn的·最大公因子。
注意的是，x在有符号数int范围内，x为负数时，p只能为奇数。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define maxn 10000003
int su[maxn/10],cnt;
bool is[maxn+3];
int cas,factor[maxn/10];
void prime()
{
    cnt=1;
    memset(is,1,sizeof(is));
    is[0]=is[1]=0;
    for(int  i=2;i<maxn-1;i++)
    {
        if(is[i])
            su[cnt++]=i;
        for(int j=1;j<cnt&&su[j]*i<=maxn;j++)
            is[su[j]*i]=0;
    }

}
void getprimefactor(long long nn)
{
    cas=0;
    memset(factor,0,sizeof(factor));
    for(int i=1;i<cnt&&su[i]*su[i]<=nn;i++)
    {
        while(nn%su[i]==0)
        {
            factor[cas]++;
            nn/=su[i];
        }
        if(factor[cas])
            cas++;
    }
    if(nn>1)
    {
        factor[cas]=1;
        cas++;
    }
}
int gcd(int a,int b)
{
    return b?gcd(b,a%b):a;
}
int main()
{
    prime();
    int t,ans1;
    long long n;
    scanf("%d",&t);
    ans1=1;
    while(t--)
    {
        int  ans=0;
        scanf("%lld",&n);
        getprimefactor(fabs(n));
        //int x=0;
        for(int i=0;i<cas;i++)
        {
            ans=gcd(ans,factor[i]);
        }
        if(n<0)
        {
            while(ans%2==0)
            ans=ans/2;
        }
        printf("Case %d: %d\n",ans1,ans);
        ans1++;
    }
    return 0;
}