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;
}