#include <iostream>
#include <string>
#include <algorithm>
#include <vector>
#include <array>
#include <queue>
#include <deque>
#include <map>
#include <set>
#include <sstream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <list>
#include <numeric>
#include <stack>
#include <iomanip>
#include <random>
#include <complex>
#include <functional>
#include <tuple>

using namespace std;

#define Rep(i,a,b) for(int i = a; i < b; ++i)
#define rep(i,b) Rep(i,0,b)
#define allof(a) (a).begin(), (a).end()

#define Yes(q) ((q) ? "Yes" : "No")
#define YES(q) ((q) ? "YES" : "NO")
#define Possible(q) ((q) ? "Possible" : "Impossible")
#define POSSIBLE(q) ((q) ? "POSSIBLE" : "IMPOSSIBLE")

using ll = long long;

using pint = std::pair<int, int>;
using pll = std::pair<ll, ll>;


constexpr int inf = 1e9 + 7;
constexpr ll infll = 1ll << 60ll;
constexpr ll mod = 1e9 + 7;
// 0~3までは右左下上 4~7までは斜め
constexpr int dx[] = { 1, 0, -1, 0, 1, 1, -1, -1 };
constexpr int dy[] = { 0, -1, 0, 1, 1, -1, -1, 1 };

namespace {
	template<typename T> void chmax(T& a, T b) { a = std::max(a, b); }
	template<typename T> void chmin(T& a, T b) { a = std::min(a, b); }
	template<typename T> void chadd(T& a, T b) { a = a + b; }

	void Cout(long long x, const char* end = "\n") { std::cout << x << end; }
	template <typename T> void Cout(const T& x, const char* end = "\n") { std::cout << x << end; }
	template <typename T> void Cout(const std::vector<T>& x, const char* sep = " ", const char* end = "\n") { for (std::size_t i = 0, sz = x.size(); i < sz; i++) { std::cout << x[i] << (i == sz - 1 ? end : sep); } }

	void CCout(long long x, const char* end = "\n") { std::cerr << x << end; }
	template <typename T> void CCout(const T& x, const char* end = "\n") { std::cerr << x << end; }
	template <typename T> void CCout(const std::vector<T>& x, const char* sep = " ", const char* end = "\n") { for (std::size_t i = 0, sz = x.size(); i < sz; i++) { std::cerr << x[i] << (i == sz - 1 ? end : sep); } }

	// 標準入出力
	struct inp {
		std::size_t sz;
		inp(std::size_t _sz = 1) : sz(_sz) {}
		template <typename T> operator T () const { T a; std::cin >> a; return a; }
		template <typename T> operator std::vector<T>() const { vector<T> a(sz); for (std::size_t i = 0; i < sz; i++) std::cin >> a[i]; return a; }
		template <typename T, typename U> operator std::pair<T, U>() const { T f; U s; std::cin >> f >> s; return std::pair<T, U>(f, s); }
	};

	inp inp1; // input one
}

int main() {

	string l = inp1;

	// 片方が0なら必ずOK
	// 繰り上げがない場合はOK (同じ桁に1がそろわない)
	size_t n = l.length();
	// AおよびBの値の上からk桁目まで確定させ，その時点ですでにA+BがL以下であることが分かっている/いない(0/1)組の数
	vector<vector<ll>> dp(n + 1, vector<ll>(2, 0));

	dp[0][1] = 1;

	rep(i, n) {
		// L以下確定
		// 1のとき
		dp[i + 1][0] += dp[i][0];
		// 0
		dp[i + 1][0] += dp[i][0] * 2;

		// 確定していない
		int d = l[i] - '0';
		if (d == 0) {
			// 0で決まり
			dp[i + 1][1] += dp[i][1];
		}
		else {
			// 0にすると確定になる (ただし，向こう側も0なので*1)
			dp[i + 1][0] += dp[i][1];
			// 1にすると不確定のまま また，0にして向こうを1にすると不確定が続く
			dp[i + 1][1] += dp[i][1] * 2;
		}

		dp[i + 1][0] %= mod;
		dp[i + 1][1] %= mod;

	}

	cerr << dp[n][0] << ", " << dp[n][1] << '\n';
	Cout((dp[n][0] + dp[n][1]) % mod);

	return 0;
}


Measured memory_bytes: 5760