指令总预算BZOJ 2809([Apio2012]dispatching-Splay启发式合并)

PS:今天上午，非常郁闷，有很多简单基础的问题搞得我有些迷茫，哎，代码几天不写就忘。目前又不当COO，还是得用心记代码哦！

    

    

2809: [Apio2012]dispatching

    Time Limit: 

    10 Sec  

    Memory Limit: 

    128 MB

    Submit: 

    606  

    Solved: 

    315

    [

    Submit

    ][

    Status

    ][

    Discuss

    ]

    

    

Description

    

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

 

1  ≤N ≤ 100,000 忍者的个数；

1  ≤M ≤ 1,000,000,000 薪水总预算； 

 

0  ≤B i < i  忍者的下级的编号；

1  ≤Ci ≤ M                     忍者的薪水；

1  ≤Li ≤ 1,000,000,000             忍者的领导力水平。

 

 

 

    

Input

    

从标准输入读入数据。

 

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

 

接下来 N行描述忍者们的下级、薪水以及领导力。其中的第 i 行包含三个整  B i , C i , L i分别表示第i个忍者的下级，薪水以及领导力。Master满足B i = 0， 并且每个忍者的老板的编号必定小于自己的编号 B i < 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      ，是可以失掉的用户满意度的最大值。

 

    

Splay启发式合并，注意备注点（重中之重）

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cstring>
#include<cctype>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define Rep(i,n) for(int i=0;i<n;i++)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define MAXN (100000+10)
#define MAXM (1000000000)
int n,m,fa[MAXN];
struct node
{
    int a,fa,lu,siz,ch[2];
    long long sum;
    node(int _a,int _lu,int _fa):a(_a),lu(_lu),siz(1),sum(_a),fa(_fa){ch[0]=ch[1]=0;}
    node():siz(0),sum(0),lu(0),a(0),fa(0){ch[0]=ch[1]=0;}
    void clear(){ch[0]=ch[1]=fa=sum=0;siz=1;sum=a;}
}a[MAXN];
/*
int edge[MAXN],next[MAXN],pre[MAXN],size=0;
void addedge(int u,int v)
{
    edge[++size]=v;
    next[size]=pre[u];
    pre[u]=size;
}*/
void update(int x)
{
    a[x].sum=a[a[x].ch[0]].sum+a[a[x].ch[1]].sum+a[x].a;
    a[x].siz=a[a[x].ch[0]].siz+a[a[x].ch[1]].siz+1;
}
void rotate(int x) //x && x.fa
{
    int y=a[x].fa,z=a[y].fa;
    int p=a[y].ch[0]==x;
    if (z)
        if (a[z].ch[0]==y) a[z].ch[0]=x;
        else a[z].ch[1]=x;
    a[x].fa=z,a[y].fa=x;
    if (a[x].ch[p]) a[a[x].ch[p]].fa=y;
    a[y].ch[p^1]=a[x].ch[p];
    a[x].ch[p]=y;
    update(y);
}
void splay(int x)
{
    while (a[x].fa)
    {
        int y=a[x].fa,z=a[y].fa;
        if (z)
            if ((a[y].ch[0]==x)^(a[z].ch[0]==y)) rotate(x);
            else rotate(y);
        rotate(x);
    }
    update(x);
}
void ins(int x,int y) //a[x] ++a[y]
{
    a[x].sum+=a[y].a,a[x].siz++;
    if (a[y].a<=a[x].a)
    {
        if (a[x].ch[0]) ins(a[x].ch[0],y);
        else a[y].fa=x,splay(a[x].ch[0]=y);
    }
    if (a[x].a<a[y].a)
    {
        if (a[x].ch[1]) ins(a[x].ch[1],y);
        else a[y].fa=x,splay(a[x].ch[1]=y);
    }
}
int q[MAXN];
void merge(int x,int y)
{
    splay(x),splay(y);
    if (a[x].siz<a[y].siz) swap(x,y);
    q[1]=y;
    int head=1,tail=1;
    while (head<=tail)
    {
        int now=q[head];
        if (a[now].ch[0]) q[++tail]=a[now].ch[0];
        if (a[now].ch[1]) q[++tail]=a[now].ch[1];
        a[now].clear();
        ins(x,now);
        x=now;// now 被旋转到根
        head++;
    }
}
long long find(int x,int maxs)
{
    if (a[x].sum<=maxs) return a[x].siz;
    int ans=0;
    if (a[x].ch[0]&&a[a[x].ch[0]].sum>=maxs) return find(a[x].ch[0],maxs);
    else if (a[a[x].ch[0]].sum+a[x].a>=maxs) return a[a[x].ch[0]].siz+bool(a[a[x].ch[0]].sum+a[x].a==maxs);
    else return a[a[x].ch[0]].siz+1+find(a[x].ch[1],maxs-a[a[x].ch[0]].sum-a[x].a);
}
int main()
{
//	freopen("bzoj2809.in","r",stdin);
//	freopen("bzoj2809.out","w",stdout);
    scanf("%d%d",&n,&m);
    For(i,n)
    {
        int b1,a1,l1;
        scanf("%d%d%d",&fa[i],&a1,&l1);
        a[i]=node(a1,l1,0);
    }
    long long ans=-1;
    ForD(x,n)
    {
        splay(x);
        ans=max(ans,(long long)a[x].lu*find(x,m));
    //  cout<<find(x,m)<<endl;
        if (fa[x]) merge(x,fa[x]);
    }
    cout<<ans<<endl;
	return 0;
}

恭喜成功转CB。

文章结束给大家分享下程序员的一些笑话语录： 3G普不普及现在已经不是看终端了，而是看应用，有好的，便宜实用的应用，花1000多买个能用的智能手机应该不是什么难事。反过来说，你200元拿一个智能手机，没有好的应用，看个电影要几十元，也是没人用3G。