#include <bits/stdc++.h>
#include <random>

using namespace std;

typedef unsigned long long ull;
typedef long long ll;
typedef double ld;
//#define ll __int128
//#define int ll
//#define int ll
//#define char ll
typedef vector<char> vc;
typedef vector<vc> vvc;
typedef vector<vvc> vvvc;
typedef pair<int, int> pii;
typedef pair<pii, pii> piii;
typedef pair<ll, ll> pll;
typedef vector<int> vi;
typedef vector<pii> vpi;
typedef vector< vi > vvi;
typedef vector< vvi > vvvi;
typedef vector<short> vs;
typedef vector<vs> vvs;
typedef vector<vvs> vvvs;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<vvl> vvvl;
typedef vector<ld> vld;
typedef vector<vld> vvld;
typedef vector<vvld> vvvld;
typedef vector<string> vst;
typedef vector<vst> vvst;
typedef pair<ld, ld> pld;

#define inmin(a, b) a = min(a, (b))
#define inmax(a, b) a = max(a, (b))
#define ALL(a) a.begin(),a.end()
#define RALL(a) a.rbegin(),a.rend()
#define sqr(x) ((x) * (x))
#define fori(i, n) for(int i = 0; i < int(n); ++i)
#define SZ(a) ((int)((a).size()))
#define triple(T) tuple<T, T, T>
#define quad(T) tuple<T, T, T, T>
#define watch(x) cerr << (#x) << " = " << (x) << endl;

#ifdef MAX_HOME
#define cerr cout
#else
#define cerr if (false) cerr
#endif

const double PI = 2 * acos(0.0);
#define rand shittttty_shit
mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
mt19937_64 rng_64(chrono::steady_clock::now().time_since_epoch().count());

const string DIGITS = "0123456789";
const string ALPH = "abcdefghijklmnopqrstuvwxyz";


template <class T0, class T1>
inline ostream & operator << (ostream &out, pair<T0, T1> &a) {
    return out << "{" << a.first << ", " << a.second << "}";
}

template <class T0, class T1>
inline istream & operator >> (istream &in, pair<T0, T1> &a) {
    return in >> a.first >> a.second;
}

template <class T0, class T1, class T2>
inline ostream & operator << (ostream &out, tuple<T0, T1, T2> &a) {
    return out << "{" << get<0>(a) << ", " << get<1>(a) << ", " << get<2>(a) << "}";
}

template <class T0, class T1, class T2, class T3>
inline ostream & operator << (ostream &out, tuple<T0, T1, T2, T3> &a) {
    return out << "{" << get<0>(a) << ", " << get<1>(a) << ", " << get<2>(a) << ", " <<  get<3>(a) << "}";
}

template<class T>
inline ostream & operator << (ostream &out, vector<T> &a) {
    out << "[";
    fori (i, a.size())
        out << a[i] << vector<string>{", ", "]  "}[i + 1 == a.size()];
    return out;
}


void smain();



signed main() {
    ios::sync_with_stdio(0);
    cin.tie(0); cout.tie(0);
#ifdef MAX_HOME
    freopen("input.txt", "r", stdin);
    clock_t start = clock();
#endif
    cout << setprecision(4) << fixed;
    smain();
#ifdef MAX_HOME
    cout << "\n\n\n\nTOTAL EXECUTION TIME: " << float( clock () - start ) /  CLOCKS_PER_SEC << endl;
#endif
    return 0;
}


const int M = 1e9 + 7;



void insum(int & a, int b) {
    a = (a + b) % M;
}

map<string, string> pattern;


string simple(string s) {
    if (s.size() <= 2)
        return s;

    if (s[0] == s[1] && s[0] == '1')
        return "11";
    for (auto pat : pattern) {
        for (int i = 0; i <= SZ(s) - SZ(pat.first); ++i) {
            if (s.substr(i, SZ(pat.first)) == pat.first) {
                s = s.substr(0, i) + pat.second + s.substr(i + SZ(pat.first), SZ(s));
                return simple(s);
            }
        }
    }
    return s;
}
const int L = 16;

pii tran[L];
vector<string> states;

void init() {


    pattern["000"] = "0";
    pattern["101"] = "1";
    pattern["011"] = "1";
    pattern["00100"] = "0";

    states = {"","0"
    ,"00"
    ,"001"
    ,"0010"
    ,"01"
    ,"010"
    ,"0100"
    ,"01001"
    ,"010010"
    ,"1"
    ,"10"
    ,"100"
    ,"1001"
    ,"10010"
    ,"11"};

    map<string, int> anti;
    fori (i, SZ(states)) {
        anti[states[i]] = i;
    }

    for (auto s : states) {
        string kek;
        kek = s + "0";
        kek = simple(kek);
        tran[anti[s]].first = anti[kek];
        kek = s + "1";
        kek = simple(kek);
        tran[anti[s]].second = anti[kek];
    }

}



void smain() {


    init();

    vi dp(L, 0);
    dp[0] = 1;

    string s;
    cin >> s;
    int n = s.size();

    int ans = 0;


    for (int i = 0; i < n; ++i) {
        vi ndp(L, 0);
        if (s[i] == '0' || s[i] == '?') {
            for (int j = 0; j < L; ++j) {
                insum(ndp[tran[j].first], dp[j]);
            }
        }
        if (s[i] == '1' || s[i] == '?') {
            for (int j = 0; j < L; ++j) {
                insum(ndp[tran[j].second], dp[j]);
            }
        }
        dp = ndp;
    }

//    for (auto x : states) {
//        cerr << x << endl;
//    }
    for (int i = 0; i < L; ++i) {
//        cerr << "dp[" << item.first << "] = " << item.second << endl;
//        watch(item.first);
//        cerr << "dp[" << i << "] = " << dp[i] << endl;
        if (count(ALL(states[i]), '1') > count(ALL(states[i]), '0'))
            insum(ans, dp[i]);
    }
    cout << ans;
}

Measured memory_bytes: 720