def v(n): r = 0 k = 1 while k <= n: r += f(n, k) r %= 10**9+7 k *= 2 return r def f(n, k): if n&k == 0: n = (n//k)*k-1 size = (n//k)//2*k + n%k+1 a = sum_range(size) b = sum_range((size//k)+1)*k*k c = (size%k)*k*((size//k)+1) return a + b + c def sum_range(n): return n*(n-1)//2 print(v(int(input())))