刚开始想直接用库函数直接水过,结果发现居然还会有负数。然后就想到因式分解,统计每个质因子的数量,求出这些数量的gcd即可。当n为负数时,结果不能为偶,需要把结果中的因子2除净。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;
const int M=1000010;
int n,c;
int vis[M];
int prime[M];
void init()
{
memset(vis,0,sizeof(vis));
for(int i=2;i<=1000;i++) if(!vis[i])
{
for(int j=i*i;j<=1000000;j+=i) vis[j]=1;
}
c=0;
for(int i=2;i<=1000000;i++) if(!vis[i])
prime[c++]=i;
}
int gcd(int a,int b)
{
return b==0?a:gcd(b,a%b);
}
int main()
{
freopen("in.txt","r",stdin);
init();
while(cin>>n)
{
if(!n) break;
int tmp=n;
vector<int> v;
for(int i=0;i<c;i++)
{
int t=0;
while(n%prime[i]==0)
{
n/=prime[i];
t++;
}
if(t>0) v.push_back(t);
}
if(n!=1&&n!=-1) v.push_back(1);
int ans=v[0];
for(int i=1;i<v.size();i++)
ans=gcd(ans,v[i]);
if(tmp<0)
{
for(int i=0;i<v.size();i++)
while(v[i]%2==0) v[i]/=2;
ans=v[0];
for(int i=1;i<v.size();i++)
ans=gcd(ans,v[i]);
}
cout<<ans<<endl;
}
return 0;
}