【二分图匹配 练习】洛谷P3386【模板】二分图最大匹配P2756飞行员配对方案问题P1129矩阵游戏P1559运动员最佳匹配问题P4014分配问题P1640连续攻击游戏

于 2021-04-16 21:51:04 发布
阅读量123 收藏 1
点赞数
分类专栏： 【算法】树，图
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/wldcmzy/article/details/115772591
版权
这些代码实现了一系列图论问题的解决方案，包括二分图最大匹配、飞行员配对、矩阵游戏和运动员最佳匹配。它们使用了深度优先搜索（DFS）策略，并通过调整权重来优化匹配过程。在洛谷P1129矩阵游戏中，检查是否存在完全匹配；在洛谷P1559运动员最佳匹配问题中，展示了两种不同方法的效率比较；洛谷P4014分配问题中，通过调整边的权重寻找最优分配。代码中还展示了如何避免因未初始化变量导致的错误。
摘要由CSDN通过智能技术生成

洛谷 P3386【模板】二分图最大匹配

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <map>
#include <stack>
#include <string>
#include <cstring>
#include <vector>
#include <cmath>
#include <set>
typedef long long int LL;
const int N = 505;
int mp[N][N],link[N],vis[N];
int n,m,e,u,v,ans;
bool dfs(int x)
{
    for(int i=1; i<=m; i++){
        if(!vis[i] && mp[x][i]){
            vis[i] = 1;
            if(!link[i] || dfs(link[i])){
                link[i] = x;
                return true;
            }
        }
    }
    return false;
}
int main()
{
    scanf("%d%d%d",&n,&m,&e);
    for(int i=0; i<e; i++){
        scanf("%d%d",&u,&v);
        mp[u][v] = 1;
    }
    for(int i=1; i<=n; i++) {
        memset(vis,0,sizeof(vis));
        ans += dfs(i);
    }
    printf("%d",ans);
    return 0;
}

洛谷 P2756飞行员配对方案问题

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <map>
#include <stack>
#include <string>
#include <cstring>
#include <vector>
#include <cmath>
#include <set>
typedef long long int LL;
const int N = 105;
int mp[N][N],vis[N],link[N];
int n,m,u,v,ans;
bool dfs(int x)
{
    for(int i=1; i<=n; i++){
        if(!vis[i] && mp[x][i]){
            vis[i] = 1;
            if(!link[i] || dfs(link[i])){
                link[i] = x;
                return true;
            }
        }
    }
    return false;
}
int main()
{
    scanf("%d%d",&m,&n);
    n-=m;
    while(scanf("%d%d",&u,&v)){
        if(u == -1 && v == -1) break;
        mp[u][v-m] = 1;
    }
    for(int i=1; i<=m; i++) {
        memset(vis,0,sizeof(vis));
        ans += dfs(i);
    }
    printf("%d\n",ans);
    for(int i=1; i<=n; i++){
        if(link[i]) printf("%d %d\n",link[i],i+m);
    }
    return 0;
}

洛谷 P1129 [ZJOI2007] 矩阵游戏

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <map>
#include <stack>
#include <string>
#include <cstring>
#include <vector>
#include <cmath>
#include <set>
typedef long long int LL;
const int N = 205;
int mp[N][N],vis[N],link[N];
int t,n,ans;
bool dfs(int x)
{
    for(int i=1; i<=n; i++){
        if(!vis[i] && mp[x][i]){
            vis[i] = 1;
            if(!link[i] || dfs(link[i])){
                link[i] = x;
                return true;
            }
        }
    }
    return false;
}
int main()
{
    scanf("%d",&t);
    while(t--){
        memset(link,0,sizeof(link));
        ans = 0;
        scanf("%d",&n);
        for(int i=1; i<=n; i++){
            for(int j=1; j<=n; j++){
                scanf("%d",&mp[i][j]);
            }
        }
        for(int i=1; i<=n; i++){
            memset(vis,0,sizeof(vis));
            ans += dfs(i);
        }
        printf(ans == n ? "Yes\n" : "No\n");
    }
    return 0;
}

洛谷P1559 运动员最佳匹配问题

