LightOJ 1220

设 n = p^k 求 k 的最大值

 

这道题就很有意思了 首先 n 可正可负，p^k 我们根据因子和公式可知道

n =  那么 我们现在要的是 把所有素数 pi 合成一个数字 从而使得 ei 合成一个最大值

举个例子 2^6 * 3^9 我们相合成一个数字，最好的方法就是 变成 4^3 * 9^3 这样k 就可以直接的出来了

那么我们稍微观察就发现，把两个合并的过程就是把系数 ei 除以一个数字使得 e1 和 e2 相等从而合并，

那么我们发现不就是在求 gcd( e1, e2 )了吗，那么我们推一下就会发现

k = gcd( e1, e2, e3, e4....en)

k解决了，但是我们别忘了负数，处理负数其实很简单，先把他变成正数，但是如果最终结果 k 是偶数 很显然答案不对，如果k是偶数我们则需要不断对 k 进行 /2 操作，必须是奇数时候才能答案p^k才能变成负数。

 

以下为 AC 代码

 

 

 

 

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 1e5+5;
ll prime[maxn], k;
bool vis[maxn];
void Prime()
{
    memset(vis, true, sizeof vis);
    k = 0;
    for(int i=2; i<maxn; i++)
    {
        if(vis[i])
        {
            prime[k ++] = i;
            for(int j=i*2; j<maxn; j+=i)
            {
                vis[j] = false;
            }
        }
    }
}
 
ll gcd(ll a, ll b)
{
    return a % b == 0 ? b : gcd(b, a % b);
}
 
ll solve(ll n)
{
    ll ans, sum;
    ans = 0;
    for(int i=0; prime[i]*prime[i] <= n; i++)
    {
        if(n%prime[i] == 0)
        {
            ll k = 0;
            while(n%prime[i] == 0)
            {
                n /= prime[i];
                k ++;
            }
            if(ans == 0) ans = k;
            else ans = gcd(ans, k);
        }
    }
    if(n>1)//如果 n 本身无法被素数分解干净，那就只能是他自己作为因子了 
        ans = 1;
    return ans;
}
 
int main()
{
    int t, cas=1;
    Prime();
    scanf("%d", &t);
    while(t --)
    {
        ll n;
        scanf("%lld", &n);
        int b = 0;
        if(n<0)
        {
            n = -n;
            b = 1;
        }
        ll sum = solve(n);
        if(b == 1 && sum%2 == 0)
        {
            while(sum%2 == 0)
                sum/=2;
        }
        printf("Case %d: %lld\n", cas++, sum);
    }
    return 0;
}