HDU 3820 Golden Eggs （最小割）

大致题意：

    给出一块n*m的区域，已知区域中的每个格子都可以放金蛋或者是银弹，已知每个格子放金蛋或者是放银弹时可以得到的收益，分别用矩阵map1[][]，map2[][]表示。如果相邻的两个格子放的蛋相同的话会扣去一定的收益，金蛋的话会扣去g，银弹s。求收益的最大值是多少。

 

大致思路：

    首先，对于某个确定的位置，只能选择金蛋或者银弹的其中一个放上去，从这里可以想到和二分图最大点权独立集有关。其次，如果相邻的位置放的蛋相同的话需要付出一定的花费这个又和hdu 3657：Game有些许相似之处。所以我们这么构图：

 

按照黑白染色把矩阵中的位置分为A,B两个集合，每个集合中的位置拆为k,k'两个点。对于A集合，从源点向k连边，容量为它放金蛋是的收获，从k再向k'连边容量设为inf，从k'连到汇点，容量为其放银弹时的收益。这样设的目的就是为了使得割必须在源点到k或者k'到汇点上。   

对于B集合，我们做一些改变，把从源点到k的容量设为放银弹的收益，k'到汇点的容量设为其放金蛋的收益,k和k'之间容量为inf的边不变。至于为什么看下一步。

然后再接下来一步加边就是关键了！我们知道如果A集合中的金蛋和B集合中的金蛋相邻的话会付出一定的话费，为了表示设A集合中的某位置a是金蛋，我们会连接源点->a。为了表示B集合中位置b是金蛋的，我们会连接b'->汇点。好的，在这里我们假设b位置和a相邻。为了表示a位置和b位置同时放金蛋需要G的花费，我们连接a->b'并把其容量设为G，对于银弹我们做同样的处理。至于为什么，好好把hdu3567刷掉就全明白了

接下来，用所有的map1[i][j]+map2[i][j]-最小割得到的就是答案。

 

 

#include<cstdio>
#include<cstring>
#include<map>
#include<vector>
#include<cmath>
#include<cstdlib>
#include<stack>
#include<queue>
#include <iomanip>
#include<iostream>
#include<algorithm>
using namespace std ;
const int N=10000 ;
const int M=N*6 ;
const int inf=0x3f3f3f3f;
struct node
{
	int u ,v,c,next;
}edge[M];
int head[N],pre[N],gap[N],cur[N],dis[N];
int g1[60][60],g2[60][60];
int top ;

void add(int u,int v,int c)
{
	edge[top].u=u;
	edge[top].v=v;
	edge[top].c=c;
	edge[top].next=head[u];
	head[u]=top++;
	edge[top].u=v;
	edge[top].v=u;
	edge[top].c=0;
	edge[top].next=head[v];
	head[v]=top++;
}

int SAP(int s,int t,int n)
{
	int u=s,v,nect,mindis,minflow,flow=0,i;
	 memset(gap,0,sizeof(gap));
	 memset(dis,0,sizeof(dis));
	 memset(pre,-1,sizeof(pre));
   for( i = 0 ; i < n ; i++)  cur[i]=head[i];	  
	  gap[0]=n;
	 while(dis[s] <n)
	 {
	 	     if(u==t)
	 	     {	 	     	
	 	     	     minflow=inf;
	 	     	     for( i = s ; i!=t ; i = edge[cur[i]].v)  
	 	     	      {	 	     	       
						     if( minflow > edge[cur[i]].c )
						     {
						     	   minflow = edge[cur[i]].c;
						     	   nect=i;
						     }						    
	 	     	      }
					 for( i = s ; i!=t ; i = edge[cur[i]].v)  
					 {
					 	    edge[cur[i]].c -=minflow;
					 	    edge[cur[i]^1].c +=minflow;
					 }  
					 u=nect;
	 	     	    flow+=minflow;
	 	     	   
	 	     }
	 	     for( i = cur[u] ; i!=-1 ; i=edge[i].next)
	 	     {
	 	         v=edge[i].v;
				 if( edge[i].c >0 && dis[u]==dis[v]+1)
				   break;
	 	     }
	 	     if(i!=-1)
	 	     {
	 	     	    cur[u]=i;
	 	     	    pre[v]=u;
	 	     	    u=v;
	 	     }
	 	     else
	 	     {   
	 	         if(--gap[dis[u]]==0) break;
	 	     	   mindis=n;
	  	   	  for( i =head[u] ; i!=-1;i=edge[i].next)
	  	       {
	  	       	    v=edge[i].v;
	  	       	     if(edge[i].c >0 && dis[v] < mindis)
	  	       	     {
	  	       	     	     mindis=dis[v];
	  	       	     	     cur[u]=i;
	  	       	     }
	  	       }
	  	      
	  	       gap[dis[u]=mindis+1]++ ;
	  	        if(u!=s) u=pre[u]; 	 	     	 	     	
	 	     }
	 }
	 return flow ;
}

int main()
{
	
	int n,m,t,G,S,x,cases=0;
	scanf("%d",&t);
	while(t--)
	{
		  scanf("%d%d%d%d",&n,&m,&G,&S);
		 top=0;
		 memset(head,-1,sizeof(head));
		 int s=0,t=2*n*m+1,sum=0;
	  	 for(int i = 1 ; i <= n; i++)
		    for(int j = 1; j <= m; j++)
		       {
		       	      scanf("%d",&g1[i][j]);
		       	      sum+=g1[i][j];
		       }
		  for(int i = 1 ; i <= n; i++)
		    for(int j = 1; j <= m; j++)
		       {
		       	      scanf("%d",&g2[i][j]);
		       	      sum+=g2[i][j];
		       }     
		     for(int i = 1 ; i<= n ; i++)
			   for(int j = 1 ; j <= m ; j++)
			      if((i+j)%2==0)
				  {
				  	     int tmp=(i-1)*m+j;
				  	       add(s,tmp,g1[i][j]);
				  	       add(tmp+n*m,t,g2[i][j]);
				  	       add(tmp,tmp+n*m,inf);
				  	   
				  	    if(i<n)   add(tmp,tmp+m+n*m,G);
				  	    if(i>1)   add(tmp,tmp-m+n*m,G);
				  	    if(j<m)   add(tmp,tmp+1+n*m,G);
				  	    if(j>1)   add(tmp,tmp-1+n*m,G);
				  }
				  else
				  {
				  	    int tmp=(i-1)*m+j;
				  	       add(s,tmp,g2[i][j]);
				  	       add(tmp+n*m,t,g1[i][j]);
				  	       add(tmp,tmp+n*m,inf);
				  	   
				  	    if(i<n)   add(tmp,tmp+m+n*m,S);
				  	    if(i>1)   add(tmp,tmp-m+n*m,S);
				  	    if(j<m)   add(tmp,tmp+1+n*m,S);
				  	    if(j>1)   add(tmp,tmp-1+n*m,S);
				  					  				  	
				  } 
			printf("Case %d: %d\n",++cases,sum-SAP(s,t,t+1));	  
	}
	  return 0;
}