网络流二十四题 ————(十七）、P3355 骑士共存问题 二分图最小割，转化为网络流最小割 -＞ 最大流求解

经典套路了，看到这种限制条件比较少的情况，我们很容易想到用最小割的思想。

先把所有点都放上骑士，然后横纵坐标为奇数的点在左部图，横纵坐标为偶数的点在右部图(二分图经典套路)。

相互不可同时出现一定不在同一部图(因为横纵坐标差3)，然后从左部图往其可达点连边，容量为inf。

然后起点S连向左部图每个点，容量为1.  右部图每个点连向终点，容量为1.

这样求最小割后，S一定不能通过中间达到右边。(中间容量为inf一定不会被割掉)。

留下来的每条边都代表一个点放骑士，且是最大数量可放置的骑士。

最小割转化为网络流最大流求解即可。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int M = 2500000+7;
const int N = 50000+7;
struct Dinic{
	//N个点，M条边 
	//add建双向边，然后D.gao ,最后输出maxflow  
	#define inf 0x3f3f3f3f
	ll maxflow;int s,t,n;
	int v[N],pre[N],d[N],now[N];
	ll incf[N];
	int head[N],cnt=1;
	struct EDGE{int to,nxt;ll w;}ee[M*2];
	inline void AD(int x,int y,ll w){ee[++cnt].nxt=head[x],ee[cnt].w=w,ee[cnt].to=y,head[x]=cnt;}
	inline void add(int x,int y,ll w){AD(x,y,w);AD(y,x,0);}
	inline bool bfs()//在残量网络上构造分层图 
	{
		for(int i=0;i<=n;i++)d[i]=0;
		queue<int>q;
		q.push(s);d[s]=1;
		while(q.size())
		{
			int x=q.front();q.pop();
			for(int i=head[x];i;i=ee[i].nxt)
			{
				int y=ee[i].to;ll w=ee[i].w;
				if(w&&!d[y])
				{
					q.push(y);
					d[y]=d[x]+1;
					if(y==t)return 1;
				}
			}
		}
		return false;
	}
	inline int dinic(int x,int flow)
	{
		if(x==t)return flow;
		ll rest = flow,k;
		for(int i=now[x];i&&rest;i=ee[i].nxt)
		{
			int y=ee[i].to;ll w= ee[i].w;
			now[x]=i;
			if(w&&d[y]==d[x]+1)
			{
				k=dinic(y,min(rest,w));
				if(!k)d[y]=0;//剪枝，去掉增广完毕的点
				ee[i].w-=k;
				ee[i^1].w+=k;
				rest-=k; 
			}
		}
		return flow - rest;
	}
	inline void gao()
	{
		int flow=0;
		while(bfs())
		{
			for(int i=0;i<=n;i++)now[i]=head[i];
			while(flow=dinic(s,inf))maxflow+=flow;
		}
	}
	inline void init(int nn,int S,int T)
	{
		cnt=1;maxflow=0;
		for(int i=0;i<=n;i++)head[i]=0;
		s=S,t=T,n=nn;
	}
}D;
//求解二分图匹配问题时的复杂度是：m*sqrt(n)的

int mp[1100][1100],n;
int id(int x,int y){
	return (x-1)*n+y;
}
int dx[8]={-1,-1,1,1,2,2,-2,-2},dy[8]={2,-2,2,-2,-1,1,-1,1};
int no[1000007];
int main()
{
	int m;
	cin>>n>>m;
	for(int i=1;i<=m;i++){
		int x,y;
		cin>>x>>y;
		no[id(x,y)]=1;
	}
	int s=n*n+1,t=n*n+2;
	D.init(n*n+2,s,t);
	int ans=n*n;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			int x=id(i,j);
			if((i+j)%2==1){
				if(!no[x])D.add(s,x,1);
				else ans--;
			}
			else{
				if(!no[x])D.add(x,t,1);
				else ans--;
			}
		}
	}
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++){
		int u=id(i,j),x,y;
		if(no[u]||(i+j)%2==0)continue;
		for(int k=0;k<8;k++){
			x=i+dx[k],y=j+dy[k];
			int v=id(x,y);
			if(x<1||x>n||y<1||y>n||no[v])continue;
			D.add(u,v,inf);
		}
	}
	D.gao();
	cout<<ans-D.maxflow<<endl;
	return 0;
}