正题

题目链接:https://www.luogu.com.cn/problem/P3288


题目大意

给出\(n\)个点\(m\)条边的一张图,没条边\(i\)流量为\(c_i\),费用是\(d_i\),然后缩小一个流量费用是\(a_i\),增加一个流量费用是\(b_i\)。

要求改动图之后最大流不减少

假设减少的费用是\(\Delta X\),改动次数是\(k\),求最大化\(\frac{\Delta X}{k}\)

\(1\leq n\leq 5000,1\leq m\leq 3000\)


解题思路

因为最大流不减少,那么显然因为初始边,最大流也不能增加,所以,每次肯定是选一条回路增流或者退流,这样就是把增流的丢到环上退流的去。

然后对于一条边增流的费用就是\(a_i-d_i\),退流的费用是\(b_i+d_i\)

然后最大化的那个显然是一个分数规划,就直接二分答案然后边权加上答案看有没有负环就好了。


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int N=5e3+10;
struct node{
int to,next;
double w;
}a[N<<1];
int n,m,tot,ls[N],cnt[N];
double f[N];bool v[N];queue<int> q;
void addl(int x,int y,double w){
a[++tot].to=y;
a[tot].next=ls[x];
a[tot].w=w;
ls[x]=tot;return;
}
bool SPFA(double w){
for(int i=1;i<=n+2;i++)f[i]=1e100,cnt[i]=0;
q.push(n+1);f[n+1]=cnt[n+1]=0;
while(!q.empty()){
int x=q.front();q.pop();v[x]=0;
for(int i=ls[x];i;i=a[i].next){
int y=a[i].to;
if(f[x]+a[i].w+w<f[y]){
f[y]=f[x]+a[i].w+w;
cnt[y]=cnt[x]+1;
if(cnt[y]>=n&&a[i].w<0)return 1;
if(!v[y])q.push(y),v[y]=1;
}
}
}
return 0;
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++){
int x,y;double A,B,C,D;
scanf("%d%d",&x,&y);
scanf("%lf%lf%lf%lf",&A,&B,&C,&D);
// if(x==n+1)A=0,B=0,D=0;
if(C>0)addl(y,x,A-D);
addl(x,y,B+D);
}
double l=0,r=1e8;
for(int i=1;i<=100;i++){
double mid=(l+r)/2.0;
if(SPFA(mid))l=mid;
else r=mid;
}
printf("%.2lf\n",(l+r)/2.0);
return 0;
}

最新文章

  1. 2016-1-28 图解HTTP(03)
  2. 精简的网站reset 和 css通用样式库
  3. [CrunchBang]中文字体美化
  4. [每日一题] OCP1z0-047 :2013-08-26 TIMESTAMP WITH LOCAL TIME ZONE....................112
  5. tarjan 算法讲解(转)
  6. PHP 8大安全函数
  7. angular.js升序降序过滤器
  8. Mac系统Git生成ssh公钥
  9. 使用JDBC连接操作数据库
  10. Perl IO:read()函数
  11. opencv学习之路(40)、人脸识别算法——EigenFace、FisherFace、LBPH
  12. PTA天梯 L3-007 天梯地图
  13. 【python练习题】程序1
  14. Entity Framework定义外键,限制通过migration命令自动更改字段名称
  15. Spring-Boot devtools项目自动重启
  16. tomcat 控制台乱码问题
  17. CSRF攻击的应对之道
  18. 单击行,自己主动选中当前行中的单选框button
  19. 关于Unity中的transform组件(二)
  20. [转]你真的了解setTimeout和setInterval吗?

热门文章

  1. mongoTemplate 条件查询
  2. 在C#中使用C++编写的类——用托管C++进行封装
  3. 梦琪小生 【转】【WPF】WPF MVVM 简单实例
  4. 【java se】java注解
  5. 一些Java知识点
  6. SpringCloud升级之路2020.0.x版-24.测试Spring Cloud LoadBalancer
  7. 基于 Mysql 实现一个简易版搜索引擎
  8. 【Azure 应用服务】Python flask 应用部署在Aure App Service 遇见的 3 个问题
  9. QT学习日记篇01(1)-QT界面初探- *.pro文件详解
  10. 动环监控系统中B接口的实现