BZOJ2809：[Apio2012]dispatching—

http://www.lydsy.com/JudgeOnline/problem.php?id=2809

题面复制于：https://www.luogu.org/problemnew/show/1552

题目背景

在一个忍者的帮派里，一些忍者们被选中派遣给顾客，然后依据自己的工作获取报偿。

题目描述

在这个帮派里，有一名忍者被称之为Master。除了Master以外，每名忍者都有且仅有一个上级。为保密，同时增强忍者们的领导力，所有与他们工作相关的指令总是由上级发送给他的直接下属，而不允许通过其他的方式发送。

现在你要招募一批忍者，并把它们派遣给顾客。你需要为每个被派遣的忍者支付一定的薪水，同时使得支付的薪水总额不超过你的预算。另外，为了发送指令，你需要选择一名忍者作为管理者，要求这个管理者可以向所有被派遣的忍者发送指令，在发送指令时，任何忍者（不管是否被派遣）都可以作为消息的传递人。管理者自己可以被派遣，也可以不被派遣。当然，如果管理者没有被排遣，你就不需要支付管理者的薪水。

你的目标是在预算内使顾客的满意度最大。这里定义顾客的满意度为派遣的忍者总数乘以管理者的领导力水平，其中每个忍者的领导力水平也是一定的。

写一个程序，给定每一个忍者i的上级Bi，薪水Ci，领导力Li，以及支付给忍者们的薪水总预算M，输出在预算内满足上述要求时顾客满意度的最大值。

输入输出格式

输入格式：

第一行包含两个整数N和M，其中N表示忍者的个数，M表示薪水的总预算。

接下来N行描述忍者们的上级、薪水以及领导力。其中的第i行包含三个整数Bi,Ci,Li分别表示第i个忍者的上级，薪水以及领导力。Master满足Bi=0，并且每一个忍者的老板的编号一定小于自己的编号Bi<i。

输出格式：

输出一个数，表示在预算内顾客的满意度的最大值。

输入输出样例

输入样例#1：

输出样例#1：

——————————————————————————————

我们将关系图建出来，发现我们需要找到最大的子树的根节点L*选取的点的个数，同时使得费用不超过m。

显然对于同一棵子树，我们取相同个数的节点，显然最优的就是取费用最小的点，去掉费用大的点。

显然我们可以用左偏树维护子树中最大费用节点，同时我们还可以顺便维护子树的节点数。

好了做完了……

但是细节很多让我们慢慢讲：

1.方便起见，我们按照建好的关系树从下往上造左偏树。

2.每当树的费用大于m的时候将根节点剔除（左右节点，权值，节点数=0），将其左右子树合并，再将原来根节点合并进来（这是为了方便我们操作，这样我们对于一个节点i，我们直接tr[find(i)].sum*l[i]就是结果，不需要别的操作）

3.那么我们就需要维护每个节点在左偏树的爸爸了，请注意使用路径压缩。

#include<cstdio>

#include<cmath>

#include<algorithm>

#include<cstring>

#include<iostream>

#include<cctype>

#include<queue>

using namespace std;

const int N=1e5+;

typedef long long ll;

inline ll read(){

    ll X=,w=;char ch=;

    while(!isdigit(ch)){w|=ch=='-';ch=getchar();}

    while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();

    return w?-X:X;

}

struct node{

    int to;

    int nxt;

}edge[N];

int cnt=,head[N];

void add(int u,int v){

    cnt++;

    edge[cnt].to=v;

    edge[cnt].nxt=head[u];

    head[u]=cnt;

    return;

}

int fa[N],root;

ll n,m,l[N];

struct tree{

    int l,r,f;

    ll dis,val,size,sum;

}tr[N];

int merge(int x,int y){

    if(x==||y==)return x+y;

    if(tr[x].val<tr[y].val)

    swap(x,y);

    tr[x].r=merge(tr[x].r,y);

    tr[tr[x].r].f=x;

    if(tr[tr[x].l].dis<tr[tr[x].r].dis)

    swap(tr[x].l,tr[x].r);

    tr[x].dis=tr[tr[x].r].dis+;

    tr[x].size=(tr[x].l?tr[tr[x].l].size:)+(tr[x].r?tr[tr[x].r].size:)+(tr[x].size?:);

    tr[x].sum=(tr[x].l?tr[tr[x].l].sum:)+(tr[x].r?tr[tr[x].r].sum:)+tr[x].val;

    return x;

}

int q[N],r;

inline int find(int x){

    if(!tr[x].f)return x;

    return tr[x].f=find(tr[x].f);

}

ll solve(){

    ll ans=;

    q[r=]=root;

    for(int left=;left<=r;left++){

    for(int k=head[q[left]];k;k=edge[k].nxt){

        int v=edge[k].to;

        r++;q[r]=v;

    }

    }

    for(int i=r;i>=;i--){

    int x=q[i],y=fa[x];

    ans=max(ans,tr[find(x)].size*l[x]);

    x=merge(find(x),find(y));

    while(tr[x].sum>m){

        tr[x].val=tr[x].sum=tr[x].size=;

        int ls=tr[x].l,rs=tr[x].r;

        tr[ls].f=tr[rs].f=tr[x].l=tr[x].r=;

        x=merge(x,merge(ls,rs));

    }

    }

    ans=max(ans,tr[find(root)].size*l[root]);

    return ans;

}

void reset(int x,ll c){

    tr[x].val=tr[x].sum=c;

    tr[x].r=tr[x].l=tr[x].dis=tr[x].f=;

    tr[x].size=;

    return;

}

int main(){

    n=read();

    m=read();

    tr[].dis=-;

    for(int i=;i<=n;i++){

    fa[i]=read();

    reset(i,read());

    l[i]=read();

    if(!fa[i])root=i;

    else add(fa[i],i);

    }

    printf("%lld\n",solve());

    return ;

}