#include <bits/stdc++.h>
#pragma GCC optimize("O3")
long mod(long a) {
  long long c = a % (long long)(1e9 + 7);
  if (c >= 0) return c;
  return c + (long long)(1e9 + 7);
}
long long po(long a, long b) {
  if (b == 0) {
    return 1;
  }
  long long z = po(a, b / 2);
  z = mod(z * z);
  if (b % 2 != 0) {
    z = mod(a * z);
  }
  return z;
}
using namespace std;
bool prime_(int n) {
  if (n == 1) {
    return false;
  } else if (n == 2) {
    return true;
  } else {
    for (int i = 2; i <= sqrt(n); i++) {
      if (n % i == 0) {
        return false;
      }
    }
    return true;
  }
}
long long gcd_(long long a, long long b) {
  if (a < b) {
    swap(a, b);
  }
  if (a % b == 0) {
    return b;
  } else {
    return gcd_(b, a % b);
  }
}
long long lcm_(long long x, long long y) { return (x / gcd_(x, y)) * y; }
map<pair<long long, long long>, long long> mp;
long long func(long long s, long long x) {
  long long res = 0;
  if (s < 0 || x < 0) {
    return 0;
  }
  if (s == 0 && x == 0) {
    return 1;
  }
  if (s == 1 && x == 0) {
    return 0;
  }
  if (s == 0 && x == 1) {
    return 0;
  }
  if (mp[pair<long long, long long>(s, x)] < 0) {
    res = 0LL - mp[pair<long long, long long>(s, x)];
    return res;
  }
  res += func(s / 2, x / 2);
  res = mod(res);
  res += func((s - 1) / 2, (x - 1) / 2);
  res = mod(res);
  res += func((s - 2) / 2, x / 2);
  res = mod(res);
  mp[pair<long long, long long>(s, x)] = 0LL - res;
  return res;
}
int main() {
  long long n;
  cin >> n;
  long long ans = func(n, n);
  cout << ans << endl;
  return 0;
}
