一开始用线段树维护区间GCD和区间和,然后每次找所有GCD不为1 的区间,区间长度是偶数则结果加上区间和,如果是奇数,枚举1,3,5,7......找出value最小值,然后区间和减去最小值。 后来发现这样不对,因为一写区间被选择之后,原来不相交但是GCD不为1的数字变为了相交,又会使结果变大。然后想每次把选出来的区间删掉,然后重新建树,重复过程,知道找不到区间GCD不为一的情况,就输出累加的和,这样无限TLE。。。。。。
后来看了题解,人家说这是DP。。。。状态转移很好写。。。。也好理解
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn = 305;
typedef long long ll;
ll key[maxn]; ll value[maxn];
ll dp[maxn][maxn];
ll sum[maxn];
ll gcd(ll a, ll b){ return a == 0 ? b : gcd(b % a, a); }
int main()
{
//freopen("in","r",stdin);
//freopen("out","w",stdout);
int cases,n;
scanf("%d",&cases);
while(cases--)
{
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%I64d",&key[i]);
for(int i=1;i<=n;i++) scanf("%I64d",&value[i]);
sum[0] = 0;
for(int i=1;i<=n;i++) sum[i] = sum[i-1] + value[i];
memset(dp,0,sizeof(dp));
for(int len=1;len<=n;len++)
{
for(int i=1;i+len<=n;i++)
{
int j = i+len;
if(gcd(key[i],key[j])!=1){
if(j-i==1) dp[i][j] = dp[i][j] + value[i] + value[j];
else if(sum[j-1]-sum[i]==dp[i+1][j-1]) dp[i][j] = max(dp[i][j],dp[i+1][j-1] + value[i] + value[j]);
}
for(int k=i;k<=j;k++) dp[i][j] = max(dp[i][j],dp[i][k]+dp[k+1][j]);
}
}
printf("%I64d\n",dp[1][n]);
}
return 0;
}