题意: 在一个区间[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;
}