#include<bits/stdc++.h>
using namespace std;
const int N=100005;
long long l,r,a[N];
int p[N], cnt;
bool flag[N];
void phi(){
    for(int i=2; i<=N; i++){
        if(!flag[i]) p[++cnt] = i;
        for(int j=1; j<=cnt; j++){
            if(i * p[j] > N)
			break;
            flag[i * p[j]] = 1;
            if(i % p[j] == 0) 
			break;
        }
    }
}
int main(){
    phi();
    cin >> l >> r;
    for(long long i=l; i<=r; i++)
        a[i - l + 1] = i; 

    for(int i=1; i<=cnt; i++){
        long long j = l / p[i];
        while(p[i] * j <= r){
            if(p[i] * j >= l)
            a[p[i]*j-l+1]=a[p[i]*j-l+1]/p[i]*(p[i]-1);
            j++; 
        }
    }
    long long ans=0;
    for(long long i=l; i<=r; i++){
        ans+=i-a[i-l+1];
        ans%=666623333;
    }
    cout<<ans<<"\n";
}