网络流，流水线模拟

 ACM Computer Factory

题意：

一共有N个机器，每个机器有P个元素，对应输入的时候输入N个机器的信息，第一个数表示这个机器可以一共能够生产多少产物，接下来2p个元素，前p个元素：其中有三种数值，1,2,0，分别表示必须有这个位子的组件，可有可没有这个位子的组件，以及不能有这个位子的组件。一个中间产物想在这个机器上进行生产接下来的产物，必须要满足这些条件才行，对应在这个机器上产生的产物的各个组件的结果是后p个元素。

solution:

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define FRE(name) freopen(#name".txt","r",stdin);
using namespace std;
const int N=1.5e4+5,M=55;
const int INF=0x3f3f3f3f;
struct edge{int u,v,cap,next;}e[N*100],last[N*100];int tot=0;
int n,p,S,T,cnt,head[N],dep[N],val[N],res[N*100][3],in[M][N],out[M][N],q[N*50];
void add(int x,int y,int z){
    e[tot].u=x;e[tot].v=y;e[tot].cap=z;e[tot].next=head[x];head[x]=tot++;
    e[tot].u=y;e[tot].v=x;e[tot].cap=0;e[tot].next=head[y];head[y]=tot++;
}
bool isNULL(int *in){//与源点相连
    for(int i=1;i<=p;i++) if(in[i]==1) return 0;
    return 1;
}
bool isFULL(int *out){//与源点相连
    for(int i=1;i<=p;i++) if(!out[i]) return 0;
    return 1;
}
bool canLink(int *out,int *in){//普通节点之间的相连
    for(int i=1;i<=p;i++) if(in[i]+out[i]==1) return 0;
    return 1;
}
void mapping(){
    S=0,T=n+1;
    for(int i=1;i<=n;i++){
        if(isNULL(in[i])) add(S,i,val[i]);
        if(isFULL(out[i])) add(i,T,val[i]);
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            if(i==j) continue;
            if(canLink(out[i],in[j])) add(i,j,val[i]);
        }
    }
}
bool bfs(){
    for(int i=S;i<=T;i++) dep[i]=-1;
    int h=0,t=1;dep[S]=0;q[t]=S;
    while(h!=t){
        int x=q[++h];
        for(int i=head[x];~i;i=e[i].next){
            if(dep[e[i].v]==-1&&e[i].cap){
                dep[e[i].v]=dep[x]+1;
                if(e[i].v==T) return 1;
                q[++t]=e[i].v;
            }
        }
    }
    return 0;
}
int dfs(int x,int f){
    if(x==T) return f;
    int used=0,t;
    for(int i=head[x];~i;i=e[i].next){
        if(e[i].cap&&dep[e[i].v]==dep[x]+1){
            t=dfs(e[i].v,min(e[i].cap,f));
            e[i].cap-=t;e[i^1].cap+=t;
            used+=t;f-=t;
            if(!f) return used;
        }
    }
    if(!used) dep[x]=-1;
    return used;
}
void dinic(){
    int ans=0;
    while(bfs()) ans+=dfs(S,INF);
    printf("%d ",ans);
}
void init(){
    tot=0;cnt=0;
    memset(head,-1,sizeof head);
}
int main(){
    //FRE(sh);
    while(scanf("%d%d",&p,&n)==2){
        init();
        for(int i=1;i<=n;i++){
            scanf("%d",&val[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]);
        }
        mapping();
        memcpy(last,e,sizeof e);
        dinic();
        for(int i=0;i<tot;i+=2){
            if(e[i].u==S||e[i].v==T) continue;
            if(last[i].cap!=e[i].cap){
                res[++cnt][0]=e[i].u;
                res[cnt][1]=e[i].v;
                res[cnt][2]=e[i^1].cap;
            }
        }
        printf("%d\n",cnt);
        for(int i=1;i<=cnt;i++){
            printf("%d %d %d\n",res[i][0],res[i][1],res[i][2]);
        }
    }
    return 0;
}