题意:n个点,m条边,每条边有容量限制 l--c,每个点满足容量平衡(流入等于流出),求可行解

无源无汇可行流问题,建立以一个超级源点和超级汇点,由于原来最大流问题时候,流量下界其实为0,

所以要转化,把边(设u-->v)的容量改为c-l,但是这样不平衡了,所以S流入v点l,u点流出到T要l,这样

保证了u,v流量平衡,用数组sumin[i]记录下i点流入下限之和,最后超级源点流入i。

最后求一次s-->t的最大流(走一遍dinic),如果添加的边都满流,说明有解(此时每条边所用流量+下限即可),

反之无解(必需要满流,否则不遵循流量平衡条件!)。(无源无汇模型和参考黑书

p366)。

#include<iostream>   //15ms
#include<cstdio>
#include<vector>
#include<queue>
using namespace std;
int n,m;const int inf=0x3f3f3f3f;
int e[90000][5];int head[210]; //链前星存边,0:to,1:pre,2,残量;3:l(下界);4,c
int sum_in[210];int sum_out[210]; //点i流入之和,流出之和
int vis[210];int level[210];
bool bfs() //dinic,小心细节!要熟练
{
for(int i=0;i<=n+1;i++)
vis[i]=level[i]=0;
queue<int>q;q.push(0);vis[0]=1;
while(!q.empty())
{
int cur=q.front();q.pop();
for(int i=head[cur];i!=-1;i=e[i][1])
{
int v=e[i][0];
if(!vis[v]&&e[i][2]>0)
{
level[v]=level[cur]+1;
if(v==n+1)return 1;
vis[v]=1;
q.push(v);
}
}
}
return vis[n+1];
}
int dfs(int u,int minf)
{
if(u==n+1||minf==0){return minf;}
int sumf=0,f;
for(int i=head[u];i!=-1&&minf;i=e[i][1])
{ int v=e[i][0];
if(level[v]==level[u]+1&&e[i][2]>0)
{
f=dfs(v,minf<e[i][2]?minf:e[i][2]);
if(f<=0)continue;
e[i][2]-=f;e[i^1][2]+=f;
sumf+=f;minf-=f;
}
}
return sumf;
}
void dinic()
{ int sumflow=0;
while(bfs())
{
sumflow+=dfs(0,inf);
}
}
bool check() //判断有无解
{
for(int i=head[0];i!=-1;i=e[i][1]) //所有从超级源点出来的流量必满,否则无解!
if(e[i][2]!=0)return 0; //满必然有解,无需再判断汇点是否满(重复了)
// int v=n; //起初多此一举判断汇点满流情况,但是要注意一点
// for(int i=head[n+1];i!=-1;i=e[i][1]) //边遍历顺序,前向星是前一条边,按添加时顺序相反
// if(e[i][2]!=sum_out[v--])return 0;//添加迟,出现早。
return 1;
}
int main()
{
while(~scanf("%d%d",&n,&m))
{
for(int i=0;i<=n+1;i++)
{
head[i]=-1;
sum_in[i]=sum_out[i]=0;
}
int a,b,l,c; int nume=0;
for( ;nume<2*m;) //读入,用每条边e[i][2]流量是残量,其他无用,只是保存起来,输出时用一下
{
scanf("%d%d%d%d",&a,&b,&l,&c);
e[nume][0]=b;e[nume][1]=head[a];head[a]=nume;
e[nume][4]=c;e[nume][3]=l;e[nume++][2]=c-l;
sum_in[b]+=l;sum_out[a]+=l;
e[nume][0]=a;e[nume][1]=head[b];head[b]=nume;
e[nume++][2]=0;
}
for(int i=1;i<=n;i++)
{
e[nume][0]=i;e[nume][1]=head[0];head[0]=nume;
e[nume++][2]=sum_in[i];
e[nume][0]=0;e[nume][1]=head[i];head[i]=nume;
e[nume++][2]=0;
e[nume][0]=n+1;e[nume][1]=head[i];head[i]=nume;
e[nume++][2]=sum_out[i];
e[nume][0]=i;e[nume][1]=head[n+1];head[n+1]=nume;
e[nume++][2]=0;
}
dinic();
if(!check())printf("NO\n");
else
{
printf("YES\n");
for(int i=0;i<2*m;i+=2)
{
printf("%d\n",e[i][4]-e[i][2]);
}
}
}
return 0;
}

