记录编号 52455 评测结果 AAAAAAAAAA
题目名称 运输问题1 最终得分 100
用户昵称 Gravatarcstdio 是否通过 通过
代码语言 C++ 运行时间 0.010 s
提交时间 2013-02-05 17:47:00 内存使用 0.39 MiB
显示代码纯文本
#include<iostream>
#include<cstdio>
using namespace std;
const int MAX=101;
int n;
int c[MAX][MAX]={0},f[MAX][MAX]={0},h[MAX]={0},e[MAX]={0};
//c为公路流量,f为前置流,h为高度,e为每个点余流,
int min(int a,int b){
	return (a<b)?a:b;
}
void initialize_preflow(void){//初始化,按算导上写的囧
	int i,s;
	s=0;
	for(i=0;i<n;i++){
		h[i]=0;
		e[i]=0;
	}
	h[s]=n;
	for(i=0;i<n;i++){
		if(c[s][i]!=0){//如果有公路
			f[s][i]=c[s][i];
			f[i][s]=-1*c[s][i];
			e[i]=c[s][i];
			e[s]-=c[s][i];
		}
	}
}
void push(int u,int v){//从u往v压
	int cf,df;//cf就是剩余流量,df就是个变量
	cf=c[u][v]-f[u][v];
	if(cf>0){//能压
		df=min(e[u],cf);
		f[u][v]+=df;
		f[v][u]=-1*f[u][v];
		e[u]-=df;
		e[v]+=df;
	}
}
void generic_push_relabel(void){
	initialize_preflow();//初始化
	bool flag=false;
	int i,j,min;
	int s=0,t=n-1;
	while(!flag){
		flag=true;
		for(i=0;i<n;i++){
			if(i==s||i==t) continue;
			if(e[i]>0){//有余流
				min=0xffffff;
				flag=false;
				for(j=0;j<n;j++){
					if(j==i) continue;
					if(c[i][j]-f[i][j]>0){//在残留图中
						if(h[j]==h[i]-1){
							push(i,j);//压
							goto NEXT;
						}
					}
					if(h[j]<min) min=h[i];
				}
				h[i]=min+1;//重标记
			}
			NEXT:;
		}
	}
}
int finalflow(void){
	int t=n-1,sum=0;
	int i;
	for(i=0;i<n;i++){
		if(i==t) continue;
		sum+=f[i][t];
	}
	return sum;
}
int main(){
	freopen("maxflowa.in","r",stdin);
	freopen("maxflowa.out","w",stdout);
	cin>>n;
	int i,j;
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			cin>>c[i][j];
		}
	}
	generic_push_relabel();
	int ans;
	ans=finalflow();
	cout<<ans<<endl;
	return 0;
}