#include <bits/stdc++.h>
using namespace std;
typedef struct point {
  int start;
  int end;
} point;
bool sort_cmp(const point& lhs, const point& rhs) {
  if (lhs.end < rhs.end)
    return true;
  else
    return false;
}
string s1;
long long powi[102];
int main() {
  long long i, j, k, l, n, m, a, t, s, x, y, ans, p;
  cin >> s1;
  ans = 0;
  p = 1;
  for (i = 0; i < 101; i++) {
    powi[i] = p;
    p *= 2;
    p %= 1000000007;
  }
  n = s1.size();
  for (i = 0; i < n; i++) {
    if (s1[n - 1 - i] == '1') {
      ans += powi[n - 1 - i] * ((powi[i] * powi[i]) % 1000000007);
      ans %= 1000000007;
    }
  }
  cout << ans << endl;
  return 0;
}
