#include<iostream>
#include<vector>
#include<algorithm>
#define MOD 1000000007
using namespace std;
int main(){
  long long int r=0,c=0;
  long long int n,m; cin >> n >> m;
  vector<int> row(n), col(m); 
  for(int i=0;i<n;i++) cin >> row[i];
  for(int i=0;i<m;i++) cin >> col[i];
  sort(row.begin(), row.end(), greater<int>());
  sort(col.begin(), col.end(), greater<int>());
  long long int ans = 1;
  int rp=0,cp=0;
  for(int i=0;i<n*m;i++){
    int num = n*m-i;
  	if(rp<n && row[rp]==num && cp<m && col[cp]==num){
    	ans = ans; r++;c++;rp++;cp++;
    }else if(rp<n && row[rp]==num){
    	ans = (c * ans)%MOD; rp++; r++;
    }else if(cp<m && col[cp]==num){
    	ans = (r * ans)%MOD; cp++; c++;
    }else{
      	ans = ((r*c - i)*ans)% MOD; 
    }
  }
  cout << ans << endl;
  return 0;
}
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <vector>
#include <math.h>
#include <iomanip>
#include <bitset>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <utility>
#include <set>
#include <queue>
using namespace std;
typedef long long int ll;
typedef long double ld;
template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; }
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; }
const long long INF = 1LL << 60;
typedef pair<ll,ll> pairs;
vector<pairs> p;
bool pairCompare(const pair<double,ll>& firstElof, const pair<double,ll>& secondElof){
    return firstElof.first < secondElof.first;
}
bool pairCompareSecond(const pair<double,ll>& firstElof, const pair<double,ll>& secondElof){
    return firstElof.second < secondElof.second;
}
#define MAX_N 100100
#define MOD 1000000007
bool x[MAX_N];
ll num[MAX_N];
ll fibl[MAX_N]={0};
// 四方向への移動ベクトル
const int dx[4] = {1, 0, -1, 0};
const int dy[4] = {0, 1, 0, -1};

ll fib(ll a){
    if (fibl[a]!=0)return fibl[a];
    if (a==0){
        return 0;
    }else if (a==1){
        return 1;
    }
    return fibl[a]=fib(a-1)+fib(a-2);
}
 
ll eratosthenes(ll n) {
    int p = 0;
    for (ll i=0; i<=n; ++i) x[i] = true;
    x[0] = x[1] = false;
    for(int i=2; i<=n; ++i) {
        if(x[i]) {
            p++;
            for(int j=2*i; j<=n; j+=i) x[j] = false;
        }
        
        num[i] = p;
    }
    return p;
}
 
ll gcd(ll a,ll b){
    if (a%b==0)return(b);
    else return(gcd(b,a%b));
}
 
ll keta(ll N){
    int tmp{};
    while( N > 0 ){
        tmp += ( N % 10 );
        N /= 10;
    }
    N = tmp;
    return N;
}

int main(){
    ld n,ave=0,c;
    cin >> n;
    ld a[(ll)n];
    for (ll i=0;i<n;i++){
        cin >> a[i];
        ave+=a[i];
    }
    ave/=n;
    c=abs(a[0]-ave);
    bool flag=true;
    for (ll i=0;i<n;i++){
        if (c!=abs(a[i]-ave)){
            flag=false;
            c=min(c,abs(a[i]-ave));
        }
    }
    if (flag){
        cout << 0 << endl;
    }else{
        for (ll i=0;i<n;i++){
            if (c==abs(a[i]-ave)){
                cout << i << endl;
                return 0;
            }
        }
    }
    return 0;
}
