[bzoj4753]最佳团体

[bzoj4753]最佳团体

---

直接跑dp没法做，你就二分个答案然后dp来确认就行了

• 代码

#include<bits/stdc++.h>
using namespace std;
typedef double db;
typedef long long ll;
const int N=2550;
int n,K;
ll f[N][N];
int fa[N];ll s[N],p[N];int sz[N];
ll v[N];
int csz[N];
bool chk(ll P){
    memset(csz,0,sizeof(csz));
    for(int i=0;i<=n;i++)
        for(int j=1;j<=K;j++)
            f[i][j]=-100000000000000000ll;
    for(int i=1;i<=n;i++){
        v[i]=1ll*p[i]-1ll*s[i]*P;
    }
    for(int i=n;i;i--){
        for(int j=min(sz[i],K);j;j--){
            f[i][j]=f[i][j-1]+v[i];
        }f[i][0]=0;
        int p=fa[i];
        for(int j=min(K,sz[i]+csz[p]);j;j--){
            for(int k=0;k<=min(j,csz[p]);k++){
                f[p][j]=max(f[p][j],f[p][k]+f[i][j-k]);
            }
        }
        csz[p]+=sz[i];
    }
    return f[0][K]>=0;
}
 
int main()
{
    scanf("%d%d",&K,&n);
    for(int i=1;i<=n;++i){
        scanf("%d%d%d",&s[i],&p[i],&fa[i]);
        p[i]*=100000;
    }
    for(int i=n;i;i--){
        sz[i]++;
        sz[fa[i]]+=sz[i];
    }
    int l=1,r=1000000000;
    while(l<r){
        int mid=(l+r+1)>>1;
        if(chk(mid))l=mid;
        else r=mid-1;
    }
    printf("%.3lf\n",(double)l/100000.00);
}
?