题目:点击打开链接
题意:给定一个区间,求出区间内所有数的最大的指数和。
分析:可以先求出1~n中有的区间指数和,再取端点相减。考虑到最大的指数不超过60(2^60>1e18),因此可以枚举所有指数,根据指数k对n开根,因为用a=pow(n,i)可能会产生不超过1的精度误差,所以可以每次判断a-1,a和a+1哪个值的k次方不超过n且更接近,并且a+1的k次方可能会超过long long表示的最大范围,要特殊处理。对于每个指数k,得到的a的值是预处理出来的小于等于a^k的个数,因为一个指数t1若是能被另一个指数t2整除,那么t2对应的a个数都在t1中出现过,相同的数会重复计算,所以在t1对应的a个数减去t2对应的a个数,所以指数从大到小枚举,会遇到上述情况,类似dp的思想,处理出所有的a后,乘上各自的指数即可。
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <stdio.h>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
#include <set>
#include <vector>
#include <map>
#define PR pair<int,int>
#define MP make_pair
#define fi first
#define se second
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define sqr(x) ((x)*(x))
#define ll long long
const ll INF = 1e18+300;
const int inf=0x3f3f3f3f;
const int M=100010;
const int N=100100;
const ll MOD=1000000007;
const double eps=1e-3;
const double pi=acos(-1.0);
using namespace std;
ll l,r;
ll dp[63];
ll qpow(ll x,int n)
{
ll ans=1;
while(n)
{
if(n&1)
{
double tmp=1.0*INF/ans;
if(x>tmp) return -1;
ans*=x;
}
n>>=1;
if(x>((ll)1<<31)&&n) return -1;
x=x*x;
}
return ans;
}
ll change(ll x,ll n)
{
ll a=(ll)pow(x,1.0/n);
ll qt=qpow(a,n);
ll ql=qpow(a-1,n);
ll qr=qpow(a+1,n);
if(qr!=-1&&qr<=x) return a+1;
else if(qt!=-1&&qt<=x) return a;
else if(ql!=-1&&ql<=x) return a-1;
}
ll calc(ll n)
{
if(n==1) return 1;
int i,j;
memset(dp,0,sizeof(dp));
dp[1]=n;
for(i=2;i<63;i++)
{
dp[i]=change(n,i)-1;
if(dp[i]==0) break;
}
int k=i;
for(i=k-1;i>0;i--)
for(j=1;j<i;j++)
if(i%j==0) dp[j]-=dp[i];
ll ans=dp[1];
for(int i=2;i<k;i++)
ans+=(i*dp[i]);
return ans;
}
int main()
{
while(~scanf("%I64d%I64d",&l,&r))
{
if(l==0&&r==0) break;
printf("%I64d\n",calc(r)-calc(l-1));
}
}