#include<bits/stdc++.h> using namespace std; long long n,a[200],c[200],sum[200][200],top,nmax,t[205]; map<long long,long long>dp; inline long long ha(long long cc[]){ long long v=0; for(int i=1;i<=n;i++){ v+=v*(t[i]+1)+cc[i]; } return v%11451419; } 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(){ // freopen("text.in","r",stdin); // freopen("text.out","w",stdout); 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]); t[i]=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]; } } cout<<dfs(c); return 0; }