//最大流(邻接阵形式)-- 求网络最大流,邻接阵形式 //求网络最大流,邻接阵形式 //返回最大流量,flow返回每条边的流量 //传入网络节点数n,容量mat,源点source,汇点sink
const int MAXN = 302; const int inf = 1000000000; int max_flow(int n,int mat[][MAXN],int source,int sink,int flow[][MAXN]) { int pre[MAXN],que[MAXN],d[MAXN],p,q,t,i,j; if (source==sink) return inf; for (i=0;i<n;i++) for (j=0;j<n;flow[i][j++]=0); for (;;) { for (i=0;i<n;pre[i++]=0); pre[t=source]=source+1,d[t]=inf; for (p=q=0;p<=q&&!pre[sink];t=que[p++]) for (i=0;i<n;i++) if (!pre[i]&&(j=mat[t][i]-flow[t][i])) pre[que[q++]=i]=t+1,d[i]=d[t]<j?d[t]:j; else if (!pre[i]&&(j=flow[i][t])) pre[que[q++]=i]=-t-1,d[i]=d[t]<j?d[t]:j; if (!pre[sink]) break; for (i=sink;i!=source;) if (pre[i]>0) flow[pre[i]-1][i]+=d[sink],i=pre[i]-1; else flow[i][-pre[i]-1]-=d[sink],i=-pre[i]-1; } for (j=i=0;i<n;j+=flow[source][i++]); return j; }