前言:
本算法是树剖DDP算法的加速版,复杂度更小(不是特意卡树剖其实看不出来),
但还是比较好写,因此在这里提一下
动态DP
动态DP听起来很黑科技,但其实并不难
首先,能动态的DP本身就很少,需要满足很多限制
而且一般动态DP的题都会给普通DP留部分分,如果有暴力写得好的,甚至可以在低风险下拿到和DDP差不多的分。
扯回正题,动态DP就是把常规DP的转移用矩阵(其实不一定是矩阵,用数据结构特有的维护方式也行,如线段树的区间标记什么的),使得转移可并,然后就能快速地用矩阵加速or树剖解决。
全局平衡二叉树
这玩意感觉没什么用,不过代码挺好写的可以学学
首先,树剖的复杂度是 O ( N l o g 2 N ) O(Nlog^2N) O(Nlog2N)的,这在一些题目中会被卡
其慢就慢在重链有 l o g N logN logN条,每条上面都是 l o g N logN logN的复杂度处理。
所以,其效率的瓶颈就在于重链的处理上。
那么,全局平衡二叉树,就是用于优化重链的处理:
对于一条重链,把它按照划分重心的方式,建一颗二叉树。
换言之,第一层的根节点,是整个树链的重心(这里计算重心要包括轻子树的大小,所以是带权重心),然后它的左右儿子是它左右两个子链的重心……
这样,可以发现,对于任意一个点,其上方所有重链的深度之和是logN的!换言之,从任意一个点出发,一定能经过logN条轻边,logN条重链上的二叉边,到达根节点。
这个很显然,就不证明了
板子题
洛谷 P4719
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define SF scanf
#define PF printf
#define MAXN 100010
#define INF 0x3FFFFFFF
using namespace std;
int n,m;
struct Matrix{
int num[2][2];
Matrix operator *(const Matrix &a) const {
Matrix tmp;
tmp.num[0][0]=tmp.num[0][1]=tmp.num[1][0]=tmp.num[1][1]=-INF;
for(int k=0;k<2;k++)
for(int i=0;i<2;i++)
for(int j=0;j<2;j++)
tmp.num[i][j]=max(tmp.num[i][j],num[