//Marcin Knapik, Potyczki Algorytmiczne 2019, Runda 5, zadanie: "Trzy kule"[B]
//złożoność O(n^2)

#pragma GCC optimize("O3")
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")


#include<bits/stdc++.h>
using namespace std;

#define ll long long
#define FOR(i, n) for(int i = 0; i < n; i++)

const int N = 10007;
const int mod = 1e9+7;

int dp1[N];
// int dp2[N][N];
int dp3[N][N];
int dp4[N][N];

int n;

struct kula{
	int r;
	int tab[N];
};

//bool cmp(kula a, kula b){
//	return a.r < b.r;
//}

kula kule[3];
int sil[N];
int odw[N];

ll fast(ll a, ll b){
	ll ret = 1;
	ll pom = a;
	while(b){
		if(b&1)
			ret = ((ll)ret*pom)%mod;
		pom = ((ll)pom*pom)%mod;
		b>>=1;
	}
	return ret;
}

void prep(){
	sil[0] = 1;
	for(int i = 1; i < N; i++)
		sil[i] = (ll)sil[i - 1] * i % mod;
	odw[N-1] = fast(sil[N - 1], mod - 2);
	for(int i = N-2; i >= 0; i--)
		odw[i] = (ll)odw[i + 1] * (i+1) % mod;
}

int newton(int a, int b){
	if(b < 0 || b > a)
		return 0;
	return ((ll)sil[a] * odw[b]%mod) * odw[a-b] % mod;
}

int wspolna(kula&A){
	int ret = 0;
	for(int i = 0; i <= A.r; i++)
		ret = (ret + newton(n, i))%mod;
	return ret;
}

int wspolna(kula&A, kula&B){
	// cout << "WSPOLNA " << endl;
	int same = 0;
	int inne = 0;
	int ret = 0;

	FOR(i, n)
		if(A.tab[i] == B.tab[i])
			same++;
		else
			inne++;

	// cout << same << ' ' << inne << endl;
	
	for(int i = 0; i < N; i++)
		dp1[i] = 0;
	
	dp1[0] = newton(same, 0);
	
	for(int i = 1; i < N; i++)
		dp1[i] = (dp1[i-1] + newton(same, i)) % mod;

	// for(int i  = 0; i < 5; i++){
		// cout << dp1[i] << ' ';
	// }
	// cout << endl;
	for(int i = 0; i <= inne; i++){
		int a = A.r - i;
		int b = B.r - (inne - i);
		int m = min(a, b);
		// cout << "M " << m << endl;
		if(m >= 0)
			ret = (ret + (ll)dp1[m]*newton(inne, i)) % mod;
	}
	// cout << "WSPOLNA DWOCH " << ret << endl;
	return ret;
}
#define ii pair<int, int>
#define f first
#define s second
ii zamien(int a, int b){
	if(a>=0 && b >= 0)
		return {a, b};
	return {0, b+a};
}

