2188. 无源汇上下界可行流(无源汇上下界最大流)

2188. 无源汇上下界可行流(无源汇上下界最大流)给定一个包含 n 个点 m 条边的有向图,每条边都有一个流量下界和流量上界。求一种可行方案使得在所有点满足流量平衡条件的前提下,所有边满足流量限制。输入格式第一行包含两个整数 n 和 m。接下来 m 行,每行包含四个整数 a,b,c,d 表示点 a 和 b 之间存在一条有向边,该边的流量下界为 c,流量上界为 d。点编号从 1 到 n。输出格式如果存在可行方案,则第一行输出 YES,接下来 m 行,每行输出一个整数,其中第 i 行的整数表示输入的第 i 条边的流量。如果不存在可行方案,直接输

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

给定一个包含 n 个点 m 条边的有向图,每条边都有一个流量下界和流量上界。

求一种可行方案使得在所有点满足流量平衡条件的前提下,所有边满足流量限制。

输入格式
第一行包含两个整数 n 和 m。

接下来 m 行,每行包含四个整数 a,b,c,d 表示点 a 和 b 之间存在一条有向边,该边的流量下界为 c,流量上界为 d。

点编号从 1 到 n。

输出格式
如果存在可行方案,则第一行输出 YES,接下来 m 行,每行输出一个整数,其中第 i 行的整数表示输入的第 i 条边的流量。

如果不存在可行方案,直接输出一行 NO。

如果可行方案不唯一,则输出任意一种方案即可。

数据范围
1≤n≤200,
1≤m≤10200,
1≤a,b≤n,
0≤c≤d≤10000

输入样例1:
4 6
1 2 1 3
2 3 1 3
3 4 1 3
4 1 1 3
1 3 1 3
4 2 1 3
输出样例1:
YES
1
2
3
2
1
1
输入样例2:
4 6
1 2 1 2
2 3 1 2
3 4 1 2
4 1 1 2
1 3 1 2
4 2 1 2
输出样例2:
NO
#include<bits/stdc++.h>
using namespace std;
const int N = 210;
const int M = (N + 10200) * 2;
const int INF = 0x3f3f3f3f;
struct Edge{ 
   
