【BZOJ 4443】 [Scoi2015]小凸玩矩阵|二分|最大匹配|匈牙利

最新推荐文章于 2019-02-24 18:53:27 发布

suxb201

最新推荐文章于 2019-02-24 18:53:27 发布

阅读量439

点赞数 1

分类专栏：下午

本文链接：https://blog.csdn.net/sxb_201/article/details/50960133

版权

下午专栏收录该内容

27 篇文章 0 订阅

订阅专栏

没YY出来

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int MAXN  = 251;
int num[MAXN][MAXN];
int n,m,k,T,cnt;
int b[MAXN],go[MAXN][MAXN],flag[MAXN];
struct H
{
	int x,y;
	int num;
}a[MAXN*MAXN];
bool cmp(H a, H b){return a.num<b.num;}
bool dfs(int x)
{
	if(b[x]==T) return false;
	b[x]=T;
	for(int i=1;i<=go[x][0];i++)
	{
		int tmp=go[x][i];
		if(!flag[tmp]||dfs(flag[tmp]))
		{
			flag[tmp]=x;
			return true;
		}
	}
	return false;
}
bool ok()
{
	memset(flag,0,sizeof(flag));
	int ans=0;
	for(int i=1;i<=n;i++)
	{
		T++;
		if(dfs(i))
		{
			ans++;
			if(ans==k)return true;
		}
	}
	return false;
}
int main()
{
	scanf("%d %d %d",&n,&m,&k);
	k=n+1-k;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			scanf("%d",&a[++cnt].num),
			a[cnt].x=i,
			a[cnt].y=j;
	sort(a+1,a+cnt+1,cmp);
	int last=0;
	int L=1,R=cnt;
	while(L!=R)
	{
		int mid=(L+R)/2;
		while(last<mid) {last++;go[a[last].x][++go[a[last].x][0]]=a[last].y;}
		while(last>mid) {go[a[last].x][0]--;last--;}
		if(ok()) R=mid;
		else L=mid+1;
	}
	printf("%d",a[L].num);
	return 0;
}