POJ 3436ACM Computer Factory 网络流

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 0x3f3f3f3f
struct 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;
}