分析:

化简一下我们可以发现,suma*ans=sumb,那么我们考虑二分ans,之后做树形背包上做剪枝。

时间复杂度证明,By GXZlegend O(nklogans)

附上代码:

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <cstring>
#include <queue>
#include <iostream>
using namespace std;
#define N 2505
#define eps 5e-4
#define max(a,b) ((a)<(b)?(b):(a))
double f[N][N],ans;
int a[N],b[N],dep[N],head[N],cnt,k,n,siz[N];
struct node
{
int to,next;
}e[N<<1];
void add(int x,int y)
{
e[cnt].to=y;
e[cnt].next=head[x];
head[x]=cnt++;
return ;
}
void dfs(int x,int from)
{
dep[x]=dep[from]+1;siz[x]=1;
if(x)f[x][1]=-ans*a[x]+b[x];
for(int i=head[x];i!=-1;i=e[i].next)
{
int to1=e[i].to;
if(to1!=from)
{
dfs(to1,x);
int o=min(k,siz[x]+siz[to1]),u=min(k,siz[to1]);
for(int j=o;j>=1;j--)
{
for(int l=max(j-siz[x],1);l<=u;l++)
{
if(l>j)break;
f[x][j]=max(f[x][j],f[x][j-l]+f[to1][l]);
}
}
siz[x]+=siz[to1];
}
}
}
bool check(double x)
{
for(int i=0;i<=n;i++)for(int j=0;j<=k;j++)f[i][j]=-1e9;
dep[0]=f[0][0]=0;
ans=x;
dfs(0,0);
if(f[0][k]>0)return 1;
return 0;
}
int main()
{
// freopen("sales.in","r",stdin);
// freopen("sales.out","w",stdout);
memset(head,-1,sizeof(head));
scanf("%d%d",&k,&n);
for(int i=1;i<=n;i++)
{
int x;
scanf("%d%d%d",&a[i],&b[i],&x);
add(x,i);
}
double l=0,r=10000;
while(l<r-eps)
{
double mid=(l+r)/2;
if(check(mid))l=mid;
else r=mid;
}
printf("%.3lf\n",l);
return 0;
}

其实这种方法就能跑的飞起,虽然是递归的,并且状态和转移比较多。

其实还有别的方法,比如说将它转化为dfs序上做背包。

我们知道,选择一个必须包含根的联通块,我们就可以这样考虑,状态和上面的差不多,f[i][j]表示dfs序上选择第i个点,在i-n中选j个的最大答案

