首页 > 编程语言 >通关搜索和图论 day_17 -- 染色法 & 匈牙利算法

通关搜索和图论 day_17 -- 染色法 & 匈牙利算法

时间:2023-01-13 11:35:54浏览次数:63  
标签:二分 return idx 17 -- ne int 染色法 color

染色法

一个图是二分图 当且仅当 她可以被2染色(不含有奇数环)

流程如下,先找到一个不在集合中的点 把他放在左边

然后遍历这个点有连接的点,把这些点放到右边,再依次遍历放到右边的点的邻点

模板

int n;      // n表示点数
int h[N], e[M], ne[M], idx;     // 邻接表存储图
int color[N];       // 表示每个点的颜色,-1表示未染色,0表示白色,1表示黑色

// 参数:u表示当前节点,c表示当前点的颜色
bool dfs(int u, int c)
{
    color[u] = c;
    for (int i = h[u]; i != -1; i = ne[i])
    {
        int j = e[i];
        if (color[j] == -1)
        {
            if (!dfs(j, !c)) return false;
        }
        else if (color[j] == c) return false;
    }

    return true;
}

bool check()
{
    memset(color, -1, sizeof color);
    bool flag = true;
    for (int i = 1; i <= n; i ++ )
        if (color[i] == -1)
            if (!dfs(i, 0))
            {
                flag = false;
                break;
            }
    return flag;
}

例题

860. 染色法判定二分图 - AcWing题库

给定一个 n 个点 m 条边的无向图,图中可能存在重边和自环。

请你判断这个图是否是二分图。

输入格式

第一行包含两个整数 n 和 m。

接下来 m 行,每行包含两个整数 u 和 v,表示点 u 和点 v 之间存在一条边。

输出格式

如果给定图是二分图,则输出 Yes,否则输出 No

数据范围

1≤n,m≤105

输入样例:

4 4
1 3
1 4
2 3
2 4

输出样例:

Yes
#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

const int N = 100010,M = 200010;
int n,m;
bool flag = true;
int h[N],e[M],ne[M],idx;
int color[N];

void add(int a,int b)
{
    e[idx] = b;
    ne[idx] = h[a];
    h[a] = idx++;
}

bool dfs(int u,int c)
{
    color[u] = c;

    for (int i = h[u]; i != -1 ; i = ne[i]) {
        int j = e[i];
        if (!color[j])
        {
            if (!dfs(j,3-c)) return false;
        }else if(color[j] == c){
            return false;
        }
    }
    return true;
}

int main()
{
    cin.tie(0);
    ios::sync_with_stdio(false);

    cin >> n >> m;

    memset(h,-1,sizeof h);

    while(m--)
    {
        int a,b;
        cin >> a >> b;
        add(a,b);
        add(b,a);
    }

    for (int i = 1; i <= n ; ++i) {
        if (!color[i])
        {
            if (!dfs(i,1))
            {
                flag = false;
                break;
            }
        }
    }

    if (flag)
    {
        cout << "Yes" << endl;
    }else{
        cout << "No" << endl;
    }
}

匈牙利算法

可以求出,左边和右边匹配成功的最大数量是多少,成功是指不存在两条边是共用的两个点

例如如下图

image-20230113104514394

我们从左一开始找,找到右二进行尝试,此时右二并没有配对,所以右二和左一配对成功

依次类推左二可以和右一配对

但是左三只能和右二配对但是右二已经和左一配对

这个时候我们找到右二配对的左一 是否有其他边,显然左一还可以和右四配对

这样最终匹配的数量就是 4

模板

int n1, n2;     // n1表示第一个集合中的点数,n2表示第二个集合中的点数
int h[N], e[M], ne[M], idx;     // 邻接表存储所有边,匈牙利算法中只会用到从第一个集合指向第二个集合的边,所以这里只用存一个方向的边
int match[N];       // 存储第二个集合中的每个点当前匹配的第一个集合中的点是哪个
bool st[N];     // 表示第二个集合中的每个点是否已经被遍历过

bool find(int x)
{
    for (int i = h[x]; i != -1; i = ne[i])
    {
        int j = e[i];
        if (!st[j])
        {
            st[j] = true;
            if (match[j] == 0 || find(match[j]))
            {
                match[j] = x;
                return true;
            }
        }
    }

    return false;
}

// 求最大匹配数,依次枚举第一个集合中的每个点能否匹配第二个集合中的点
int res = 0;
for (int i = 1; i <= n1; i ++ )
{
    memset(st, false, sizeof st);
    if (find(i)) res ++ ;
}

例题

861. 二分图的最大匹配 - AcWing题库

给定一个二分图,其中左半部包含 n1 个点(编号 1∼n1),右半部包含 n2 个点(编号 1∼n2),二分图共包含 m 条边。

数据保证任意一条边的两个端点都不可能在同一部分中。

