#include<bits/stdc++.h> using namespace std; long long n,a[200],c[200],sum[200][200],top,nmax; long long dp[11451419]; inline long long ha(long long cc[]){ long long v=0; for(int i=1;i<=n;i++){ v+=v*(i+1)+cc[i]; } return v; } long long dfs(long long c[]){ long long hi=ha(c),deep=0; if(dp[hi]!=0)return dp[hi]; for(int i=1;i<=n;i++){ deep+=i*c[i]; } for(int i=1;i<=n;i++){ if(c[i]!=0){ c[i]--; dp[hi]=max(dp[hi],dfs(c)+sum[deep-i+1][deep]); c[i]++; } } return dp[hi]; } int main(){ scanf("%lld",&n); for(int i=1;i<=n;i++){ scanf("%lld",&a[i]); } for(int i=1;i<=n;i++){ scanf("%lld",&c[i]); } for(int i=1;i<=n;i++){ for(int j=i;j<=n;j++){ sum[i][j]=sum[i][j-1]; for(int k=i;k<=j;k++){ sum[i][j]+=2*a[k]*a[j]/(__gcd(a[k],a[j])); } sum[i][j]-=a[j]; } } dfs(c); long long nmax=0; for(int i=1;i<=11451418;i++){ nmax=max(nmax,dp[i]); } printf("%lld",nmax); return 0; }