HDU3820 最小割

完全没思绪，看了大牛博客的构图才慢慢懂了。。太想不到啊

#include<cstdio>
#include<cstring>
const int inf=99999999;
int t,n,m,s,cnt,G,S,sum,T;
int a[60][60],b[60][60],head[10000],gap[10000],d[10000];
struct EDGE
{
	int to,f,nxt;
}edge[500000];
int min(int x,int y)
{
	return x<y?x:y;
}
void init()
{
	cnt=sum=0;
	memset(head,-1,sizeof(head));
	memset(gap,0,sizeof(gap));
	memset(d,0,sizeof(d));
}
void add(int x,int y,int c)
{
	edge[cnt].to=y;
	edge[cnt].f=c;
	edge[cnt].nxt=head[x];
	head[x]=cnt++;
}
void gao(int x,int y,int c)
{
	add(x,y,c);
	add(y,x,0);
}
int dfs(int x,int flow)
{
	if(x==t)
		return flow;
	int temp=flow;
	int  pos=t;
	int j;
	for(j=head[x];j!=-1;j=edge[j].nxt)
	{
		int y=edge[j].to;
		int c=edge[j].f;
		if(c>0&&d[x]==d[y]+1)
		{
			int temp_flow=dfs(y,min(c,temp));
			temp-=temp_flow;
			edge[j].f-=temp_flow;
			edge[j^1].f+=temp_flow;
			if(temp==0||d[s]==t)
				return flow-temp;
		}
		if(c>0&&pos>d[y])
			pos=d[y];
	}
	if(temp==flow)
	{
		if(!(--gap[d[x]]))
			d[s]=t;
		else
			gap[d[x]=pos+1]++;
	}
	return flow-temp;
}
int sap()
{
	int maxflow=0;
	gap[0]=t;
	while(d[s]<t)
	{
		maxflow+=dfs(s,inf);
	}
	return sum-maxflow;
}
int main()
{
	int pro=0;
	scanf("%d",&T);
	while(T--)
	{
		init();
		scanf("%d%d%d%d",&n,&m,&G,&S);
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++)
			{
				scanf("%d",&a[i][j]);
				sum+=a[i][j];
			}
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++)
			{
				scanf("%d",&b[i][j]);
				sum+=b[i][j];
			}
		s=2*n*m+1,t=s+1;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++)
			{
				int num=(i-1)*m+j;
				gao(num,num+n*m,inf);
			}
		for(int i=1;i<=n;i++)	
		for(int j=1;j<=m;j++)
			{
				int num=(i-1)*m+j;
				if((i+j)%2==0)
				{
					gao(s,num,a[i][j]);
					gao(num+n*m,t,b[i][j]);
				}
				else
				{
					gao(s,num,b[i][j]);
					gao(num+n*m,t,a[i][j]);
				}
			}
			for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++)
			{
				int num=(i-1)*m+j;
					if((i+j)%2==0)
					{
						if(i>1)
						gao(num,num-m+n*m,G);
						if(i<n)
						gao(num,num+m+n*m,G);
						if(j>1)
						gao(num,num-1+n*m,G);
						if(j<m)
						gao(num,num+1+n*m,G);
					}
					else
					{
						if(i>1)
						gao(num,num-m+n*m,S);
						if(i<n)
						gao(num,num+m+n*m,S);
						if(j>1)
						gao(num,num-1+n*m,S);
						if(j<m)
						gao(num,num+1+n*m,S);
					}	
			}
			printf("Case %d: %d\n",++pro,sap());
	}
	return 0;
}

，实在弱啊。