#include<bits/stdc++.h>
#define LL long long
#define LLL __int128
#define uint unsigned
#define ldb long double
#define uLL unsigned long
using namespace std;
const int N=20,P=1e9+7,M=4005;
int n,all,pn;vector<pair<int,int>>vec[N],ved[N][N];vector<int>pwk[M];
int a[N],b[N+1],ans[N+1],d[N],C[N+1][N+1],p[M],invp[M],mx[M],cur[M];
inline int qpow(int x,int y,int z=1){
    for(;y;(y>>=1)&&(x=1ll*x*x%P))if(y&1)z=1ll*z*x%P;return z;
}
inline void split(int x,vector<pair<int,int>>&vec){
    for(int i=2;i*i<=x;++i)if(x%i==0){
        int c=0;while(x%i==0)x/=i,++c;vec.emplace_back(i,c),p[pn++]=i;
    }
    if(x>1)vec.emplace_back(x,1),p[pn++]=x;
}
inline int fpow(int x,int y){
    if(pwk[x].empty())pwk[x].emplace_back(1);
    while(pwk[x].size()<=y)pwk[x].emplace_back(1ll*pwk[x].back()*invp[x]%P);
    return pwk[x][y];
}
signed main(){
    cin.tie(0)->sync_with_stdio(0);
    cin>>n;
    for(int i=0;i<n;++i)cin>>a[i],split(a[i],vec[i]);
    for(int i=0;i<n;++i)for(int j=i+1;j<n;++j)split(a[j]-a[i],ved[i][j]);
    sort(p,p+pn),pn=unique(p,p+pn)-p;
	for(int i=0;i<pn;++i)invp[i]=qpow(p[i],P-2);
    for(int i=0;i<n;++i)for(auto&[x,c]:vec[i])x=lower_bound(p,p+pn,x)-p;
    for(int i=0;i<n;++i)for(int j=i+1;j<n;++j)for(auto&[x,c]:ved[i][j])x=lower_bound(p,p+pn,x)-p;
    for(int i=2;i<=n;++i)cin>>b[i];
    for(int S=1;S<(1<<n);++S){
        int m=0;
        for(int j=0;j<n;++j)if(S>>j&1)d[m++]=j;
        if(m<2)continue;
        int r=2*(a[d[1]]-a[d[0]])%P;
        const auto ins=[&](const vector<pair<int,int>>&vec){
            for(const auto&[p,c]:vec)if((cur[p]+=c)>mx[p])mx[p]=cur[p];
        };
        const auto del=[&](const vector<pair<int,int>>&vec){
            for(const auto&[p,c]:vec)cur[p]-=c;
        };
        const auto clr=[&](const vector<pair<int,int>>&vec){
            for(const auto&[p,c]:vec)if(mx[p])r=1ll*r*fpow(p,mx[p])%P,cur[p]=0,mx[p]=0;
        };
        for(int j=2;j<m;++j)r=1ll*r*(a[d[j]]-a[d[j-1]])%P,ins(ved[d[j-1]][d[j]]);
        for(int j=1;j<m;++j){
			ins(vec[d[j-1]]),ins(vec[d[j]]),
            del(vec[d[j-1]]),del(vec[d[j]]);
			if(j+1<m)
				del(ved[d[j]][d[j+1]]),
            	ins(ved[d[j-1]][d[j]]);
		}
        for(int j=1;j<m;++j)clr(ved[d[j-1]][d[j]]);
        for(int j=0;j<m;++j)clr(vec[d[j]]);
        (ans[m]+=r)>=P&&(ans[m]-=P);
    }
    for(int i=0;i<=n;++i)
        for(int j=C[i][0]=1;j<=i;++j)
            C[i][j]=(C[i-1][j]+C[i-1][j-1])%P;
    for(int i=n;i>=2;--i)for(int j=2;j<i;++j)ans[j]=(ans[j]-1ll*C[i][j]*ans[i])%P;
    int tot=0;
    for(int i=2;i<=n;++i)tot=(tot+1ll*ans[i]*b[i])%P;
    cout<<(tot<0?tot+P:tot)<<'\n';
    return 0;
}
/*
*/