#include <iostream>
#include <algorithm>
#include <vector>
#include <iomanip>
#include <map>
#include <queue>

static const int MOD = 1000000007;
using ll = long long;
using u32 = unsigned;
using namespace std;

template<class T>
constexpr T INF = ::numeric_limits<T>::max() / 32 * 15 + 208;

map<ll, ll> m;

ll f(ll n){
    if(m[n] != 0) return m[n];
    if(n%2 == 0) {
        return (m[n] = (f(n/2) + 2*f(n/2-1)) % MOD);
    }
    else return (m[n] = (2*f(n/2) + f(n/2-1)) % MOD);
}

int main() {
    ll n;
    cin >> n;
    m[0] = 1;
    m[1] = 2;
    cout << f(n) << "\n";
    return 0;
}
