zoj 2124 || poj 1730 Perfect Pth Powers

 

数论，求使得 x = b^p的最大的p。

 

这题之前用pow神马的，一直过不了。昨天看到了，有个想法，想试试。

 

先求质因数，质因数每个的个数肯定和p有关。

 

求出之后，想了想，呀，这不就是求这些质因数重复个数的最大公约数嘛！！！激动~！

 

写了 = =。。WA。但是数据都能过。还用了long long也过不了。搜了下，说这题有负的数据，FT！！！

 

改了改。如果是负数的话，质因数重复个数如果是偶数一直除2一直到奇数，因为如果不这么弄的话，求出来的最大公约可能就是2的倍数了，而这个对于负数来说是不可能的。然后再求最大公约数即可。用int就好。

 

#include <cstdio> #include <cstdlib> #include <iostream> #include <string.h> #include <math.h> #include <limits.h> #define MAX 50000 using namespace std; int f[MAX]; int num[MAX]; int pri[MAX]; int cou; void init() { int i,j; cou = 0; memset(num,0,sizeof(num)); for(i=2; i<MAX; i++) for(j=2; j*i<MAX; j++) if( num[i] == 0 ) num[j*i] = 1; for(i=2; i<MAX; i++) if( num[i] == 0 ) pri[cou++] = i; } int gcd(int n,int m) { return m == 0 ? n : gcd(m,n%m); } int Analy(int n) { int fac = 0; int i = 0 ,num; int x = n; while( n != 1 ) { num = 0; if( abs(pri[i]) > sqrt(fabs((double)n)) ) break; if( n % pri[i] == 0 ) { while( n % pri[i] == 0 ) { num++; n /= pri[i]; } f[fac] = num; fac++; } i++; } if( fac == 0 ) return 1; if( x < 0 ) for(i=0; i<fac; i++) while( f[i] % 2 == 0 ) f[i] /= 2; if( fac == 1 ) return f[0]; int tmp = gcd(f[0],f[1]); for(i=2; i<fac; i++) tmp = gcd(f[i],tmp); return tmp; } int main() { init(); int n; int ans; while( scanf("%d",&n) != EOF && n ) { memset(f,0,sizeof(f)); ans = Analy(n); printf("%d/n",ans); } return 0; }