J. Spy（多重匹配KM算法）

最新推荐文章于 2021-02-16 20:25:48 发布
最新推荐文章于 2021-02-16 20:25:48 发布
阅读量743 收藏 1
点赞数 1
分类专栏： 二分图匹配 文章标签： KM 完全匹配
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Fighting_Peter/article/details/112503634
版权

J. Spy

随机打，最后答案乘n，因为我方等概率的遇见敌人，相当于与n个敌人都打一遍，然后贡献累加作为匹配边权

bfs版本的KM
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=510;
ll w[N][N],lx[N],ly[N];//顶标
bool visx[N],visy[N];//记录是否在交错树上
int match[N];//匹配
int n,p[N];
ll delta,c[N];//delta和更新后的delta
void bfs(int x)
{
    int a,y=0,y1=0;
    for(int i=1;i<=n;i++) p[i]=0,c[i]=1e18;
    match[y]=x;
    do{
        a=match[y],delta=1e18,visy[y]=1;
        for(int b=1;b<=n;b++)
            if(!visy[b])
            {
                if(c[b]>lx[a]+ly[b]-w[a][b])
                    c[b]=lx[a]+ly[b]-w[a][b],p[b]=y;
                if(c[b]<delta)//Δ还是取最小的
                    delta=c[b],y1=b;
            }
        for(int b=0;b<=n;b++)
            if(visy[b])
                lx[match[b]]-=delta,ly[b]+=delta;
            else
                c[b]-=delta;
        y=y1;
    }while(match[y]);
    while(y) match[y]=match[p[y]],y=p[y];
}
ll KM()
{
    for(int i=1;i<=n;i++)
        match[i]=lx[i]=ly[i]=0;
    for(int i=1;i<=n;i++)
    {
        memset(visy,0,sizeof visy);
        bfs(i);
    }
    ll res=0;
    for(int i=1;i<=n;i++)
        res+=w[match[i]][i];
    return res;
}
ll A[N],P[N],B[N],C[N];
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++) scanf("%lld",&A[i]);
    for(int i=1;i<=n;i++) scanf("%lld",&P[i]);
    for(int i=1;i<=n;i++) scanf("%lld",&B[i]);
    for(int i=1;i<=n;i++) scanf("%lld",&C[i]);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        {
            ll sum=0;
            for(int k=1;k<=n;k++)
                if(A[k]<B[i]+C[j])
                    sum+=P[k];
            w[i][j]=sum;
        }
    printf("%lld\n",KM());
    return 0;
}
费用流

时间复杂度 Θ ( n m f ) \Theta(nmf) Θ(nmf)
此题满流，边级别是 n 2 n^2 n2，于是费用流时间复杂度 Θ ( n 4 ) \Theta(n^4) Θ(n4)
卡费用流就离谱，不得不多看一个KM算法www

#define IO ios::sync_with_stdio(false);cin.tie();cout.tie(0)
#pragma GCC optimize(2)
#include<iostream>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;
const int N=810,M=330010,INF=0x3f3f3f3f;
int n,S,T;
int h[N],e[M],ne[M],w[M],f[M],idx;
int d[N],flow[N],pre[N];
bool st[N];
queue<int> q;
void add(int a,int b,int c,int d)
{
    e[idx]=b,f[idx]=c,w[idx]=d,ne[idx]=h[a],h[a]=idx++;
    e[idx]=a,f[idx]=0,w[idx]=-d,ne[idx]=h[b],h[b]=idx++;
}
bool spfa()
{
    memset(d,-0x3f,sizeof d);
    memset(flow,0,sizeof flow);
    d[S]=0,flow[S]=INF;
    q.push(S);
    while(q.size())
    {
        int t=q.front();q.pop();
        st[t]=0;
        for(int i=h[t];i!=-1;i=ne[i])
        {
            int j=e[i];
            if(f[i]&&d[j]<d[t]+w[i])
            {
                d[j]=d[t]+w[i];
                pre[j]=i;
                flow[j]=min(f[i],flow[t]);
                if(!st[j])
                {
                    q.push(j);
                    st[j]=1;
                }
            }
        }
    }
    return flow[T];
}
int EK()
{
    int cost=0;
    while(spfa())
    {
        int t=flow[T];
        cost+=d[T]*t;
        for(int i=T;i!=S;i=e[pre[i]^1])
        {
            f[pre[i]]-=t;
            f[pre[i]^1]+=t;
        }
    }
    return cost;
}
ll A[N],P[N],B[N],C[N];
int main()
{
    IO;
    cin>>n;
    S=0,T=2*n+1;
    memset(h,-1,sizeof h);
    for(int i=1;i<=n;i++) cin>>A[i];
    for(int i=1;i<=n;i++) cin>>P[i];
    for(int i=1;i<=n;i++) cin>>B[i];
    for(int i=1;i<=n;i++) cin>>C[i];
    for(int i=1;i<=n;i++) 
        add(S,i,1,0),add(i+n,T,1,0);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        {
            int sum=0;
            for(int k=1;k<=n;k++)
                if(A[k]<B[i]+C[j])
                    sum+=P[k];
            if(sum) add(i,j+n,INF,sum);
        }
    cout<<EK()<<'\n';
    return 0;
}

