zoj 2314 Reactor Cooling

 

 

好郁闷。最大流。但是没见过这种类型的，想了种方法，但是答案不对 = =。。。

 

搜了下，这个是另一种类型，无源汇（其实我觉得汇源比较顺口。。）的有上下界的最大流。

 

引用http://blog.imzzl.com/2010/07/630.html的一段话

 

 

---

总的来说就是新建一个原点汇点 s,t，对原来的每个点i，设m(i)=sum{b<j,i>}-sum{b<i,j>}其中b<>为流量下界，若m(i)>0，连一条<s,i>容量为m(i)的边；若m(i)<0，连一条<i,t>容量为|m(i)|的边。然后将原来边的容量变为c<i,j>-b<i,j>。求一次最大流。如果与s和t关联的边全部满流，则可行流存在，且每条边的流量为现在的流量+流量的下界。否则不存在可行流。

 

---

 

 

下了个《一种简易的方法求解流量有上下界的网络中网络流问题》明天印出来看看。。

 

#include <cstdio> #include <cstdlib> #include <iostream> #include <string.h> #include <queue> #include <limits.h> #define MAX 210 using namespace std; int n,m; int cap[MAX][MAX]; int flow[MAX][MAX]; int EKarp(int s,int t) { int a[MAX],u,v,f,pre[MAX]; f = 0; queue<int> Q; memset(flow,0,sizeof(flow)); while(1) { Q.push(s); memset(a,0,sizeof(a)); a[s] = INT_MAX; while( !Q.empty() ) { u = Q.front(); Q.pop(); for(v=1; v<=n+1; v++) if( !a[v] && cap[u][v] > flow[u][v] ) { Q.push(v); a[v] = a[u] < cap[u][v] - flow[u][v] ? a[u] : cap[u][v] - flow[u][v]; pre[v] = u; } } if( a[t] == 0 ) break; for(u=t; u!=s; u=pre[u]) { flow[pre[u]][u] += a[t]; flow[u][pre[u]] -= a[t]; } f += a[t]; } return f; } int f[MAX*MAX],t[MAX*MAX]; int l[MAX*MAX]; int main() { int i,k,c,ncases; scanf("%d",&ncases); while( ncases-- ) { scanf("%d%d",&n,&m); memset(cap,0,sizeof(cap)); for(i=1; i<=m; i++) { scanf("%d%d%d%d",&f[i],&t[i],&l[i],&c); cap[f[i]][t[i]] = c - l[i]; cap[0][t[i]] += l[i]; cap[f[i]][n+1] += l[i]; } EKarp(0,n+1); int flag = 1; for(i=1; i<=m; i++) if( cap[0][t[i]] - flow[0][t[i]] != 0 || cap[f[i]][n+1] - flow[f[i]][n+1] != 0 ) { flag = 0; break; } if( flag ) { printf("YES/n"); for(i=1; i<=m; i++) printf("%d/n",flow[f[i]][t[i]]+l[i]); } else printf("NO/n"); } return 0; }