标签:col view continue 网络流 stdin map txt while 必须
一共有N个机器,每个机器有P个元素,对应输入的时候输入N个机器的信息,第一个数表示这个机器可以一共能够生产多少产物,接下来2p个元素,前p个元素:其中有三种数值,1,2,0,分别表示必须有这个位子的组件,可有可没有这个位子的组件,以及不能有这个位子的组件。一个中间产物想在这个机器上进行生产接下来的产物,必须要满足这些条件才行,对应在这个机器上产生的产物的各个组件的结果是后p个元素。
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 }
标签:col view continue 网络流 stdin map txt while 必须
原文地址:https://www.cnblogs.com/klaycf/p/9693680.html