BZOJ 4443 [Scoi2015]小凸玩矩阵(二分答案+二分图匹配)
2024-08-26 00:52:25
【题目链接】http://www.lydsy.com/JudgeOnline/problem.php?id=4443
【题目大意】
从矩阵中选出N个数,其中任意两个数字不能在同一行或同一列
求选出来的N个数中第K大的数字的最小值是多少。
【题解】
我们二分这个第k大数字的大小,将其以上的数字全部删除,
在剩余的部分按行列连边,如果二分图匹配的数量大于n-k那么说明该答案可行。
【代码】
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>
using namespace std;
const int MAX_V=1000;
const int INF=0x3f3f3f3f;
int V,match[MAX_V];
vector<int> G[MAX_V];
bool used[MAX_V];
void add_edge(int u,int v){
G[u].push_back(v);
G[v].push_back(u);
}
bool dfs(int v){
used[v]=1;
for(int i=0;i<G[v].size();i++){
int u=G[v][i],w=match[u];
if(w<0||!used[w]&&dfs(w)){
match[v]=u;
match[u]=v;
return 1;
}
}return 0;
}
int bipartite_matching(){
int res=0;
memset(match,-1,sizeof(match));
for(int v=0;v<V;v++){
if(match[v]<0){
memset(used,0,sizeof(used));
if(dfs(v))res++;
}
}return res;
}
const int N=300;
int n,m,k,a[N][N];
bool check(int x){
V=n+m+1;
for(int i=0;i<V;i++)G[i].clear();
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++)if(a[i][j]<=x){
add_edge(i,n+j);
}
}return bipartite_matching()>=n-k+1;
}
int main(){
while(~scanf("%d%d%d",&n,&m,&k)){
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++)scanf("%d",&a[i][j]);
}
int l=1,r=1000000000,ans=r;
while(l<=r){
int mid=(l+r)>>1;
if(check(mid))r=mid-1,ans=mid;
else l=mid+1;
}printf("%d\n",ans);
}return 0;
}
最新文章
- php 快速排序
- json写入和读取代码
- 用 perl 统计 fasta 文件序列的总长
- Google地图接口API之地图类型(六)
- dom4j-full.jar 解析 XML
- DebugViewHierarchy
- Eclipse中通过Android模拟器调用OpenGL ES2.0函数操作步骤
- git 使用过程(二、创建版本库)
- Chapter 1 First Sight——9
- mysql 查询 最大值,最小值,第二大,第三大 一共四个值
- Tensorflow计算加速
- Android 常用RGB值以及中英文名称
- pkg-config 用法
- ViewController之间的转场和数据传递
- python函数定义语法总结
- C++设计模式(第一周)
- APP之红点提醒三个阶段
- umount nfs文件系统 显示 umount.nfs: device is busy
- adb shell 运行时报错";adb server version (26) doesn&#39;t match this client (39); killing...";的解决方案
- Deal with relational data using libFM with blocks