#include <bits/stdc++.h>
using namespace std;
const double PI = 3.14159265358979323846;
const double eps = (1e-8);
const int MOD = 1000000007;
long long mod(long long x) { return (x % MOD + MOD) % MOD; }
long long mult(long long x, long long y) { return (x % MOD * y % MOD) % MOD; }
long long add(long long x, long long y) { return (x % MOD + y % MOD) % MOD; }
long long power(long long base, int exp) {
  long long ret = 1;
  while (exp) {
    if (exp & 1) ret = mult(ret, base);
    base = mult(base, base);
    exp /= 2;
  }
  return ret;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  cout.tie(0);
  string s;
  cin >> s;
  long long ans = 0;
  for (int i = 0; i < (int)(s).size(); ++i) {
    long long temp = 0;
    if (s[i] == '1')
      temp = mult(power(2, (int)(s).size() - i - 1),
                  power(2, (int)(s).size() - i - 1));
    ans = add(ans, mult(temp, power(2, i)));
  }
  cout << ans;
  return 0;
}
