#include <bits/stdc++.h>
using namespace std;
using lint = long long int;
constexpr lint MOD = 1000000007;
lint N;

int main()
{
    cin >> N;
    int digit = 0;
    lint tmp = N;
    while (tmp)
    {
        digit++;
        tmp >>= 1;
    }
    vector<vector<lint>> dp(digit+1, vector<lint>(3));
    dp[0][0] = 1;
    for (int i = 1; i <= digit; i++) {
        if (N >> (digit-i) & 1) {
            dp[i][0] = dp[i - 1][0];
            dp[i][1] = dp[i - 1][0] + dp[i - 1][1];
            dp[i][2] = 2 * dp[i - 1][1] + 3 * dp[i - 1][2];
        }
        else {
            dp[i][0] = dp[i - 1][0] + dp[i - 1][1];
            dp[i][1] = dp[i - 1][1];
            dp[i][2] = dp[i - 1][1] + 3 * dp[i - 1][2];
        }
        for (auto &v : dp[i]) v %= MOD;
    }
    cout << accumulate(dp[digit].begin(), dp[digit].end(), 0LL) % MOD << endl;
}
