poj 1730 Perfect Pth Powers

题目：

对于整数x，求满足p^q = x条件的q的最大值

数论题

对整数x分解得：x = a1^b1 * a2^b2 *...*an^bn 其中a1、a2、、、an为素数

由x = p^q，则x的素因子和p的素因子相同

则有：a1^b1 * a2^b2 *...*an^bn = p^q

易得，q的最大值即各个素因子个数的最大公约数

本题有负整数，对于负整数如果最大公约数为偶数，需要额外处理

例如：-64 = (-4)^3 而不是 -64 = (-2)^6

设最大公约数为m，则这时所有素因子个数的奇数因子一定是m的因子，如果不是，假如有奇数因子y

则m*y也为公因子，这与m是最大公因子矛盾

所以，可以通过对m不断除以2直至得到奇数，即为解

// POJ 1730
/*
 * re wa ac 324K 16MS
 */

#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;

#define MAX 50000

bool prime[MAX + 10];
int x;
int len;
int p[MAX + 10];
int num[MAX + 10];

void create() {
    int i, j;

    prime[2] = true;
    len = 0;
    for(i = 3; i < MAX; i += 2)
        prime[i] = true;
    for(i = 2; i <= MAX / 2; ++i) {
        if(prime[i])
            for(j = 1; i + i * j < MAX; ++j)
                prime[i + i * j] = false;
    }
    for(i = 2; i < MAX; ++i)
        if(prime[i])
            p[len++] = i;
}

int gcd(int a, int b) {
    return (b == 0) ? a : gcd(b, a % b);
}

int proc() {
    int i;
    int tem_len = 0;
    int num_len = 0;
    bool sign = false;
    if(x < 0) {
        x = -x;
        sign = true;
    }
    while(x != 1) {
        int tem_num = 0;
        if(tem_len == len)
            break;
        while(x % p[tem_len] == 0) {
            x /= p[tem_len];
            ++tem_num;
        }
        if(tem_num > 0)
            num[num_len++] = tem_num;
        ++tem_len;
    }
    if(num_len == 0)
        return 1;
    int a = num[0];
    for(i = 1; i < num_len; ++i)
        if(num[i] > a)
            a = gcd(num[i], a);
        else
            a = gcd(a, num[i]);
    if(sign)
        while(!(a & 1))
            a /= 2;
    return a;
}

int main()
{
    create();
    while(cin>>x, x) {
        cout<<proc()<<endl;
    }
    return 0;
}