bzoj4443 [Scoi2015]小凸玩矩阵

37 篇文章 0 订阅
30 篇文章 0 订阅

Description


小凸和小方是好朋友,小方给小凸一个N*M(N<=M)的矩阵A,要求小秃从其中选出N个数,其中任意两个数字不能在同一行或同一列,现小凸想知道选出来的N个数中第K大的数字的最小值是多少。

1<=K<=N<=M<=250,1<=矩阵元素<=10^9

来自 https://www.lydsy.com/JudgeOnline/problem.php?id=4443

Solution


第k大可以是第n-k+1小,那么二分答案把不大于mid的点行列连边,若最大匹配数>n-k+1就可行,否则扩大二分下界

Code


#include <stdio.h>
#include <string.h>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))

const int N=200005;
const int E=800005;
const int L=505;

struct edge {int x,y,next;} e[E];

int rc[L][L];
int ls[N],edCnt;
int vis[N],link[N];

void add_edge(int x,int y) {
    e[++edCnt]=(edge) {x,y,ls[x]}; ls[x]=edCnt;
    e[++edCnt]=(edge) {y,x,ls[y]}; ls[y]=edCnt;
}

bool find(int x,int id) {
    for (int i=ls[x];i;i=e[i].next) {
        if (vis[e[i].y]==id) continue;
        vis[e[i].y]=id;
        if (find(link[e[i].y],id)||!link[e[i].y]) {
            link[e[i].y]=x;
            return true;
        }
    }
    return false;
}

int hungary(int n) {
    fill(link,0);
    fill(vis,0);
    int ret=0;
    rep(i,1,n) ret+=find(i,i);
    return ret;
}

void build_graph(int n,int m,int lim) {
    fill(ls,0); edCnt=0;
    rep(i,1,n) rep(j,1,m) if (rc[i][j]<=lim) add_edge(i,j+n);
}

int main(void) {
    int n,m,k; scanf("%d%d%d",&n,&m,&k); k=n-k+1;
    rep(i,1,n) rep(j,1,m) scanf("%d",&rc[i][j]);
    int ans=0;
    for (int l=0,r=1000000000;l<=r;) {
        int mid=(l+r)>>1;
        build_graph(n,m,mid);
        int ret=hungary(n);
        if (ret>=k) r=mid-1,ans=mid;
        else l=mid+1;
    }
    printf("%d\n", ans);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值