#include <iostream>

using namespace std;

const int MOD = 1000000007;

int main(){
	long long N;
	while(cin >> N){
		long long dp[2][3];
		for(int i=0;i<3;i++) dp[0][i] = dp[1][i] = 0;
		dp[0][0] = 1;
		for(int i=60;i>=0;i--){
			long long* cur = dp[i%2];
			long long* next = dp[1-i%2];
			for(int j=0;j<3;j++) next[j] = 0;
			for(int j=0;j<3;j++){
				int r = 2*j + (N>>i)%2;
				for(int k=0;k<=2;k++){
					if(r-k < 0) break;
					next[min(2, r-k)] += cur[j];
				}
			}
			for(int j=0;j<3;j++) next[j] %= MOD;
		}
		int res = 0;
		for(int j=0;j<3;j++) res = (res + dp[1][j])%MOD;
		cout << res << endl;
	}
}
