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

int64_t f(int a, int b, int c, int n) {
    if (n < 0 || a == 0) return 0;
    if (a >= c || b >= c) 
        return (a / c) * (1LL * n * (n + 1) / 2) + 1LL * (b / c) * (n + 1) + f(a % c, b % c, c, n);
    int64_t m = (1LL * a * n + b) / c;
    return 1LL * n * m - f(c, c - b - 1, a, m - 1);
}

int main(){
  cin.tie(nullptr);
  ios_base::sync_with_stdio(false);
  int n,m;cin>>n>>m;
  vector<int>a(n),b(m);
  const long long mod=1e9+7;
  for(int i=0;i<n;i++)cin>>a[i];
  for(int i=0;i<m;i++)cin>>b[i];
  long long ret=0;
  for(int i=0;i<n;i++){
    for(int j=0;j<m;j++){
      long long add=f(a[i],0,b[j],b[j]);
      ret=(ret+add*2)%mod;
    }
  }
  cout<<ret<<endl;
}