#include<bits/stdc++.h>
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rof(i,a,b) for(int i=(a);i>=(b);--i)
using namespace std;
const int Maxn=1e6,Mod=1e9+7;

inline int Pow(int x,int y)
{
    int res=1;
    while(y)
    {
        if(y&1) res=1ll*res*x%Mod;
        x=1ll*x*x%Mod,y>>=1;
    }
    return res;
}

int n; string s;

inline void Solve()
{
    n=s.size(),s=' '+s;
    int f[2]={1,0};
    For(i,1,n)
    {
        int k=((n-1)&(i-1))==(i-1);
        if(k==0)
        {
            f[0]=(1ll+(s[i]=='?'))*f[0]%Mod;
            f[1]=(1ll+(s[i]=='?'))*f[1]%Mod;
        }
        else
        {
            int g0=0,g1=0;
            if(s[i]!='1') (g0+=f[0])%=Mod,(g1+=f[1])%=Mod;
            if(s[i]!='0') (g1+=f[0])%=Mod,(g0+=f[1])%=Mod;
            f[0]=g0,f[1]=g1;
        }
    } cout<<f[1]<<endl;
}

int main()
{
    // freopen("1.in","r",stdin);

    while(cin>>s) Solve();
    return 0;
}