题目链接:https://www.luogu.org/problem/P2774
思路:根据题意可以将图中的点分成两种点,一种(i+j)为奇数的点(红点),一种(i+j)为偶数的点(黑点),i为横坐标,j为纵坐标,设一个超级起点和一个超级终点,起点到红点建正向边和反向边,正向边权值为该红点的值,红点到黑点建边,黑点到终点建边,权值为黑点的值,然后跑网络流,输出为所有的点的总和减去最大流
AC代码:
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
#include<cmath>
#include<stack>
#include<map>
#include<set>
using namespace std;
#define LL long long
const int mod=1e9+7;
const int inf=0x3f3f3f3f;
const LL inff=0x3f3f3f3f3f3f3f3f;
const LL N=105;
const LL M=105;
#define MEF(x) memset(x,-1,sizeof(x))
#define ME0(x) memset(x,0,sizeof(x))
#define MEI(x) memset(x,inf,sizeof(x))
using namespace std;
LL poww(LL a,LL b)
{
LL ans=1;
while(b)
{
if(b&1)
ans=(ans*a)%mod;
a=(a*a)%mod;
b/=2;
}
return ans;
}
struct Edge
{
int v,w,next;
}edge[N*M*16];
int cnt,first[N*M],dis[N*M];
void init()
{
cnt=0;
MEF(first);
}
void add(int u,int v,int w)
{
edge[cnt].v=v;
edge[cnt].w=w;
edge[cnt].next=first[u];
first[u]=cnt++;
}
int n,m,al[N][M],sal;
int fx[]={0,1,0,-1,1,0,-1,0};
int bfs(int s,int t)//建层次图
{
queue<int> q;
MEF(dis);
q.push(s);
dis[s]=0;
while(!q.empty())
{
int x=q.front();
q.pop();
if(x==t)
{
return 1;
}
for(int i=first[x];i!=-1;i=edge[i].next)
{
int v=edge[i].v;
if(edge[i].w&&dis[v]==-1)
{
dis[v]=dis[x]+1;
q.push(v);
}
}
}
if(dis[t]==-1)
{
return 0;
}
else
{
return 1;
}
}
int dfs(int s,int flow)//flow 到大某一节点的最大流
{
if(s==n*m+1)//s==终点
{
return flow;
}
int ans=0;
for(int i=first[s];i!=-1;i=edge[i].next)
{
int v=edge[i].v;
if(edge[i].w&&dis[v]==dis[s]+1)
{
int f=dfs(v,min(flow-ans,edge[i].w));//(flow-ans) 某一节点的最大流-流过该节点的流量
edge[i].w-=f;
edge[i^1].w+=f;
ans+=f;
if(ans==flow)//某一节点流量达到最大流
{
return ans;
}
}
}
return ans;//某一节点的最终流量
}
int dinc(int s,int t)
{
int flow=0;
while(bfs(s,t))
{
flow+=dfs(s,inf);
}
return flow;
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);
cout.tie(0);
cin>>n>>m;
sal=0;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
cin>>al[i][j];
sal+=al[i][j];
}
}
init();
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
if((i+j)%2==0)
{
add(0,m*(i-1)+j,al[i][j]);
add(m*(i-1)+j,0,0);
for(int k=0;k<8;k+=2)
{
int lx=i+fx[k],ly=j+fx[k+1];
if(lx>=1&&lx<=n&&ly>=1&&ly<=m)
{
add(m*(i-1)+j,m*(lx-1)+ly,inf);
add(m*(lx-1)+ly,m*(i-1)+j,0);
}
}
}
else
{
add(m*(i-1)+j,n*m+1,al[i][j]);
add(n*m+1,m*(i-1)+j,0);
}
}
}
cout<<sal-dinc(0,n*m+1)<<endl;
return 0;
}