记录编号 52506 评测结果 AAAAAAAAA
题目名称 [USACO 4.2] 完美的牛栏 最终得分 100
用户昵称 Gravatarcstdio 是否通过 通过
代码语言 C++ 运行时间 0.678 s
提交时间 2013-02-06 17:22:05 内存使用 1.55 MiB
显示代码纯文本
#include<iostream>
#include<cstdio>
#include<iomanip>
using namespace std;
const int MAX=403;
int n;
int c[MAX][MAX]={0},f[MAX][MAX]={0},h[MAX]={0},e[MAX]={0};
//c为公路流量,f为前置流,h为高度,e为每个点余流,
int s,t;
int min(int a,int b){
	return (a<b)?a:b;
}
void initialize_preflow(void){//初始化,按算导上写的囧
	int i;
	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){
	s=0,t=n-1;
	initialize_preflow();//初始化
	bool flag=false;
	int i,j,min;
	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 sum=0;
	int i;
	for(i=0;i<n;i++){
		if(i==t) continue;
		sum+=f[i][t];
	}
	return sum;
}
int main(){
	freopen("stall4.in","r",stdin);
	freopen("stall4.out","w",stdout);
	int cow,pen;
	scanf("%d%d",&cow,&pen);//牛,场
	n=cow+pen+2;
	int i,j,sum,temp;
	for(i=1;i<=cow;i++){
		scanf("%d",&sum);
		for(j=0;j<sum;j++){
			scanf("%d",&temp);
			c[i][cow+temp]=1;
		}
	}
	for(i=1;i<=cow;i++) c[0][i]=c[i][0]=1;
	for(i=cow+1;i<=cow+pen;i++) c[i][n-1]=c[n-1][i]=1;
	generic_push_relabel();
	int ans;
	ans=finalflow();
	cout<<ans<<endl;
	return 0;
}