poj 1730 Perfect Pth Powers

这个有2种方法。

一种是通过枚举p的值（p的范围是从1-32），这样不会超时，再就是注意下精度用1e-8就可以了，还有要注意负数的处理……

 #include<iostream>
 #include<stdio.h>
 #include<algorithm>
 #include<cmath>
 #include<iomanip>
 #include<string>
 using namespace std;
 int fun(double n)
 {
     int p=,f=;
     if(n<)
     {
         p=;
         f=;
         n=-n;
     }
     for(;p>=;p-=f)
     {
         double a=pow(n,1.0/p);
         long t1=a;
         long t2=a+;
         if(a-t1<=1e-||t2-a<=1e-)
             break;
     }
     return p;
 }
 int main()
 {
     int x,p;
     while(cin>>x&&x)
     {
         p=fun(x);
         cout<<p<<endl;
     }
     return ;
 }

第二种方法就是分解素数因子，求出因子的指数的最大公约数就是答案了，但是我这测试没错，提交超时……

 #include<iostream>
 #include<stdio.h>
 #include<algorithm>
 #include<vector>
 #include<cmath>
 #include<iomanip>
 #include<string>
 using namespace std;
 int prime[],a[];
 void init()
 {
     int i,j,k=;
     for(i=;i<=;i++)
     {
         if(!a[i])
         {
             prime[k++]=i;
             for(j=i+i;j<=;j+=i)
                 a[j]=;
         }
     }
 }
 int gcd(int a,int b)
 {
     int c;
     if(a<b) swap(a,b);
     while(b)
     {
         c=a;
         a=b;
         b=c%b;
     }
     return a;
 }
 int main()
 {
     int n,i,j;
     bool flag;
     init();
     while(scanf("%d",&n),n)
     {
         flag=;
         if(n<)
         {
             flag=;
             n=-n;
         }
         vector<int> key;
         for(i=;prime[i]*prime[i]<=n;i++)
         {
             if(n%prime[i]==)
             {
                 int c=;
                 for(;n%prime[i]==;n/=prime[i],c++);
                 key.push_back(c);
             }
         }
         if(n>) key.push_back();
         int sum=key[];
         for(i=;i<key.size();i++)
             sum=gcd(sum,key[i]);
         if(flag)
         {
             while(sum%==)
                 sum/=;
         }
         printf("%d\n",sum);
     }
     return ;
 }