[bzoj2809][主席树]dispatching

Description

在一个忍者的帮派里，一些忍者们被选中派遣给顾客，然后依据自己的工作获取报偿。在这个帮派里，有一名忍者被称之为 Master。除了
Master以外，每名忍者都有且仅有一个上级。为保密，同时增强忍者们的领导力，所有与他们工作相关的指令总是由上级发送给他的直接下属，而不允许通过其他的方式发送。现在你要招募一批忍者，并把它们派遣给顾客。你需要为每个被派遣的忍者
支付一定的薪水，同时使得支付的薪水总额不超过你的预算。另外，为了发送指令，你需要选择一名忍者作为管理者，要求这个管理者可以向所有被派遣的忍者
发送指令，在发送指令时，任何忍者（不管是否被派遣）都可以作为消息的传递
人。管理者自己可以被派遣，也可以不被派遣。当然，如果管理者没有被排遣，就不需要支付管理者的薪水。你的目标是在预算内使顾客的满意度最大。这里定义顾客的满意度为派遣的忍者总数乘以管理者的领导力水平，其中每个忍者的领导力水平也是一定的。写一个程序，给定每一个忍者
i的上级 Bi，薪水Ci，领导力L i，以及支付给忍者们的薪水总预算 M，输出在预算内满足上述要求时顾客满意度的最大值。

1 ≤N ≤ 100,000 忍者的个数； 1 ≤M ≤ 1,000,000,000 薪水总预算； 0 ≤Bi < i
忍者的上级的编号； 1 ≤Ci ≤ M 忍者的薪水； 1 ≤Li ≤ 1,000,000,000
忍者的领导力水平。

Input

从标准输入读入数据。 第一行包含两个整数 N和 M，其中 N表示忍者的个数，M表示薪水的总预算。 接下来
N行描述忍者们的上级、薪水以及领导力。其中的第 i 行包含三个整 Bi , C i , L
i分别表示第i个忍者的上级，薪水以及领导力。Master满足B i = 0，并且每一个忍者的老板的编号一定小于自己的编号 Bi < i。

Output

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

Sample Input

5 4

0 3 3

1 3 5

2 2 2

1 2 4

2 3 1

Sample Output

6

HINT

如果我们选择编号为 1的忍者作为管理者并且派遣第三个和第四个忍者，薪水总和为 4，没有超过总预算
4。因为派遣了 2 个忍者并且管理者的领导力为 3，

用户的满意度为 2 ，是可以得到的用户满意度的最大值。

题解

很容易发现忍者构成了一棵树的结构
那么当一个忍者为管理者时，能派出的只有他这棵子树的忍者
由于答案为派出的忍者数*管理者领导力，所以派出的肯定是这棵子树里总花费最少的忍者
在主席树上二分即可

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
struct edge{int x,y,next;}a[210000];int len,last[110000];
void ins(int x,int y){len++;a[len].x=x;a[len].y=y;a[len].next=last[x];last[x]=len;}
struct LSnode{LL y;int p;}w[110000];int Rank[110000],gg[110000],S[110000];
bool cmp(LSnode n1,LSnode n2){return n1.y<n2.y;}
struct node
{
    int lc,rc,c;LL s;
}tr[4110000];int trlen,ROOT;
int rt[110000];
void add(int &now,int l,int r,int p,LL ad)
{
    if(now==0)now=++trlen;
    tr[now].c++;tr[now].s+=ad;
    if(l==r)return ;
    int mid=(l+r)/2;
    if(p<=mid)add(tr[now].lc,l,mid,p,ad);
    else add(tr[now].rc,mid+1,r,p,ad);
}
void merge(int &x,int y)
{
    if(x==0){x=y;return ;}
    if(y==0)return ;
    tr[x].c+=tr[y].c;tr[x].s+=tr[y].s;
    merge(tr[x].lc,tr[y].lc);
    merge(tr[x].rc,tr[y].rc);
}
LL ans,m;int n;
LL col[110000],sum;
int gtans(int x,int l,int r)
{
    if(tr[x].s+sum<=m)return tr[x].c;
    if(l==r)return (m-sum)/S[l];
    int lc=tr[x].lc,rc=tr[x].rc;
    int mid=(l+r)/2;
    if(tr[lc].s+sum<=m){sum+=tr[lc].s;return gtans(rc,mid+1,r)+tr[lc].c;}
    else return gtans(lc,l,mid);
}
int cnt;
void sol(int x,int fa)
{
    add(rt[x],1,n,Rank[x],w[gg[x]].y);
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(y!=fa)
        {
            sol(y,x);
            merge(rt[x],rt[y]);
        }
    }
    sum=0;
    int tmp=gtans(rt[x],1,n);
    ans=max(ans,(LL)col[x]*tmp);
}
int main()
{
    scanf("%d%lld",&n,&m);
    len=0;memset(last,0,sizeof(last));
    for(int i=1;i<=n;i++)
    {
        int op;
        scanf("%d%lld%lld",&op,&w[i].y,&col[i]);w[i].p=i;
        if(op==0)ROOT=i;
        else ins(op,i),ins(i,op);
    }
    sort(w+1,w+1+n,cmp);
    int tt=1;Rank[w[1].p]=1;gg[w[1].p]=1;
    for(int i=2;i<=n;i++)
    {
        if(w[i].y!=w[i-1].y)tt++;
        Rank[w[i].p]=tt;gg[w[i].p]=i;S[tt]=w[i].y;
    }
    ans=0;sol(ROOT,0);
    printf("%lld\n",ans);
    return 0;
}