左偏树

左偏树的定义：左偏树是一种可并堆的实现。左偏树是一棵二叉树，它的节点除了和二叉树的节点一样具有左右子树外，还有两个属性：键值和距离。节点 i 的距离(dis(i)) 是节点 i 到它的后代中，最近的外节点所经过的边数（节点i称为外节点(external node)，当且仅当节点 i 的左子树或右子树为空）。特别的，如果节点 i 本身是外节点，则它的距离为0；而空节点的距离规为-1 。

左偏树的性质：

【性质1】节点的键值满足堆性质。（小顶堆或者大顶堆）

【性质2】节点的左子节点的距离不小于右子节点的距离。即 dis[ls[i]] >= dis[rs[i]]。

【性质3】节点的距离等于它的右子节点的距离加1。即 dis[i] = dis[rs[i]] + 1;

【引理1】若左偏树的距离为一定值，则节点数最少的左偏树为完全二叉树。

【定理1】若一棵左偏树的距离为 k，则这棵左偏树至少有2^(k+1)-1个节点。

【性质4】一棵N个节点的左偏树距离最多为log(N+1) -1。

左偏树的操作：

合并：

1、要合并的两棵树中其中有一棵为空，则返回另一棵树。

2、若A和B都非空，我们假设A的根节点小于等于B的根节点（否则交换A,B），把A的根节点作为新树C的根节点，剩下的事就是合并A的右子树 right(A) 和B了。

3、合并了right(A) 和B之后，right(A) 的距离可能会变大，当right(A) 的距离大于left(A) 的距离时，左偏树的性质2会被破坏。在这种情况下，我们只须要交换left(A) 和right(A)。

4、最后，由于right(A) 的距离可能发生改变，我们必须更新A的距离：dis(A)  = dis(right(A)) + 1。

int Merge(int A, int B){
    if(!A || !B) return A+B;
    if(C[A] < C[B]) swap(A, B);
    rs[A] = Merge(rs[A], B);
    if(dis[ls[A]] < dis[rs[A]]) swap(ls[A], rs[A]);
    dis[A] = dis[rs[A]] + 1;
    return A;
}

删除：合并左右子树就好了。

int Delete(int A){
    return Merge(ls[A], rs[A]);
}

模板：洛谷P1552

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define ll long long
using namespace std;
const int maxn = 100010;
const int maxm = 1000010;
int n, Master, ls[maxn], rs[maxn],dis[maxn];
int head[maxn], tot;
int ver[maxm], Next[maxm];
ll m, C[maxn], L[maxn], sum[maxn], sz[maxn], ans;
void add(int x, int y){
    ver[++tot] = y, Next[tot] = head[x];
    head[x] = tot;
}
int Merge(int A, int B){
    if(!A || !B) return A+B;
    if(C[A] < C[B]) swap(A, B);
    rs[A] = Merge(rs[A], B);
    if(dis[ls[A]] < dis[rs[A]]) swap(ls[A], rs[A]);
    dis[A] = dis[rs[A]] + 1;
    return A;
}
int Delete(int A){
    return Merge(ls[A], rs[A]);
}
int dfs(int x){
    int A = x, B;
    sum[x] = C[x], sz[x] = 1;
    for(int i = head[x]; i; i = Next[i]){
        int y = ver[i];
        B = dfs(y);
        A = Merge(A, B);
        sum[x] += sum[y];
        sz[x] += sz[y];
    }
    while(sum[x] > m){
        sum[x] -= C[A];
        sz[x]--;
        A = Delete(A);
    }
    ans = max(ans, L[x]*sz[x]);
    return A;
}
int main()
{
    int B;
    cin >> n >> m;
    for(int i = 1; i <= n; i++){
        cin >> B >> C[i] >> L[i];
        if(!B) Master = i;
        else add(B, i);
    }
    dfs(Master);
    cout << ans << endl;
    return 0;
}