#include <bits/stdc++.h>
using namespace std;
#define ll long long
const ll mod=1e9+7;
const int N=1e4+5;
ll qpow(ll x,ll y)
{
    ll res=1;while(y){if(y&1ll) res=res*x%mod; x=x*x%mod,y>>=1;} return res;
}
int mul[N],inv[N];
void init()
{
    mul[0]=1;for(int i=1;i<N;i++) mul[i]=1ll*mul[i-1]*i%mod;
    inv[N-1]=qpow(mul[N-1],mod-2);
    for(int i=N-2;i>=0;i--) inv[i]=1ll*inv[i+1]*(i+1ll)%mod;
}
ll C(int x,int y){return x<y?0:1ll*mul[x]*inv[y]%mod*inv[x-y]%mod;}
int S[N],f[N][N],g[N][N],a,b,c,d,r,s,t,n;
int lim1,lim2,dn;
char ch[3][N];
inline int sum(int l,int r){return l>r?0:(S[r]-(l>0?S[l-1]:0)+mod)%mod;}
int main()
{
    //freopen("hypercube.in","r",stdin);
    //freopen("hypercube.out","w",stdout);
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);cin>>n;
    //cout<<qpow(2,n);
    init();
    cin>>r>>(ch[0]+1)>>s>>(ch[1]+1)>>t>>(ch[2]+1);
    ++r,++s,++t;
    for(int i=1;i<=n;i++)
    {
        int x=ch[0][i]-'0',y=ch[1][i]-'0',z=ch[2][i]-'0';
        y^=x,z^=x;if(!y&&!z) ++a; if(y&&!z) ++b;
        if(!y&&z) ++c; if(y&&z) ++d;
    }lim1=max(t,r),lim2=s-b,dn=s-b-d-c;
    S[0]=1;for(int i=1;i<=a;i++) S[i]=(S[i-1]+C(a,i))%mod;
    for(int A=0;A<=lim1;A++)
    {
        f[0][A]=sum(A,a);
        for(int i=1;i<=b;i++)
            f[i][A]=(f[i-1][A]+C(b,i)*sum(max(A-i,0),a)%mod)%mod;
    }
    for(int B=dn;B<=lim2;B++)
    {
        g[0][B-dn]=sum(max(B,0),a);
        for(int i=1;i<=b;i++)
            g[i][B-dn]=(g[i-1][B-dn]+C(b,i)*sum(max(B+i,0),a)%mod)%mod;
    }ll ans=0;
    for(int z=0;z<=c;z++)
        for(int w=0;w<=d;w++)
        {
            int A=max(max(t-c-d+z+w,r-z-w),0),B=s-b-d+w-z;
            int mid=min((A-B<0?-1:(A-B)/2),b); ll res=(mid>-1?f[mid][A]:0);
            if(mid<b) res=(res+g[b][B-dn]-(mid>-1?g[mid][B-dn]+mod:0))%mod;
            ans=(ans+C(c,z)*C(d,w)%mod*res%mod)%mod;
        }
    cout<<(qpow(2,n)+mod-ans)%mod;
    return 0;
}