#include using namespace std; typedef long long int ll; const int N=17705360; bool p[N]; int mob[N]; void mobius(int n) { int i,j; for(i=0;i<=n;i++) {p[i]=false; mob[i]=1;} for(i=2;i<=n;i++) if(!p[i]) { mob[i]=-1; for(j=2;j<=n/i;j++) { p[i*j]=true; mob[i*j]*=-1; } ll sq=1LL*i*i; for(j=1;j<=n/sq;j++) mob[j*i*i]=0; } } int maxpower(ll n,int p) { if(n==0) return 1; int l=1,r=(int)1e9,m,ans,i; if(p==5) r=3190; ans=r; while(l<=r) { m=(l+r)>>1; ll v=1; for(i=0;i mxi; int sum=0,i; for(i=imax-1;i>0;i--) { int cur=1; int root=maxpower(n/i,2),sqd=maxpower(root,2),j; for(j=1;j<=root/(sqd+1);j++) cur-=1LL*(root/j-root/(j+1))*mob[j]; for(j=2;j<=sqd;j++) if(root/j<=d) cur-=mob[root/j]; else cur-=mxi[imax-j*j*i-1]; mxi.push_back(cur); sum+=cur; } return sum-(imax-1)*mob[d]; } int main() { ll a,b; cin>>a>>b; a--; int imax=maxpower(b,5),ima=maxpower(a,5); int d=maxpower(b/imax,2),d2=maxpower(a/ima,2); int dm=max(d,d2); mobius(dm); ll s1=0,s2=0; int i; mob[0]=0; for(i=1;i<=dm;i++) { if(i<=d) s1+=(b/(1LL*i*i))*mob[i]; if(i<=d2) s2+=(a/(1LL*i*i))*mob[i]; mob[i]+=mob[i-1]; } s1+=squarefree(b,imax,d); s2+=squarefree(a,ima,d2); cout<<(s1-s2)<<"\n"; return 0; }