版权声明:本文为博主原创文章,未经博主允许必须转载。 https://blog.csdn.net/qq_35950004/article/details/88423801
题目
这个题,求最大值,那么先把所有值加起来-最小割。
那么我们的二元关系应该是选同样有
的代价。
相异代价-相同代价 < 0。
需要反转,中间代价为C之和。
具体看:pty《浅析一类最小割问题》
AC Code:
#include<bits/stdc++.h>
#define maxn 30010
#define maxm maxn * 100
#define inf 0x3f3f3f3f
using namespace std;
int n,m,mark[205][105],S,T,tot;
int a[205][205],b[205][205],c[205][205],d[205][205];
int dir[2][2]={{-1,0},{0,-1}};
int dis[maxn];
int info[maxn],Prev[maxm],to[maxm],cap[maxm],cnt_e=1;
inline void Node(int u,int v,int c){ Prev[++cnt_e]=info[u],info[u]=cnt_e,to[cnt_e]=v,cap[cnt_e]=c; }
inline void Line(int u,int v,int c,int d=0){ Node(u,v,c),Node(v,u,d); }
int aug(int now,int Max)
{
if(now == T) return Max;
int inc , st = Max;
for(int i=info[now];i;i=Prev[i])
if(cap[i] && dis[to[i]]+1 == dis[now])
{
inc = aug(to[i],min(cap[i] , st));
if(inc) st -= inc , cap[i] -= inc , cap[i^1] += inc;
else dis[to[i]] = 0;
if(!st) break;
}
return Max - st;
}
bool BFS()
{
static queue<int>q;
memset(dis,-1,sizeof dis);
q.push(T),dis[T]=0;
for(int now;!q.empty();)
{
now = q.front() , q.pop();
for(int i=info[now];i;i=Prev[i])
if(cap[i^1] && dis[to[i]]==-1)
{
dis[to[i]] = dis[now] + 1;
q.push(to[i]);
}
}
return dis[S] != -1;
}
int main()
{
scanf("%d%d",&n,&m);
S = ++tot , T = ++tot;
int ans = 0;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
mark[i][j] = ++tot;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&a[i][j]);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&b[i][j]);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&c[i][j]);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
ans += a[i][j] + b[i][j];
if((i+j)&1) Line(S,mark[i][j],a[i][j]),Line(mark[i][j],T,b[i][j]);
else Line(mark[i][j],T,a[i][j]),Line(S,mark[i][j],b[i][j]);
for(int k=0,x,y;k<2;k++)
if((x=i+dir[k][0])>=1 && x<=n && (y=j+dir[k][1])>=1 && y<=m)
{
int tmp;
ans += (tmp = c[x][y] + c[i][j]);
Line(mark[i][j],mark[x][y],tmp,tmp);
}
}
for(;BFS();)
ans -= aug(S,inf);
printf("%d\n",ans);
}