#include<stdio.h>
#include<string.h>
#include<iostream>
long long int MOD=1000000007;
const int LIMIT=2004;
long long int dp[2][2][2][2];
long long int nextdp[2][2][2][2];
//dp[R][L][old][now]
//next[R][L][now][next]
int as[LIMIT];

int f(int a,int b,int c){
	if((a==0)&&(b==0)&&(c==1)){
		return 1;
	}
	if((a==0)&&(b==1)&&(c==0)){
		return 1;
	}
	if((a==0)&&(b==1)&&(c==1)){
		return 1;
	}
	if((a==1)&&(b==0)&&(c==1)){
		return 1;
	}
	if((a==1)&&(b==1)&&(c==0)){
		return 1;
	}
	return 0;
}

int main(){
	int n;
	scanf("%d",&n);
	for(int i=0;i<n;i++){
		scanf("%d",&as[i]);
	}
	
	memset(dp,0,sizeof(dp));
	for(int i=0;i<2;i++){
		for(int j=0;j<2;j++){
			for(int k=0;k<2;k++){
				if(f(i,j,k)==as[0]){
					dp[i][j][j][k]=1;	
				}
			}
		}
	}
	long long int ans=0;
	for(int c=1;c<n;c++){
		memset(nextdp,0,sizeof(nextdp));
		for(int r=0;r<2;r++){
			for(int l=0;l<2;l++){
				for(int old=0;old<2;old++){
					for(int now=0;now<2;now++){
						int b1=as[c];
						if(c==n-2){
							int b=f(old,now,r);
							if(b==b1){
								nextdp[r][l][now][r]=(nextdp[r][l][now][r]+dp[r][l][old][now])%MOD;
							}
						}else if(c==n-1){
							if(r!=now)continue;
							int b=f(old,r,l);
							if(b==b1){
								ans=(ans+dp[r][l][old][now])%MOD;
							}
						}else{
							int b=f(old,now,0);
							if(b==b1){
								nextdp[r][l][now][0]=(nextdp[r][l][now][0]+dp[r][l][old][now])%MOD;
							}
							b=f(old,now,1);
							if(b==b1){
								nextdp[r][l][now][1]=(nextdp[r][l][now][1]+dp[r][l][old][now])%MOD;
							}
						}
					}
				}
			}
		}
		memcpy(dp,nextdp,sizeof(nextdp));
	}
	std::cout<<ans<<"\n";
}