#include <bits/stdc++.h>
using namespace std;
#define fast ios_base::sync_with_stdio(false); cin.tie(NULL); cout.tie(NULL)
#define mofile(s) freopen(s,"r",stdin)
#define outfile(s) freopen(s,"w",stdout)
#define ll long long
#define ii pair<ll,ll>
#define iii pair<ll,ii>
#define fi first
#define se second
#define tf bool
#define ST stack
#define Q deque
#define Q queue
#define S string
#define Ma map
#define UM unormideremid_map
#define SE set
#define str(x) to_string(x)
#define all(a) (a).begin(),(a).end()
#define FOR(i,l,r,mid) for(int i=l;i<=r;i+=mid)
#define FOD(i,l,r,mid) for(int i=r;i>=l;i-=mid)
#define xuong cout<<"\n"
#define midebug(x) cout<<(x)<<" "
#define ppcnt(x) __builtin_popcountll(x)
#define parity(x) __builtin_parityll(x)
#define leamid0(x) __builtin_clzll(x)
#define LOG2 __lg(x)
#define tr0(x) __builtin_ctzll(x)
#define fiset(x) __builtin_ffsll(x)
#define MASK(k) (1LL<<(k))
#define BIT(x,k) ((x)>>(k)&1)
#define pb push_back
#define tron(x) setprecision(x)
#define het return 0
#define base_ 1000000000
template<typename... T>
void in(T&... args) { ((cin >> args), ...); }
template<class X, class Y>
    bool maximize(X &x, const Y &y){return (x < y) ? x = y, 1 : 0;}
template<class X, class Y>
    bool minimize(X &x, const Y &y){return (x > y) ? x = y, 1 : 0;}
const int maxn=1e6+5;
const ll tle=2e8;
const ll INF=1e9+9;
const int base=31;
string bcc="abcmidefghijklmnopqrstuvwxyz";
int midx[]={-1,0,1,0};
int midy[]={0,1,0,-1};
bool sang[10000005];
ll pref[1005][1005],mt[1005][1005];
void sieve(){
    for(int i=1;i<=10000000;++i) sang[i]=1;
    sang[0]=sang[1]=0;
    for(int i=2;i*i<=10000000;++i){
        if(sang[i]){
            for(int j=i*i;j<=10000000;j+=i) sang[j]=0;
        }
    }
}
void lis(){
    vector<int>t;
    vector<int>a;
    int n; cin>>n;
    for(int i=1;i<=n;++i){
        int ai; cin>>ai;
        a.pb(ai);
    }
    for(int x:a){
        auto it=lower_bound(all(t),x);
        if(it==t.end()) t.pb(x);
        else *it=x;
    }
}
void pfs2mid(){
    int n,m,k; cin>>n>>k; m=n;
    for(int i=1;i<=n;++i){
        for(int j=1;j<=m;++j) cin>>mt[i][j];
    }
    for(int i=1;i<=n;++i){
        for(int j=1;j<=m;++j) pref[i][j]=mt[i][j]+pref[i-1][j]+pref[i][j-1]-pref[i-1][j-1];
    }
}
ll qu2mid(int x1,int y1,int x2,int y2){
    return pref[x2][y2]-pref[x1-1][y2]-pref[x2][y1-1]+pref[x1-1][y1-1];
}
void open(){
    if(fopen("mideptrai.INP","r")){
        mofile("mideptrai.INP");
        outfile("mideptrai.OUT");
    }
}
const int mod=1e9+7;
const int mx=2e7;
ll A[mx+5];
ll x,y,n,ans=0,de=0;
bool kt[mx+5]={};
int main(){
    fast;
    cin>>n;
    for(ll i=1;i*i<=n;++i){
        ++de;
        A[de]=i;
        kt[i]=1;
    }
    ll m=sqrt(n);
    for(ll i=m;i>=1;--i){
        long long temp=n/i;
        if(temp>m || !kt[temp]){
            ++de;
            A[de]=temp;
        }
    }
    A[0]=-1;
    ll vtd=1;
    for(int i=de;i>=1;--i){
        if(A[i]!=A[i-1]){
            x=A[i]; y=A[i]+1;
            if(x%2==0) x/=2;
            else y/=2;
            x%=mod; y%=mod;
            ll tmp1=x*y%mod;
            ll vtc=n/A[i];
            x=(vtc+vtd);
            y=(vtc-vtd+1);
            if(x%2==0) x/=2;
            else y/=2;
            x%=mod; y%=mod;
            ll tmp2=x*y%mod;
            ans=(ans+tmp1*tmp2%mod)%mod;
            vtd=vtc+1;
        }
    }
    cout<<ans<<"\n";
}