6287. 2019.08.09【NOIP提高组A】扭动的树

Description

Input

Output

Sample Input

Sample 1:
4
2 3
6 4
9 8
12 1

Sample 2:
20
64978574415886122 263411
40589037247202745 239844
19724737874528206 167360
49216095485959384 760606
65063121727264647 659450
16572376111094320 726552
72014092598616298 133699
52843699826658793 427487
43374492289647376 552030
22047612465142862 605387
92386136280598953 718860
6436388687842008 368771
87727847161227820 880866
43622103777719758 352810
36870904328895185 322737
48993192459657624 456880
93250693206986868 619976
77407991580158822 861256
974508361120026 344635
77136053229840400 465474

Sample Output

Sample1:
51

Sample2:
101007480

Data Constraint

Solution

dp。
原来二叉搜索树是用key来排序的。。。
我们先以key为关键字将a从小到大排序，之后我们对于每个i它的左边比它小的都是它的左二子，右边比它大的都是它的右儿子，那么我们就相当于做一个区间dp，设f[i][j][k]表示在i-j这段区间它们的父亲是k的最大sum和，但是我们发现因为[i,j]这段区间的父亲只可能是i-1或j+1,所以k只用开0/1表示是前面的还是后面的就可以了，对于转移，我们先预处理出两两之间的gcd，val的前缀和以及区间长度为1的情况，f[i][i][0/1]条件是i和i-1或i+1的gcd不为1。然后转移让区间长度从小到大枚举，f[i][j][0/1]可以从f[i][k-1][1]+f[k+1][j][0]+sum[j]-sum[i-1]，注意边界情况要特殊处理。最后答案就是f[1][n][0/1]。注意最先将f赋值成负数。

Code

#include<cstdio>
#include<algorithm>
#include<cstring>
#define I int
#define F(i,a,b) for(I i=a;i<=b;i++)
#define Fd(i,a,b) for(I i=a;i>=b;i--)
#define N 500
#define ll long long
#define mem(a,b) memset(a,b,sizeof(a))
#define open(x) freopen(x".in","r",stdin);freopen(x".out","w",stdout);
using namespace std;
ll rd(){
	ll x=0,w=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-') w=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*w;
}
ll gcd(ll x,ll y){
	if(x%y==0) return y;
	return gcd(y,x%y);
}
ll n,f[N][N][2],g[N][N],s[N],l,r,ans,st;
struct node{ll k,v;}a[N];
I cmp(node x,node y){return x.k<y.k;}
I main(){
	open("tree");
	n=rd();
	F(i,1,n) a[i]=node{rd(),rd()};
	sort(a+1,a+1+n,cmp);
	mem(f,-0x3f);
	F(i,1,n)
		F(j,1,n) g[i][j]=gcd(a[i].k,a[j].k);
	F(i,1,n){
		s[i]=s[i-1]+a[i].v;
		if(i!=1&&g[i][i-1]!=1) f[i][i][0]=a[i].v;
		if(i!=n&&g[i][i+1]!=1) f[i][i][1]=a[i].v;
	}
	F(i,2,n){
		F(l,1,n-i+1){
			r=l+i-1;
			F(k,l,r){
				if(k==l) st=f[l+1][r][0]+s[r]-s[l-1];
				else if(k==r) st=f[l][r-1][1]+s[r]-s[l-1];
				else st=f[l][k-1][1]+f[k+1][r][0]+s[r]-s[l-1];
				if(l!=1&&g[k][l-1]!=1) f[l][r][0]=max(f[l][r][0],st);
				if(r!=n&&g[k][r+1]!=1) f[l][r][1]=max(f[l][r][1],st);
				if(i==n) ans=max(ans,st);
			}
		}
	}
	if(!ans) printf("-1\n");
	else printf("%lld\n",ans);
	return 0;
}

作者：zsjzliziyang 
QQ:1634151125 
转载及修改请注明 
本文地址:https://blog.csdn.net/zsjzliziyang/article/details/98997948