BZOJ4443: [Scoi2015]小凸玩矩阵 解题报告

日常吐槽，可跳过这部分

这道题第一眼就觉得要二分一下答案，然后…然后我发现我不会……..冷静下来，想了一会，感觉可以二分后可以贪心一下，然后就欢快的打了起来

打完过了样例，然后WA了（不对呀我这么正确的 正解怎么会WA呢）

打了个对拍死活拍不出来，好吧我弱我看题解唉呀我怎么没想到二分匹配…….（但我真的觉得我原来的做法没错的）

没救了

吐槽结束下面是题解

二分一下第k大，然后二分匹配一下，不大于二分的第k大的都将这列和这行连边，

如果最大匹配数≥n−k+1就可行

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<ctime>
#include<cmath>
#include<string>
#include<vector>
#include<string>
#include<cstdio>
#include<cstring>
#include<climits>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 270;
int mat[maxn][maxn],d[maxn*maxn];
int match[maxn];
bool v[maxn];
int n,m,k,num;

struct edge
{
int y,next;
}a[maxn*maxn]; int len,first[maxn];
void insert(int x,int y)
{
len++;
a[len].y=y;
a[len].next=first[x];first[x]=len;
}
bool find_( int x )
{
if( v[x] ) return false;
v[x] = true;
for( int k=first[x];k;k=a[k].next )
{
int y=a[k].y;
if(!match[y]||find_(match[y]))
{
match[y]=x;
return true;
}
}
return false;
}
bool judge(int now)
{
memset(first,0,sizeof first); len=0;
memset(match,0,sizeof match);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
if(mat[i][j]<=now) insert(i,j);
int ret=0;
for(int i=1;i<=n;i++)
{
memset(v,false,sizeof v);
if( find_(i) ) ret++;
}
return ret>n-k;
}
int solve()
{
int l=1,r=num;
while(l<=r)
{
int mid=(l+r)>>1;
if(judge(d[mid])) r=mid-1;
else l=mid+1;
}
return r+1;
}

int main()
{
//freopen("tmp.in","r",stdin);
//freopen("tmp.out","w",stdout);

scanf("%d%d%d",&n,&m,&k);
num=0;
for( int i=1;i<=n;i++ )
for( int j=1;j<=m;j++ )
{
scanf("%d",&mat[i][j]);
d[++num]=mat[i][j];
}
sort(d+1,d+num+1);
printf("%d\n",d[solve()]);

return 0;
}