luogu P2774 方格取数问题

关于构图:

首先,这是一道最小割的题目。

我们可以想到一个性质,最大点独立集=点数-最小点覆盖,通俗来说最大和=全局和-舍弃和,而舍弃和=最小割=最大流。


我们可以看到上面的那一幅图,就会有所启发。

当一个点(i,j)满足(i+j)mod 2=1时,为白色的点,否则为黑色的点。

对于任意一个点,我们不能取它上下左右的点,于是就可以得到上面的那一幅图,白色的点为一个集合,黑色的点为一个集合,两个集合中的点都满足这一规定。

你或许会问,这可以用贪心(ans=max(sigma白,sigma黑))解决啊。

但是,一个白色的点可以与一个不与它相临的黑色的点连边,这时,就得用到最小割。一个白色的点与与它相临的黑色的点建边,但这些便是需要“割”掉的。有了这些边,就可以从一个白色的点流向一个不与它相临的黑色的点了。

构图:

1.首先,需要一个超级源点和汇点。

2.源点连向每一个黑色的点,流量为这个点的值;同理,每一个白色的点都连向汇点,流量为这个点的值。

3.每一个黑色的点都连向与它相临的白色的点,流量为INF,这些边需要“割”掉一些。

最后就可以根据最小割=最大流这个定理&最大和=全局和-舍弃和求解了。


关于代码:

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
#define INF 2147483647
using namespace std; 
queue<int> f;
	const int fx[4]={0,0,1,-1};
	const int fy[4]={1,-1,0,0};
	struct node{int x,y,z,next;} a[200000];
	int last[1000];
	int n,m,len=-1,ans=0,st,ed;
void ins(int x,int y,int z)
{
	a[++len].x=x;a[len].y=y;a[len].z=z;a[len].next=last[x];last[x]=len;
}
int h[1000];
bool bfs()
{
	memset(h,0,sizeof(h));
	h[st]=1;
	f.push(st);
	while(!f.empty())
	{
		int x=f.front();
		for(int i=last[x];i>=0;i=a[i].next)
		{
			int y=a[i].y;
			if(a[i].z>0&&h[y]==0)
			{
				h[y]=h[x]+1;
				f.push(y);
			}
		}
		f.pop();
	}
	if(h[ed]) return true; else return false;
}
int dfs(int x,int f)
{
	int s=0,t;
	if(x==ed) return f;
	for(int i=last[x];i>=0;i=a[i].next)
	{
		int y=a[i].y;
		if(a[i].z>0&&h[y]==h[x]+1&&f>s)
		{
			s+=(t=(dfs(y,min(f-s,a[i].z))));
			a[i].z-=t;
			a[i^1].z+=t;
		}
	}
	if(!s) h[x]=0;
	return s;
}
int dinic()
{
	int sum=0;
	while(bfs())
		sum+=dfs(st,INF);
	return sum;
}
int main()
{
	int x,id;
	scanf("%d %d",&n,&m);
	st=0,ed=n*m+1;
	memset(last,-1,sizeof(last));
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&x);
			ans+=x;
			id=(i-1)*m+j;
			if((i+j)&1) ins(st,id,x),ins(id,st,0); else ins(id,ed,x),ins(ed,id,0);
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			if((i+j)&1)
			{
				for(int k=0;k<=3;k++)
				{
					int t1=i+fx[k],t2=j+fy[k];
					if(t1<=0||t1>n||t2<=0||t2>m) continue;
					id=(i-1)*m+j;
					ins(id,(t1-1)*m+t2,INF),ins((t1-1)*m+t2,id,0);
				}
			}
	printf("%d",ans-dinic());
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值