博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ 4753 [Jsoi2016]最佳团体 ——01分数规划 树形DP
阅读量:6701 次
发布时间:2019-06-25

本文共 1250 字,大约阅读时间需要 4 分钟。

要求比值最大,当然用分数规划。

二分答案,转化为选取一个最大的联通块使得它们的和大于0

然后我们直接DP。

复杂度$O(n^2\log {n})$

#include #include 
#include
#include
#include
#include
#include
using namespace std;#define F(i,j,k) for (int i=j;i<=k;++i)#define D(i,j,k) for (int i=j;i>=k;--i)#define ll long long#define inf 1e15#define eps 1e-6#define mp make_pair#define maxn 2505int k,n;double s[maxn],p[maxn],l,r,mid,dp[maxn][maxn],tmp[maxn];int h[maxn],to[maxn],ne[maxn],en=0,siz[maxn]; void add(int a,int b){to[en]=b;ne[en]=h[a];h[a]=en++;} void Tree_DP(int o){ dp[o][1]=p[o]-mid*s[o];siz[o]=1; for (int i=h[o];i>=0;i=ne[i]) { Tree_DP(to[i]); F(j,0,siz[o]+siz[to[i]]) tmp[j]=-inf; F(j,1,siz[o]) F(k,0,siz[to[i]]) tmp[j+k]=max(tmp[j+k],dp[o][j]+dp[to[i]][k]); F(j,0,siz[o]+siz[to[i]]) dp[o][j]=max(dp[o][j],tmp[j]); siz[o]+=siz[to[i]]; }} int main(){ memset(h,-1,sizeof h); scanf("%d%d",&k,&n);k++; F(i,1,n) { int x; scanf("%lf%lf%d",&s[i],&p[i],&x); add(x,i); } l=0;r=1e10; while (fabs(r-l)>eps) { mid=(l+r)/2; F(i,0,n) F(j,0,k) dp[i][j]=-inf; s[0]=1;p[0]=mid; Tree_DP(0); int flag=0; if (dp[0][k]>=0) l=mid; else r=mid; } printf("%.3f\n",l);}

  

转载于:https://www.cnblogs.com/SfailSth/p/6736890.html

你可能感兴趣的文章
去除vue项目中的#及其ie9兼容性
查看>>
linux实例 批量修改图片文件名
查看>>
day15(mysql 的多表查询,事务)
查看>>
IOS
查看>>
beta冲刺第三天
查看>>
beta第二天
查看>>
工作中常用的但是又容易忽略的问题
查看>>
css特殊情况
查看>>
洛谷 P2735 电网 Electric Fences Label:计算几何--皮克定理
查看>>
Vue + webpack 项目配置化、接口请求统一管理
查看>>
关于电脑的基础单词笔记
查看>>
安卓App设计博文
查看>>
11.8 开课二个月零四天 (Jquery)
查看>>
ZEN CART 在LINUX系统下设置邮箱方法---用GMAIL设置,方法选择SMTPAUTH
查看>>
ofstream的使用方法--超级精细。C++文件写入、读出函数(转)
查看>>
DOM剪切板
查看>>
10.高效分布
查看>>
装饰器概念及运用
查看>>
CCF 201312-4 有趣的数
查看>>
screenX clientX pageX的区别
查看>>