from math import gcd
MOD=10**9+7
N=int(input())
M=int(input())
A=list(map(int,input().split()))
B=list(map(int,input().split()))
ans=0
for a in A:
    for b in B:
        d=b//gcd(a,b)
        W=b
        H=W*a//b
        ans+=(H+1)*(W+1)-2*W+(W-1)//d
        ans%=MOD
print(ans)