[NOI2009]二叉查找树
题目大意：
给定一棵严格的treap，父亲节点的优先级必然小于儿子节点的。权值按照二叉树的定义，左儿子小于父亲小于右儿子。
深度从1开始定义，每个点除优先级、数值之外，还有一个访问频度。
访问频度所产生的代价是：访问频度*该点深度（这和事实相符）
可以用给定的k的代价，修改任意个点的优先级为任意实数（当然，修改优先级，树的形态，各点深度就可能变化了）
最终的总代价为：频度产生代价+修改代价。
最小化这个总代价。
N<=70，1<=K<=30000000
分析：
平衡树是一个动态的数据结构，难以抓住形态的变化，也不方便记录深度之类。所以必须抓住不变的量当做突破口。
不管平衡树怎么转，根据二叉树的定义，它的中序遍历一定是不变的。
所以我们可以找到这棵树的中序遍历，就把这棵树变成了一个静态的区间，只不过每个区间所代表的点的优先级可能会变。
发现，每一个连续的子区间，都对应treap的连续一部分。可以把小的区间先建树，再把大的区间用小的区间合并。我们合并的时候枚举的划分点，就是这部分treap的树根
区间DP顺理成章。
除了f[l][r]之外，为了维护优先级的关系，必然要再记录一维。
发现，只要根节点的优先级确定，子树的优先级的范围就确定了。
所以考虑记录根节点优先级。（这里优先级只考虑相对大小，而且范围又大，所以要离散化为1~n）
但是，朴素的f[l][r][w]中，w单单记录根节点优先级的话，由于子树所有大于w的都可以转移，还要for一遍。状态n^3,转移n^2，会爆。
所以，我们令f[l][r][w]表示，将l~r这段区间建成treap，其中根节点优先级大于等于w的最小代价。
根据枚举的根节点是否修改，可以设计转移方程是：
修改：
f[l][r][w]=min(f[l][r][w],f[l][k-1][w]+f[k+1][r][w]+K+sum[r]-sum[l-1]) ————其中，sum[i]表示，区间中，1~i的访问频度和
当划分点的优先级ch[k]大于w时，可以不修改。
f[l][r][w]=min(f[l][r][w],f[l][k-1][ch[k]]+f[k+1][r][ch[k]]+sum[r]-suim[l-1])
最后答案就是：f[1][n][1];
注意，o循环的时候，必须倒序！！因为ch[k]>=o时候，要从o更大的地方获取最小值，必须先把o较大的处理完。
代码1（未简化）：
#include<bits/stdc++.h> using namespace std; typedef long long ll; const int N=77; const ll inf=2e18; ll f[N][N][N]; int n; ll m; int a[N]; int tot; int w[N],p[N],d[N]; int prio[N]; ll sum[N]; int ch[N];//离散化后的优先级 bool cmp(int a,int b) { return w[a]<w[b]; } int main() { scanf("%d%lld",&n,&m); for(int i=1;i<=n;i++) scanf("%d",&w[i]),a[i]=i; for(int i=1;i<=n;i++) scanf("%d",&p[i]),prio[i]=p[i]; for(int i=1;i<=n;i++) scanf("%d",&d[i]); sort(a+1,a+n+1,cmp); sort(prio+1,prio+n+1); for(int i=1;i<=n;i++) ch[i]=lower_bound(prio+1,prio+n+1,p[a[i]])-prio;//离散化 for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) for(int k=1;k<=n;k++) f[i][j][k]=inf; for(int i=1;i<=n;i++) for(int k=n;k>=1;k--) { if(k<=ch[i]) f[i][i][k]=d[a[i]]; else f[i][i][k]=d[a[i]]+m;//注意，k>ch的时候，不一定是+oo，可以通过修改改变 }//l=1的初值 for(int i=1;i<=n;i++) sum[i]=sum[i-1]+d[a[i]];//前缀和 for(int l=2;l<=n;l++) for(int i=1;i<=n;i++) { int j=l+i-1; if(j>n) break; if(l==2)//长度为二的时候，只能二并一 { for(int o=n;o>=1;o--) { f[i][j][o]=min(f[i][j][o],f[i][i][o]+m+sum[j]-sum[i-1]); f[i][j][o]=min(f[i][j][o],f[j][j][o]+m+sum[j]-sum[i-1]); if(ch[i]>=o) f[i][j][o]=min(f[i][j][o],f[i][i][ch[i]]+f[j][j][ch[i]]+sum[j]-sum[i-1]-d[a[i]]); if(ch[j]>=o) f[i][j][o]=min(f[i][j][o],f[i][i][ch[j]]+f[j][j][ch[j]]+sum[j]-sum[i-1]-d[a[j]]); } } else{ for(int o=n;o>=1;o--) for(int k=i;k<=j;k++) { if(k==i)//k在端点处，只能用端点和右边所有部分合并 { f[i][j][o]=min(f[i][j][o],f[i+1][j][o]+m+sum[j]-sum[i-1]); if(ch[i]>=o) f[i][j][o]=min(f[i][j][o],f[i][i][ch[i]]+f[i+1][j][ch[i]]+sum[j]-sum[i-1]-d[a[i]]); } else if(k==j)//同理 { f[i][j][o]=min(f[i][j][o],f[i][j-1][o]+m+sum[j]-sum[i-1]); if(ch[j]>=o) f[i][j][o]=min(f[i][j][o],f[i][j-1][ch[j]]+f[j][j][ch[j]]+sum[j]-sum[i-1]-d[a[j]]); } else{//正宗转移方程 f[i][j][o]=min(f[i][j][o],f[i][k-1][o]+f[k+1][j][o]+m+sum[j]-sum[i-1]); if(ch[k]>=o) f[i][j][o]=min(f[i][j][o],f[i][k-1][ch[k]]+f[k+1][j][ch[k]]+sum[j]-sum[i-1]); } } } } printf("%lld",f[1][n][1]); return 0; }
太恶心了。为了保证l<=r，做出了巨大的讨论。
其实不用这么麻烦，只要让l>r的时候，赋值为0就好，相当于不存在。根本不影响答案。
代码2（化简）
#include<bits/stdc++.h> using namespace std; typedef long long ll; const int N=77; const ll inf=2e18; ll f[N][N][N]; int n; ll m; int a[N]; int tot; int w[N],p[N],d[N]; int prio[N]; ll sum[N]; int ch[N];//离散化后的优先级 bool cmp(int a,int b) { return w[a]<w[b]; } int main() { scanf("%d%lld",&n,&m); for(int i=1;i<=n;i++) scanf("%d",&w[i]),a[i]=i; for(int i=1;i<=n;i++) scanf("%d",&p[i]),prio[i]=p[i]; for(int i=1;i<=n;i++) scanf("%d",&d[i]); sort(a+1,a+n+1,cmp); sort(prio+1,prio+n+1); for(int i=1;i<=n;i++) ch[i]=lower_bound(prio+1,prio+n+1,p[a[i]])-prio; for(int i=1;i<=n;i++) for(int j=i;j<=n;j++) for(int o=0;o<=n;o++) f[i][j][o]=inf; for(int i=1;i<=n;i++) for(int o=0;o<=n;o++) f[i][i-1][o]=0;//其实这步不需要，因为上面就没有给它赋值，只是在这里强调一下。 for(int i=1;i<=n;i++) sum[i]=sum[i-1]+d[a[i]]; for(int o=n;o>=1;o--) for(int i=n;i>=1;i--) for(int j=i;j<=n;j++) for(int k=i;k<=j;k++) { f[i][j][o]=min(f[i][j][o],f[i][k-1][o]+f[k+1][j][o]+m+sum[j]-sum[i-1]); if(ch[k]>=o) f[i][j][o]=min(f[i][j][o],f[i][k-1][ch[k]]+f[k+1][j][ch[k]]+sum[j]-sum[i-1]); }//不放心，可以考虑代入长度小于等于2的情况。0的作用就出来了。 //连初始化l=1都省了。 printf("%lld",f[1][n][1]); return 0; }
总结：
1.对于琢磨不透的变化，一定有不变的东西。一定要抓住其中的不变量，作为突破口。
2.循环顺序要注意，一个是不能有后效性，一个是要保证能影响到这个状态的所有状态都处理完了。
3.注意考虑清楚所有可能转移的方式。
