思路:提前处理最小质因子。对每个格子i有:dp[i]=max(dp[i],dp[j]+a[i]),j为可以转移到i的格子;
#include<bits/stdc++.h>
#define ll long long
using namespace std;
int n;
int a[10010],dp[10010];
int mip[10010];
bool pri[10010];
void prime()
{
memset(pri,true,sizeof(pri));
mip[1]=1;
for(int i=2;i<=10000;i++)
{
if(pri[i]==false)continue;
mip[i]=i;
for(int j=i*2;j<=10000;j+=i)
{
if(pri[j])
{
mip[j]=i;
pri[j]=false;
}
}
}
}
int main()
{
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
prime();
cin>>n;
for(int i=1;i<=n;i++)
{
cin>>a[i];
dp[i]=-1e9-10;
}
dp[1]=a[1];
for(int i=2;i<=n;i++)
{
for(int j=1;j<i;j++)
{
if(i<=j+mip[n-j])
{
dp[i]=max(dp[i],dp[j]+a[i]);
}
}
}
cout<<dp[n]<<"\n";
return 0;
}