http://www.elijahqi.win/archives/3664
Description
JSOI信息学代表队一共有N名候选人,这些候选人从1到N编号。方便起见,JYY的编号是0号。每个候选人都由一位
编号比他小的候选人Ri推荐。如果Ri=0则说明这个候选人是JYY自己看上的。为了保证团队的和谐,JYY需要保证,
如果招募了候选人i,那么候选人Ri”也一定需要在团队中。当然了,JYY自己总是在团队里的。每一个候选人都有
一个战斗值Pi”,也有一个招募费用Si”。JYY希望招募K个候选人(JYY自己不算),组成一个性价比最高的团队。
也就是,这K个被JYY选择的候选人的总战斗值与总招募总费用的比值最大。
Input
输入一行包含两个正整数K和N。
接下来N行,其中第i行包含3个整数Si,Pi,Ri表示候选人i的招募费用,战斗值和推荐人编号。
对于100%的数据满足1≤K≤N≤2500,0< ” Si,Pi”≤10^4,0≤Ri< i
Output
输出一行一个实数,表示最佳比值。答案保留三位小数。
Sample Input
1 2
1000 1 0
1 1000 1
Sample Output
0.001
HINT
2017.9.12新加数据一组 By GXZlegend
Source
分数规划+dfs序 背包 常数大 差点被卡T
二分答案 然后结果就是sigma p[x]-s[x]*mid
设dp[i][j]表示i子树选了j个点的最大代价
具体可上网搜索 依赖型背包
只要看选k个能否>0即可
#include<cmath>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define eps 1e-9
using namespace std;
inline char gc(){
static char now[1<<16],*S,*T;
if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
return *S++;
}
inline int read(){
int x=0,f=1;char ch=gc();
while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
return x*f;
}
const int N=2550;
struct node{
int y,next;
}data[N<<1];
double dp[N][N],mid;int num,h[N],k,n,s[N],p[N];
inline void insert1(int x,int y){
data[++num].y=y;data[num].next=h[x];h[x]=num;
}
inline void dfs(int x,int fa){
if (x) for (int i=k-1;~i;--i) dp[x][i+1]=dp[x][i]+p[x]-mid*s[x];
for (int i=h[x];i;i=data[i].next){
int y=data[i].y;if (x==0) dp[y][0]=0;
for (int j=1;j<=k;++j) dp[y][j]=dp[x][j];dfs(y,x);
}
for (int i=1;i<=k;++i) dp[fa][i]=max(dp[fa][i],dp[x][i]);
}
inline bool check(double md){
memset(dp,0xc2,sizeof(dp));
mid=md;dp[0][0]=0;
dfs(0,0);return dp[0][k]>eps;
}
int main(){
freopen("bzoj4753.in","r",stdin);
k=read();n=read();
for (int i=1;i<=n;++i){
s[i]=read();p[i]=read();int r=read();
insert1(r,i);
}
double l=0,r=2e7;
while(fabs(r-l)>eps){
double md=(l+r)/2;
if (check(md)) l=mid;else r=mid;
}printf("%.3f\n",r);
return 0;
}