题意:给定一个网格,每个网格有选取代价和占据收益。每个点被占据,需要满足以下两个条件至少一个条件:1.被选取 2.邻近方格都被选取(有公共边被称为邻近) 不一定要占据所有方格,求最大收益。
第一直观感受和文理分科那道题很像,这类肯定用最小割,这种题一般都这样搞,但是建图是个大问题。这道题建出来的图要满足,如果一个点要保留收益,那么要么自己的花费边被割,要么邻近的被割,怎么建呢?
考虑先黑白染色,拆点,然后我们S连向黑色,容量为花费,黑色向自己的分身连收益边,并且黑色向相邻的白色点的分身连INF,黑色分身向白色连INF,白色分身向自己连收益,白色向T连花费。
仔细观察我们发现,确实能满足要求。
这种题一般都是套路,我也不知道怎么就要这么建边,不过可以总结出一些东西,比如一旦有关系,两者之间都会连INF以确保能彼此影响又不会被最小割割中,然后花费和收益,拆点怎么安排就要看具体的题目了。
上代码(与原题的输入不一样,是自己写的)
#include<bits/stdc++.h>
using namespace std;
#define N 25
#define INF 1e9
#define id(x,y) ((x-1)*m+y)
inline int read(){
int x=,f=; char a=getchar();
while(a<'' || a>'') {if(a=='-') f=-; a=getchar();}
while(a>='' && a<='') x=x*+a-'',a=getchar();
return x*f;
}
const int dir[][]={,,-,,,,,-};
int n,m,be[N][N],co[N][N],S,T,P,ans,d[],head[],cur[],cnt;
bool vis[];
queue<int>q;
struct edges{
int to,cap,flow,next;
}e[];
inline void insert(int u,int v,int c){
e[cnt]=(edges){v,c,,head[u]};head[u]=cnt++;
e[cnt]=(edges){u,,,head[v]};head[v]=cnt++;
}
inline bool bfs(){
memset(vis,,sizeof(vis));
vis[S]=; d[S]=; q.push(S);
while(!q.empty()){
int x=q.front(); q.pop();
for(int i=head[x];i>=;i=e[i].next){
if(!vis[e[i].to] && e[i].cap>e[i].flow)
d[e[i].to]=d[x]+,q.push(e[i].to),vis[e[i].to]=;
}
}
return vis[T];
}
int dfs(int x,int a){
if(x==T || !a) return a;
int f,flow=;
for(int& i=cur[x];i>=;i=e[i].next){
if(d[e[i].to]==d[x]+ && (f=dfs(e[i].to,min(a,e[i].cap-e[i].flow)))>)
e[i].flow+=f,flow+=f,e[i^].flow-=f,a-=f;
if(!a) break;
}
return flow;
}
inline int maxflow(){
int flow=;
while(bfs()){
for(int i=S;i<=T;i++) cur[i]=head[i];
flow+=dfs(S,INF);
}
return flow;
}
int main(){
memset(head,-,sizeof(head));
n=read(); m=read();
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
be[i][j]=read(),ans+=be[i][j];
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
co[i][j]=read();
S=; T=*n*m+; P=n*m;
for(int i=;i<=n;i++)
for(int j=;j<=m;j++){
int x,y;
if((i+j)%){
insert(S,id(i,j),co[i][j]),insert(id(i,j),id(i,j)+P,be[i][j]);
for(int k=;k<;k++){
x=i+dir[k][],y=j+dir[k][];
if(x< || x>n || y< || y>m) continue;
insert(id(i,j),id(x,y)+P,INF);
insert(id(i,j)+P,id(x,y),INF);
}
}
else insert(id(i,j)+P,id(i,j),be[i][j]),insert(id(i,j),T,co[i][j]);
}
ans-=maxflow();
printf("%d\n",ans);
return ;
}