zoukankan      html  css  js  c++  java
  • KM算法 带权二分匹配 O(n^3)

    #include<cstdio>
    #include<cstdlib>
    #include<cstring>
    #include<string>
    #include<algorithm>
    #include<iostream>
    #include<queue>
    #include<map>
    #include<cmath>
    #include<set>
    #include<stack>
    #define LL long long
    using namespace std;
    const LL MAXN = 550;
    const LL INF = 1000000009;
    LL n, m, e;
    LL love[MAXN][MAXN];   // 记录每个妹子和每个男生的好感度
    LL ex_girl[MAXN];      // 每个妹子的期望值
    LL ex_boy[MAXN];       // 每个男生的期望值
    bool vis_girl[MAXN];    // 记录每一轮匹配匹配过的女生
    bool vis_boy[MAXN];     // 记录每一轮匹配匹配过的男生
    LL match[MAXN];        // 记录每个男生匹配到的妹子 如果没有则为-1
    LL slack[MAXN];        // 记录每个汉子如果能被妹子倾心最少还需要多少期望值
    
    
    bool dfs(int girl)
    {
        vis_girl[girl] = true;
    
        for (int boy = 0; boy < m; ++boy) {
    
            if (vis_boy[boy]) continue; // 每一轮匹配 每个男生只尝试一次
    
            LL gap = ex_girl[girl] + ex_boy[boy] - love[girl][boy];
    
            if (gap == 0) {  // 如果符合要求
                vis_boy[boy] = true;
                if (match[boy] == -1 || dfs(match[boy])) {    // 找到一个没有匹配的男生 或者该男生的妹子可以找到其他人
                    match[boy] = girl;
                    return true;
                }
            }
            else {
                slack[boy] = min(slack[boy], gap);  // slack 可以理解为该男生要得到女生的倾心 还需多少期望值 取最小值 备胎的样子【捂脸
            }
        }
    
        return false;
    }
    
    LL KM()
    {
        if (n > m)return -1;
        memset(match, -1, sizeof match);    // 初始每个男生都没有匹配的女生
        memset(ex_boy, 0, sizeof ex_boy);   // 初始每个男生的期望值为0
    
                                            // 每个女生的初始期望值是与她相连的男生最大的好感度
        for (int i = 0; i < n; ++i) {
            ex_girl[i] = love[i][0];
            for (int j = 1; j < m; ++j) {
                ex_girl[i] = max(ex_girl[i], love[i][j]);
            }
        }
    
        // 尝试为每一个女生解决归宿问题
        for (int i = 0; i < n; ++i) {
    
            fill(slack, slack + m, INF);    // 因为要取最小值 初始化为无穷大
    
            while (1) {
                // 为每个女生解决归宿问题的方法是 :如果找不到就降低期望值,直到找到为止
    
                // 记录每轮匹配中男生女生是否被尝试匹配过
                memset(vis_girl, false, sizeof vis_girl);
                memset(vis_boy, false, sizeof vis_boy);
    
                if (dfs(i)) break;  // 找到归宿 退出
    
                                    // 如果不能找到 就降低期望值
                                    // 最小可降低的期望值
                LL d = INF;
    
                for (int j = 0; j < m; ++j)
                    if (!vis_boy[j]) d = min(d, slack[j]);
                for (int j = 0; j < n; j++)
                {
                    if (vis_girl[j])
                        for (int z = 0; z < m; z++)
                            if (!vis_boy[z] && d > ex_girl[j] + ex_boy[z] - love[j][z])
                                d = ex_girl[j] + ex_boy[z] - love[j][z];
                }
                for (int j = 0; j < n; ++j) {
                    // 所有访问过的女生降低期望值
                    if (vis_girl[j]) ex_girl[j] -= d;
                }
                for (int j = 0; j < m; j++)
                {
                    // 所有访问过的男生增加期望值
                    if (vis_boy[j]) ex_boy[j] += d;
                    // 没有访问过的boy 因为girl们的期望值降低,距离得到女生倾心又进了一步!
                    else slack[j] -= d;
                }
    
    
            }
        }
    
        // 匹配完成 求出所有配对的好感度的和
        LL res = 0, cnt = 0;
        for (int i = 0; i < m; ++i)
        {
            if (match[i] != -1)
            {
                if (love[match[i]][i] < 0) return -1;
                res += love[match[i]][i];
                cnt++;
            }
        }
        if (cnt != n)return -1;
    
        return res;
    }
    
    int main()
    {
        cin.sync_with_stdio(false);
        int cas = 1;
        while (cin >> n >> m >> e)
        {
            for (int i = 0; i < n; i++)for (int j = 0; j < m; j++) love[i][j] = -INF;
            for (int i = 0; i < e; i++)
            {
                LL fr, to, vv;
                cin >> fr >> to >> vv;
                if (vv < 0) continue;
                love[fr][to] = vv;
            }
            cout << "Case " << cas++ << ": ";
            cout << KM() << endl;
        }
        return 0;
    }
  • 相关阅读:
    SpringMVC 返回 json 字符串中文乱码
    两个HTML页面之间传值
    NOIP2018 保卫王国
    Arch Linux使用cisco anyconnect
    CSP模拟赛 巨神兵
    计蒜客The Fake Fake Friends
    Manjaro安装,配置,美化指南
    [国家集训队] Crash 的文明世界
    ZJOI2011 营救皮卡丘
    CF1198E Rectangle Painting 2
  • 原文地址:https://www.cnblogs.com/LukeStepByStep/p/7451794.html
Copyright © 2011-2022 走看看