#include<cstdio>
#include<algorithm>
using namespace std;
long long floor_sum(long long N,long long M,long long A,long long B)
//Sum[floor((A*i+B)/M),{i,0,N-1}]
{
	long long ans=0;
	if(B>=M)
	{
		long long BM=B/M;
		ans+=N*BM;
		B-=BM*M;
	}
	while(true)
	{
		long long AM=A/M;
		ans+=N*(N-1)/2*AM;
		A-=AM*M;
		if(A*N+B<M)break;
		long long Ym=(A*N+B)/M,Xm=Ym*M-B;
		long long TX=(Xm+A-1)/A;
		ans+=(N-TX)*Ym;
		N=Ym;
		B=TX*A-Xm;
		swap(A,M);
	}
	return ans;
}
int N,M;
int A[2000],B[2000];
main()
{
	scanf("%d%d",&N,&M);
	for(int i=0;i<N;i++)scanf("%d",&A[i]);
	for(int i=0;i<M;i++)scanf("%d",&B[i]);
	long ans=0;
	for(int i=0;i<N;i++)
	{
		for(int j=0;j<M;j++)
		{
			ans+=2*floor_sum(B[j],B[j],A[i],A[i]);
		}
		ans%=(int)1e9+7;
	}
	printf("%ld\n",ans);
}