ICode9
精准搜索请尝试： 精确搜索
首页
编程语言>
数据库
系统相关
互联网
其他分享
Python
Java
JavaScript
android
PHP
首页 > 其他分享> 文章详细

网络流棋盘模型 | P3355 骑士共存问题 P4304 [TJOI2013]攻击装置

2022-09-09 13:33:36  阅读：173  来源： 互联网

标签：P3355 head return val int gap TJOI2013 P4304 now


题面（骑士共存问题）

在一个 \(n \times n\) 个方格的国际象棋棋盘上，马（骑士）可以攻击的棋盘方格如图所示。棋盘上某些方格设置了障碍，骑士不得进入。

image

对于给定的 \(n \times n\) 个方格的国际象棋棋盘和 \(m\) 个障碍标志，计算棋盘上最多可以放置多少个骑士，使得它们彼此互不攻击。

对于全部的测试点，保证 \(1 \leq n \leq 200\)，\(0 \leq m \lt n^2\)。

思路

二分图的最大独立集。

先对这个棋盘进行黑白染色，然后发现，同颜色的格子的马永远无法互相攻击。

首先，先创立源点 \(s\)，汇点 \(t\)。

然后对于白点 \((i,j)\)，连边 \((s,(i,j),1)\)。对于黑点 \((k,l)\)，连边 \(((k,l),t,1)\)，对于可以互相攻击的点 \((a,b),(c,d)\)，连边 \(((a,b),(c,d),+\infty)\)。

然后求最大独立集。

代码

P3355 骑士共存问题

#include <bits/stdc++.h>
#define int long long
using namespace std;

int n,m,s,t;

struct edge{
	int from,to,val;
}e[2000001];int head[2000001],cur[2000001],siz=1;
void addedge(int x,int y,int z){
	e[++siz].to=y,e[siz].val=z;
	e[siz].from=head[x],head[x]=siz;
}
void add(int x,int y,int z){
	addedge(x,y,z);
	addedge(y,x,0);
}

int gap[2000001];
bool bfs(){
	memset(gap,0,sizeof(gap));
	fill(gap+1,gap+1+n,0);
	queue<int> q;
	q.push(s);
	gap[s]=1;
	while(!q.empty()){
		int now=q.front();
		q.pop();
		for(int i=head[now];i;i=e[i].from){
			int u=e[i].to;
			if(e[i].val&&!gap[u]){
				gap[u]=gap[now]+1;
				q.push(u);
			}
		}
	}
	return (gap[t]);
}
int dfs(int now,int val){
	if(now==t) return val;
	for(int &i=cur[now];i;i=e[i].from){
		int u=e[i].to;
		if(e[i].val&&gap[now]+1==gap[u]){
			int F=dfs(u,min(e[i].val,val));
			if(F){
				e[i].val-=F;
				e[i^1].val+=F;
				return F;
			}
		}
	}
	return 0;
}
int dinic(){
	int ret=0;
	while(bfs()){
		copy(head,head+1+n,cur);
		int F=0;
		while(F=dfs(s,10000000000000)){
			ret+=F;
		}
	}
	return ret;
}

int knight[100000][1000];

int mk(int x,int y){
	x--;
	return (x*n+y);	
}

int delta[][2]={{1,2},{1,-2},{-1,2},{-1,-2},{2,1},{2,-1},{-2,1},{-2,-1}};

bool valid(int x){
	return (x>=1)&&(x<=n);
}

signed main(){
	cin>>n>>m;
	for(int i=1;i<=m;i++){
		int x,y;
		cin>>x>>y;
		knight[x][y]=1;
	}
	s=0,t=n*n+1;
	
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if((i+j)%2&&!knight[i][j]){
				add(s,mk(i,j),1);
			}
			if(!((i+j)%2)&&!knight[i][j]){
				add(mk(i,j),t,1);
			}
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if((i+j)%2==0){
				continue;
			}
			for(int k=0;k<8;k++){
				int tx=i+delta[k][0];
				int ty=j+delta[k][1];
				if(valid(tx)&&valid(ty)){
					if(!knight[tx][ty]){
						add(mk(i,j),mk(tx,ty),INT_MAX);
					}
				}
			}
		}
	}
	n=n*n;
	cout<<n-m-dinic()<<'\n';
	return 0;
}

P4304 [TJOI2013]攻击装置

#include <bits/stdc++.h>
#define int long long
using namespace std;

int n,m,s,t;

struct edge{
	int from,to,val;
}e[200001];int head[200001],cur[200001],siz=1;
void addedge(int x,int y,int z){
	e[++siz].to=y,e[siz].val=z;
	e[siz].from=head[x],head[x]=siz;
}
void add(int x,int y,int z){
	addedge(x,y,z);
	addedge(y,x,0);
}

int gap[200001];
bool bfs(){
	memset(gap,0,sizeof(gap));
	fill(gap+1,gap+1+n,0);
	queue<int> q;
	q.push(s);
	gap[s]=1;
	while(!q.empty()){
		int now=q.front();
		q.pop();
		for(int i=head[now];i;i=e[i].from){
			int u=e[i].to;
			if(e[i].val&&!gap[u]){
				gap[u]=gap[now]+1;
				q.push(u);
			}
		}
	}
	return (gap[t]);
}
int dfs(int now,int val){
	if(now==t) return val;
	for(int &i=cur[now];i;i=e[i].from){
		int u=e[i].to;
		if(e[i].val&&gap[now]+1==gap[u]){
			int F=dfs(u,min(e[i].val,val));
			if(F){
				e[i].val-=F;
				e[i^1].val+=F;
				return F;
			}
		}
	}
	return 0;
}
int dinic(){
	int ret=0;
	while(bfs()){
		copy(head,head+1+n,cur);
		int F=0;
		while(F=dfs(s,10000000000000)){
			ret+=F;
		}
	}
	return ret;
}

