#include <iostream>
using namespace std;
typedef long long ll;

const ll M=60,mod=1e9+7;

ll N;

ll res=0,dp[4],b=0;

int main(){
	cin>>N;
	for(int i=M;i>=0;i--){
		ll DP[4];
		if(N&1ll<<i){
			DP[0]=(dp[0]+!b)%mod;
			DP[1]=(2*dp[0]+3*dp[1]+b)%mod;
			DP[2]=(dp[0]+2*dp[2])%mod;
			DP[3]=(dp[1]+2*dp[2]+4*dp[3])%mod;
			b=1;
		}
		else{
			DP[0]=2*dp[0]%mod;
			DP[1]=(3*dp[1]+b)%mod;
			DP[2]=2*dp[2]%mod;
			DP[3]=(dp[1]+4*dp[3])%mod;
		}
		for(int j=0;j<4;j++) dp[j]=DP[j];
	}
	cout<<(dp[2]+dp[3])%mod<<endl;
}