KM算法

P6577 【模板】二分图最大权完美匹配
KM算法——二分图最大带权匹配
【算法笔记】二分图最大权匹配 - KM算法（dfs版O(n4) + bfs版O(n3)）

  • 顶标：两边点都有的标记（左 x i x_i xi​右 y i y_i yi​）满足 x i + y i ≥ w i , j x_i+y_i\ge w_{i,j} xi​+yi​≥wi,j​，不唯一

  • 相等边： x i + y i = w i , j x_i+y_i=w_{i,j} xi​+yi​=wi,j​的边 相等子图：相等边构成的子图。

  • 交错树：增广路径形成的树。

  • KM算法核心定理 :对于某组可行顶标，如果其相等子图存在完美匹配，那么，该匹配就是原二分图的最大权完美匹配。

dfs时间复杂度 O ( n 4 ) O(n^4) O(n4)
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=510;
ll w[N][N],lx[N],ly[N];//顶标
bool visx[N],visy[N];//记录是否在交错树上
int match[N];//匹配
int n,m;
ll delta;//delta和更新后的delta
bool dfs(int x)
{
    visx[x]=1;
    for(int y=1;y<=n;y++)
        if(!visy[y])
        {
            if(lx[x]+ly[y]==w[x][y])
            {
                visy[y]=1;
                if(!match[y]||dfs(match[y])) 
                    return match[y]=x,1;
            }
            else delta=min(delta,lx[x]+ly[y]-w[x][y]);// visx[x]=1,visy[y]=0更新delta
        }
    return 0;
}
ll KM()
{
    memset(lx,-0x3f,sizeof lx);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            lx[i]=max(lx[i],w[i][j]);
    for(int i=1;i<=n;i++)
        while(1)
        {
            memset(visx,0,sizeof visx);
            memset(visy,0,sizeof visy);
            delta=1e18;
            if(dfs(i)) break;
            for(int j=1;j<=n;j++)
            {//把所有在交错树上的左右部点都左加右减 Δ
                if(visx[j]) lx[j]-=delta;
                if(visy[j]) ly[j]+=delta;
            }
        }
    ll res=0;
    for(int i=1;i<=n;i++)
        res+=w[match[i]][i];
    return res;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++) 
            w[i][j]=-1e18;
    for(int i=1;i<=m;i++)
    {
        int a,b;ll c;
        scanf("%d%d%lld",&a,&b,&c);
        w[a][b]=max(w[a][b],c);
    }
    printf("%lld\n",KM());
    for(int i=1;i<=n;i++)
        printf("%d ",match[i]);
    puts("");
    return 0;
}

KM算法详解和其更快优化方法(吊打传统的slack优化)

bfs 时间复杂度 O ( n 3 ) O(n^3) O(n3)

把DFS换成 BFS。本质上没有什么区别
但是每个左边的点只会进队、搜索一次。
用p数组记录的是增广交错树。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=510;
ll w[N][N],lx[N],ly[N];//顶标
bool visx[N],visy[N];//记录是否在交错树上
int match[N];//匹配
int n,m,p[N];
ll delta,c[N];//delta和更新后的delta
void bfs(int x)
{
    int a,y=0,y1=0;
    for(int i=1;i<=n;i++) p[i]=0,c[i]=1e18;
    match[y]=x;
    do{
        a=match[y],delta=1e18,visy[y]=1;
        for(int b=1;b<=n;b++)
            if(!visy[b])
            {
                if(c[b]>lx[a]+ly[b]-w[a][b])
                    c[b]=lx[a]+ly[b]-w[a][b],p[b]=y;
                if(c[b]<delta)//Δ还是取最小的
                    delta=c[b],y1=b;
            }
        for(int b=0;b<=n;b++)
            if(visy[b])
                lx[match[b]]-=delta,ly[b]+=delta;
            else
                c[b]-=delta;
        y=y1;
    }while(match[y]);
    while(y) match[y]=match[p[y]],y=p[y];
}
ll KM()
{
    for(int i=1;i<=n;i++)
        match[i]=lx[i]=ly[i]=0;
    for(int i=1;i<=n;i++)
    {
        memset(visy,0,sizeof visy);
        bfs(i);
    }
    
    ll res=0;
    for(int i=1;i<=n;i++)
        res+=w[match[i]][i];
    return res;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++) 
            w[i][j]=-1e18;
    for(int i=1;i<=m;i++)
    {
        int a,b;ll c;
        scanf("%d%d%lld",&a,&b,&c);
        w[a][b]=max(w[a][b],c);
    }
    printf("%lld\n",KM());
    for(int i=1;i<=n;i++)
        printf("%d ",match[i]);
    puts("");
    return 0;
}

