#pragma GCC target("avx2")

#define rd() ({long _v=0,_c;while(_c=*rp++-48,_c>=0)_v=_v*10+_c;_v;})
#define wt(v) {ulong _z=v;do*--wp=_z%10+48;while(_z/=10);}
#define rep(v,e) for(long v=0;v<e;++v)

typedef unsigned long ulong;

typedef struct {
	union { ulong i; double f; } u;
	struct { unsigned a,b; } v;
} S;

void radix_sort_aux(S*a,S*b,int n){
	int c[256];
	for(int i=0;i<256;++i){
		c[i]=0;
	}
	for(int i=0;i<n;++i){
		++c[a[i].u.i&255];
	}
	int t=0;
	for(int i=256;i--;){
		int u=c[i];
		c[i]=t;
		t+=u;
	}
	for(int i=0;i<n;++i){
		int j=c[a[i].u.i&255]++;
		b[j].u.i=a[i].u.i>>8|a[i].u.i<<56;
		b[j].v=a[i].v;
	}
}

void radix_sort(S*a,int n){
	static S b[200000];
	radix_sort_aux(a,b,n);
	radix_sort_aux(b,a,n);
	radix_sort_aux(a,b,n);
	radix_sort_aux(b,a,n);
	radix_sort_aux(a,b,n);
	radix_sort_aux(b,a,n);
	radix_sort_aux(a,b,n);
	radix_sort_aux(b,a,n);
}

S a[200000];

#define MD 1000000007

int main(){
	char*mmap();
	char*rp=mmap(0l,1l<<25,1,2,0,0ll);
	long n=rd();
	rep(i,n){
		a[i].v.a=rd();
	}
	rep(i,n){
		a[i].v.b=rd();
		a[i].u.f=(a[i].v.b-1.)/a[i].v.a;
	}
	radix_sort(a,n);
	ulong z=0,x=1;
	rep(i,n){
		z=(z+a[i].v.a*x)%MD;
		x=x*a[i].v.b%MD;
	}

	char wbuf[64],*wp=wbuf+sizeof wbuf;
	wt(z);
	write(1,wp,wbuf+sizeof wbuf-wp);
	_exit(0);
}