码迷,mamicode.com
首页 > 其他好文 > 详细

网络流,流水线模拟

时间:2018-09-23 22:35:49      阅读:211      评论:0      收藏:0      [点我收藏+]

标签:col   view   continue   网络流   stdin   map   txt   while   必须   

 ACM Computer Factory

题意:

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

solution:

技术分享图片
  1 #include<cstdio>
  2 #include<cstdlib>
  3 #include<cstring>
  4 #include<iostream>
  5 #define FRE(name) freopen(#name".txt","r",stdin);
  6 using namespace std;
  7 const int N=1.5e4+5,M=55;
  8 const int INF=0x3f3f3f3f;
  9 struct edge{int u,v,cap,next;}e[N*100],last[N*100];int tot=0;
 10 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];
 11 void add(int x,int y,int z){
 12     e[tot].u=x;e[tot].v=y;e[tot].cap=z;e[tot].next=head[x];head[x]=tot++;
 13     e[tot].u=y;e[tot].v=x;e[tot].cap=0;e[tot].next=head[y];head[y]=tot++;
 14 }
 15 bool isNULL(int *in){//与源点相连
 16     for(int i=1;i<=p;i++) if(in[i]==1) return 0;
 17     return 1;
 18 }
 19 bool isFULL(int *out){//与源点相连
 20     for(int i=1;i<=p;i++) if(!out[i]) return 0;
 21     return 1;
 22 }
 23 bool canLink(int *out,int *in){//普通节点之间的相连
 24     for(int i=1;i<=p;i++) if(in[i]+out[i]==1) return 0;
 25     return 1;
 26 }
 27 void mapping(){
 28     S=0,T=n+1;
 29     for(int i=1;i<=n;i++){
 30         if(isNULL(in[i])) add(S,i,val[i]);
 31         if(isFULL(out[i])) add(i,T,val[i]);
 32     }
 33     for(int i=1;i<=n;i++){
 34         for(int j=1;j<=n;j++){
 35             if(i==j) continue;
 36             if(canLink(out[i],in[j])) add(i,j,val[i]);
 37         }
 38     }
 39 }
 40 bool bfs(){
 41     for(int i=S;i<=T;i++) dep[i]=-1;
 42     int h=0,t=1;dep[S]=0;q[t]=S;
 43     while(h!=t){
 44         int x=q[++h];
 45         for(int i=head[x];~i;i=e[i].next){
 46             if(dep[e[i].v]==-1&&e[i].cap){
 47                 dep[e[i].v]=dep[x]+1;
 48                 if(e[i].v==T) return 1;
 49                 q[++t]=e[i].v;
 50             }
 51         }
 52     }
 53     return 0;
 54 }
 55 int dfs(int x,int f){
 56     if(x==T) return f;
 57     int used=0,t;
 58     for(int i=head[x];~i;i=e[i].next){
 59         if(e[i].cap&&dep[e[i].v]==dep[x]+1){
 60             t=dfs(e[i].v,min(e[i].cap,f));
 61             e[i].cap-=t;e[i^1].cap+=t;
 62             used+=t;f-=t;
 63             if(!f) return used;
 64         }
 65     }
 66     if(!used) dep[x]=-1;
 67     return used;
 68 }
 69 void dinic(){
 70     int ans=0;
 71     while(bfs()) ans+=dfs(S,INF);
 72     printf("%d ",ans);
 73 }
 74 void init(){
 75     tot=0;cnt=0;
 76     memset(head,-1,sizeof head);
 77 }
 78 int main(){
 79     //FRE(sh);
 80     while(scanf("%d%d",&p,&n)==2){
 81         init();
 82         for(int i=1;i<=n;i++){
 83             scanf("%d",&val[i]);
 84             for(int j=1;j<=p;j++) scanf("%d",&in[i][j]);
 85             for(int j=1;j<=p;j++) scanf("%d",&out[i][j]);
 86         }
 87         mapping();
 88         memcpy(last,e,sizeof e);
 89         dinic();
 90         for(int i=0;i<tot;i+=2){
 91             if(e[i].u==S||e[i].v==T) continue;
 92             if(last[i].cap!=e[i].cap){
 93                 res[++cnt][0]=e[i].u;
 94                 res[cnt][1]=e[i].v;
 95                 res[cnt][2]=e[i^1].cap;
 96             }
 97         }
 98         printf("%d\n",cnt);
 99         for(int i=1;i<=cnt;i++){
100             printf("%d %d %d\n",res[i][0],res[i][1],res[i][2]);
101         }
102     }
103     return 0;
104 }
View Code

 

网络流,流水线模拟

标签:col   view   continue   网络流   stdin   map   txt   while   必须   

原文地址:https://www.cnblogs.com/klaycf/p/9693680.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!