#include<bits/stdc++.h>
using namespace std;
const int N=3e5+5,P=1e9+7;
char a[N]; int n,f[N][3][3];
int main()
{
	//freopen(".in","r",stdin);
	//freopen(".out","w",stdout);
	scanf("%s",a+1); n=strlen(a+1); f[0][0][0]=1;
	for(int i=0;i<n;i++) 
		for(int j=0;j<3;j++) for(int k=0;k<3;k++) if(f[i][j][k]){
			int v=f[i][j][k];
			for(int x=0;x<2;x++) if(x==a[i+1]-'0' || a[i+1]=='?'){
				if(j){
					int y;
					if(!x){
						y=j+1; if(y==3) y=1;
					}
					else y=j-1;
					(f[i+1][y][k]+=v)%=P;
				}
				else{
					if(k==2) (f[i+1][0][2]+=v)%=P;
					else{
						if(!x) (f[i+1][1][k]+=v)%=P;
						else (f[i+1][0][k+1]+=v)%=P;
					}
				}
			}
		}
	int ans=f[n][0][1];
	for(int i=0;i<3;i++) (ans+=f[n][i][2])%=P;
	cout<<ans<<endl;
	return 0;
}
