hdu3280 Golden Eggs(网络流最小割)

这题难点在于建图。参考了网上大牛的建图，膜拜下orz

首先这类题，分成两个集合A，B

集合内部的点都是不相交的。

然后对于每一个点，有三种选择，金，银，无。

每一个格拆为两个点，金和银，之间为正无穷，因为金银只能选一种。

对于A集合 为源-金-银-汇

对于B集合 为源-银-金-汇。

A、B集合之间为相邻的点连权值对应的边。

最小割割去的都是不要的，所以答案为Sum-最小割。

这个题会了，黑白格染色就简单许多了。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#include<queue>

using namespace std;
int T,N,M,G,S;
const int maxn = 55,INF= 100000;
int g[maxn][maxn],s[maxn][maxn];
struct EDGE{
    int u,v,w,next;
}edge[100010];

int cnt;
int head[maxn*maxn*2],mark[maxn*maxn*2];
int src,des;
int dx[4]={0,0,1,-1};
int dy[4]={1,-1,0,0};
int Sum=0;
#define MIN(a,b) a<b?a:b

void addedge(int u,int v,int w){
    edge[cnt].u=u;edge[cnt].v=v;edge[cnt].w=w;
    edge[cnt].next=head[u];head[u]=cnt++;
    edge[cnt].u=v;edge[cnt].v=u;edge[cnt].w=0;
    edge[cnt].next=head[v];head[v]=cnt++;
}
bool check(int x,int y){
    if(x<0||x>=N||y<0||y>=M)return 0;
    return 1;
}
void build_map(){
    memset(head,-1,sizeof(head));
    cnt=0;
    src=N*M*2,des=N*M*2+1;
    for(int i=0;i<N;i++){
        for(int j=0;j<M;j++){
            if((i+j)%2==0){
                addedge(src,i*M+j,g[i][j]);
                addedge(i*M+j+N*M,des,s[i][j]);
                addedge(i*M+j,i*M+j+N*M,INF);
            }
            else{
                addedge(src,i*M+j,s[i][j]);
                addedge(i*M+j+N*M,des,g[i][j]);
                addedge(i*M+j,i*M+j+N*M,INF);
            }
            for(int k=0;k<4;k++){
                int xx=i+dx[k],yy=j+dy[k];
                if(check(xx,yy)){
                    if((i+j)%2==0){
                        addedge(i*M+j,xx*M+yy+N*M,G);
                    }
                    else{
                        addedge(i*M+j,xx*M+yy+N*M,S);
                    }
                }
            }
        }
    }
}
bool BFS(){
    queue<int>q;
    q.push(src);
    memset(mark,-1,sizeof(mark));
    mark[src]=1;
    while(!q.empty()){
        int k=q.front();q.pop();
        for(int i=head[k];i!=-1;i=edge[i].next){
            int v = edge[i].v;
            if(edge[i].w>0&&mark[v]==-1){
                mark[v]=mark[k]+1;
                q.push(v);
            }
        }
    }
    return mark[des]!=-1;
}
int dfs(int x,int delta){
    int cost=0,min;
    if(x==des)return delta;
    for(int i=head[x];i!=-1;i=edge[i].next){
        int v=edge[i].v;
        if(mark[v]==mark[x]+1&&edge[i].w>0){
            min=dfs(v,MIN(edge[i].w,delta-cost));
            if(min>0){
                edge[i].w-=min;
                edge[i^1].w+=min;
                cost+=min;
                if(cost==delta)break;
            }
            else{
                mark[v]=-1;
            }
        }
    }
    return cost;
}
int ica;
void Dinic(){
    int flow=0;
    while(BFS()){
        flow+=dfs(src,INF);
    }
    printf("Case %d: %d\n",ica++,Sum-flow);
}
int main(){
    
    scanf("%d",&T);
    ica=1;
    while(T--){
        Sum=0;
        scanf("%d%d%d%d",&N,&M,&G,&S);
        for(int i=0;i<N;i++){
            for(int j=0;j<M;j++){
                scanf("%d",&g[i][j]);
                Sum+=g[i][j];
            }
        }
        for(int i=0;i<N;i++){
            for(int j=0;j<M;j++){
                scanf("%d",&s[i][j]);
                Sum+=s[i][j];
            }
        }
        build_map();
        Dinic();
    }
    
    return 0;
}