两种敲代码方法，一种dfs完算tmp，一种在dfs中算tmp。不确定是否一种比另一种更优，这道题后者时间占优，然而下一道题洛谷 P4014 分配问题前者时间占优。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <map>
#include <stack>
#include <string>
#include <cstring>
#include <vector>
#include <cmath>
#include <set>
typedef long long int LL;
const int N = 205,INF = 0x7fffffff;
int mp[N][N],lx[N],ly[N],visx[N],visy[N],link[N];
int n,tmp,v;
bool dfs(int x)
{
    visx[x] = 1;
    for(int i=1; i<=n; i++){
        if(!visy[i] && lx[x]+ly[i]==mp[x][i]){
            visy[i] = 1;
            if(!link[i] || dfs(link[i])){
                link[i] = x;
                return true;
            }
        }
    }
    return false;
}
int one()
{
    int ans = 0;
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++){
        lx[i] = std::max(lx[i],mp[i][j]);
    }
    for(int i=1; i<=n; i++){
        while(true){
            memset(visx,0,sizeof(visx));
            memset(visy,0,sizeof(visy));
            tmp = INF;
            if(dfs(i)) break;
            for(int i=1; i<=n; i++) if(visx[i]){
                for(int j=1; j<=n; j++) if(!visy[j]){
                    tmp = std::min(tmp,lx[i]+ly[j]-mp[i][j]);
                }
            }
            for(int j=1; j<=n; j++){
                if(visx[j]) lx[j] -= tmp;
                if(visy[j]) ly[j] += tmp;
            }
        }
    }
    for(int i=1; i<=n; i++){
        ans += ly[i] + lx[link[i]];
        //std::cout << "ly[" << i << "] = " << ly[i] << " lx[" << link[i] << "] = " << lx[link[i]] << " ans = " << ans << std::endl;
    }
    return ans;
}
int main()
{
    //freopen("D:\\EdgeDownloadPlace\\P4014_1.in","r",stdin);
    scanf("%d",&n);
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++){
        scanf("%d",&v);
        mp[i][j] = -v;
    }
    for(int i=0; i<N; i++) lx[i] = -INF;
    printf("%d\n",-one());
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++){
        mp[i][j] = -mp[i][j];
    }
    memset(lx,0,sizeof(lx));
    memset(ly,0,sizeof(ly));
    memset(link,0,sizeof(link));
    printf("%d\n",one());
    return 0;
}

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <map>
#include <stack>
#include <string>
#include <cstring>
#include <vector>
#include <cmath>
#include <set>
typedef long long int LL;
const int N = 23,INF = 0x7fffffff;
int mp[N][N],vm[N],vf[N],vism[N],visf[N],link[N];
int n,v,ans,tmp;
bool dfs(int x)
{
    vism[x] = 1;
    for(int i=1; i<=n; i++){
        if(!visf[i]){
            int t = vm[x]+vf[i]-mp[x][i];
            if(t==0){
                visf[i] = 1;
                if(!link[i] || dfs(link[i])){
                    link[i] = x;
                    return true;
                }

            }
            else if(t>0){
                tmp = std::min(t,tmp);
            }
        }
    }
    return false;
}
int main()
{
    scanf("%d",&n);
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++) {
        scanf("%d",&mp[i][j]);
    }
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++) {
        scanf("%d",&v);
        mp[j][i] *= v;
    }
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++) {
        vm[i] = std::max(vm[i],mp[i][j]);
    }
    for(int i=1; i<=n; i++){
        while(true){
            memset(visf,0,sizeof(visf));
            memset(vism,0,sizeof(vism));
            tmp = INF;
            if(dfs(i))break;
            for(int i=1; i<=n; i++){
                if(visf[i]) vf[i] += tmp;
                if(vism[i]) vm[i] -= tmp;
            }
        }
    }
    for(int i=1; i<=n; i++){
        ans += vf[i] + vm[link[i]];
    }
    printf("%d",ans);
    return 0;
}

洛谷 P4014 分配问题

