#include using namespace std; typedef signed long long ll; #undef _P #define _P(...) (void)printf(__VA_ARGS__) #define FOR(x,to) for(x=0;x<(to);x++) #define FORR(x,arr) for(auto& x:arr) #define FORR2(x,y,arr) for(auto& [x,y]:arr) #define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++) #define ALL(a) (a.begin()),(a.end()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) //------------------------------------------------------- int N; const ll mo=1000000007; ll slow(int v) { ll ret=0; int i; for(i=2;i<=v;i++) { ll p=v/i; if(p<1) break; ret+=p-1; } return ret%mo; } ll hoge(int v) { if(v<=100) return slow(v); ll ret=0; int sq=sqrt(v)+2; map M; // M[x]= |N/i|がxになる数 int i; for(i=2;i<=sq;i++) ret+=v/i-1; ret%=mo; for(i=1;i<=sq;i++) { int L=v/i; int R=max(sq+1,v/(i+1)+1); if(L>=R) (ret+=1LL*(i-1)*(L-R+1))%=mo; else break; } return ret%mo; } void solve() { int i,j,k,l,r,x,y; string s; //for(i=1;i<=20;i++) cout<>N; ll ret=0; map M; // M[x]= |N/i|がxになる数 for(i=1;i<=min(20000,N);i++) M[N/i]++; for(i=1;i<=20000;i++) { int L=N/i; int R=max(20001,N/(i+1)+1); if(L>=R) M[i]+=L-R+1; } // ? x aX FORR(m,M) { //cout<1) (ret+=(N-2)*b%mo*m.second)%=mo; // ax x bx (a>=1, b>1, a!=b)を引く ll a=(m.first-2); (ret-=a*b%mo*m.second)%=mo; } // X aX abX (1