请你求出二分图的最大匹配数。

二分图的匹配:给定一个二分图 G,在 G 的一个子图 M 中,M 的边集 {E} 中的任意两条边都不依附于同一个顶点,则称 M 是一个匹配。

二分图的最大匹配:所有匹配中包含边数最多的一组匹配被称为二分图的最大匹配,其边数即为最大匹配数。

输入格式

第一行包含三个整数 n1、 n2 和 m。

接下来 m 行,每行包含两个整数 u 和 v,表示左半部点集中的点 u 和右半部点集中的点 v 之间存在一条边。

输出格式

输出一个整数,表示二分图的最大匹配数。

数据范围

1≤n1,n2≤500,
1≤u≤n1,
1≤v≤n2,
1≤m≤105

输入样例:

2 2 4
1 1
1 2
2 1
2 2

输出样例:

2
#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

const int N = 100010,M = 200010;
int n1,n2,m;
int h[N],e[M],ne[M],idx;
int match[N];   // 右边的点对应的点
bool st[N]; //判重,每次不要重复搜一个点



void add(int a,int b)
{
    e[idx] = b;
    ne[idx] = h[a];
    h[a] = idx++;
}

bool find(int x)
{
    // 枚举这个点所有可能的点
    for (int i = h[x]; i != -1 ; i = ne[i]) {
        int j = e[i];
        if (!st[j])
        {
            st[j] = true;
            if (match[j] == 0 || find(match[j]))
            {
                match[j] = x;
                return true;
            }
        }
    }
    return false;
}


int main()
{
    cin.tie(0);
    ios::sync_with_stdio(false);

    cin >> n1 >> n2 >> m;

    memset(h,-1,sizeof h);

    while(m--)
    {
        int a,b;
        cin >> a >> b;
        add(a,b);
    }

    int res = 0;
    for (int i = 1; i <= n1; ++i) {
        memset(st,false,sizeof st);
        if (find(i))
        {
            res++;
        }
    }

    cout << res << endl;
    return 0;
}

标签:二分,return,idx,17,--,ne,int,染色法,color
From: https://www.cnblogs.com/ShibuyaKanon/p/17049076.html

相关文章

  • Thanos 配置Prometheus的高可用 (一)
    标签(空格分隔):Prometheus系列一:Prometheus的介绍与架构1.1Prometheus的概述1.prometheus的介绍Prometheus是一个开源的系统监控和告警工具包,最初由SoundCloud开......
  • (六)elasticsearch 源码之选主流程分析
    1.概述es(elasticsearch,下同)的选举流程相对来说比较简单,使用的bully算法,简而言之,就是谁强谁就是老大,待会儿看下怎么判定谁更强。2.选主流程在启动篇中我们讲解了节点启动......
  • 新一代云原生日志架构 - Loggie的设计与实践
    Loggie萌芽于网易严选业务的实际需求,成长于严选与数帆的长期共建,持续发展于网易数帆与网易传媒、中国工商银行的紧密协作。广泛的生态,使得项目能够基于业务需求不断完善、......
  • shiro登录认证过程讲解
    先粘出登录的代码​@RequestMapping(value="/submitLogin",method=RequestMethod.POST)@ResponseBodypublicMap<String,Object>submitLogin(Stringusername,String......
  • MySQL必知必会第七章-数据过滤
    数据过滤组合WHERE子句操作符(operator)用来联结或改变WHERE子句中的子句的关键字。也称为逻辑操作符(logicaloperator)。AND操作符为了通过不止一个列进行过滤,可使用AND......
  • java中关于继承,多态及方法调用的底层细节
    java中关于继承,多态及方法调用的底层细节一、继承继承已存在的类就是复用(继承)这些类的方法和域。在此基础上,还可以添加一些新的方法和域,以满足新的需求。子类会拥有......
  • 河北稳控科技振弦采集模块的各种参数操作
    河北稳控科技振弦采集模块的各种参数操作固件版本读取点击指令区【读取版本】按钮,读取当前连接模块的固件版本信息,读取到的版本信息显示于按钮右侧。VMTool会根据读取......
  • CKS 考试题整理 (01)-NetworkPolicy
       Task创建一个名为pod-restriction的NetworkPolicy来限制对在namespacedev-team中运行的Podproducts-service的访问。 只运行以下Pod连接到Pod......
  • Visual Studio Code快捷键的使用-RK3568开发
    在VisualStudioCode软件中使用快捷键可以增加自身的效率,一些常用的快捷键如下所示:F1打开可以输入命令F2重命名变量,方便重构F5运行和调试代码F12去到定义的地方shift......
  • Solution Set - 寒假睡眠记录
    目录「2018集训队互测」「LOJ#2504」小H爱染色  前面有一些题解还没补啊,但有人急着看就先发出来啦.「2018集训队互测」「LOJ#2504」小H爱染色  Link&Su......