题意:就是给你一个等式x=b^p,现在给你x,b未知,求p的最大值;
这道题我一读,没感觉呢?但是应该想得到唯一分解定理:
p1,p2…pn是素数,这个可以打表解决;
但是怎么求p最大呢?
可以多写几个例子比如:
其实就可以大概推测出来p=gcd(x1,x2,x3…xn);
所以很明显素数表扫一遍就OK;但是注意了:
x可以是负数,所以它的p只能是奇数次幂;还有一点需要注意一个数论结论:如果一个数能被1e5以内的素数完全分解为1,那么这个数一定是素数;还有一点需要注意gcd(0,x)==x;因为0除以任何数都为0,还有就是注意开ll;
注意这些就OK了;
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
//素数筛选
const int N = 1e6 + 10;
int pri[700010], k;//k记录素数的个数这里的素数是从下标为0开始的
bool Isprime[N];//这个除了辅助筛素数,也可以用来判断是不是素数if(Isprime[n])) n为素数;
void prime()//打表
{
k = 0;
memset(Isprime, true, sizeof(Isprime));
Isprime[1] = false;
for(int i = 2 ; i < N ; i++)
{
if(Isprime[i])
{
pri[k++] = i;
for(int j = 2 ; i * j < N ;j++)
Isprime[i * j] = false;
}
}
}
ll gcd(ll a,ll b){
return b?gcd(b,a%b):a;
}
int main(){
prime();
ll T,g=1;
scanf("%lld",&T);
while(T--){
ll x,cnt=0,ans=0;
scanf("%lld",&x);
ll t=x;
if(x<0)x=-x;
for(int i=0;i<k&&pri[i]*pri[i]<=x;i++){//分解x
cnt=0;
if(x%pri[i]==0){
while(x%pri[i]==0){
cnt++;
x/=pri[i];
}
ans=gcd(ans,cnt);
}
}
if(x>1) ans=1;//因为不能完全分解所以pmax=1;
if(t<0){
while(ans%2==0)ans/=2;//如果是负数就只能是奇数次幂
}
printf("Case %lld: %lld\n",g++,ans);
}
return 0;
}