int knight[10000][1000];

int mk(int x,int y){
	x--;
	return (x*n+y);	
}

int delta[][2]={{1,2},{1,-2},{-1,2},{-1,-2},{2,1},{2,-1},{-2,1},{-2,-1}};

bool valid(int x){
	return (x>=1)&&(x<=n);
}

int gc(){
	char ch;
	do{
		ch=getchar();	
	}while(!isdigit(ch));
	return ch-'0';
}

signed main(){
	cin>>n;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			knight[i][j]=gc();
			m+=knight[i][j];
		}
	}
	s=0,t=n*n+1;
	
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if((i+j)%2&&!knight[i][j]){
				add(s,mk(i,j),1);
			}
			if(!((i+j)%2)&&!knight[i][j]){
				add(mk(i,j),t,1);
			}
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if((i+j)%2==0){
				continue;
			}
			for(int k=0;k<8;k++){
				int tx=i+delta[k][0];
				int ty=j+delta[k][1];
				if(valid(tx)&&valid(ty)){
					if(!knight[tx][ty]){
						add(mk(i,j),mk(tx,ty),INT_MAX);
					}
				}
			}
		}
	}
	n=n*n;
	cout<<n-m-dinic()<<'\n';
	return 0;
}

习题：P5030 长脖子鹿放置

题面

众周所知，在西洋棋中，我们有城堡、骑士、皇后、主教和长脖子鹿。

如图所示，西洋棋的“长脖子鹿”，类似于中国象棋的马，但按照“目”字攻击，且没有中国象棋“别马腿”的规则。（因为长脖子鹿没有马腿）

image

给定一个\(N * M\),的棋盘，有 \(k\) 个格子禁止放棋子。问棋盘上最多能放多少个不能互相攻击的长脖子鹿。

对于\(100\)%的数据，\(1 ≤ N,M ≤ 200\)

代码

#include <bits/stdc++.h>
#define int long long
using namespace std;

int n,m,s,t;

struct edge{
	int from,to,val;
}e[2000001];int head[2000001],cur[2000001],siz=1;
void addedge(int x,int y,int z){
	e[++siz].to=y,e[siz].val=z;
	e[siz].from=head[x],head[x]=siz;
}
void add(int x,int y,int z){
	addedge(x,y,z);
	addedge(y,x,0);
}

int gap[2000001];
bool bfs(){
	memset(gap,0,sizeof(gap));
	fill(gap+1,gap+1+n,0);
	queue<int> q;
	q.push(s);
	gap[s]=1;
	while(!q.empty()){
		int now=q.front();
		q.pop();
		for(int i=head[now];i;i=e[i].from){
			int u=e[i].to;
			if(e[i].val&&!gap[u]){
				gap[u]=gap[now]+1;
				q.push(u);
			}
		}
	}
	return (gap[t]);
}
int dfs(int now,int val){
	if(now==t) return val;
	for(int &i=cur[now];i;i=e[i].from){
		int u=e[i].to;
		if(e[i].val&&gap[now]+1==gap[u]){
			int F=dfs(u,min(e[i].val,val));
			if(F){
				e[i].val-=F;
				e[i^1].val+=F;
				return F;
			}
		}
	}
	return 0;
}
int dinic(){
	int ret=0;
	while(bfs()){
		copy(head,head+1+n,cur);
		int F=0;
		while(F=dfs(s,10000000000000)){
			ret+=F;
		}
	}
	return ret;
}

int knight[100000][1000];

int mk(int x,int y){
	x--;
	return (x*m+y);	
}

int delta[][2]={{1,3},{1,-3},{-1,3},{-1,-3},{3,1},{3,-1},{-3,1},{-3,-1}};

bool valid(int x){
	return (x>=1)&&(x<=n);
}
bool validy(int x){
	return (x>=1)&&(x<=m);
}

int kkk,k;

signed main(){
	cin>>n>>m>>kkk;
	for(int i=1;i<=kkk;i++){
		int x,y;
		cin>>x>>y;
		if(!knight[x][y]){
			k++;
		}
		knight[x][y]=1;
	}
	s=0,t=n*m+1;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			if((i)%2){
				add(s,mk(i,j),1);
			}
			if(!((i)%2)){
				add(mk(i,j),t,1);
			}
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			if(knight[i][j]){
				continue;
			}
			if((i)%2==0){
				continue;
			}
			for(int k=0;k<8;k++){
				int tx=i+delta[k][0];
				int ty=j+delta[k][1];
				if(valid(tx)&&validy(ty)){
					if(!knight[tx][ty]){
						add(mk(i,j),mk(tx,ty),1);
					}
				}
			}
		}
	}
	n=n*m;
	cout<<n-k-dinic()<<'\n';
	return 0;
}

标签：P3355,head,return,val,int,gap,TJOI2013,P4304,now
来源： https://www.cnblogs.com/zheyuanxie/p/network-flow-chessboard-model.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

关于我们 | 联系我们 | 留言反馈

专注分享技术，共同学习，共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有