#include <bits/stdc++.h>
using namespace std;
template <class T>
T cross(complex<T> a, complex<T> b) {
  return imag(a * conj(b));
}
template <class T>
T min(T a, T b, T c) {
  return min(a, min(b, c));
}
template <class T>
T max(T a, T b, T c) {
  return max(a, max(b, c));
}
template <class T>
void setmin(T &a, T b) {
  if (b < a) a = b;
}
template <class T>
void setmax(T &a, T b) {
  if (b > a) a = b;
}
const int MAXN = 200;
const long long MOD = 1000000007;
int N;
string s;
long long powMod(long long p, long long e) {
  long long out = 1;
  for (int i = 0; i < e; i++) out *= p, out %= MOD;
  return out;
}
void io(istream &in, ostream &out) {
  in >> s;
  N = s.size();
  long long ans = 0;
  for (int i = 0; i < N; i++) {
    if (s[i] == '1') ans += powMod(2, i) * powMod(4, N - i - 1);
    ans %= MOD;
  }
  out << ans << '\n';
}
int main() {
  ifstream fin("A.in");
  if (fin.good()) {
    ofstream fout("A.out");
    io(fin, fout);
    fin.close();
    fout.close();
  } else
    io(cin, cout);
  return 0;
}
