容斥原理——hdu3208

和hdu2204有点像

这题要特别注意精度问题，如pow的精度需要自己搞一下，然后最大的longlong可以设为1<<31

/*
只要求[1,n]范围内的sum即可
那么先枚举幂次k[1,63]
然后求 x^k<=n，x的最大取值sum[i]（这里要扩大一下pow的精度）
然后进行容斥，j%i==0,si-=sj
*/
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define INF (ll)1<<31
#define esp 1e-18
ll sum[];
const double inf=1e18+;
/*
ll Pow(ll a,ll b){//防越界的快速幂
    ll res=1;
    while(b){
        if(b%2){
            double judge=1.0*inf/res;
            if(a>judge)return -1;//必须先判断res*a是否会越界
            res=res*a;
        }
        b>>=1;
        if(a>INF && b>0)return -1;//判断a是否越界了
        a=a*a;
    }
    return res;
} 

ll calc(ll k,ll n){
    ll x=(ll)pow(n,1.0/k),t,p;
    p=Pow(x,k);
    if(p==n)return x;//刚好相等
    if(p>n || p==-1) x--;//越界了
    else {//判断x是否能再+1
        t=Pow(x+1,k);
        if(t!=-1 && t<=n)x++;
    }
    return x;
} */

ll multi(ll a,ll b)//快速乘
{
    ll ans=;
    while(b){
        if(b&){
            double judge=1.0*inf/ans;
            if(a>judge) return -;
            ans*=a;
        }
        b>>=;
        if(a>INF&&b>) return -;
        a=a*a;
    }
    return ans;
}
ll Find(ll x,ll k)//手动扩大pow精度
{
    ll r=(ll)pow(x,1.0/k);
    ll t,p;
    p=multi(r,k);
    if(p==x) return r;
    if(p>x||p==-) r--;
    else
    {
        t=multi(r+,k);
        if(t!=-&&t<=x) r++;
    }
    return r;
}
ll solve(ll n){
    if(n==)return ;
    ll res=,Max=;
    memset(sum,,sizeof sum);
    sum[]=n;
    for(ll k=;k<=;k++){
        sum[k]=Find(n,k);//求扩精度后的最大的x
        if(sum[k])sum[k]--;//1特判
        if(sum[k]==){Max=k;break;}
    }
    //容斥
    for(int j=Max;j>=;j--)
        for(int i=;i<j;i++)
            if(j%i==)sum[i]-=sum[j];
    for(int i=;i<=Max;i++)res+=sum[i]*i;
    return res;
}

int main(){
    ll L,R;
    while(cin>>L>>R && R)
        cout<<solve(R)-solve(L-)<<endl;
}