#include<iostream>
using namespace std;
long mod=1e9+7,ans;
int T[8]={0,1,1,1,0,1,1,0};
int N,A[2000];
long dp[2020][4];
main()
{
	cin>>N;
	for(int i=0;i<N;i++)cin>>A[i];
	for(int x=0;x<4;x++)
	{
		for(int i=0;i<=N;i++)for(int j=0;j<4;j++)dp[i][j]=0;
		dp[0][x]=1;
		for(int i=1;i<N-1;i++)
		{
			for(int j=0;j<4;j++)
			{
				for(int k=0;k<2;k++)
				{
					if(A[i]==T[j*2+k])
					{
						(dp[i][j%2*2+k]+=dp[i-1][j])%=mod;
					}
				}
			}
		}
		for(int j=0;j<4;j++)
		{
			if(A[N-1]==T[j*2+x/2]&&A[0]==T[j%2*4+x])(ans+=dp[N-2][j])%=mod;
		}
	}
	cout<<ans<<endl;
}