HDU 3208 Integer’s Power 【容斥】

题目：点击打开链接

题意：给定一个区间，求出区间内所有数的最大的指数和。

分析：可以先求出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));
    }
}