传送门
果然和dp有关的东西我绝对做不出来啊……
设\(dp[i][j]\)表示消完区间\([i,j]\)中的数之后能得到的最大分数，如果消不完则为\(-inf\)，否则枚举断点。顺便如果\(a[i],a[j]\)不互质可以用\(dp[i+1][j-1]+b[i]+b[j]\)来更新答案
然后设\(f[i]\)为前缀的答案，直接普通的dp即可
//minamoto
#include
#define R register
#define ll long long
#define inf 1e18
#define fp(i,a,b) for(R int i=a,I=b+1;iI;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
templateinline bool cmax(T&a,const T&b){return a'9'||ch<'0')(ch=='-')&&(f=-1);
for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
return res*f;
}
char sr[1<<21],z[20];int C=-1,Z=0;
inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;}
void print(R int x){
if(C>1<<20)Ot();if(x<0)sr[++C]='-',x=-x;
while(z[++Z]=x%10+48,x/=10);
while(sr[++C]=z[Z],--Z);sr[++C]='\n';
}
const int N=1005;
int gcd(int x,int y){return y?gcd(y,x%y):x;}
ll dp[N][N],f[N];int a[N],b[N];int n;bool Gcd[N][N];
int main(){
// freopen("testdata.in","r",stdin);
n=read();
fp(i,1,n)a[i]=read();
fp(i,1,n)b[i]=read();
fp(i,1,n)fp(j,i+1,n)Gcd[i][j]=gcd(a[i],a[j])>1?1:0;
fp(i,1,n-1)dp[i][i+1]=Gcd[i][i+1]?b[i]+b[i+1]:-inf;
fp(j,4,n)if(~j&1)fp(i,1,n-j+1){
dp[i][i+j-1]=-inf;
if(Gcd[i][i+j-1])dp[i][i+j-1]=dp[i+1][i+j-2]+b[i]+b[i+j-1];
fp(k,i+1,i+j-2)if((k-i)&1)
cmax(dp[i][i+j-1],dp[i][k]+dp[k+1][i+j-1]);
}fp(i,2,n){
f[i]=f[i-1];
fp(j,1,i-1)cmax(f[i],f[j-1]+dp[j][i]);
}printf("%lld\n",f[n]);
return 0;
}