#pragma GCC optimize("Ofast") #pragma GCC target("avx2") char*mmap(); #define rd(v) long v=0;{int _c;while(_c=*rp++-48,_c>=0)v=v*10+_c;} #define wt(v) {long _z=v;do*--wp=_z%10+48;while(_z/=10);} #define MD 1000000007 #define LIM 100000 #define A1 1, 2, 3, 5, 6,10, 1, 5, 3, 7, 1,13, #define A2 1, 2, 8, 5, 1,10, 1,10, 3, 2, 1,13, #define A3 6, 2, 3, 5, 1,15, 1, 5, 3, 2, 6,13, #define A4 1, 2, 3,10, 1,10, 1, 5, 8, 2, 1,13, #define A5 1, 7, 3, 5, 1,10, 6, 5, 3, 2, 1,18, #define A A1 A2 A3 A4 A5 #define B A A A A A A A A A A #define C B B B B B B B B B B #define D C C C C C C C C C C int dp[]={0,D}; main(){ char buf[64]; char*rp=buf; read(0,rp,sizeof buf); rd(m); long*p=dp+1; long r=0; r+=1*(m/1-LIM/2/1); r+=1*(m/2-LIM/2/2); r+=2*(m/3-LIM/2/3); r+=3*(m/4-LIM/2/4); r+=5*(m/5-LIM/2/5); r+=6*(m/6-LIM/2/6); for(long i=7;i<=LIM/2;++i){ long a=dp[i-1]; for(long j=i;j<=LIM/2;j+=i){ long b=a+dp[j]; dp[j]=b+(b>=MD?-MD:0); } r+=a*((m+0)/i-(LIM/2+0)/i); } if(m<=LIM/2){ r=dp[m]; }else{ r=(r+dp[LIM/2])%MD; } char*wp=buf+sizeof buf; wt(r); write(1,wp,buf+sizeof buf-wp); _exit(0); }