全WA，我的输出莫名其妙出现了"-"，洛谷判题机的锅？
破案了，one()里的ans没初始化，Windows自动给个0，Linux就挂了，初始化成0就好了

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <map>
#include <stack>
#include <string>
#include <cstring>
#include <vector>
#include <cmath>
#include <set>	
typedef long long int LL;
const int N = 105,INF = 0x7fffffff;
int mp[N][N],lx[N],ly[N],visx[N],visy[N],link[N];
int n,tmp,v;
bool dfs(int x)
{
    visx[x] = 1;
    for(int i=1; i<=n; i++){
        if(!visy[i] && lx[x]+ly[i]==mp[x][i]){
            visy[i] = 1;
            if(!link[i] || dfs(link[i])){
                link[i] = x;
                return true;
            }
        }
    }
    return false;
}
int one()
{
    int ans;
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++){
        lx[i] = std::max(lx[i],mp[i][j]);
    }
    for(int i=1; i<=n; i++){
        while(true){
            memset(visx,0,sizeof(visx));
            memset(visy,0,sizeof(visy));
            tmp = INF;
            if(dfs(i)) break;
            for(int i=1; i<=n; i++){
                if(visx[i]){
                    for(int j=1; j<=n; j++){
                        if(!visy[j]){
                            tmp = std::min(tmp,lx[i]+ly[j]-mp[i][j]);
                        }
                    }
                }
            }
            for(int j=1; j<=n; j++){
                if(visx[j]) lx[j] -= tmp;
                if(visy[j]) ly[j] += tmp;
            }
        }
    }
    for(int i=1; i<=n; i++){
        ans += ly[i] + lx[link[i]];
    }
    return ans;
}
int main()
{
    //freopen("D:\\EdgeDownloadPlace\\P4014_1.in","r",stdin);
    scanf("%d",&n);
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++){
        scanf("%d",&v);
        mp[i][j] = -v;
    }
    for(int i=0; i<N; i++) lx[i] = -INF;
    printf("%d\n",-one());
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++){
        mp[i][j] = -mp[i][j];
    }
    memset(lx,0,sizeof(lx));
    memset(ly,0,sizeof(ly));
    memset(link,0,sizeof(link));
    printf("%d\n",one());
    return 0;
}

洛谷 P1640连续攻击游戏