要加油哦~
板子还是记住吧，具体啥的也没啥

秋天俯身采种子
关注
专栏目录
二分图判定，二分图匹配（匈牙利算法），多重匹配，最小顶点覆盖，带权二分图匹配（KM算法）
wa的一声就哭了
01-03 565
目录 二分图判定： 最大匹配 最小点覆盖 二分图多重匹配 带权二分图最大匹配 二分图判定： 两组顶点，每组顶点内没有相交。 如何判断二分图？ 递归，染色。 随意将一顶点染色为1，将与其相连的染色为-1，再将与-1相连的点染为1，如此重复直到遍历完整个图 如果有矛盾（两个相连点染色相同）则不为二分图 代码 int dfs(int u,int sum) { f...
poj 3189(多重匹配算法求解)
08-27 634
多重匹配的流解法：http://www.cppblog.com/MatoNo1/archive/2011/03/26/142766.aspx 多重匹配问题：二分图的最大匹配中左部图和右部图的顶点是一一对应的，然而在多重匹配中，左部图的的点可以和右部图的多个点匹配，但是右
J.Spy
spnooyseed的博客
12-02 328
Over time, those coaches in BDN Collegiate Progranuning Contest split into two camps.The big danger is that, while optimists and pessimists battle it out, the environment of this area becomes ever mor...
二分图带权匹配、最佳匹配与KM算法
weixin_30780221的博客
08-05 109
---------------------以上转自ByVoid神牛博客，并有所省略。 [二分图带权匹配与最佳匹配] 什么是二分图的带权匹配？二分图的带权匹配就是求出一个匹配集合，使得集合中边的权值之和最大或最小。而二分图的最佳匹配则一定为完备匹配，在此基础上，才要求匹配的边权值之和最大或最小。二分图的带权匹配与最佳匹配不等价，也不互相包含。 我们可以使用KM算法实现求二分图的最佳匹...
2019ICPC南京J - Spy
EIP_silly的博客
12-13 659
2019ICPC南京J - Spy 有两个队伍进行比拼 两边各有n个人 第一行表示Amy队伍每个人的战斗力 第二行表示击败每个人能获得的声望值 第三行表示Bob第一只队伍的战斗力 第三行表示Bob第二支队伍的战斗力 Bob的两支队伍的队员两两配，战斗力相加，然后和Amy的队员进行战斗。 每个人只能战斗一次。 如果Bob方的战斗力大于Amy的就获得声望值，如果战斗力≤Amy队员战斗力，就无法获...
2019 ICPC Asia Nanjing Regional J.Spy（KM算法O(n^3)板子题）
繁凡さん的博客
10-09 591
整理的算法模板合集： ACM模板 前面好几段又在讲故事… 题目大意： 　　a[i]表示对手的每个队伍战斗力 　　p[i]表示打败对手后获得的分数 　　b[i]表示我方第一种人的战斗力 　　c[i]表示我方第二种人的战斗力 　　定义我方一组选手的战斗力为b[i]+c[j]，第一种选手与第二种选手某种顺序两两组队后，与对方进行pk，共有 n!n!n! 种pk顺序，求最大期望×n 期望就是加权平均，期望*n实际上就是所有方案中的最大权值，也就是二分图最大加权匹配 所以这是一个板子题，但是数据卡O(n4)O(n
Card.Spy-darwin-x64_zip_
10-01
标题中的"Card.Spy-darwin-x64_zip_"表明这是一个与操作系统安全和恶意软件相关的主题，特别是针对苹果的Darwin/x64平台。Darwin是Apple macOS的基础，基于UNIX，而x64指的是64位处理器架构。"zip"标签则提示我们这...
microsoft.spy
05-03
"microsoft.spy" 是一款专用于Windows窗口事件跟踪查看的工具，主要针对开发者和系统管理员，帮助他们深入了解操作系统中窗口的交互过程。该工具能够记录并显示系统中的各种窗口事件，如窗口创建、销毁、大小调整、...
org.eclipse.tools.layout.spy.nl_zh_4.20.0.v20210630020001.jar
最新发布
04-26
Eclipse-plugins 插件，eclipse 各种插件 .jar 包，免费下载 eclipse202106 -plugins 各种插件 .jar 包，免费下载 如果下载不了，关注我，评论区联系我， Eclipse-plugins eclipse, plugins, eclipse202106, ...
mybatis-plus p6spy 报错，java.lang.NoClassDefFoundError: Could not initialize class com.p6spy.engine.sp
HI,我是小瑞！
12-02 2616
按照 mybatis-plus的文档来加入 p6spy 后，执行单元测试出现两个错误， https://baomidou.com/guide/p6spy.html 执行 SQL 分析打印 该功能依赖p6spy组件，完美的输出打印 SQL 及执行时长3.1.0以上版本 错误1： java.lang.NoClassDefFoundError: Could not initialize class com.p6spy.engine.spy.P6ModuleManager at com.p...
二分图匹配问题（匈牙利及KM算法）
05-15
资源为ppt格式，主要讲解匈牙利算法和KM算法等二分图匹配的问题
匹配优化--KM算法 详解+模板
蓝星潮
11-17 1372
KM算法 详解+模板   (2012-07-11 11:20:00) 转载▼ 标签：  km算法   详解   模板   讲解   slack 分类： 图论 先说KM算法求二分图的最佳匹配思想，再详讲KM的实现。 【KM算法求二分图的最佳匹配思想】 对于具有二部划分( V1, V2 )的加权完
[KM算法]uoj#80 二分图最大权匹配
weixin_34128501的博客
02-03 144
题意：给定二分图，有边权，求最大边权匹配。边权非负。 解：KM算法求解最大权完备匹配。 完备匹配就是点数少的那一边每个点都有匹配。 为了让完备匹配与最大权匹配等价，我们添加若干条0边使之成为完全二分图(自造名词别在意......) 为了让左边成为点数较少的一边，我们还要添加一些虚点，m = max(n,m) 然后求解完备匹配。 KM的DFS写法会被卡成n4，如果你不在意可以写........
Going Home（最大匹配km算法）
weixin_33753845的博客
11-12 65
Going Home Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 20115   Accepted: 10189 Description On a grid map there are n little men and n houses. In each unit ...
二分图-多重匹配
zxy160的博客
07-16 1208
poj2289 二分图匹配一个点只能匹配一个，但是多重匹配能匹配多个，只要改变限制条件就可以了，先判断这个点匹配的数量有没有达到它自身或二分mid的限制，如果达到限制可以从匹配的点中选取能转移的匹配点。 //感觉kuangbin板子上的代码太费时了，还是自己写的好 #include&lt;stdio.h&gt; #include&lt;sstream&gt; #include&lt;str...
poj 4468Spy(kmp算法)
爱情魔法师的专栏
07-23 916
Spy Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 204    Accepted Submission(s): 96 Problem Description “Be subtle! Be subtle! An
序列匹配(五)——重复匹配问题的动态规划算法
生信了
04-30 3461
原创： hxj7 前言： 蛋白质序列中常有重复的功能域（domain）或模体（motif）拷贝，由此衍生出一个抽象的序列多重匹配的问题，即如何从一个序列中找出另一个序列的某部分（如功能域或模体）的多个无交叠（non-overlapping)拷贝。本文给出了该问题的示例、关键计算公式以及C语言实现代码。 问题及算法描述 更具体地描述上面的问题：有序列x和y，其中y是包含结构域的序列，x是要从中找...
二分图最大权值匹配KM算法
smsmn的专栏
11-06 848
<br />#include <cstdio> #include <climits> #include <iostream> using namespace std; const int MAX = 1024; int n; // X 的大小int weight [MAX] [MAX]; // X 到 Y 的映射（权重） int lx [MAX], ly [MAX]; // 标号 bool sx [MAX], sy [MAX]; // 是否被搜索过 int match [MAX];
java实现多模匹配算法_多模字符串匹配算法-Aho–Corasick
weixin_42512903的博客
02-16 360
背景在做实际工作中，最简单也最常用的一种自然语言处理方法就是关键词匹配，例如我们要对n条文本进行过滤，那本身是一个过滤词表的，通常进行过滤的代码如下for (String document : documents) {for (String filterWord : filterWords) {if (document.contains(filterWord)) {//process ...}}}...
net.spy.spymemcached net.spy.spymemcached
05-11
net.spy.spymemcached是一个用于Java的开源Memcached客户端库，用于与Memcached NoSQL数据库进行交互。它提供了简单的API，允许Java应用程序通过Memcached服务器存储和检索数据。net.spy.spymemcached支持Memcached...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值