- 思路:二分答案+动态规划(结合dfs序)
- 类型:选/不选:最大比值
- 代码:
#include<stdio.h>
#include<bits/stdc++.h>
using namespace std;
const int N=2505;
int rt=0,k,n,f[N];
double eps=1e-6,v[N],s[N],p[N],dp[N][N],sum;
int tot,head[N],nxt[N*2],to[N*2],size[N],Time,od[N];
void add_edge(int u,int v) {
tot++; nxt[tot]=head[u]; to[tot]=v; head[u]=tot;
}
void dfs(int u) {
size[u]=1;
od[Time++]=u;
for(int i=head[u];i;i=nxt[i]) {
int v=to[i];
dfs(v);
size[u]+=size[v];
}
}
bool check(double ans) {
for(int i=1;i<=n;i++) v[i]=p[i]-s[i]*ans;
for(int i=1;i<=n+1;i++) for(int j=0;j<=k;j++) dp[i][j]=-sum;
dp[1][0]=0.0;
for(int i=1;i<=n;i++) {
for(int j=0;j<=min(i,k);j++) {
dp[i+1][j+1]=max(dp[i+1][j+1],dp[i][j]+v[od[i]]);
dp[i+size[od[i]]][j]=max(dp[i+size[od[i]]][j],dp[i][j]);
}
}
if(dp[1+n][k]>=0) return true;
return false;
}
double solve(double l,double r) {
double ans,mid;
while(r-l>=eps) {
mid=(l+r)/2;
if(check(mid)) {
ans=mid; l=mid;
}
else r=mid;
}
return ans;
}
int main() {
scanf("%d%d",&k,&n);
for(int i=1;i<=n;i++) {
scanf("%lf%lf%d",&s[i],&p[i],&f[i]);
sum+=p[i];
add_edge(f[i],i);
}
dfs(0);
double ans=solve(0,sum);
printf("%.3lf",ans);
return 0;
}
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
· 基于Microsoft.Extensions.AI核心库实现RAG应用
· Linux系列:如何用heaptrack跟踪.NET程序的非托管内存泄露
· 开发者必知的日志记录最佳实践
· SQL Server 2025 AI相关能力初探
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· winform 绘制太阳,地球,月球 运作规律
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 上周热点回顾(3.3-3.9)
· 超详细:普通电脑也行Windows部署deepseek R1训练数据并当服务器共享给他人