P1768-天路【负环,SPFA,01分数规划,二分答案】
正题
题目链接:https://www.luogu.org/problemnew/show/P1768
题目大意
求一条回路使得路上∑vi∑pi\frac{\sum v_i}{\sum p_i}∑pi∑vi最大。
解题思路
考虑01分数规划
∑vi∑pi=ans\frac{\sum v_i}{\sum p_i}=ans∑pi∑vi=ans
∑vi=ans∗∑pi\sum v_i=ans*\sum p_i∑vi=ans∗∑pi
转为二分判定
∑vi−ans∗∑pi>0\sum v_i-ans*\sum p_i>0∑vi−ans∗∑pi>0
ans∗∑pi−∑vi<0ans*\sum p_i-\sum v_i<0ans∗∑pi−∑vi<0
然后对于每次二分,计算边权,然后求负环
codecodecode
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
const int N=7100,M=21000;
struct node{int to,next;double v,p,w;
}a[M*2];
queue<int> q;
double f[N];
int n,m,tot,ls[N],cnt[N];
double l,r;
bool v[N];
void addl(int x,int y,double v,double p)
{a[++tot].to=y;a[tot].next=ls[x];a[tot].v=v;a[tot].p=p;ls[x]=tot;
}
bool spfa()
{for(int i=1;i<=n;i++)f[i]=2147483647/3;memset(cnt,0,sizeof(cnt));q.push(n);cnt[n]=1;f[n]=0;v[n]=1;while(!q.empty()){int x=q.front();v[x]=0;q.pop();for(int i=ls[x];i;i=a[i].next){int y=a[i].to;if(f[x]+a[i].w<f[y]){f[y]=f[x]+a[i].w;cnt[y]=cnt[x]+1;if(cnt[y]>=n&&a[i].w<0)return true;if(!v[y]){v[y]=1;q.push(y);}}} }return false;
}
bool check(double ans){for(int i=1;i<=tot;i++)a[i].w=a[i].p*ans-a[i].v;if(spfa()) return true;else return false;
}
int main()
{scanf("%d%d",&n,&m);for(int i=1;i<=m;i++){int x,y;double v,p;scanf("%d%d%lf%lf",&x,&y,&v,&p);addl(x,y,v,p);}n++;for(int i=1;i<n;i++)addl(n,i,0,0);l=0;r=500;while(r-l>1e-2){double mid=(l+r)/2;if(check(mid)) l=mid;else r=mid;}if(l==0) printf("-1");else printf("%.1lf",l);
}
P1768-天路【负环,SPFA,01分数规划,二分答案】相关推荐
- 算法提高课-图论-负环-AcWing 361. 观光奶牛:spfa判正环、负环、01分数规划、二分
文章目录 题目分析 题目链接 题目分析 来源:acwing 分析: 题目要求ΣfiΣgi\frac{\Sigma{f_i}}{\Sigma{g_i}}ΣgiΣfi的最大值,这种问题称为01分数规 ...
- 解题报告:AcWing 1165. 单词环(01分数规划、hash、经验优化)
本题的关键在于: 建图 01分数规划 本题的数据过大,如果直接spfa判断会TLE,因此我们使用经验优化,就是如果所有的点入队的次数过多,比如大于100000,那么我们直接认为它是存在正环的.(免去T ...
- POJ-3621 Sightseeing Cows 最优比率环、01分数规划
题目链接:http://poj.org/problem?id=3621 这题是01分数规划问题,详细资料点这里.有了01分数规划的基础后,这个题目就很简单了.构建函数f(mid)=Σ(w[i]-mid ...
- [BJOI2019]奥术神杖——AC自动机+DP+分数规划+二分答案
题目链接: [BJOI2019]奥术神杖 答案是$ans=\sqrt[c]{\prod_{i=1}^{c}v_{i}}=(\prod_{i=1}^{c}v_{i})^{\frac{1}{c}}$. 这 ...
- poj 2728 Desert King(最小比率生成树 / 0-1分数规划 / 二分)
二分答案,我们要找最小的答案,如果有更小的答案说明 ∑W−Z∗∑L<=0∑W−Z∗∑L <= 0∑W−Z∗∑L<=0. #include<cstdio> #include ...
- [jzoj 4230] 淬炼神体{ 0/1分数规划+二分答案}
题目 Description 王仙女将你提供的答案填在<葵花宝典>上,突然,宝典发出耀眼的白光,一股强大的吸力瞬间将仙女吸入宝典中. 一阵眩晕过后,仙女发现自己来到了一个浮岛上,四周的半空 ...
- 【图论专题】负环与01分数规划
整理的算法模板合集: ACM模板 题目列表: 题目 算法 AcWing 904. 虫洞 spfa判负环 AcWing 361. 观光奶牛 最优比率环.01分数规划 AcWing 1165. 单词环 0 ...
- P1768 天路(spfa+二分答案+01分数规划)
题目 题目传送门 输入输出样例 输入样例 5 6 1 2 1 1 4 1 6 2 5 4 8 1 2 3 2 2 5 2 4 1 3 5 6 4 输出样例 2.3 题解 本题我们要求max∑vi∑p ...
- 算法提高课-图论-负环-AcWing 1165. 单词环:spfa判正环、二分、01分数规划
文章目录 题目分析 题目链接 题目分析 来源:acwing 分析: 如何建图? 这样建图.以样例举例.起点是前两个字母,终点是末尾两个字母,边权是字符串的长度. 我们求的是什么呢? 题目要求Σ边权Σ1 ...
最新文章
- php编译工具 知乎,关于知乎回答问题编辑框用Ctrl+V 粘贴图片是如何实现的详解...
- /bin/bash: jar: command not found
- mysql约束sex_MySQL笔记--约束
- Thrift入门及Java实例演示
- 数据结构--队列Queue--循环顺序队列
- java list 合并去重复_java 实现多个list 合并成一个去掉重复的案例
- 计算1970距今的秒数
- 黑马博客——详细步骤(四)项目功能的实现之数据分页
- SparkSQL架构
- 获取照片经纬度(wgs84)转高德经纬度(gcj02)
- 广州艺术博物院走进春睡画院旧址“云赏画”
- windows pagefile.sys配置调整
- 第五章 社会存在发展的基础和基本结构
- HDU6411 带劲的and和(2018百度之星复赛,并查集,位运算,思路)
- C++调用PHP程序
- 全国地图省份联动效果
- 函数打桩原理_难重现问题定位“函数打桩”
- Automatic Inference of Search Patterns for Taint-Style Vulnerabilities
- apidoc使用教程-编写漂亮的api文档
- DOSLinux命令大全
热门文章
- java如何用键盘输入_java中如何从键盘输入(附代码)
- matlab提示未定义wc,WooCommerce 教程:修复致命错误调用未定义的函数wc_get_order() - WooCommerce 微站...
- python hexdigest,Java方法,可以为Hex中的HMAC-SHA256提供与Python方法相同的输出
- java原生的编译软件_Java 转原生平台代码 RoboVM
- 两台思科交换机vlan划分_Cisco交换机Vlan划分及ACL配置详细步骤 | 吴文辉博客
- java 1..0 openjdk_java-1.7.0-openjdk-i386和java-7-openjdk-i386有什么区别
- leetcode406. 根据身高重建队列
- [蓝桥杯][2013年第四届真题]剪格子-dfs
- AcWing 523. 组合数问题
- php fpm 调试模式,调试 – nginx php-fpm xdebug netbeans只能启动一个调试会话