【bzoj4443】[Scoi2015]小凸玩矩阵

第K大也就是第n-K+1小，所以就可以的二分答案了 (江哥讲过一道类似题)

二分答案找出比当前答案小的数的位置的坐标，判断一下是否可以选出满足不在同一行同一列的n-K+1个数，然后就可以跑匈牙利了，对于一个坐标(x,y)如果满足a[x][y]≤a[x][y]当前答案，就把第x行向第y列连边，然后跑匈牙利判断最大匹配是否大于n-K+1

 #include<algorithm>
 #include<iostream>
 #include<cstdlib>
 #include<cstring>
 #include<cstdio>
 #include<cmath>
 #include<queue>
 #include<set>
 using namespace std;

 #define M 300
 #define N 1000001

 struct edge
 {
     int v,next;
 }e[N<<];
 int head[N<<];
 int cnt;

 int map[M][M];
 int ly[N],f[N];

 int n,m,k;
 int l,r;
 int maxn,tot,ans;

 void link(int x,int y)
 {
     e[++cnt]=(edge){y,head[x]};
     head[x]=cnt;
 }

 bool find(int d)
 {
     for (int i=head[d];i;i=e[i].next)
     {
         int t=e[i].v;
         if (ly[t]!=tot)
         {
             ly[t]=tot;
             if (!f[t] || find(f[t]))
             {
                 f[t]=d;
                 return true;
             }
         }
     }
     return false;
 }

 int work(int x)
 {
     ans=cnt=;
     memset(head,,sizeof(head));
     memset(f,,sizeof(f));
     for (int i=;i<=n;i++)
         for (int j=;j<=m;j++)
             if (map[i][j]<=x)
                 link(i,j);
     for (int i=;i<=n;i++)
     {
         tot++;
         ans+=find(i);
     }
     return ans>=k ? : ;
 }

 int main()
 {
     scanf("%d%d%d",&n,&m,&k);
     k=n-k+;
     for (int i=;i<=n;i++)
         for (int j=;j<=m;j++)
             scanf("%d",&map[i][j]),maxn=max(maxn,map[i][j]);
     l=;
     r=maxn;
     while (l<r)
     {
         int m=l+r>>;
         if (work(m))
             r=m;
         else
             l=m+;
     }
     printf("%d\n",l);
     return ;
 }