记录编号	346388	评测结果	AAAAAAAAAA
题目名称	[网络流24题] 骑士共存	最终得分	100
用户昵称	AntiLeaf	是否通过	通过
代码语言	C++	运行时间	0.643 s
提交时间	2016-11-12 08:55:38	内存使用	69.93 MiB

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=50010,maxe=3000010,dx[]={1,1,-1,-1,2,2,-2,-2},dy[]={-2,2,-2,2,-1,1,-1,1};
struct edge{int to,cap,prev;}e[maxe<<1];
void addedge(int,int,int);
void AddEdge(int,int,int);
int Dinic();
void bfs();
int dfs(int,int);
int n,m,last[maxn],len=0,id[210][210],cnt=0,s,t,d[maxn],now[maxn],q[maxn],head,tail,x,y;
bool ban[210][210]={{false}};
int main(){
#define MINE
#ifdef MINE
	freopen("knight.in","r",stdin);
	freopen("knight.out","w",stdout);
#endif
	memset(last,-1,sizeof(last));
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)id[i][j]=++cnt;
	s=cnt+1;t=s+1;
	for(int i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		ban[x][y]=true;
	}
	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++){
		if(ban[i][j])continue;
		if((i+j)&1){
			addedge(s,id[i][j],1);
			for(int k=0;k<8;k++){
				x=i+dx[k];y=j+dy[k];
				if(x<1||y<1||x>n||y>n||ban[x][y])continue;
				addedge(id[i][j],id[x][y],1);
			}
		}
		else addedge(id[i][j],t,1);
	}
	printf("%d",cnt-m-Dinic());
#ifndef MINE
	printf("\n-------------------------DONE-------------------------\n");
	for(;;);
#endif
	return 0;
}
void addedge(int x,int y,int z){
	AddEdge(x,y,z);
	AddEdge(y,x,0);
}
void AddEdge(int x,int y,int z){
	e[len].to=y;
	e[len].cap=z;
	e[len].prev=last[x];
	last[x]=len++;
}
int Dinic(){
	int ans=0;
	for(;;){
		bfs();
		if(d[t]==0x3f3f3f3f)break;
		memset(now,0,sizeof(now));
		ans+=dfs(s,~(1<<31));
	}
	return ans;
}
void bfs(){
	int x;
	memset(d,63,sizeof(d));
	head=tail=0;
	q[tail++]=s;
	d[s]=0;
	while(head!=tail){
		x=q[head++];
		for(int i=last[x];i!=-1;i=e[i].prev)if(e[i].cap>0&&d[e[i].to]>d[x]+1){
			d[e[i].to]=d[x]+1;
			q[tail++]=e[i].to;
		}
	}
}
int dfs(int x,int delta){
	if(x==t||!delta)return delta;
	int flow=0,f;
	for(int i=(now[x]?now[x]:last[x]);i!=-1;i=e[i].prev)if(e[i].cap>0&&d[e[i].to]==d[x]+1){
		now[x]=i;
		f=dfs(e[i].to,min(delta,e[i].cap));
		if(!f)continue;
		e[i].cap-=f;
		e[i^1].cap+=f;
		flow+=f;
		delta-=f;
		if(!delta)break;
	}
	return flow;
}