网络流
网络流未处理文件的代码:
/*656.网络流-矩阵计算 (10分) C时间限制:3000?毫秒?|? C内存限制:3000?Kb 题目内容: 有一个n行m列的整数矩阵A, 知道每行的和以及每列的和,还知道一些矩阵元素的约束如A[i][j]<x, 或者A[i][j]>y等, 判断该是否存在满足上述条件的可行矩阵。 输入描述 第一行是测试用例的数目c 每个测试用例的第一行是n,m 表示行和列 接下来一行是n个行和 接下来一行是m个列和 然后一行是约束个数k 接下来k行是约束,每个约束如 a b c d, 其中a,b是某个元素的行列坐标,c是一个字符(>,=,<), d是一个整数, 2 3 > 4 表示的意思是A[2][3]>4。矩阵左上角坐标规定为(1,1),所以一个约束的a为0,则表示b列所有的元素, 而如果b为0,则表示a行所有的元素。 输出描述 如果存在,则输出这个矩阵;否则输出“不存在” */ //#include<fstream> // 文件读写头文件 #include<iostream> #include<stdio.h> #include<string.h> #include <time.h> #include <stdlib.h> using namespace std; #define maxM 50000 #define maxN 500 #define inf 1<<30 struct node{ int u,v,f,next; }edge[maxM]; int head[maxN],p,lev[maxN],cur[maxN]; int que[maxM],tre[maxN],up[maxN][50],low[maxN][50]; bool bfs(int s,int t){ int qin=0,qout=0,u,i,v; memset(lev,0,sizeof(lev)); lev[s]=1,que[qin++]=s; while(qout!=qin){ u=que[qout++]; for(i=head[u];i!=-1;i=edge[i].next){ if(edge[i].f>0 && lev[v=edge[i].v]==0){ lev[v]=lev[u]+1,que[qin++]=v; if(v==t) return 1; } } } return lev[t]; } int dinic(int s,int t){ int qin,u,i,k,f; int flow=0; while(bfs(s,t)){ memcpy(cur,head,sizeof(head)); u=s,qin=0; while(1){ if(u==t){ for(k=0,f=inf;k<qin;k++) if(edge[que[k]].f<f) f=edge[que[i=k]].f; for(k=0;k<qin;k++) edge[que[k]].f-=f,edge[que[k]^1].f+=f; flow+=f,u=edge[que[qin=i]].u; } for(i=cur[u];cur[u]!=-1;i=cur[u]=edge[cur[u]].next) if(edge[i].f>0 && lev[u]+1==lev[edge[i].v]) break; if(cur[u]!=-1) que[qin++]=cur[u],u=edge[cur[u]].v; else{ if(qin==0) break; lev[u]=-1,u=edge[que[--qin]].u; } } } return flow; } void addedge(int u,int v,int f){ edge[p].u=u,edge[p].v=v,edge[p].f=f,edge[p].next=head[u],head[u]=p++; edge[p].u=v,edge[p].v=u,edge[p].f=0,edge[p].next=head[v],head[v]=p++; } void init1(int n,int m){ p=0,memset(head,-1,sizeof(head)),memset(tre,0,sizeof(tre)); for(int i=0;i<=n;i++) for(int j=0;j<=m;j++) up[i][j]=inf,low[i][j]=0; } bool buit(int n,int m){ for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) if(low[i][j]>up[i][j]) return 0; else{ tre[i]-=low[i][j],tre[j+n]+=low[i][j]; addedge(i,j+n,up[i][j]-low[i][j]); } return 1; } void limitflow(int s,int t,int n,int m){ int i,j,x,y; x=t+1,y=t+2; for(i=0;i<=t;i++){ if(tre[i]>0) addedge(x,i,tre[i]); else if(tre[i]<0) addedge(i,y,-tre[i]); } addedge(t,s,inf); dinic(x,y); for(i=head[x];i!=-1;i=edge[i].next) if(edge[i].f){ printf("不可能\n"); return ; } for(i=head[t];i!=-1;i=edge[i].next) if(edge[i].v==s) break; if(i<0){ printf("不可能\n"); return; } for(i=1;i<=n;i++){ for(j=1;j<m;j++) printf("%d ",edge[((i-1)*m+j)*2-1].f+low[i][j]); printf("%d\n",edge[i*m*2-1].f+low[i][j]); } // printf("\n"); } int main(){ // ifstream fin("file.txt");// 打开文件, "fin"可以自己命名 double start,finish; start=(double) clock(); //获得开始 时间 int cas,cas1,n,m,i,j,sum1,sum2; int u,v,d,f1,t1,f2,t2,s,t; char c[5]; scanf("%d",&cas); for(int tt=1;tt<=cas;tt++){ scanf("%d%d",&n,&m); s=0,t=n+m+1,sum1=0,sum2=0; init1(n,m); for(i=1;i<=n;i++) scanf("%d",&u),tre[s]-=u,tre[i]+=u,sum1+=u;; for(i=n+1;i<=n+m;i++) scanf("%d",&u),tre[i]-=u,tre[t]+=u,sum2+=u; scanf("%d",&cas1); while(cas1--){ scanf("%d%d%s%d",&u,&v,c,&d); f1=t1=u,f2=t2=v; if(u==0) f1=1,t1=n; if(v==0) f2=1,t2=m; for(i=f1;i<=t1;i++) for(j=f2;j<=t2;j++){ if(c[0]=='='){ low[i][j]=max(d,low[i][j]),up[i][j]=min(d,up[i][j]); }else if(c[0]=='>'){ low[i][j]=max(d+1,low[i][j]); }else if(c[0]=='<'){ up[i][j]=min(d-1,up[i][j]); } } } if(sum1==sum2 && buit(n,m)) limitflow(s,t,n,m); else printf("不可能\n"); } //finish=(double)clock(); //获得结束 时间 // 输出 时间差 毫秒 数: // printf("%.4f ms",finish-start); // system( "wmic process list MEMORY >>log.txt"); //用system函数执行wmic命令,生成一个log.txt文件, return 0; }