import math
MOD=1000000007
n=int(input())
m=int(input())
a=list(map(int,input().split()))
b=list(map(int,input().split()))
ans=0
for i in a:
    for j in b:
        ans=(ans+(i+1)*(j+1)+(math.gcd(i,j)+1)-2*(j+1))%MOD
print(ans)