【质因数分解 && n = x^p求最大p】LightOJ - 1220 Mysterious Bacteria

#include<bits/stdc++.h>
using namespace std;
#define nn 100005
bool vis[nn];
int prime[nn/10], m;
int a[nn], b[nn], t;
void get_Prime()//素数筛法
{
m = 0;
memset(vis, 0, sizeof(vis));
for(long long i = 2; i <= nn - 5; i++)
{
if(!vis[i])
{
prime[m++] = i;
for(long long j = 2*i; j <= nn - 5; j = j + i)
vis[j] = 1;
}
}
}
void factor(long long n)//质因数分解
{
int mm = sqrt(n);
if(mm * mm > n) mm--;
t = 0;
for(int i = 0; i < m && prime[i] <= mm; i++)
{
if(n%prime[i] == 0)
{
a[t] = prime[i];
b[t] = 0;
while(n%prime[i] == 0)
{
b[t]++;
n /= prime[i];
}
t++;
}
}
if(n != 1){
a[t] = n;
b[t++] = 1;
}
}
int gcd(int x, int y)//求x,y的最大公约数
{
if(!y) return x;
else return gcd(y, x%y);
}
int main()
{
int T, i, Case = 1;
long long n;
scanf("%d", &T);
get_Prime();
while(T--)
{
scanf("%lld", &n);
int flag = 0;
if(n < 0) {//如果是负数标记一下，变正
n = -1 * n;
flag = 1;
}
factor(n);
int ans = b[0];
for(i = 1; i < t; i++)
{
ans = gcd(ans, b[i]);
}
if(flag && ans % 2 == 0) {//n为负数，ans为偶。那么需要变为奇
while(ans%2 == 0)
{
ans >>= 1;
}
printf("Case %d: %d\n", Case++, ans);
}
else
printf("Case %d: %d\n", Case++, ans);
}
}