    int v,next,low,high,w;
}edge[M];
int A[N];
int head[N],cnt;
void add(int u,int v,int l,int h){ 
   
    //正向边
    edge[cnt].v = v;
    edge[cnt].low = l;
    edge[cnt].high = h;
    edge[cnt].w = h - l;
    edge[cnt].next = head[u];
    head[u] = cnt ++;
}
int n,m,s,e;
int d[N],cur[N],q[N],tt = 0,hh = 0;
bool bfs(){ 
   
    memset(d,-1,sizeof d);
    hh = 0,tt = 0;
    d[s] = 0,q[tt ++] = s,cur[s] = head[s];
    while(hh < tt){ 
   
        int t = q[hh ++];
        for(int i = head[t];~i;i = edge[i].next){ 
   
            int v = edge[i].v,w = edge[i].w;
            if(d[v] == -1 && w){ 
   
                d[v] = d[t] + 1;
                cur[v] = head[v];
                if(v == e)return true;
                q[tt ++] = v;
            }
        }
    }
    return false;
}
int dfs(int u,int limit){ 
   
    if(u == e)return limit;
    int flow = 0;
    for(int i = cur[u];~i && flow < limit;i = edge[i].next){ 
   
        int v = edge[i].v,w = edge[i].w;
        cur[u] = i;
        if(d[v] == d[u] + 1 && w > 0){ 
   
            int t = dfs(v,min(limit - flow,w));
            if(!t)d[v] = -1;
            flow += t,edge[i].w -= t,edge[i ^ 1].w += t;
        }
    }
    return flow;
}
int dinic(){ 
   
    int maxflow = 0,flow;
    while(bfs())while(flow = dfs(s,INF))maxflow += flow;
    return maxflow;
}
int main(){ 
   
    memset(head,-1,sizeof head);
    cin>>n>>m;
    int x,y,l,h;
    s = 0,e = n + 1;
    for(int i = 1;i <= m;i ++){ 
   
        cin>>x>>y>>l>>h;
        add(x,y,l,h);
        add(y,x,0,0);               
        A[y] += l,A[x] -= l;        //如果一个点流入流量需要减去的流量大于流出流量需要减去的值,那么需要从S连一个边到这个点
    }                               //如果一个点流入流量需要减去的流量小于流出流量需要减去的值,那么需要从这个点连一个边到E
    int tot = 0;
    for(int i = 1;i <= n;i ++){ 
   
        if(A[i] > 0)add(s,i,0,A[i]),add(i,s,0,0),tot += A[i];       //
        else if(A[i] < 0)add(i,e,0,-A[i]),add(e,i,0,0);
    }
    if(tot != dinic())cout<<"NO"<<endl;     //跑一遍Dinic之后S的出流全满,E的流入流也全满。因为需要满足流量守恒
    else{ 
   
        cout<<"YES"<<endl;
        for(int i = 0;i < 2 * m;i += 2){ 
   
            cout<<edge[i ^ 1].w + edge[i].low<<endl;    //当前边的流量等于残差网络反向边的流量
        }
    }
    return 0;
}
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/168591.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 消息队列 rabbitmq面试题(中间件面试题)

    文章目录为什么使用MQ?MQ的优点消息队列有什么优缺点?RabbitMQ有什么优缺点?你们公司生产环境用的是什么消息中间件?Kafka、ActiveMQ、RabbitMQ、RocketMQ有什么优缺点?MQ有哪些常见问题?如何解决这些问题?什么是RabbitMQ?rabbitmq的使用场景RabbitMQ基本概念RabbitMQ的工作模式如何保证RabbitMQ消息的顺序性?消息如何分发?消…

    2022年4月14日
    61
  • 深度残差网络(ResNet)之ResNet34的实现和个人浅见[通俗易懂]

    深度残差网络(ResNet)之ResNet34的实现和个人浅见[通俗易懂]残差网络是由来自MicrosoftResearch的4位学者提出的卷积神经网络,在2015年的ImageNet大规模视觉识别竞赛(ImageNetLargeScaleVisualRecognitionChallenge,ILSVRC)中获得了图像分类和物体识别的优胜。**残差网络的特点是容易优化,并且能够通过增加相当的深度来提高准确率。其内部的残差块使用了跳跃连接(shortcut),缓解了在深度神经网络中增加深度带来的梯度消失问题**。残差网络(ResNet)的网络结构图举例如下:

    2022年10月6日
    3
  • matlab 实现 garch 模型波动率估计

    matlab 实现 garch 模型波动率估计matlab 实现 garch 模型波动率估计 matlab 实现 garch 模型波动率估计代码高亮问题数据获取数据处理描述性统计时间序列平稳性检验相关和偏自相关 arch 效应检验建立 garch 模型波动率估计代码及文档地址代码高亮问题这边首先说一个问题 你们看到的 matlab 代码没有高亮 看上去很不舒服 但这个锅是 csdn 的 真垃圾 连个语法高亮

    2025年6月10日
    4
  • centos7安装pycharm_pycharm配置环境变量

    centos7安装pycharm_pycharm配置环境变量Centos下pycharm的安装与配置1下载安装pycharm首先在下面的网址下载安装包:https://www.jetbrains.com/pycharm/download/#section=linux然后使用下列指令将安装包放入合适的目录下(本文将安装包放入了/usr/local目录下):cd/home/yue/Downloadsmvpycharm-community-2020.1.tar.gz/usr/local然后进行安装包所在目录,进行解压:tar-zxvfpyc

    2022年8月29日
    3
  • 源码网_python源码从哪下载

    源码网_python源码从哪下载源码目录结构我们首先来看下models.py的代码结构我们可以看到这个模块中定义了12个属性和22个模型类,我们依次来看属性源码分析importosfromenumimportEnu

    2022年7月30日
    7

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号