题目:
对于整数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;
}