POJ 3436ACM Computer Factory  网络流

一、题意

ACM的电脑由p个部分组成，只有p个部分都存在才能正常工作，现有n台机器，每台机器有三个参数，即：生产效率q，输入标准，输出标准，为了提高生产效率，对这些机器进行重新组网，求最大的生产效率。

二、解题思路

这题可以这样建图：(1)将每台机器看成两个点，入点和出点，对于每台机器的生产效率可以看作，入点和出点之间有一条容量（cap）为q的变。
（3）分别设置源点s，和汇点t，若某机器的输入标准中能接受初始状态（即：没有任何配件的状态）则s和该机器的入点有一条INF的边，同理，若某机器的输出为可接受状态（即：全为1的状态），则，该机器的出点与汇点t有一条INF的边。

（2）枚举两台不同机器的有序组合（i，j），若i的输出标准满足j的输入标准，则在i的出点与j的入点之间建一条cap为无穷大（INF）的边。

图建好后，s-t的最大流即为所求，用dinic算法进行处理。这一题要求输出答案的边的流量，用dinic算法处理后的残余网络，与原图进行对比即可。

三、代码

#include<iostream>#include<cstdio>#include<cstring>#include<cstdlib>#include<cmath>#include<algorithm>#include<map>#include<string>#include<queue>#include<vector>#include<list>//#pragma comment(linker,"/STACK:1024000000,1024000000")using namespace std;const int MAX_V=205;#define INF 0x3f3f3f3fstruct edge{int to,cap,rev;};vector<edge>G[MAX_V];int level[MAX_V];int iter[MAX_V];void add_edge(int from,int to,int cap){    G[from].push_back((edge){to,cap,G[to].size()});    G[to].push_back((edge){from,0,G[from].size()-1});}void bfs(int s){    memset(level,-1,sizeof level);    queue<int> que;    level[s]=0;    que.push(s);    while(!que.empty())    {        int v=que.front();que.pop();        for(int i=0;i<G[v].size();i++)        {            edge &e=G[v][i];            if(e.cap>0&&level[e.to]<0)            {                level[e.to]=level[v]+1;                que.push(e.to);            }        }    }}int dfs(int v,int t,int f){    if(v==t) return f;    for(int &i=iter[v];i<G[v].size();i++)    {        edge& e=G[v][i];        if(e.cap>0&&level[v]<level[e.to])        {            int d=dfs(e.to,t,min(f,e.cap));            if(d>0)            {                e.cap-=d;                G[e.to][e.rev].cap+=d;                return d;            }        }    }    return 0;}int max_flow(int s,int t){    int flow=0;    for(;;)    {        bfs(s);        if(level[t]<0) return flow;        memset(iter,0,sizeof iter);        int f;        while((f=dfs(s,t,INF))>0)        {            flow+=f;        }    }}int p,n;int mac[55][15][2],cap[55];vector<edge> res;int main(){    while(~scanf("%d %d",&p,&n)&&p+n)    {        res.clear();        for(int i=0;i<MAX_V;i++) G[i].clear();        for(int i=1;i<=n;i++)        {            scanf("%d",&cap[i]);            for(int j=1;j<=p;j++) scanf("%d",&mac[i][j][0]);            for(int j=1;j<=p;j++) scanf("%d",&mac[i][j][1]);            add_edge(2*i-1,2*i,cap[i]);            int flag=0;for(int j=1;j<=p;j++) if(mac[i][j][0]==1) {flag=1;break;} if(!flag) add_edge(0,2*i-1,INF);            flag=0;for(int j=1;j<=p;j++) if(mac[i][j][1]==0) {flag=1;break;} if(!flag) add_edge(2*i,2*n+1,INF);        }        for(int i=1;i<=n;i++)            for(int j=1;j<=n;j++)        {            if(i==j) continue;            int flag=1;            for(int k=1;k<=p;k++)            {                if(mac[i][k][1]==0&&mac[j][k][0]==1) {flag=0;break;}                if(mac[i][k][1]==1&&mac[j][k][0]==0) {flag=0;break;}            }            if(flag) add_edge(2*i,2*j-1,INF);        }        int ans=max_flow(0,2*n+1);        for(int i=2;i<=2*n;i+=2)        {            for(int j=0;j<G[i].size();j++)            {                edge e=G[i][j];                if(e.to!=2*n+1&&e.to!=i-1&&e.cap!=INF) res.push_back((edge){i/2,e.to/2+1,INF-e.cap});            }        }        printf("%d %d\n",ans,res.size());        for(int i=0;i<res.size();i++)        {            printf("%d %d %d\n",res[i].to,res[i].cap,res[i].rev);        }    }    return 0;}