#include using namespace std; using ll = long long; constexpr char newl = '\n'; // https://noshi91.hatenablog.com/entry/2019/03/31/174006 template struct ModInt { using u64 = std::uint_fast64_t; static constexpr u64 MOD = Modulus; u64 val; constexpr ModInt(const u64 x = 0) noexcept : val(x % MOD) {} constexpr ModInt operator+() const noexcept { return ModInt(*this); } constexpr ModInt operator-() const noexcept { ModInt res(*this); if (res.val != 0) res.val = MOD - res.val; return res; } constexpr bool operator==(const ModInt& rhs) const noexcept { return val == rhs.val; } constexpr bool operator!=(const ModInt& rhs) const noexcept { return val != rhs.val; } // prefix increment/decrement constexpr ModInt& operator++() noexcept { return *this += ModInt(1); } constexpr ModInt& operator--() noexcept { return *this -= ModInt(1); } // postfix increment/decrement constexpr ModInt operator++(int) noexcept { ModInt tmp(*this); ++*this; return tmp; } constexpr ModInt operator--(int) noexcept { ModInt tmp(*this); --*this; return tmp; } constexpr ModInt operator+(const ModInt& rhs) const noexcept { return ModInt(*this) += rhs; } constexpr ModInt operator-(const ModInt& rhs) const noexcept { return ModInt(*this) -= rhs; } constexpr ModInt operator*(const ModInt& rhs) const noexcept { return ModInt(*this) *= rhs; } constexpr ModInt operator/(const ModInt& rhs) const noexcept { return ModInt(*this) /= rhs; } constexpr ModInt& operator+=(const ModInt& rhs) noexcept { val += rhs.val; if (val >= MOD) val -= MOD; return *this; } constexpr ModInt& operator-=(const ModInt& rhs) noexcept { if (val < rhs.val) val += MOD; val -= rhs.val; return *this; } constexpr ModInt& operator*=(const ModInt& rhs) noexcept { val = val * rhs.val % MOD; return *this; } // prime Modulus only constexpr ModInt& operator/=(const ModInt& rhs) noexcept { return *this *= rhs.inv(); } // prime Modulus only constexpr ModInt inv() const noexcept { return pow(*this, MOD - 2); } }; template constexpr ModInt pow(ModInt x, std::uint_fast64_t n) { ModInt res(1); while (n) { if (n & 1) res *= x; x *= x; n >>= 1; } return res; } template istream& operator>>(istream& is, ModInt& x) { std::uint_fast64_t val; is >> val; x = ModInt(val); return is; } template ostream& operator<<(ostream& os, const ModInt& x) { return os << x.val; } using mint = ModInt<1000000007>; int main() { cin.tie(nullptr); ios::sync_with_stdio(false); ll n; cin >> n; mint ans = 0; for (ll i = 1; i * i <= n; i++) { if (i * i < n) { // i * x + j == n (0 <= j < x, i < x) ll num_x = (n / i) - max((n / (i + 1)), i); mint xsum = mint(n / i * 2 - num_x + 1) * mint(num_x) / mint(2); mint jsum = mint(n) * mint(num_x) - xsum * i; ans += jsum + mint(num_x) * i; } if (i == 1) continue; ll nn = n; while (nn > 0) { ans += nn % i; nn /= i; } } cout << ans << newl; return 0; }