bzoj4753 [Jsoi2016]最佳团体

Description

---

JSOI信息学代表队一共有N名候选人，这些候选人从1到N编号。方便起见，JYY的编号是0号。每个候选人都由一位
编号比他小的候选人Ri推荐。如果Ri=0则说明这个候选人是JYY自己看上的。为了保证团队的和谐，JYY需要保证，
如果招募了候选人i，那么候选人Ri”也一定需要在团队中。当然了，JYY自己总是在团队里的。每一个候选人都有
一个战斗值Pi”，也有一个招募费用Si”。JYY希望招募K个候选人（JYY自己不算），组成一个性价比最高的团队。
也就是，这K个被JYY选择的候选人的总战斗值与总招募总费用的比值最大。

对于100%的数据满足1≤K≤N≤2500,0<”Si,Pi”≤10^4,0≤Ri

Solution

---

居然没人写题解?!
看到分数的形式容易想到二分答案变成判定性问题
若$mid<\frac{\sum a_i}{\sum b_i}$，等价于$\sum a_i-mid\sum b_i>0$，此时mid小于所求的答案

那么可以把mid套进柿子里做树形dp求最大值，设f[i,j]表示以i为根的子树中选了j个的最大贡献，转移的时候强制让i选自己转移就能保证所选结果是合法的（具体可以看代码）

需要注意数组的大小，开太大容易T

如果发现和洛谷的题解撞了，那是我写的呀

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))

const double EPS=1e-4;
const int N=2505;
const int E=5005;

struct edge {int x,y,next;} e[E];

double f[N][N],tmp[N],v[N];
int size[N],a[N],b[N];
int ls[N],edCnt;

void add_edge(int x,int y) {
    e[++edCnt]=(edge) {x,y,ls[x]}; ls[x]=edCnt;
    e[++edCnt]=(edge) {y,x,ls[y]}; ls[y]=edCnt;
}

void dfs(int now,int fa) {
    size[now]=1; f[now][0]=0; f[now][1]=v[now];
    for (int i=ls[now];i;i=e[i].next) {
        if (e[i].y==fa) continue;
        dfs(e[i].y,now);
        rep(j,1,size[now]+size[e[i].y]) tmp[j]=f[0][N-1];
        rep(j,1,size[now]) rep(k,0,size[e[i].y]) {
            tmp[j+k]=std:: max(tmp[j+k],f[now][j]+f[e[i].y][k]);
        }
        size[now]+=size[e[i].y];
        rep(j,1,size[now]) f[now][j]=tmp[j];
    }
}

int main(void) {
    int n,m; scanf("%d%d",&m,&n);
    rep(i,1,n) {
        int fa; scanf("%d%d%d",&b[i],&a[i],&fa);
        add_edge(fa,i);
    }
    double l,r;
    for (l=0,r=10000;r-l>=EPS;) {
        double mid=(l+r)*0.5;
        fill(f,0xc2);
        rep(i,1,n) v[i]=(double)a[i]-mid*(double)b[i];
        dfs(0,0);
        if (f[0][m+1]>0) l=mid;
        else r=mid;
    }
    printf("%.3lf\n", l);
    return 0;
}