int wspolna(kula&A, kula&B, kula&C){
	int same = 0;
	int _A = 0;
	int _B = 0;
	int _C = 0;

	FOR(i, n)
		if(A.tab[i] == B.tab[i] && A.tab[i] == C.tab[i])
			same++;
		else if(A.tab[i] == B.tab[i])
			_C++;
		else if(A.tab[i] == C.tab[i])
			_B++;
		else
			_A++;

	// cout << same << ' ' << _A << ' ' << _B << ' ' << _C << endl;
	// cout << endl;

	// for(int i = 0; i <= 3; i++){
	// 	for(int j = 0; j <= 3; j++){
	// 		cout << dp2[i][j] << ' ';
	// 	}
	// 	cout << endl;
	// }
		
	for(int i = 0; i <= _C; i++)
		for(int j = 0; j <= same; j++)
			if(j + (_C - i) <= C.r)
				dp3[i][j] = dp4[i][j] = (ll)newton(same, j) * newton(_C, i) % mod;

	// cout << endl;
	
	// for(int i = 0; i <= 3; i++){
	// 	for(int j = 0; j <= 3; j++){
	// 		cout << dp3[i][j] << ' ';
	// 	}
	// 	cout << endl;
	// }

	for(int j = n; j >= 1; j--)
		for(int i = n; i >= 0; i--)
			dp3[i][j] = (dp3[i][j] + dp3[i][j-1]) % mod;

	// cout << endl;

	// for(int i = 0; i <= 3; i++){
	// 	for(int j = 0; j <= 3; j++){
	// 		cout << dp3[i][j] << ' ';
	// 	}
	// 	cout << endl;
	// }

	for(int j = 1; j <= n; j++)
		for(int i = 0; i <= n; i++){
			if(i > 0)
				dp3[i][j] = (dp3[i][j] + dp3[i-1][j-1]) % mod;
			if(j > 1 && i > 0) // TU może być bug
				dp3[i][j] = (dp3[i][j] - dp3[i][j-2] + mod) % mod; 
			dp3[i][j] = (dp3[i][j] + dp3[i+1][j-1]) % mod;
			dp4[i][j] = (dp4[i][j] + dp4[i+1][j-1]) % mod;
		}

	// cout << endl;

	// for(int i = 0; i <= 3; i++){
	// 	for(int j = 0; j <= 3; j++){
	// 		cout << dp3[i][j] << ' ';
	// 	}
	// 	cout << endl;
	// }

	// cout << endl;

	// for(int i = 0; i <= 3; i++){
	// 	for(int j = 0; j <= 3; j++){
	// 		cout << dp4[i][j] << ' ';
	// 	}
	// 	cout << endl;
	// }

	// cout << endl;

	int wynik = 0;
	for(int i = 0; i <= _A; i++){
		for(int j = 0; j <= _B; j++){
			int teraz = (ll)newton(_A, i) * newton(_B, j) % mod;
			int val = min(A.r - i - (_B - j), B.r - j - (_A - i)); // ograniczenie drugie
			int inc = (_B - j) + (_A - i);
			// cout << val << ' ' << C.r << ' ' << _C << ' ' << inc << endl;
			int pom = val + (C.r - _C - inc);
			// inc + (_C - c) + d <= C.r
			// inc == 2 więc

			int d = pom / 2;
			int c = val - d;
			if(pom & 1)
				c--;

			// cout << i << ' ' << j << ' ' << c << ' ' << d << ' ' << pom << endl; 
			if(pom < 0){
				continue;
			}
			ii pom2 = zamien(c, d);
			if(pom2.f>=0 && pom2.s>=0 && pom2.f < N && pom2.s < N){
				// ii pom = {c, d};
				wynik = (wynik + (ll)teraz * dp3[pom2.f][pom2.s]) % mod;	
			}
			if(pom & 1){
				pom2 = zamien(c+1, d); 
				if(pom2.f>=0 && pom2.s>=0 && pom2.f < N && pom2.s < N){
					wynik = (wynik + (ll)teraz * dp4[pom2.f][pom2.s]) % mod;						
				}
			}
		}
	}
	// cout << "ZWRACAM " << wynik << endl;
	return wynik;
}

int main(){
	ios::sync_with_stdio(0);
	cin.tie(0);

	cin >> n;
	prep();
	FOR(i, 3){
		cin >> kule[i].r;
		FOR(j, n){
			char znak;
			cin >> znak;
			if(znak == '1')
				kule[i].tab[j] = 1;
			else
				kule[i].tab[j] = 0;
		}
	}
//	sort(kule, kule+3);
	
	int odpowiedz = 0;
	// cout << "COOO? " << endl;
	FOR(i, 3)
		odpowiedz = (odpowiedz + wspolna(kule[i]))%mod;
	FOR(i, 3)
		odpowiedz = (odpowiedz - wspolna(kule[i], kule[(i+1)%3]) + mod)%mod;
	odpowiedz = (odpowiedz + wspolna(kule[0], kule[1], kule[2]))%mod;

	cout << odpowiedz << '\n';
}