/*

只要求[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;

}