无源汇上下界可行流_格怎么找最小上界和最大下界

无源汇上下界可行流_格怎么找最小上界和最大下界给定一个包含 n 个点 m 条边的有向图,每条边都有一个流量下界和流量上界。给定源点 S 和汇点 T,求源点到汇点的最小流。输入格式第一行包含四个整数 n,m,S,T。接下来 m 行,每行包含四个整数 a,b,c,d 表示点 a 和 b 之间存在一条有向边,该边的流量下界为 c,流量上界为 d。点编号从 1 到 n。输出格式输出一个整数表示最小流。如果无解,则输出 No Solution。数据范围1≤n≤50003,1≤m≤125003,1≤a,b≤n,0≤c≤d≤21474836

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

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

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

给定源点 S 和汇点 T,求源点到汇点的最小流。

输入格式
第一行包含四个整数 n,m,S,T。

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

点编号从 1 到 n。

输出格式
输出一个整数表示最小流。

如果无解,则输出 No Solution。

数据范围
1≤n≤50003,
1≤m≤125003,
1≤a,b≤n,
0≤c≤d≤2147483647,
数据保证答案不超过 int 范围。

输入样例:
7 12 6 7
6 1 0 2147483647
1 7 0 2147483647
6 2 0 2147483647
2 7 0 2147483647
6 3 0 2147483647
3 7 0 2147483647
6 4 0 2147483647
4 7 0 2147483647
6 5 0 2147483647
5 7 0 2147483647
5 1 1 2147483647
3 4 1 2147483647
输出样例:
2
#include<bits/stdc++.h>
using namespace std;
const int N = 5e4 + 10;
const int M = 2*(125e3 + 10 + N);
const int INF = 2147483647;
typedef long long ll;
struct Edge{ 
   
    int v,next,w,l,h;
}edge[M];
int head[N],cnt;
void add(int u,int v,int l,int h){ 
   
    edge[cnt].v = v;
    edge[cnt].h = h;
    edge[cnt].l = l;
    edge[cnt].w = h - l;
    edge[cnt].next = head[u];
    head[u] = cnt ++;
}
int d[N],cur[N],q[N],tt = 0,hh = 0;
int A[N];
int n,m,s,e;
bool bfs(){ 
   
    hh = 0,tt = 0;
    memset(d,-1,sizeof d);
    d[s] = 0,cur[s] = head[s],q[tt ++] = 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 > 0){ 
   
                q[tt ++] = v;
                cur[v] = head[v];
                d[v] = d[t] + 1;
                if(v == e)return true;
            }
        }
    }
    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(w,limit - flow));
            if(!t)d[v] = -1;
            flow += t,edge[i].w -= t,edge[i ^ 1].w += t;
        }
    }
    return flow;
}
ll dinic(){ 
   
    int maxflow = 0,flow;
    while(bfs())while(flow = dfs(s,INF))maxflow += flow;
    return maxflow;
}
int main(){ 
   
    memset(head,-1,sizeof head);
    int ts,te,S,E;      //ts te代表源点和汇点 S E代表虚拟源点和虚拟汇点
    cin>>n>>m>>ts>>te;
    S = 0,E = n + 1;
    int x,y,l,h;
    for(int i = 0;i < m;i ++){ 
   
        scanf("%d %d %d %d",&x,&y,&l,&h);
        add(x,y,l,h),add(y,x,0,0);
        A[y] += l,A[x] -= l;
    }
    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);
    }
    add(te,ts,0,INF),add(ts,te,0,0);
    s = S,e = E;
    if(tot != dinic()){ 
   
        cout<<"No Solution"<<endl;
    }
    else{ 
   
        int res = edge[cnt - 1].w;
        edge[cnt - 1].w = edge[cnt - 2].w = 0;
        s = te,e = ts;
        cout<<res - dinic()<<endl;
    }
    return 0;
}
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

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

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


相关推荐

  • 方程自己解(1)——物理信息神经网络(PINN)

    方程自己解(1)——物理信息神经网络(PINN)文章目录前言(一)物理神经网络(PINN)解读1.1PINN基本背景1.2PINN算法描述前言  最近正在看利用“深度学习”(大概吧,其实只是利用了neuralnetwork的自动微分特性(AD)),在看一些文章的同时,将文章中提到的开源代码用起来和复现一些基本方程求解能够加快我学习的进度,这里将持续贴出一些方程和代码的求解过程。当然非常希望同学和朋友们可以给我指出错误,最后如果能够坚持下去,这个系列的工作会分享到github,同时如果对各个开源程序有所帮助那就更好了!  由于会参考到很多文章

    2022年10月23日
    0
  • Python文件名后缀_python获取目录下所有文件的文件名

    Python文件名后缀_python获取目录下所有文件的文件名method使用os.path.splitext(file)[1]可获得以.开头的文件后缀名。codeimportosfile="Hello.py"assertos.path.splitext(file)[1]==".py"assertos.path.splitext(file)[1][1:]=="py"

    2022年9月22日
    0
  • 11条javascript知识

    1.局部变量和全局变量var操作符定义的变量将成为定义该变量作用域中的局部变量。这个局部变量会在函数退出后销毁。不同于其他语言,javaScript不存在块级作用域。全局变量就是window对象的属性

    2021年12月20日
    49
  • c++int转换成char_字符数字转为int型

    c++int转换成char_字符数字转为int型1.通过ascii码:chara=’0′;intia=(int)a;/*notethattheintcastisnotnecessary–intia=awouldsuffice*/cout&lt;&lt;ia&lt;&lt;endl;结果如下:可以看出这种方法得到的其实是char对应的ascii码。因为ascii码…

    2022年10月2日
    0
  • Android开发-API指南- Calendar Provider

    Android开发-API指南- Calendar ProviderAndroid官方文档原创译文《CalendarProvider》。CalendarProvider是用来存放用户日历事件(event)的数据库。通过CalendarProvider的

    2022年7月1日
    24
  • kmp的最小循环节

    kmp的最小循环节

    2021年9月27日
    47

发表回复

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

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