转移:f[i][j]=f[i+1][j-1]+x*a[idx[i]]+b[idx[i]]和f[i][j]=max(f[i][j],f[i+siz[idx[i]][j]);

最后判断f[1][k]是否大于0

附上代码:

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <cstring>
#include <queue>
#include <iostream>
using namespace std;
#define N 2505
#define eps 5e-4
#define max(a,b) ((a)<(b)?(b):(a))
double f[N][N],ans;
int a[N],b[N],head[N],cnt,k,n,idx[N],x,siz[N],tims;
struct node
{
int to,next;
}e[N<<1];
void add(int x,int y)
{
e[cnt].to=y;
e[cnt].next=head[x];
head[x]=cnt++;
return ;
}
void dfs(int x,int from)
{
idx[++tims]=x;siz[x]=1;
for(int i=head[x];i!=-1;i=e[i].next)
{
int to1=e[i].to;
if(to1!=from)
{
dfs(to1,x);
siz[x]+=siz[to1];
}
}
}
bool check(double x)
{
for(int i=1;i<=n+2;i++)
{
for(int j=0;j<=k;j++)
{
f[i][j]=-1e9;
}
}
f[n+2][0]=0;
for(int i=n+1;i>=1;i--)
{
int t=idx[i];
int u=min(k,n+2-i);
if(i!=1)
{
for(int j=u;j>=1;j--)
{
f[i][j]=f[i+1][j-1]-x*a[t]+b[t];
}
}else
{
for(int j=k;j>=1;j--)
{
f[i][j]=f[i+1][j];
}
}
for(int j=u;j>=0;j--)
{
f[i][j]=max(f[i][j],f[i+siz[t]][j]);
}
}
if(f[1][k]>0)return 1;
return 0;
}
int main()
{
// freopen("sales.in","r",stdin);
// freopen("sales.out","w",stdout);
memset(head,-1,sizeof(head));
scanf("%d%d",&k,&n);
for(int i=1;i<=n;i++)
{
int x;
scanf("%d%d%d",&a[i],&b[i],&x);
add(x,i);
}
dfs(0,0);
double l=0,r=10000;
while(l<r-eps)
{
double mid=(l+r)/2;
if(check(mid))l=mid;
else r=mid;
}
printf("%.3lf\n",l);
return 0;
}

  

[Jsoi2016]最佳团体 BZOJ4753 01分数规划+树形背包/dfs序的更多相关文章

  1. LUOGU P4322 [JSOI2016]最佳团体(0/1分数规划+树形背包)

    传送门 解题思路 一道0/1分数规划+树上背包,两个应该都挺裸的,话说我常数为何如此之大..不吸氧洛谷过不了啊. 代码 #include<iostream> #include<cst ...

  2. bzoj 4753: [Jsoi2016]最佳团体【01分数规划+二分+树上背包】

    01分数规划,二分答案然后把判别式变成Σp[i]-Σs[i]*mid>=0,然后树上背包判断,设f[i][j]为在i点子树里选j个的最大收益,随便背包一下就好 最丧病的是神卡常--转移的时候要另 ...

  3. BZOJ.4753.[JSOI2016]最佳团体(01分数规划 树形背包DP)

    题目链接 \(Description\) 每个点有费用si与价值pi,要求选一些带根的连通块,总大小为k,使得 \(\frac{∑pi}{∑si}\) 最大 \(Solution\) 01分数规划,然 ...

  4. bzoj 4753 [Jsoi2016]最佳团体——0/1分数规划

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4753 0/1分数规划裸题. #include<iostream> #includ ...

  5. bzoj 4753 最佳团体 —— 01分数规划+树形背包

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4753 注意赋初值为 -inf: eps 设为 1e-3 会 WA ... 代码如下: #in ...

  6. 【bzoj4753】[Jsoi2016]最佳团体 分数规划+树形背包dp

    题目描述 JSOI信息学代表队一共有N名候选人,这些候选人从1到N编号.方便起见,JYY的编号是0号.每个候选人都由一位编号比他小的候选人Ri推荐.如果Ri=0则说明这个候选人是JYY自己看上的.为了 ...

  7. 【BZOJ】4753: [Jsoi2016]最佳团体 01分数规划+树上背包

    [题意]n个人,每个人有价值ai和代价bi和一个依赖对象ri<i,选择 i 时 ri 也必须选择(ri=0时不依赖),求选择k个人使得Σai/Σbi最大.n<=2500,ai,bi< ...

  8. P1642 规划 01分数规划+树形DP

    $ \color{#0066ff}{ 题目描述 }$ 某地方有N个工厂,有N-1条路连接它们,且它们两两都可达.每个工厂都有一个产量值和一个污染值.现在工厂要进行规划,拆除其中的M个工厂,使得剩下的工 ...

  9. BZOJ 4753 [Jsoi2016]最佳团体 ——01分数规划 树形DP

    要求比值最大,当然用分数规划. 二分答案,转化为选取一个最大的联通块使得它们的和大于0 然后我们直接DP. 复杂度$O(n^2\log {n})$ #include <map> #incl ...

随机推荐

  1. vue 结合 animate.css

    这里说的是vue2.0和animate.css的结合使用.其实用过就知道用法是比较简单的.但是就是刚开始使用的时候,难免有的会遇到各种问题.简单的说说我所用过并且遇过的坑. 首先是transition ...

  2. JS获取填报扩展单元格控件的值

    1. 问题描述 填报预览时,我们想获取到某个控件的值相对来说较容易.但如果控件是扩展的,就只能获取到第一个值,无法根据扩展一行行获取对应的值. 例:本意是想获取到袁成洁,结果还是获取到第一个单元格值孙 ...

  3. CVE-2018-18820 icecast 栈缓冲区越界写漏洞分析

    前言 icecast 是一款开源的流媒体服务器 , 当服务器配置了 url 认证时,服务器在处理 HTTP 头部字段时错误的使用了 snprintf 导致栈缓冲区的越界写漏洞( CVE-2018-18 ...

  4. nginx www解析失败问题解决

    nginx www解析失败: nginx代理IIS下域名时 xxxx.xxx可以解析 但www.xxxx.xxx解析失败 IIS增加ip解析:配置下127.0.0.1就可以解析了.

  5. Android 经典欧美小游戏 guess who

    本来是要做iOS开发的,因为一些世事无常和机缘巧合与测试工作还有安卓系统结下了不解之缘,前不久找到了guess who 源码,又加入了一些自己的元素最终完成了这个简单的小游戏. <?xml ve ...

  6. beego快速入门

    beego的官方网址:https://beego.me 参考文档:https://beego.me/quickstart 1:安装 您需要安装 Go 1.1+ 以确保所有功能的正常使用. 需要已经设置 ...

  7. windows环境下 nginx+iis 反向代理解决跨域问题

    项目基本完成,是时候花点时间整理一下最近的姿势了 1 什么是跨域? 网上对于跨域的概念会有大篇幅的文章去解释,似乎有点玄乎,初学者很容易对这个概念产生恐惧,跨域其实很简单,其实只要知道一点,无法跨域访 ...

  8. windows10下安装source insight 4.0(破解版)

    1.从官网下载source insight4.0版本(不用下载,在后面已经把所有需要的文件都准备好了); 2.安装source insightt4.0; 3.使用下载好的sourceinsight4. ...

  9. sysbench压力测试工具安装及使用

    使用sysbench指定连接的时候不需要写上mysql-socket参数 如果自己使用sysbench来测试mysql性能的话,建议使用源码编译比较好,版本为0.4版本. 步骤是: .yum inst ...

  10. ansible 增加本机/etc/hosts 下hostsname 与IP

    --- - hosts: all vars: IP: "{{ ansible_eth0['ipv4']['address'] }}" tasks: - name: 将原有的host ...