[luoguP3355] 骑士共存问题(二分图最大独立集)
模型
二分图最大独立集,转化为二分图最大匹配,从而用最大流解决。
实现
首先把棋盘黑白染色,使相邻格子颜色不同。
把所有可用的黑色格子看做二分图X集合中顶点,可用的白色格子看做Y集合顶点。
建立附加源S汇T,从S向X集合中每个顶点连接一条容量为1的有向边,从Y集合中每个顶点向T连接一条容量为1的有向边。
从每个可用的黑色格子向骑士一步能攻击到的可用的白色格子连接一条容量为无穷大的有向边。
求出网络最大流,要求的结果就是可用格子的数量减去最大流量。
分析
用网络流的方法解决棋盘上的问题,一般都要对棋盘黑白染色,使之成为一个二分图。放尽可能多的不能互相攻击的骑士,就是一个二分图最大独立集问题。有关二分图最大独立集问题,更多讨论见《最小割模型在信息学竞赛中的应用》作者胡伯涛。
该题规模比较大,需要用效率较高的网络最大流算法解决。(使用Dinic+当前弧优化)
——代码
- #include <queue>
- #include <cstdio>
- #include <cstring>
- #include <iostream>
- #define INF 1e9
- #define N 1000001
- #define min(x, y) ((x) < (y) ? (x) : (y))
- int n, m, cnt, tot, sum, s, t;
- int map[][], head[N], to[N], val[N], next[N], dis[N], cur[N];
- int dx[] = {-, -, , , , , -, -},
- dy[] = {, , , , -, -, -, -};
- inline int read()
- {
- int x = , f = ;
- char ch = getchar();
- for(; !isdigit(ch); ch = getchar()) if(ch == '-') f = -;
- for(; isdigit(ch); ch = getchar()) x = (x << ) + (x << ) + ch - '';
- return x * f;
- }
- inline void add2(int x, int y, int z)
- {
- to[cnt] = y;
- val[cnt] = z;
- next[cnt] = head[x];
- head[x] = cnt++;
- }
- inline void add(int x, int y, int z)
- {
- add2(x, y, z);
- add2(y, x, );
- }
- inline bool bfs()
- {
- int i, u, v;
- std::queue <int> q;
- memset(dis, -, sizeof(dis));
- q.push(s);
- dis[s] = ;
- while(!q.empty())
- {
- u = q.front(), q.pop();
- for(i = head[u]; i ^ -; i = next[i])
- {
- v = to[i];
- if(val[i] && dis[v] == -)
- {
- dis[v] = dis[u] + ;
- if(v == t) return ;
- q.push(v);
- }
- }
- }
- return ;
- }
- inline int dfs(int u, int maxflow)
- {
- if(u == t) return maxflow;
- int v, d, ret = ;
- for(int &i = cur[u]; i ^ -; i = next[i])
- {
- v = to[i];
- if(val[i] && dis[v] == dis[u] + )
- {
- d = dfs(v, min(val[i], maxflow - ret));
- ret += d;
- val[i] -= d;
- val[i ^ ] += d;
- if(ret == maxflow) return ret;
- }
- }
- return ret;
- }
- int main()
- {
- int i, j, k, x, y;
- n = read();
- m = read();
- s = , t = n * n + ;
- memset(head, -, sizeof(head));
- for(i = ; i <= m; i++)
- {
- x = read();
- y = read();
- map[x][y] = -;
- }
- for(i = ; i <= n; i++)
- for(j = ; j <= n; j++)
- if(!map[i][j])
- {
- map[i][j] = ++tot;
- if((i + j) & ) add(s, tot, );
- else add(tot, t, );
- }
- for(i = ; i <= n; i++)
- for(j = ; j <= n; j++)
- if(map[i][j] ^ - && (i + j) & )
- for(k = ; k < ; k++)
- {
- x = i + dx[k];
- y = j + dy[k];
- if(x >= && x <= n && y >= && y <= n && map[x][y] ^ -) add(map[i][j], map[x][y], INF);
- }
- while(bfs())
- {
- for(i = s; i <= t; i++) cur[i] = head[i];
- sum += dfs(s, INF);
- }
- printf("%d\n", tot - sum);
- return ;
- }
[luoguP3355] 骑士共存问题(二分图最大独立集)的更多相关文章
- 洛谷P3355 骑士共存问题 二分图_网络流
Code: #include<cstdio> #include<cstring> #include<queue> #include<vector> #i ...
- 【Codevs1922】骑士共存问题(最小割,二分图最大独立集转最大匹配)
题意: 在一个n*n个方格的国际象棋棋盘上,马(骑士)可以攻击的棋盘方格如图所示.棋盘上某些方格设置了障碍,骑士不得进入. 对于给定的n*n个方格的国际象棋棋盘和障碍标志,计算棋盘上最多可以放置多少个 ...
- 「CODVES 1922 」骑士共存问题(二分图的最大独立集|网络流)&dinic
首先是题目链接 http://codevs.cn/problem/1922/ 结果发现题目没图(心情复杂 然后去网上扒了一张图 大概就是这样了. 如果把每个点和它可以攻击的点连一条边,那问题就变成了 ...
- Cogs 746. [网络流24题] 骑士共存(最大独立集)
[网络流24题] 骑士共存 ★★☆ 输入文件:knight.in 输出文件:knight.out 简单对比 时间限制:1 s 内存限制:128 MB 骑士共存问题 «问题描述: 在一个n*n个方格的国 ...
- COGS746. [网络流24题] 骑士共存
骑士共存问题«问题描述:在一个n*n个方格的国际象棋棋盘上,马(骑士)可以攻击的棋盘方格如图所示.棋盘 上某些方格设置了障碍,骑士不得进入. «编程任务:对于给定的n*n个方格的国际象棋棋盘和障碍标志 ...
- 【刷题】LOJ 6226 「网络流 24 题」骑士共存问题
题目描述 在一个 \(\text{n} \times \text{n}\) 个方格的国际象棋棋盘上,马(骑士)可以攻击的棋盘方格如图所示.棋盘上某些方格设置了障碍,骑士不得进入. 对于给定的 \(\t ...
- P3355 骑士共存问题
P3355 骑士共存问题 题目描述 在一个 n*n (n <= 200)个方格的国际象棋棋盘上,马(骑士)可以攻击的棋盘方格如图所示.棋盘上某些方格设置了障碍,骑士不得进入 对于给定的 n*n ...
- BZOJ3175:[TJOI2013]攻击装置(二分图最大独立集)
Description 给定一个01矩阵,其中你可以在0的位置放置攻击装置.每一个攻击装置(x,y)都可以按照“日”字攻击其周围的 8个位置(x-1,y-2),(x-2,y-1),(x+1,y-2), ...
- [网络流24题] 骑士共存(cogs 746)
骑士共存问题«问题描述:在一个n*n个方格的国际象棋棋盘上,马(骑士)可以攻击的棋盘方格如图所示.棋盘 上某些方格设置了障碍,骑士不得进入. «编程任务:对于给定的n*n个方格的国际象棋棋盘和障碍标志 ...
随机推荐
- [dp uestc oj] G - 邱老师玩游戏
G - 邱老师玩游戏 Time Limit: 3000/1000MS (Java/Others) Memory Limit: 65535/65535KB (Java/Others) Submi ...
- codeforces Gym 100286J Javanese Cryptoanalysis (二染色)
每一单词相邻两个字母,不能同时为元音或者辅音... 各种姿势都可以过:7个for,dp,黑白染色,dfs,并查集.... 最主要的思路就是相邻字母连边,把元音和辅音看成两个集合,那么有连边的两个字母一 ...
- 2018.4.28 基于java的聊天系统(带完善)
Java聊天系统 1.Socket类 Socket(InetAddress address, int port) 创建一个流套接字并将其连接到指定 IP 地址的指定端口号. Socket(String ...
- centos7 python3 Saltstack配置
Python安装完毕后,提示找不到ssl模块 pip is configured with locations that require TLS/SSL, however the ssl module ...
- java基础—static关键字
一.static关键字
- 服务器上搭建flowvisor平台
之前全是在virtualbox上的Ubuntu虚拟机上测试的ovs以及pox, 现在我们开始在服务器上开始了 两台服务器上的ovs均是1.4.6版本 遇到一个问题:之前装的ovs down了 然后什么 ...
- ios sinaweibo 客户端(三)
这个页面要讲述的是用户的粉丝列表,下面是效果图: 可以看到这个视图明显也是一个tableview,在每一个cell中包含的有三个部分的内容:粉丝头像image,粉丝昵称label,我和粉丝之间的相互关 ...
- poj3335 Rotating Scoreboard
题目描述: vjudge POJ 题解: 半平面交判核的存在性. 重点在于一个点的核也算核. 这样的话普通的求多边形的版本就要加一个特判. 就是把剩下的一个节点暴力带回所有直线重判,这时判叉积是否$\ ...
- 201621123080《JAVA程序设计》第八周学习总结
作业08-集合 1. 本周学习总结 2. 书面作业 1. ArrayList代码分析 1.1 解释ArrayList的contains源代码 根据代码,首先在源数组里找到下标,若下标符合>=0 ...
- LeetCode935
问题:935. 骑士拨号器 国际象棋中的骑士可以按下图所示进行移动: . 这一次,我们将 “骑士” 放在电话拨号盘的任意数字键(如上图所示)上,接下来,骑士将会跳 N-1 步. ...