http://poj.org/problem?id=3436

題意好難懂，有n個工廠加工電腦半成品，半成品電腦需要p個部件。

輸入p(需要的部件數) n（工廠數）

Q_i? (工廠最大生產量) ?S_i,1?Si_,2...S_{i,P ?（s表示在工廠i聲場必須滿足的條件1必須已經有了，0可以沒有，2可以有，也可以沒有）}?D_i,1?D_i,2...D_i,P,(經過i廠生產之后，所滿足的條件1有，0木有)

題目建圖不好了理解，建完圖后就是套最大流模板的問題了。。。

首先構建網絡圖,一臺電腦在進入生產線之前各部分的零件狀態肯定全為0,離開生產線時,各部分的零件狀態肯定全為1. 所以源點與輸入狀態全為0的機器相連,輸出狀態全為1機器與匯點相連,并且設定容量為INF,然后對每臺機器的輸入,輸出 拆為兩個點,容量為Q;最后判斷兩臺機器之間,某臺機器的輸出狀態是否可以為另一臺機器的輸入狀態,是,則相連,容量為 INF,至此網絡流圖建立完畢,下面就是求最大流問題了. 但是為了輸出相連的機器信息,我們在構建網絡流圖完畢后,應該備份一下殘留矩陣residual為cresidual,當EK()完以后, 通過比較cresidual與residual的值才可以確定那兩個機器之間被利用了. View Code #include <cstdio>
#include <cstring>
#include <queue>
#include <iostream>
#define maxn 200
#define inf 0x7fffffff
using namespace std;

int g[maxn][maxn],flow[maxn],tmp[maxn][maxn];
int pre[maxn];
struct node
{
int in[20],out[20];
int num;
}tp[maxn];
int s,e,p,n,ans;
void input()
{
memset(g,0,sizeof(g));
for (int i = 1; i <= n; ++i)
{
cin>>tp[i].num;
g[i][i+n] = tp[i].num;
for (int j = 0; j < p; ++j)
cin>>tp[i].in[j];
for (int j = 0; j < p; ++j)
cin>>tp[i].out[j];
}
}
bool isok(node a,node b)
{
for (int i = 0; i < p; ++i)
{
if (a.out[i] + b.in[i] == 1)
return false;
}
return true;
}
int getf(int *l)
{
int t0 = 0;
int t1 = 0;
for (int i = 0; i < p; ++i)
{
if (l[i] == 1) t1++;
else t0++;
}
if (t1 == p) return 1;
if (t0 == p) return 0;
return -1;
}
void build()
{
s = 0;
e = 2*n + 1;
int i,j;
for (i = 1; i <= n; ++i)
{
for (j = 1; j <= n; ++j)
{
if (i != j && isok(tp[i],tp[j]))
{
g[i + n][j] = inf;
}
}
}
for (i = 1; i <= n; ++i)
{
int flag1 = getf(tp[i].in);
int flag2 = getf(tp[i].out);
if (flag1 == 0) g[s][i] = inf;
if (flag2 == 1) g[i + n][e] = inf;
}
memcpy(tmp,g,sizeof(g));
}
int bfs()
{
int i;
memset(pre,-1,sizeof(pre));
for (i = s; i <= e; ++i)
flow[i] = inf;
queue<int>q;
q.push(s);
while (!q.empty())
{
int pos = q.front(); q.pop();
if (pos == e) return flow[e];
for (i = s; i <= e; ++i)
{
if (i == s || g[pos][i] == 0 || pre[i] != -1) continue;
pre[i] = pos;
flow[i] = min(flow[pos],g[pos][i]);
q.push(i);
}
}
return -1;
}
void EK()
{
ans = 0;
while (1)
{
int d = bfs();
if (d == -1) break;
ans += d;
for (int i = e; i != s; i = pre[i])
{
g[pre[i]][i] -= d;
g[i][pre[i]] += d;
}
}
}
void print()
{
int ct = 0;
int i,j;
for (i = n + 1; i < e; ++i)
{
for (j = 1; j <= n; ++j)
{
if (i != j && g[i][j] < tmp[i][j])
ct++;
}
}
printf("%d %d\n",ans,ct);
for (i = n + 1; i < e; ++i)
{
for (j = 1; j <= n; ++j)
{
if (i != j && g[i][j] < tmp[i][j])
printf("%d %d %d\n",i - n,j,tmp[i][j] - g[i][j]);
}
}
}
int main()
{
//freopen("d.txt","r",stdin);
while (~scanf("%d%d",&p,&n))
{
input();
build();
EK();
print();
}
return 0;
}

轉載于:https://www.cnblogs.com/E-star/archive/2012/03/20/2408422.html

總結

以上是生活随笔為你收集整理的pku 3436 ACM Computer Factory ——最大流 EK的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。