import std.stdio, std.string, std.conv;
import std.algorithm, std.math;
import std.range;
import std.container.rbtree, std.container.dlist;
import std.container.binaryheap, std.container.array;
import std.typecons;

alias mstring = char[];
const long INF = 1L << 60L;
const long mod = 1_000_000_000 + 7;

void chmin (T)(ref T x, T y) {
  x = min(x, y);
}

void chmax (T)(ref T x, T y) {
  x = max(x, y);
}

// 単一の数値を取得
// readln.chomp.to!int;
//  or
// int a;
// readf("%s\n", a);

// 複数の数値を取得(可変数個の場合推奨)
// readln.chomp.split.map!(to!long).array

// 複数の数値を取得(固定数個の場合、推奨)
// int a, b;
// readf("%s %s\n", &a, &b);

// インデントを一個ずらして複数の数値を取得(累積和とかで1-indexedの方が良い時がある)
// long[] vs = readln.chomp.split.map!(to!long).array;
// vs = [0L] ~ vs;

// 小数点は以下で指定
// double ret = 10.0;
// writefln("%.12f", ret);

// 配列の最後の要素は arr[$-1] でアクセスできる

void main () {
  long N = readln.chomp.to!long;
  const long maxK = 62L;

  auto dp = new long[][][][][](maxK + 1, 2, 2, 3, 2);
  dp[maxK][0][0][0][0] = 1L;
  for (long i = maxK; i > 0L; i--) {
    int ni = (N >> (i - 1L)) & 1L;

    foreach (j; 0 .. 2)
      foreach (k; 0 .. 2)
        foreach (l; 0 .. 3)
          foreach (m; 0 .. 2)
            foreach (xi; 0 .. 2)
              foreach (yi; 0 .. 2) {
                if (j == 0 && xi > yi)
                  continue;
                if (k == 0 && yi > ni)
                  continue;

                int nj = cast(int)(j == 1 || xi < yi);
                int nk = cast(int)(k == 1 || yi < ni);
                int nl;
                if (l == 0) {
                  if (xi == yi && xi == 0)
                    nl = 0;
                  else if (xi == yi && xi == 1)
                    nl = 2;
                  else
                    nl = 1;
                } else
                  nl = l;
                int nm = cast(int)(m == 1 || (xi == yi && xi == 1));

                dp[i - 1][nj][nk][nl][nm] += dp[i][j][k][l][m];
                dp[i - 1][nj][nk][nl][nm] %= mod;
              }
  }

  long ret = 0L;
  foreach (j; 0 .. 2)
    foreach (k; 0 .. 2) {
      ret += dp[0][j][k][1][1];
      ret %= mod;
    }

  writeln(ret);
}