一个最大流问题.
首先给出有几个点,每个点的容量是多少.
然后给出有几条边,给出每条边并给出边的最大流量.
然后给出有几个源点,几个汇点.源点有哪些,汇点有哪些.
比起标准的最大流,这题多了两个地方,一个是点也具有容量,另一个是源点和汇点都不只一个.
但是我们可以通过建图时的预处理,把它变成标准模板题.
首先多个源点的问题,我们可以另外加一个点(超级源点),把它和所有的源点相连,并且边的最大流就是那个源点的容量.这样就只剩一个源点,并且不影响.多个汇点用的是同样的办法.
还有就是点的容量问题,我们可以把点的容量转化为边的容量,办法就是拆点,把一个点变成两个.
例如,总共有7个点.
现在我们的图中有一段是 2-1-3 ,2流到1在流到3.这时1这个点的容量是10
那么我们可以把它变成2-1-8-3.我们把点1,拆成1 和 1 + n(总共n个点)并且把1和1+n的边的最大流量设成10.
这样我们把一个n个点的图.变成2n个点,但是成功把点的容量转化为边的容量.这样就可以直接套用最大流模板算法Edmonds-Karp来求解.
AC代码:
#include<stdio.h>
#include<string.h>
#include<queue>
using namespace std;
const int N = 300;
const int INF = 0x3f3f3f3f;
int cap[N][N];
int flow[N][N];
int a[N];
int p[N];
int f;
int n,m;
int ek() {
queue<int> q;
memset(flow , 0 ,sizeof(flow));
f = 0;
while(1) {
memset(a, 0 ,sizeof(a));
a[0] = INF;
q.push(0);
while(!q.empty()) {
int u = q.front() ;
q.pop();
for (int v = 0 ; v <= 2 * n + 1; v++) {
if (!a[v] && cap[u][v] - flow[u][v] > 0) {
p[v] = u ;
q.push(v);
a[v] = a[u] < cap[u][v] - flow[u][v] ? a[u] : cap[u][v] - flow[u][v];
}
}
}
if (a[2 * n + 1] == 0)
break;
for (int u = 2 * n + 1 ; u != 0 ;u = p[u]) {
flow[p[u]][u] += a[2 * n + 1];
flow[u][p[u]] -= a[2 * n + 1];
}
f += a[2 * n + 1];
}
return f;
}
int main () {
while(scanf("%d",&n) != EOF) {
memset(cap , 0 ,sizeof(cap));
int a,b,temp,sn,en;
for (int i = 1 ; i <= n ; i++) {
scanf("%d",&temp);
cap[i][i + n] = temp;
}
scanf("%d",&m);
for (int i = 0 ; i < m ;i++) {
scanf("%d%d%d",&a,&b,&temp);
cap[a + n][b] = temp;
}
scanf("%d%d",&sn,&en);
for (int i = 0 ; i < sn ;i++) {
scanf("%d",&temp);
cap[0][temp] = cap[temp][temp + n];
}
for (int i = 0 ; i < en ;i++) {
scanf("%d",&temp);
cap[temp + n][2 * n + 1] = cap[temp][temp + n];
}
printf("%d\n",ek()) ;
}
}