分析(引用 BYvoid大牛的分析):
二分图点权最大独立集,转化为最小割模型,从而用最大流解决。
建模方法:
首先把棋盘黑白染色,使相邻格子颜色不同,所有黑色格子看做二分图 X集合中顶点,白色格子看做Y集合顶点,建立附加
源 S汇T。
1、从 S向 X 集合中每个顶点连接一条容量为格子中数值的有向边。
2、从 Y集合中每个顶点向 T 连接一条容量为格子中数值的有向边。
3、相邻黑白格子 Xi,Yj 之间从 Xi向 Yj连接一条容量为无穷大的有向边。
求出网络最大流,要求的结果就是所有格子中数值之和减去最大流量。
建模分析:
这是一个二分图最大点权独立集问题,就是找出图中一些点,使得这些点之间没有边相连,这些点的权值之和最大。独立集
与覆盖集是互补的,求最大点权独立集可以转化为求最小点权覆盖集(最小点权支配集)。最小点权覆盖集问题可以转化为最小割
问题解决。结论:最大点权独立集 = 所有点权 - 最小点权覆盖集 = 所有点权 - 最小割集 = 所有点权 - 网络最大流。
对于一个网络,除去冗余点(不存在一条 ST路径经过的点),每个顶点都在一个从 S到T的路径上。割的性质就是不存在从
S到 T的路径,简单割可以认为割边关联的非 ST节点为割点,而在二分图网络流模型中每个点必关联到一个割点(否则一定还有
增广路,当前割不成立),所以一个割集对应了一个覆盖集(支配集)。最小点权覆盖集就是最小简单割,求最小简单割的建模方
法就是把 XY集合之间的变容量设为无穷大,此时的最小割就是最小简单割了。
有关二分图最大点权独立集问题,更多讨论见《最小割模型在信息学竞赛中的应用》作者胡伯涛。
#include<stdio.h>
#include<iostream>
#include <string.h>
using namespace std;
const int oo=1e9;
/**oo 表示无穷大*/
const int mm=111111111;
/**mm 表示边的最大数量,记住要是原图的两倍,在加边的时候都是双向的*/
const int mn=999;
/**mn 表示点的最大数量*/
int node,src,dest,edge;
/**node 表示节点数,src 表示源点,dest 表示汇点,edge 统计边数*/
int ver[mm],flow[mm],next[mm];
/**ver 边指向的节点,flow 边的容量 ,next 链表的下一条边*/
int head[mn],work[mn],dis[mn],q[mn];
void prepare(int _node, int _src,int _dest)
{
node=_node,src=_src,dest=_dest;
for(int i=0; i<node; ++i)head[i]=-1;
edge=0;
}
/**增加一条 u 到 v 容量为 c 的边*/
void addedge( int u, int v, int c)
{
ver[edge]=v,flow[edge]=c,next[edge]=head[u],head[u]=edge++;
ver[edge]=u,flow[edge]=0,next[edge]=head[v],head[v]=edge++;
}
/**广搜计算出每个点与源点的最短距离,如果不能到达汇点说明算法结束*/
bool Dinic_bfs()
{
int i,u,v,l,r=0;
for(i=0; i<node; ++i)dis[i]=-1;
dis[q[r++]=src]=0;
for(l=0; l<r; ++l)
for(i=head[u=q[l]]; i>=0; i=next[i])
if(flow[i]&&dis[v=ver[i]]<0)
{
/**这条边必须有剩余容量*/
dis[q[r++]=v]=dis[u]+1;
if(v==dest) return 1;
}
return 0;
}
/**寻找可行流的增广路算法,按节点的距离来找,加快速度*/
int Dinic_dfs( int u, int exp)
{
if(u==dest) return exp;
/**work 是临时链表头,这里用 i 引用它,这样寻找过的边不再寻找*/
for( int &i=work[u],v,tmp; i>=0; i=next[i])
if(flow[i]&&dis[v=ver[i]]==dis[u]+1&&(tmp=Dinic_dfs(v,min(exp,flow[i])))>0)
{
flow[i]-=tmp;
flow[i^1]+=tmp;
/**正反向边容量改变*/
return tmp;
}
return 0;
}
int Dicnic_flow()
{
int i,ret=0,delta;
while(Dinic_bfs())
{
for(i=0; i<node; ++i)work[i]=head[i];
while(delta=Dinic_dfs(src,oo))
ret+=delta;
}
return ret;
}
int d[4][2]= {{0,-1},{0,1},{-1,0},{1,0}};
struct sa
{
int val,num;
bool col;
} a[205][205];
int main()
{
ios::sync_with_stdio(false);
int m,n;
while(cin>>m>>n)
{
memset(a,0,sizeof(a));
int cnt=0,sum=0;
for(int i=1; i<=m; i++)
{
for(int j=1; j<=n; j++)
{
cin>>a[i][j].val;
sum+=a[i][j].val;
if((i+j)%2==1)
{
a[i][j].col=1;
}
else a[i][j].col=0;
a[i][j].num=++cnt;
}
prepare(cnt+2,0,cnt+1);
for(int i=1; i<=m; i++)
{
for(int j=1; j<=n; j++)
{
if(a[i][j].col==1)
{
addedge(0,a[i][j].num,a[i][j].val);
for(int k=0; k<4; k++)
{
int x,y;
x=i+d[k][0];
y=j+d[k][1];
if(x>=1&&x<=m&&y>=1&&y<=n)
{
addedge(a[i][j].num,a[x][y].num,oo);
}
}
}
else if(a[i][j].col==0)
{
addedge(a[i][j].num,cnt+1,a[i][j].val);
}
}
}
}
cout<<sum-Dicnic_flow()<<endl;
}
return 0;
}