题解 【SCOI2015】小凸玩矩阵

题面

解析

这题其实也是网络流建图。。

首先，转换下思路，

求第k大的数的最小值，

其实就是求一个最小的值，

使选取的点中能有（n-k+1)个的值比它小。

因此，可以采用二分答案，

每次判断一个值，

将比它小的点加到图中跑最大流，

看流量是否大于（n-k+1）。

那么，怎么连边呢？

其实，我们可以每一行连源点，流量为1，

每一列连汇点，流量为1，

中间源点与汇点连INF。

最后判断就能AC了！

上AC代码：

#include<bits/stdc++.h>
using namespace std;

inline int read(){
    int sum=,f=;char ch=getchar();
    while(ch>'' || ch<''){if(ch=='-')f=-;ch=getchar();}
    while(ch>='' && ch<=''){sum=sum*+ch-'';ch=getchar();}
    return f*sum;
}

const int INF=0x3f3f3f3f;
struct node{
    int next,to,v;
}e[];
int n,k,m,ans,s,t;
int a[][];
int maxn=;
int head[],cnt;
int d[],cur[];

void add(int x,int y,int v){
    e[++cnt].to=head[x];e[cnt].next=y;
    e[cnt].v=v;head[x]=cnt;
    e[++cnt].to=head[y];e[cnt].next=x;
    e[cnt].v=;head[y]=cnt;
}

bool bfs(){
    memset(d,,sizeof(d));
    memcpy(cur,head,sizeof(cur));
    queue <int>    que;
    que.push(s);
    d[s]=;
    while(!que.empty()){
        int x=que.front();
        que.pop();
        for(int i=head[x];i;i=e[i].to){
            int k=e[i].next;
            if(!e[i].v||d[k]) continue;
            d[k]=d[x]+;
            que.push(k);
        }
    }
    return d[t];
}

int dfs(int x,int mi){
    if(x==t||!mi) return mi;
    int flow=;
    for(int &i=cur[x];i;i=e[i].to){
        int k=e[i].next;
        if(!e[i].v||d[k]!=d[x]+) continue;
        int ret=dfs(k,min(mi,e[i].v));
        flow+=ret;mi-=ret;
        e[i].v-=ret;e[i^].v+=ret;
        if(!mi) break;
    }
    if(mi) d[x]=-;
    return flow;
}

bool DINIC(int x){
    int ans=;
    while(bfs()){
        ans+=dfs(s,INF);
    }
    if(ans>=n-k+) return ;
    return ;
}

bool check(int x){
    memset(e,,sizeof(e));
    memset(head,,sizeof(head));
    cnt=;
    for(int i=;i<=n;i++) add(s,i,);
    for(int j=;j<=m;j++) add(j+n,t,);
    for(int i=;i<=n;i++){
        for(int j=;j<=m;j++){
            if(a[i][j]<=x) add(i,j+n,INF);
        }
    }
    if(DINIC(x)) return ;
    return ;
}

int main(){
    //freopen("matrix.in","r",stdin);
    //freopen("matrix.out","w",stdout);
    n=read();m=read();k=read();
    s=m+n+;t=s+;
    for(int i=;i<=n;i++){
        for(int j=;j<=m;j++){
            a[i][j]=read();
            maxn=max(maxn,a[i][j]);
        }
    }
    int l=,r=maxn;
    while(l<=r){
        int mid=(l+r)>>;
        if(check(mid)){
            r=mid-;ans=mid;
        }
        else l=mid+;
    }
    printf("%d\n",ans);
    return ;
}