POJ3426 ACM Computer Factory——最大流（EK+输出路径）

点这里

题意： 题目大意就是问把000加工到111（多少个零视输入的P而定）的最大零件个数，并且输出路径。
  但是具体的加工规则真的是把我看晕了 然后我们来说说具体的加工规则：N台机器，每台机器能同时加工Q台电脑，并且都有输入、输出两端，两端都有P个参数。输入端有0、1、2三种状态，输出端都有0、1两种状态。例如第一个样例中的第一行15 0 0 0 0 1 0表示这台机器能同时加工15台电脑，并且能将状态为000的电脑加工成010。这里我尝试理解了很久才看懂 规则就是每台机器只能加工其输入端对应状态的电脑，并加工成其输出端的状态。特殊点在于输入端的状态2表示0或1都能处理。
  求最多同时能加工多少电脑，并且输出所有的加工路径。

---

题解： 因为P和N的范围都不大，所以用EK算法就能解决。实际上是我一开始用Dinic没写出来。 至于输出路径 ，只要复制一下一开始的图，跑完模板之后对比前后的图，容量减少的边都要输出。

---

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
const int inf = 0x3f3f3f3f;
const int N = 60;

int p, n, w[N];
int in[N][12], out[N][12];
int pre[N], G[N][N], ori[N][N];	//记录前驱节点、存图
int bfs(int s, int t){
	int flow[N];		//到节点i的流量
	memset(pre, -1, sizeof pre);	//重置前驱节点
	flow[s] = inf;	pre[s] = 0;		//源点的容量无限
	queue<int> Q; Q.push(s);
	while(!Q.empty()){
		int u = Q.front(); Q.pop();
		if(u == t)	break;
		for(int i = 1; i <= t; i++)
			if(i != s && G[u][i] > 0 && pre[i] == -1){	//保证i不是起点，i没有去过，i还有残留容量
				pre[i] = u;
				Q.push(i);
				flow[i] = min(flow[u], G[u][i]);
			}
	}
	return pre[t] == -1 ? -1 : flow[t];
}
int maxflow(int s, int t){
	int Maxflow = 0;
	while(1){
		int flow = bfs(s, t);
		if(flow == -1)	break;
		int cur = t;
		while(cur != s){
			int fa = pre[cur];
			G[fa][cur] -= flow;
			G[cur][fa] += flow;	//更新残留容量
			cur = fa;
		}
		Maxflow += flow;	//这条增广路的流量汇入汇点
	}
	return Maxflow;
}
void print_flow(){
	memcpy(ori, G, sizeof G);
	int res = maxflow(0, n + 1);
	int cnt = 0, ans[110][3];
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= n; j++)
			if(ori[i][j] > G[i][j])
				ans[cnt][0] = i, ans[cnt][1] = j, ans[cnt++][2] = ori[i][j] - G[i][j];
	printf("%d %d\n", res, cnt);
	for(int i = 0; i < cnt; i++)	printf("%d %d %d\n", ans[i][0], ans[i][1], ans[i][2]);
}
int main(){
	while(~scanf("%d%d", &p, &n)){
		memset(G, 0, sizeof G);
		for(int i = 1; i <= n; i++){
			scanf("%d", w + i);
			for(int j = 1; j <= p; j++)	scanf("%d", in[i] + j);
			for(int j = 1; j <= p; j++)	scanf("%d", out[i] + j);
		}
		for(int i = 1; i <= n; i++){
			int f1 = 1, f2 = 1;
			for(int j = 1; j <= p; j++){
				if(in[i][j] == 1)	f1 = 0;
				if(out[i][j] == 0)	f2 = 0;}
			if(f1)	G[0][i] = w[i];
			if(f2) 	G[i][n + 1] = w[i];
			for(int j = 1; j <= n; j++){
				if(i == j)	continue;
				int flag = 1;
				for(int k = 1; k <= p; k++)	
					if(out[i][k] + in[j][k] == 1){	flag = 0;	break;}
				if(flag)	G[i][j] = min(w[i], w[j]);
			}
		}
		print_flow();
	}
	return 0;
}