lightoj 1220 唯一分解定理

#include<bits/stdc++.h>
using namespace std;
#define maxn 1000005
#define ll long long

int v[maxn],prime[maxn],m;
void init(int n){
    memset(v,,sizeof v);
    memset(prime,,sizeof prime);
    m=;
    for(int i=;i<=n;i++){
        if(v[i]==){
            v[i]=i;
            prime[++m]=i;
        }
        for(int j=;j<=m;j++){
            if(prime[j]>v[i] || prime[j]*i>n)break;
            v[i*prime[j]]=prime[j];
        }
    }
}
ll p[],c[],mm;
void divide(ll n){
    memset(p,,sizeof p);
    memset(c,,sizeof c);
    mm=;
    for(int i=;i<=m;i++){
        if(prime[i]>n)break;
        if(n%prime[i]==){
            p[++mm]=i,c[mm]=;
            while(n%prime[i]==)n/=prime[i],c[mm]++;
        }
    }
    if(n>)p[++mm]=n,c[mm]++;
}

int main(){
    init();
    int T;
    ll n;
    cin>>T;

    for(int tt=;tt<=T;tt++){
        cin>>n;
        int flag=;
        if(n<)n=-n,flag=;
        divide(n);
        ll Max=;
        for(int i=;i<=mm;i++)
            Max=max(Max,c[i]);

        ll ans=;
        for(int i=;i<=mm;i++)
            ans=__gcd(ans,c[i]);

        if(flag){
            while(ans%==)ans/=;
        }
        printf("Case %d: %d\n",tt,ans);
    }
}