题意:
x==B^p,对于每个给定的x找到最大的p。
思路:
注意x是所有32位有符号整数,所以用long long存安全,并且有负数。对于每个x可以根据唯一分解定理分为不同素因子ai,B一定是不同素因子ai的乘积,每个B都包含一定数量的每个ai,数量就是所有bi的最大公约数,正好平均把所有素因子分成几分。所以只要求所有素因子数量的gcd就行了,注意500ms所以要姿势正确。
然后就是负数要保证每个B中含有的,每样素因子都一定是奇数个,否则会变成正数,所以特判下负数把所有bi变成奇数。
代码:
#define N 112345
long long n,m;
int flag,sum,ave,ans,res,len,ans1,ans2;
int gcd(int x,int y)
{
if(x<y)x^=y,y^=x,x^=y;
int t;
while(y)
t=y,y=x%y,x=t;
return x;
}
bool mark[N];
long long pri[N],cnt;
void SP()
{
cnt=0;
memset(mark,true,sizeof(mark));
mark[0]=mark[1]=false;
for(int i=2;i<N;i++)
{
if(mark[i])
pri[cnt++]=i;
for (int j=0;(j<cnt)&&(i*pri[j]<N);j++)
{
mark[i*pri[j]]=false;
if (i%pri[j]==0)
break;
}
}
}
int getFactor(long long x)
{
int sum=0;
for(int i=0;i<cnt&&pri[i]*pri[i]<=x;i++)
{
int res=0;
while(x%pri[i]==0&&x)x/=pri[i],res++;
if(flag)while(res%2==0&&res)res/=2;
if(res&&!sum)sum=res;
else if(res)sum=gcd(sum,res);
}
if(x>1)sum=1;
return sum;
}
int main()
{
int i,j,k,kk,t,x,y,z;
SP();
scanf("%d",&k);
kk=0;
while(k--)
{
scanf("%lld",&n);
flag=0;
if(n<0)n=-n,flag=1;
printf("Case %d: %d\n",++kk,getFactor(n));
}
return 0;
}