ACM Computer Factory

ACM Computer Factory

题目链接：http://poj.org/problem?id=3436

网络流Dinic

将一个机器拆分成两个点，这两个点之间的容量为机器的加工量；建立一个超级源点和超级汇点，并将满足条件的点连接在一起，容量均为inf。之后跑Dinic算法即可。

注意输出的时候不要求输出为路径顺序（这里坑了一晚上才知道不需要=。=）

代码如下：

#include <cstdio>
#include <vector>
#include <cstring>
#include <queue>
#define pb(x) push_back(x)
#define N 200
using namespace std;
int const inf=1<<30;
struct edge{
    int to,rev,cap;
    edge(int _to=0,int _cap=0,int _rev=0){
        to=_to;rev=_rev;cap=_cap;
    }
};
struct node{
    int from,to,c;
    node(int _from=0,int _to=0,int _c=0){
        from=_from;to=_to;c=_c;
    }
};
struct dinic{
    vector<edge>g[N];
    int level[N],iter[N];
    void init(int s,int t){
        for(int i=s;i<=t;++i)g[i].clear();
    }
    void add(int u,int v,int c){
        g[u].pb(edge(v,c,g[v].size()));
        g[v].pb(edge(u,0,g[u].size()-1));
    }
    bool bfs(int s,int t){
        memset(level,-1,sizeof(level));
        queue<int>q;
        q.push(s);
        level[s]=0;
        while(!q.empty()){
            int u=q.front();q.pop();
            for(int i=0;i<g[u].size();++i){
                edge e=g[u][i];
                if(level[e.to]<0&&e.cap>0){
                    level[e.to]=level[u]+1;
                    q.push(e.to);
                }
            }
        }
        return level[t]>=0;
    }
    int dfs(int u,int t,int f){
        if(u==t)return f;
        for(int &i=iter[u];i<g[u].size();++i){
            edge &e=g[u][i];
            if(level[u]<level[e.to]&&e.cap>0){
                int d=dfs(e.to,t,f>e.cap?e.cap:f);
                if(d>0){
                    e.cap-=d;
                    g[e.to][e.rev].cap+=d;
                    return d;
                }
            }
        }
        return 0;
    }
    int maxflow(int s,int t){
        int flow=0,f;
        while(bfs(s,t)){
            memset(iter,0,sizeof(iter));
            while((f=dfs(s,t,inf))>0)flow+=f;
        }
        return flow;
    }
    void print(int n){
        printf("%d ",maxflow(0,1));
        int ans=0;
        for(int i=1;i<=2*n;++i){
            for(int j=0;j<g[2*i+1].size();++j){
                edge e=g[2*i+1][j];
                if(e.to!=1&&e.cap>0&&e.cap<inf&&e.to!=2*i)
                    ans++;
            }
        }
        printf("%d
",ans);
        for(int i=1;i<=2*n;++i){
            for(int j=0;j<g[2*i+1].size();++j){
                edge e=g[2*i+1][j];
                if(e.to!=1&&e.cap>0&&e.cap<inf&&e.to!=2*i)
                    printf("%d %d %d
",i,e.to/2,inf-e.cap);
            }
        }
    }
}G;
int p,n;
struct nod{
    int c,in[15],out[15];
    void add(){
        scanf("%d",&c);
        for(int i=0;i<p;++i)
            scanf("%d",&in[i]);
        for(int i=0;i<p;++i)
            scanf("%d",&out[i]);
    }
    bool match(nod t){
        for(int i=0;i<p;++i)
            if(out[i]+t.in[i]==1)return 0;
        return 1;
    }
    bool ifs(){
        for(int i=0;i<p;++i)
            if(in[i]==1)return 0;
        return 1;
    }
    bool ift(){
        for(int i=0;i<p;++i)
            if(out[i]==0)return 0;
        return 1;
    }
}a[55];
void init(){
    G.init(0,2*n+2);
    for(int i=1;i<=n;++i){
        a[i].add();//s=0,t=1
        G.add(2*i,2*i+1,a[i].c);//将i拆成2*i和2*i+1
        if(a[i].ifs())G.add(0,2*i,inf);
        if(a[i].ift())G.add(2*i+1,1,inf);
    }
    for(int i=1;i<=n;++i)
    for(int j=1;j<=n;++j)
    if(i!=j&&a[i].match(a[j]))
        G.add(2*i+1,2*j,inf);
}
int main(void){
    while(~scanf("%d%d",&p,&n)){
        init();
        G.print(n);
    }
}