匈牙利优化 每次memset置零vis数组太慢了，TLE，直接开一个变量timmer当作时间计数，当timmer-vis[x]>0就说明在这次dfs中x节点没有访问过。
.
初始timmer为1， 每次dfs后timmer自增1

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <string>
#include <queue>
#include <vector>
#include <stack>
#include <map>
#include <set>
typedef long long int LL;
const int INF = 0x7fffffff, SCF = 0x3fffffff;
const int N = 1e6+5, M = 220, K = 1e6+4;
int head[N],link[N],vis[N];
int n,idx,a,b,timmer;
struct Edge
{
    int to,nxt;
}edg[N<<1];
void addE(int fr,int to)
{
    edg[idx].to = to;
    edg[idx].nxt = head[fr];
    head[fr] = idx++;
}
bool dfs(int x)
{
    for(int e=head[x]; e; e=edg[e].nxt){
        int y = edg[e].to;
        if(timmer-vis[y]>0){
            vis[y] = timmer;
            if(!link[y] || dfs(link[y])){
                link[y] = x;
                return true;
            }
        }
    }
    return false;
}
int main()
{
    idx = 2;
    scanf("%d",&n);
    for(int i=1; i<=n; i++){
        scanf("%d%d",&a,&b);
        addE(a,i);
        addE(b,i);
    }
    int i;
    for(i=1; i<=n; i++){
        timmer = i;
        if(!dfs(i)) break;
    }
    printf("%d",i-1);
    return 0;
}
StarsWhisper
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
广告商收益化最大问题_M?n_KM算法_二分图最大匹配_
10-01
n-KM算法与二分图最大匹配是解决广告商收益化最大问题的关键工具。通过对广告商和广告位进行高效的匹配，可以实现双方的共赢，提高整个系统的效益。在处理实际问题时，还需要结合其他策略如预算约束、点击预测等，以...
ZJPP.zip_graph matching_二分图_二分图 matlab_最佳匹配_最大权匹配
07-14
在本讨论中，我们将深入探讨如何利用MATLAB来解决二分图的最大权匹配问题。 最大权匹配，也被称为最佳匹配，是指在二分图中找到一个匹配，使得所有匹配的边的权重之和最大。这个问题在很多实际问题中都有应用，如...
二分图最优匹配matlab代码.zip
07-11
二分图最优匹配是图论中的一个重要概念，它在许多实际问题中有着广泛的应用，比如分配问题、网络流问题、任务调度等。在本压缩包中，提供的是一组使用MATLAB语言实现的二分图最优匹配算法代码。MATLAB作为一种强大的...
xqq.rar_matlab二分图_二分 最大 匹配 matlab_二分图_匹配 二分_最大完美匹配
09-24
例如，在匹配学生与导师的问题中，可以构建一个二分图，学生和导师分别位于两个集合，然后寻找最大完美匹配，以确定最佳的配对方案。 总的来说，理解和掌握二分图最大完美匹配的理论与实现，对于IT从业者尤其是算法...
Adafruit Feather 328P资料 包含原理图/pcb/用户手册-电路方案
04-21
有了这款Feather，我们对ATmega328P（经典的'Arduino'芯片）有了一点怀旧之情 - 这款Adafruit Feather 328P运行3.3V和8 MHz。 羽毛是Adafruit的新开发板，就像它的名字一样，它轻薄，让你飞翔！ 我们将Feather设计为...
大数据处理框架：Flink：FlinkSQL入门与实践.docx
09-02
大数据处理框架：Flink：FlinkSQL入门与实践.docx
数据集成工具：Talend：Talend数据集成最佳实践与案例分析.docx
09-02
数据集成工具：Talend：Talend数据集成最佳实践与案例分析.docx
实时计算：Azure Stream Analytics：数据流窗口与聚合操作.docx
09-02
实时计算：Azure Stream Analytics：数据流窗口与聚合操作.docx
利于python爬虫爬取淘宝司法拍卖和京东司法拍卖土地的每日信息.zip
09-02
利于python爬虫爬取淘宝司法拍卖和京东司法拍卖土地的每日信息
51单片机心型流水灯 多种样式
09-02
51单片机 多种花样流水灯
#-ssm-024-mysql-物流管理系统vue-(1实用1).zip
09-02
摘 要 企业的发展离不开物流的运输，在一个大型的企业中，商品的生产和建设，推广只是前期的一些工作，在后期的商品销售和物流方面的建立，才能让一个企业得到大力的发展。 企业在物流方面的竞争随之日益激烈，企业开始逐渐的认识到物流管理信息化的重要性。对于企业来说，加强物流信息化管理，不仅有助于提高企业员工的工作效率，同时还助于降低企业物流成本，并且还有助于提高企业物流管理的质量。因此，针对物流管理系统的问题及措施展开探究十分的有必要。因此，开发这款物流管理系统是至关重要的。 我们根据需求调研分析后得出结论，为了使用方便，该平台必须采用BS架构的模式进行开发，利用java编程语言和现今成熟的SSM框架进行整合开发。系统数据库采用方便安装的mysql进行数据管理。系统前端主要采用VUE框架技术设计开发，该框架在网页前端设计具有很强大的功能和成熟的设计方式。该系统是一款多用户的后台管理系统，主要包括物流的一系列流程，从采购管理，供应商信息管理以及配送运输，入库出库等进行综合管理。最终达到系统的应用目标。 关键词：物流信息化；配送管理；SSM框架平台技术；mysql小型数据库
1024_+_深度强化学习（Deep_Reinforcement_Learning_+_1024_G_DI-1024.zip
09-02
1024_+_深度强化学习（Deep_Reinforcement_Learning_+_1024_G_DI-1024
Av1an-5.6.post1-py3-none-any.whl.zip
09-02
Av1an-5.6.post1-py3-none-any.whl.zip
CentOS-Base.repo
09-02
CentOS-Base.repo
消息队列：RabbitMQ：RabbitMQ高级特性：死信队列技术教程.docx
09-02
消息队列：RabbitMQ：RabbitMQ高级特性：死信队列技术教程.docx
数据集成工具：Azure数据工厂：控制流活动：构建复杂工作流.docx
09-02
数据集成工具：Azure数据工厂：控制流活动：构建复杂工作流.docx
C语言小作业：含有附加功能的倒计时器
最新发布
09-02
可输入1-999的数字并开始计时 可选择三种图形显示数字的方式：“ . ” “ * ” “ |-/ ” 整十秒、倒数最后五秒和倒计时结束时会有音效和部分提示信息 可以结束后选择是否再次计时，也可以退出 中途想中断可直接关闭命令行窗口，想再次计时可以直接再次运行 输入错误指令后会显示错误，重新输入          -倒计时结束可以设置彩蛋 -如果修改部分代码后图形错版，请对图形函数进行检查 -环境应该是visual studio 2012，vs2015后gets()无法使用。 -如果由于硬件原因，不同电脑运行时候倒计时速度有小误差，可在Sleep(978)处修改值，越大越慢，越小越快 -如果音效不奏效，请检查电脑设置或音量
外呼神器V1.0 免费下载
09-02
外呼神器V1.0 免费下载
Av1an-1.13.post6-py3-none-any.whl.zip
09-02
Av1an-1.13.post6-py3-none-any.whl.zip
二分图匹配问题与匈牙利算法解析
"二分图匹配的理论及应用实例" 二分图匹配是图论中的一个...二分图匹配是一个强大的工具，广泛应用于各种优化问题，如任务分配、网络路由、配对问题等。理解并掌握二分图匹配的理论和算法对于解决实际问题至关重要。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值