#include<set>
#include<map>
#include<queue>
#include<vector>
#include<algorithm>
#include<bits/stdc++.h>
#define pr pair
#define f first
#define s second
#define ll long long
#define mp make_pair
#define pll pr<ll,ll>
#define pii pr<int,int>
#define piii pr<int,pii>
using namespace std;
const ll m=1e9+7;
ll jc[10004];
ll jy[10004];
int pf[20004][10004];
ll fd[10004],fa[10004],fb[10004],fc[10004];
ll C(int a,int b)
{
	return jc[a]*jy[b]%m*jy[a-b]%m;
}
ll pw(ll a,ll b=m-2)
{
	ll r=1;
	while(b)
	{
		if(b&1) r=r*a%m;
		a=a*a%m;
		b>>=1;
	}
	return r;
}
int n;
ll calc(int a,int b,string sa,string sb)
{
	int td=0,ta=0;
	for(int i=0;i<n;i++)
	{
		int fa=sa[i]-'0',fb=sb[i]-'0';
		if(fb>=1) fa=1-fa,fb=1-fb;
		if(fa) ta++;
		else td++;
	}
	for(int i=0;i<=ta;i++) fa[i]=C(ta,i);
	for(int i=0;i<=td;i++) fd[i]=C(td,i);
	for(int i=1;i<=td;i++) fd[i]=(fd[i-1]+fd[i])%m;
	ll rt=0;
	for(int i=0;i<=ta;i++)
	{
		int g=min(a-ta+i,b-i);
		if(g<0) continue;
		rt+=fa[i]*fd[g]%m;
	}
	return rt%m;
}
int main()
{
	jc[0]=1;
	for(int i=1;i<=10000;i++) jc[i]=jc[i-1]*i%m;
	jy[10000]=pw(jc[10000]);
	for(int i=9999;i>=0;i--) jy[i]=jy[i+1]*(i+1)%m;
	ios_base::sync_with_stdio(0);
	cin>>n;
	int a,b,c;
	string sa,sb,sc;
	cin>>a>>sa>>b>>sb>>c>>sc;
	ll ans=0;
	for(int i=0;i<=a;i++) ans+=C(n,i);
	for(int i=0;i<=b;i++) ans+=C(n,i);
	for(int i=0;i<=c;i++) ans+=C(n,i);
	int td=0,ta=0,tb=0,tc=0;
	for(int i=0;i<n;i++)
	{
		int fa=sa[i]-'0',fb=sb[i]-'0',fc=sc[i]-'0';
		if(fa+fb+fc>=2) fa=1-fa,fb=1-fb,fc=1-fc;
		if(fa) ta++;
		else if(fb) tb++;
		else if(fc) tc++;
		else td++;
	}
	cerr<<"F "<<ans<<endl;
	cerr<<"G "<<td<<' '<<ta<<' '<<tb<<' '<<tc<<endl;
	for(int i=0;i<=ta;i++) fa[i]=C(ta,i);
	for(int i=0;i<=tb;i++) fb[i]=C(tb,i);
	for(int i=0;i<=tc;i++) fc[i]=C(tc,i);
	for(int i=0;i<=td;i++) fd[i]=C(td,i);
	for(int i=0;i<=tc;i++) for(int j=0;j<=td;j++)
	{
		pf[n+j-i][i+j]=fc[i]*fd[j]%m;
	}
	for(int i=1;i<=n*2;i++) for(int j=0;j<=n;j++) pf[i][j]=(pf[i][j]+pf[i-1][j])%m;
	for(int i=0;i<=n*2;i++) for(int j=1;j<=n;j++) pf[i][j]=(pf[i][j]+pf[i][j-1])%m;
	for(int i=0;i<=ta;i++) for(int j=0;j<=tb;j++)
	{
		int ab=min(a-ta+i-j,b-tb+j-i),gc=c-tc-i-j+n;
		if(ab<0||gc<0) continue;
		ans+=fa[i]*fb[j]%m*pf[gc][ab]%m;
	}
	cerr<<"F "<<ans<<endl;
	ans+=m-calc(a,b,sa,sb);
	ans+=m-calc(b,c,sb,sc);
	ans+=m-calc(c,a,sc,sa);
	cout<<ans%m<<endl;
	return 0;
}