[Jsoi2016]最佳团体 BZOJ4753 01分数规划+树形背包/dfs序
分析:
化简一下我们可以发现,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序的更多相关文章
- LUOGU P4322 [JSOI2016]最佳团体(0/1分数规划+树形背包)
传送门 解题思路 一道0/1分数规划+树上背包,两个应该都挺裸的,话说我常数为何如此之大..不吸氧洛谷过不了啊. 代码 #include<iostream> #include<cst ...
- bzoj 4753: [Jsoi2016]最佳团体【01分数规划+二分+树上背包】
01分数规划,二分答案然后把判别式变成Σp[i]-Σs[i]*mid>=0,然后树上背包判断,设f[i][j]为在i点子树里选j个的最大收益,随便背包一下就好 最丧病的是神卡常--转移的时候要另 ...
- BZOJ.4753.[JSOI2016]最佳团体(01分数规划 树形背包DP)
题目链接 \(Description\) 每个点有费用si与价值pi,要求选一些带根的连通块,总大小为k,使得 \(\frac{∑pi}{∑si}\) 最大 \(Solution\) 01分数规划,然 ...
- bzoj 4753 [Jsoi2016]最佳团体——0/1分数规划
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4753 0/1分数规划裸题. #include<iostream> #includ ...
- bzoj 4753 最佳团体 —— 01分数规划+树形背包
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4753 注意赋初值为 -inf: eps 设为 1e-3 会 WA ... 代码如下: #in ...
- 【bzoj4753】[Jsoi2016]最佳团体 分数规划+树形背包dp
题目描述 JSOI信息学代表队一共有N名候选人,这些候选人从1到N编号.方便起见,JYY的编号是0号.每个候选人都由一位编号比他小的候选人Ri推荐.如果Ri=0则说明这个候选人是JYY自己看上的.为了 ...
- 【BZOJ】4753: [Jsoi2016]最佳团体 01分数规划+树上背包
[题意]n个人,每个人有价值ai和代价bi和一个依赖对象ri<i,选择 i 时 ri 也必须选择(ri=0时不依赖),求选择k个人使得Σai/Σbi最大.n<=2500,ai,bi< ...
- P1642 规划 01分数规划+树形DP
$ \color{#0066ff}{ 题目描述 }$ 某地方有N个工厂,有N-1条路连接它们,且它们两两都可达.每个工厂都有一个产量值和一个污染值.现在工厂要进行规划,拆除其中的M个工厂,使得剩下的工 ...
- BZOJ 4753 [Jsoi2016]最佳团体 ——01分数规划 树形DP
要求比值最大,当然用分数规划. 二分答案,转化为选取一个最大的联通块使得它们的和大于0 然后我们直接DP. 复杂度$O(n^2\log {n})$ #include <map> #incl ...
随机推荐
- 【代码笔记】iOS-TableViewOfTwoSecton
一,效果图. 二,工程图. 三,代码. RootViewController.h #import <UIKit/UIKit.h> @interface RootViewController ...
- AngularJS学习之 登录表单 清爽验证(边学边更新)
注册过程的确好多需要验证的,但是注册成功之后的登录就简单多了, 只要用户 输入 用户名和密码, ajax向后台提交登录请求, 根据返回的结果确定用户名或者密码是否正确即可登录. 所以这个登录表单的验证 ...
- 高德地图 JS API - 根据地名实现标记定位
德地图 JS API 使用前的准备工作请参考官方网站说明: https://lbs.amap.com/api/javascript-api/guide/abc/prepare 根据地名实现地图标记定位 ...
- MySql UNIX_TIMESTAMP和FROM_UNIXTIME函数讲解
MySql UNIX_TIMESTAMP和FROM_UNIXTIME函数讲解 by:授客 QQ:1033553122 1. unix_timestamp(date)将时间转换为时间戳,如果参数为空,则 ...
- 树莓派上启动nfs server
1. nfs 是什么 (略)http://vbird.dic.ksu.edu.tw/linux_server/linux_redhat9/0330nfs.php 2. 安装 nfs-kernel-se ...
- textbox只允许输入数字
private void txtUserId_KeyPress(object sender, KeyPressEventArgs e) { //如果输入的不是数字键,也不是回车键.Backspace键 ...
- TCP的socket资源被耗尽的问题
一. 故障现象 部分机顶盒用户出现大面积登录APP时,界面停留在登陆页面,无反应. 二. 现象初步分析 本次问题出现时,所有AAA出现了异常流量波动,在AAA异常流量段期间接到用户故障报障.此时主要表 ...
- 【转】Redis学习---NoSQL和SQL的区别及使用场景
什么是NoSQL NoSQL,指的是非关系型的数据库.NoSQL有时也称作Not Only SQL的缩写,是对不同于传统的关系型数据库的数据库管理系统的统称,它具有非关系型.分布式.不提供ACID的数 ...
- MFC 常用功能总结
目录 ◆ 获取启动参数◆ 调整窗体大小◆ 关闭Dialog窗体◆ Dialog启动时指定控件为焦点◆ Dialog中禁止ESC回车关闭窗体◆ 回车后焦点自动跳到下一个控件◆ 窗体中某个控件捕获右键菜单 ...
- Geth命令用法-参数详解 and 以太坊源码文件目录
本文是对以太坊客户端geth命令的解析 命令用法 geth [选项] 命令 [命令选项] [参数-] 版本 1.7.3-stable 命令 account 管理账户 attach 启动交互式JavaS ...