SGU 194 无源无汇可行流求解的更多相关文章

  1. sgu 194 Reactor Cooling(有容量上下界的无源无汇可行流)

    [题目链接] http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=20757 [题意] 求有容量上下界的无源无汇可行流. [思路] ...

  2. sgu 194 上下界网络流可行流判定+输出可行流

    #include <cstdio> #include <cstring> #define min(a,b) ((a)<(b)?(a):(b)) #define oo 0x ...

  3. 【HDU 4940】Destroy Transportation system(无源无汇带上下界可行流)

    Description Tom is a commander, his task is destroying his enemy’s transportation system. Let’s repr ...

  4. [ACdream 1211 Reactor Cooling]无源无汇有上下界的可行流

    题意:无源无汇有上下界的可行流 模型 思路:首先将所有边的容量设为上界减去下界,然后对一个点i,设i的所有入边的下界和为to[i],所有出边的下界和为from[i],令它们的差为dif[i]=to[i ...

  5. ZOJ 1314 Reactor Cooling | 上下界无源汇可行流

    ZOJ 1314 Reactor Cooling | 上下界无源汇可行流 题意 有一个网络,每条边有流量的上界和下界,求一种方案,让里面的流可以循环往复地流动起来. 题解 上下界无源汇可行流的模型: ...

  6. ZOJ 2314 Reactor Cooling | 无源汇可行流

    题目: 无源汇可行流例题 http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1314 题解: 证明什么的就算了,下面给出一种建图方式 ...

  7. 算法复习——无源汇可行流(zoj2314)

    题目: The terrorist group leaded by a well known international terrorist Ben Bladen is buliding a nucl ...

  8. 【BZOJ-2055】80人环游世界 上下界费用流 (无源无汇最小费用最大流)

    2055: 80人环游世界 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 321  Solved: 201[Submit][Status][Discus ...

  9. SGU 194 Reactor Cooling Dinic求解 无源无汇有上下界的可行流

    题目链接 题意:有向图中有n(1 <= n <= 200)个点,无自环或者环的节点个数至少为3.给定每条边的最小流量和最大流量,问每条边的可行流量为多少? 思路:一般求解的网络流并不考虑下 ...

随机推荐

  1. 【学习笔记】C++文件操作详解(ifstream、ofstream、fstream)

    C++ 通过以下几个类支持文件的输入输出: ofstream: 写操作(输出)的文件类 (由ostream引申而来) ifstream: 读操作(输入)的文件类(由istream引申而来) fstre ...

  2. objectbox基础

    objectbox基础 参考链接 官网地址 http://objectbox.io github地址 https://github.com/objectbox/objectbox-java https ...

  3. B树、B+树、红黑树、AVL树

    定义及概念 B树 二叉树的深度较大,在查找时会造成I/O读写频繁,查询效率低下,所以引入了多叉树的结构,也就是B树.阶为M的B树具有以下性质: 1.根节点在不为叶子节点的情况下儿子数为 2 ~ M2. ...

  4. 第8章 应用协议 图解TCP/IP 详解

    第8章 应用协议 图解TCP/IP 详解 8.1 应用层协议概要 应用层协议的定义 TCP和IP等下层协议是不依赖上层应用类型.实用性非常广的协议.而应用协议则是为了实现某种应用而设计和创造的协议. ...

  5. 秦晓波著的编写高质量代码-改善Java程序的151个建议一书中的线程解释错误.

    位置: 建议127: Lock与synchronized是不一样的 首先在概念上纠正这一篇内容: 援引Java源码中关于ReentrantLock的开篇说明: * A reentrant mutual ...

  6. lavarel功能总结

    详细可参见笔记:laraval学习笔记(二) 路由 route 绑定模型,绑定参数 模版 blade .blade.php后缀,有laravel自己的模版语法 模型 model 如果用create创建 ...

  7. 解决python pip安装提示"not a supported wheel on this platform"

    python下载.whl安装文件后使用pip安装有时候会提示报错"xxxxx not a supported wheel on this platform",应该是下载的库文件版本 ...

  8. 7-Java-C(小题答案)

    1:58497 2:171700 3:145 4:i + j+2 == k+1 || i + k+2 == j+1 || k + j+2 == i+1 5:s + " " + (c ...

  9. wkhtmltopdf导出html到pdf

    1.使用背景     最近公司需要把html页面的内容生成pdf并下载,试过很多方法都没有满意的效果,后来找到wkhtmltopdf这款软件,终于解决了这个问题. wkhtmltopdf是exe文件, ...

  10. 谷歌 滚动条样式 scrollbar

    谷歌  滚动条样式 scrollbar main.less :: 是css3 最新的伪类 /*滚动条整体样式*/ ::-webkit-scrollbar { width: 10px; height: ...