#include<bits/stdc++.h>
#define mem(i,a) memset(i,a,sizeof(i))
#define inf 0x3f3f3f3f3f3f3f3f
#define PII pair<int,int>
#define int  long long
#define endl '\n'
using namespace std;
const int maxn = 1e6 + 5,mod=1e9+7;
int f[maxn];
void solve() {
	string s;
	while(cin>>s){
		string s2;
		int n=s.size();
		if(n==1){
			if(s=="1"||s=="?")cout<<1<<endl;
			else cout<<0<<endl;
			continue;
		}
		// if(n%2==0){
			// for(int i=0;i<n-1;i++){
				// if(s[i]=='?'||s[i+1]=='?')s2+='?';
				// else{
					// int a=s[i]-'0',b=s[i+1]-'0';
					// a=abs(a-b);
					// s2+=(a+'0');
				// }
			// }
			// s=s2;
			// cout<<s2<<endl;
		// }
		int sum=0;
		for(int i=1;i<(n-1);i++){
			sum+=(s[i]=='?');
		}
		int ans=0;
		if((s[0]=='1'||s[0]=='?')&&(s.back()=='0'||s.back()=='?')){
			if(n%2==1)ans=(ans+f[sum]);
			else ans=(ans+f[max(0ll,sum-1)]);
		}
		if((s[0]=='0'||s[0]=='?')&&(s.back()=='1'||s.back()=='?')){
			if(n%2==1)ans=(ans+f[sum])%mod;
			else ans=(ans+f[max(0ll,sum-1)])%mod;
		}
		cout<<ans<<endl;
	}
	
}
signed main() {
	f[0]=1;
	for(int i=1;i<=1000000;i++){
		f[i]=f[i-1]*2%mod;
	}
	cin.tie(0);
	cout.tie(0);
	ios::sync_with_stdio(0);
		solve();
	return 0;
}