import sys MOD = 10**9 + 7 def main(): M = int(sys.stdin.readline()) if M == 0: print(0) return # Maximum possible value of s is M, so s-1 can be up to M-1 max_n = M - 1 # Precompute divisors for all numbers up to max_n (M-1) divisors = [[] for _ in range(max_n + 1)] for d in range(1, max_n + 1): for multiple in range(d, max_n + 1, d): divisors[multiple].append(d) # Compute f(s) for all s from 1 to M max_s = M f = [0] * (max_s + 1) f[1] = 1 # Base case for s in range(2, max_s + 1): n = s - 1 total = 0 for d in divisors[n]: s_new = n // d if s_new >= 1: total += f[s_new] f[s] = total % MOD # Find all divisors of M M_divisors = set() for d in range(1, int(M**0.5) + 1): if M % d == 0: M_divisors.add(d) M_divisors.add(M // d) # Sum f(M/d) for all divisors d of M ans = 0 for d in M_divisors: s = M // d ans = (ans + f[s]) % MOD print(ans % MOD) if __name__ == '__main__': main()