POJ3436

题意

第一行输入数字N为机器个数，P为零件的种类数。

接下来的N行中，第一个数字为此工厂生产产品个数。

前P个数字中是需要的种类，0表示不需要，1表示必须要，2表示可要可不要。

后P个数字是产品的零件种类，0表示有，1表示没有。

建图：

建立一个超级源一个超级汇。

把一个机器分成两个节点i和i+N，流量限制为该机器的总产量。

如果p的产品可以作为q的生产原料则把p+n到q连一条流量限制无限大的边。

匹配的条件是p的产品不是2的地方和q的生产原料1和0的情况相同。

超级源和原料建立边的条件是原料中没有1存在。

产品和超级汇建立边的条件是产品中所有的数字都是1。

这道题的目的就是求这个图的最大流量。用ISAP算法实现（同一个系列中有文章介绍）

输出的生产线的条件是产品到加工品的边的流量不为0。

#include<iostream>    using namespace std;   struct Node{int cap, flow;};int p, n, size;int dp1[50][100], dp2[50][100];Node edges[1000][1000];int visit[200], s, t;int p1[1000000], num[1000000];int q[1000000], d[1000000];int Count;int BFS(){int i, front, tail,temp;for(i = 1; i <= size; i++){visit[i] = 0;}front = 0; tail = 1;q[front] = t;visit[t] = 1; d[t] = 0;num[d[t]]++;while(front < tail){temp = q[front];front++;for(i = 1; i <= size; i++){if(!visit[i] && edges[i][temp].cap > edges[i][temp].flow){q[tail] = i;d[i] = d[temp] + 1;num[d[i]]++;visit[i] = 1;tail++;}}}return visit[t];}int augment(){int a = 1000000, temp;for(int i = 0; i < Count - 1; i++){temp = edges[p1[i]][p1[i + 1]].cap - edges[p1[i]][p1[i + 1]].flow;if(temp < a)a = temp;}for(int i = 0; i < Count - 1; i++){edges[p1[i]][p1[i + 1]].flow += a;edges[p1[i + 1]][p1[i]].flow -= a;}//for(int i = 0; i < Count ; i++)//{//cout << "aaa   " << p[i] << " " << a <<endl;//}//cout <<endl;return a;}int maxFlow(){BFS();int x = s,i, flow = 0, ok, min, m;Count = 0; p1[Count++] = s;while(d[s] < size){if(x == t){flow += augment();x = s;  Count = 0;p1[Count++] =s;}ok = 0;for(i = 1; i <= size ; i++){if(d[x] == d[i] + 1 && edges[x][i].cap > edges[x][i].flow){ok = 1;x = i; p1[Count++] = i;}}if(ok == 0){min = size - 1;for(i = 1; i <= size ; i++){if(edges[x][i].cap > edges[x][i].flow){if(d[i] < min){min = d[i];m = i;}}}if(--num[d[x]] == 0)break;    d[x] = min + 1; num[d[x]]++;//p[Count++] = m;//x = m;if(x != s){x = p1[Count - 2];Count --;}}}return flow;}void build(){int i, j, k;bool flag;s = 2 * n + 1; t = 2 * n + 2;for(i = 1; i <= n; i++){flag = true;for(j = 1; j <= p; j++){if(dp1[i][j] == 1){flag = false;break;}}if(flag){edges[2 * n + 1][i].cap = 100000000;}}for(i = 1; i <= n; i++){flag = true;for(j = 1; j <= p; j++){if(dp2[i][j] == 0){flag = false;break;}}if(flag){edges[i + n][2 * n + 2].cap = 100000000;}}for(i = 1; i <= n; i++){for(j = 1; j <= n; j++){flag = true;if(i != j){    for(k = 1; k <= p; k++)    {    if(dp1[i][k] == dp2[j][k] || dp1[i][k] == 2)    continue;    else    {    flag = false;break;    }     }     if(flag)    {    edges[j + n][i].cap = 100000000;    }}}}}int main()    {   int i, j, temp, num;while(scanf("%d", &p) != EOF){    scanf("%d", &n);    size = 2 * n + 2;   for(i = 1; i <= size; i++)   {   for(j = 1; j <= size; j++)   {   edges[i][j].cap = 0; edges[i][j].flow = 0;   }   }   for(i = 1; i <= n; i++)   {   scanf("%d", &temp);   edges[i][i + n].cap = temp;    for(j = 1; j <= p; j++)   scanf("%d", &dp1[i][j]);    for(j = 1; j <= p; j++)   scanf("%d", &dp2[i][j]);   }   build();   cout << maxFlow();num = 0;for(i = n + 1; i <= 2 * n; i++){for(j = 1; j <= n; j++){if(edges[i][j].flow > 0){num ++;}}}cout << " " << num << endl;for(i = n + 1; i <= 2 * n; i++){for(j = 1; j <= n; j++){if(edges[i][j].flow > 0)cout << (i - n) << " " << j << " " << edges[i][j].flow << endl;}}}}