hdu 3208 容斥原理+高精度

题意： 在一个区间[l,r]，其中每一个数必定可以写成写成a^p的形式，现在设p的和为res，问你在[l,r]区间内最大的res为多少。

思路：
其实和hdu2204很相似，都是枚举幂次，2^63>1e18，所以我们只用枚举62以内的幂次有多少数即可，从大到小枚举，然后找到能整除他的幂次，从中删去即可，但这里pow精度不够，需要手动在p-1,p,p+1中用快速幂处理选取，越界的时候记得返回-1

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll sum[63];
const double eps=1e-18;
const ll inf = (ll)(1e18) + 500;
const ll INF = (ll)1 << 31;
 
ll mypow(ll a,ll b)
{
    ll ans=1;
    while(b){
        if(b&1){
            double judge=1.0*inf/ans;
            if(a>judge) return -1;
            ans*=a;
        }
        b>>=1;
        if(a>INF&&b>0) return -1;
        a=a*a;
    }
    return ans;
}
ll Find(ll n,ll k)//手动扩大pow精度
{
    ll a=1ll*pow(n,1.0/k);
    ll q1=mypow(a,k);
    ll q2=mypow(a-1,k);
    ll q3=mypow(a+1,k);
    if(q3!=-1&&q3<=n)return a+1;
    else if(q1!=-1&&q1<=n)return a;
    else if(q2!=-1&&q2<=n)return a-1;
}
 
ll solve(ll n){
    sum[1]=n;
    int maxx=0;
    for(int i=2;i<=62;i++)//枚举63位幂
    { 
        sum[i]=Find(n,1ll*i)-1;
        if(!sum[i]){
            maxx=i;//记录达到的最大的幂是多少
            break;
        }
    }
    for(int i=maxx;i>=2;i--)//从大到小枚举容斥
        for(int j=1;j<i;j++)
            if(i%j==0) sum[j]-=sum[i];
    ll res=0;
    for(int i=1;i<=maxx;i++)
        res+=sum[i]*i;
    return res;
}
int main()
{
    ll a,b;
    while(~scanf("%lld%lld",&a,&b)&&a)
        printf("%lld\n",solve(b)-solve(a-1));
    return 0;
}