#include #include #include #include using namespace std; using i64 = int64_t; constexpr int mod = static_cast(powl(10, 9)) + 7; string to_bin(const i64 a) { if(a == 0) { return "0"; } string s = bitset<64>(a).to_string(); int p = static_cast(s.find('1')); return s.substr(p); } int cmp(int x, int y) { if(x < y) { return 0; } if(x == y) { return 1; } return 2; } int mycmp(int x, int y, int val) { int res = cmp(x, y); if(res == 1) { res = val; } return res; } i64 f(i64 x) { string s = to_bin(x); int n = static_cast(s.size()); // dp[yがNより未満/丁度/超過][xがyより未満/丁度/超過][ANDがXORより未満/丁度/超過][XORがORより未満/丁度/超過] := パターン数 vector>>> dp(3, vector>>(3, vector>(3, vector(3, 0)))); dp[1][1][1][1] = 1; for(int i=n-1; i>=0; --i) { vector>>> ndp(3, vector>>(3, vector>(3, vector(3, 0)))); for(int flag0=0; flag0<3; ++flag0) { for(int flag1=0; flag1<3; ++flag1) { for(int and_xor=0; and_xor<3; ++and_xor) { for(int xor_or=0; xor_or<3; ++xor_or) { if(!dp[flag0][flag1][and_xor][xor_or]) { continue; } for(int yi=0; yi<2; ++yi) { int nflag0 = mycmp(yi, s[i]-'0', flag0); for(int xi=0; xi<2; ++xi) { int nflag1 = mycmp(xi, yi, flag1), n_and_xor = mycmp(xi & yi, xi ^ yi, and_xor), n_xor_or = mycmp(xi ^ yi, xi | yi, xor_or); ndp[nflag0][nflag1][n_and_xor][n_xor_or] += dp[flag0][flag1][and_xor][xor_or]; ndp[nflag0][nflag1][n_and_xor][n_xor_or] %= mod; } } } } } } dp = ndp; } i64 res = 0; for(int flag0=0; flag0<2; ++flag0) { for(int flag1=0; flag1<2; ++flag1) { res += dp[flag0][flag1][0][0]; res %= mod; } } return res; } int main(void) { i64 x; scanf("%ld", &x); i64 res = f(x); printf("%ld\n", res); return 0; }