最小生成树|二分图

最小生成树跟边的正负没有任何关系。

最小生成树

朴素Prime

该算法和Dijkstr算法很像。

1. 先把所有距离初始化为正无穷

2. 进行n次迭代

1. 找到不在集合（集合指当前的生成树）当中的点，s数组表示当前已经在连通块（生成树）中的所有点。

2. 找到集合外距离最近的点赋给t，用t更新其它点到集合的距离

3. s加入到t当中去。

首先让所有点到集合的距离都是正无穷，然后找到一个集合外到集合距离最近的点，这四个点距离都是正无穷，我们则可以随便挑一个点

这里我们挑1号点，之后找到1号点到集合的最短距离，若该点没有边连到集合，则该点到集合的距离还是正无穷。

若找到了1号点到集合的最短距离，则用1号点去更新所有其它点到集合的距离

之后把1号点（绿点）加到集合当中去。

第二次迭代，从剩下的点当中选一个到集合距离最近的点，这里2号点最近，选中2号点之后，用2号点更新其它点到集合的距离（这个例子中更新前后结果不变，因为3到1和3到2的距离都是2，而且2和4没有连接），更新完之后，把2号点加到集合里面去。

重复进行上述操作 ，在3和4中找距离最近的点，这里3最近，用3更新4到集合的距离时，3到4的距离是4，比1到4的距离3要大，所以这里不更新。我们的目的是将距离更新的更短，距离如果变长，我们就不更新。

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N= 510,INF=0x3f3f3f3f;
int n, m;
int g[N][N];//某俩点间的距离
int dist[N];//某点到集合的距离
bool st[N];
int prim()
{
    memset(dist, 0x3f, sizeof dist);
    int res = 0;//存储最小生成树所有边的长度之和
    for (int i = 0; i < n; ++i)
    {
        int t = -1;
        for (int j = 1; j <= n; j++)//找到集合外所有距离最小的点
        {
            if (!st[j] && (t == -1 || dist[t] > dist[j]))//t=-1表示当前还没有找到距离最小的点
                //如果t代表的点到集合的距离大于j代表的点到集合的距离，也表示当前没找到
                t = j;
        }
            if (i && dist[t] == INF)//说明当前距离最近的点到集合的距离是正无穷，即所有点没有连通
            {
                return INF;
            }
            if (i) res += dist[t];//dist[t]表示当前点和另一个点（该点已与集合连好）的距离
            for (int j = 1; j <=n; ++j)
            {
                dist[j] = min(dist[j], g[t][j]);//更新其它点
            }
        
            st[t] = true;//修改状态
        
    }
    return res;
}
int main()
{
    scanf("%d %d", &n, &m);
    memset(g, 0x3f, sizeof g);
    while (m--)//读入所有边
    {
        int a, b, c;
        scanf("%d %d %d", &a, &b, &c);
        g[a][b] = g[b][a]=min(g[a][b],c);//由于不用考虑边的正负，a到b的距离和b到a的距离都一样,因为可能有重边，我们取一个最小值
    }
    int t = prim();
    if (t == INF)//不存在生成树
        puts("impossible");
    else//存在生成树
        printf("%d", t);
    return 0;
}

克鲁斯卡尔算法

1. 先将所有边按权重从小到大进行排序

2. 枚举每条边a，b权重为c

1. 如果当前a，b不连通，就把a，b这条边加到集合里面来

第二步其实是前面所学并查集的应用

#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 200010;
int n, m;
int p[N];

struct Edge
{
    int a, b, w;//a,b是边的俩条点,w是权重
    bool operator<(const Edge& W) const//重载小于号方便排序
    {
        return w < W.w;
    }
}edges[N];
int find(int x)
{
    if (p[x] != x)//如果x不是祖宗节点
        p[x] = find(p[x]);
    return p[x];
}
int main()
{
    scanf("%d %d", &n, &m);
    for (int i = 0; i < m; ++i)
    {
        int a, b, w;
        scanf("%d %d %d", &a, &b, &w);//把所有的边读进来
        edges[i] = { a,b,w };
    }
    sort(edges, edges + m);//把所有的边排序
    for (int i = 1; i <=n; ++i)//初始化并查集
        p[i] = i;
    int res = 0, cnt = 0;//res存的是最小生成树当中所有树边的权重之和，cnt是当前加了多少条边
    for (int i = 0; i < m; ++i)//从小到大枚举所有边
    {
        int a = edges[i].a, b = edges[i].b, w = edges[i].w;
        a = find(a), b = find(b);//a=a的祖宗节点，b=b的祖宗节点，find函数是并查集的find函数
        if (a != b)//如果俩个祖宗节点不连通
        {
            p[a] = b;//这条边加进来之后，要把俩个点各自所在的集合进行合并
            res += w;
            cnt++;
        }
    }
    if (cnt < n - 1)//如果连的边数<n-1，说明这个树是不连通的
        puts("impossible");
    else
        printf("%d\n",res);
    return 0;
}

二分图

二分图定义：所有点划分到俩边区，使所有的边都是在集合之间的，集合内部没有边。如果一个图中有奇数环，则一定不是二分图。把图分为俩个区域，一个是左边，另一个是右边。

染色法

如果一个图中不含有奇数环则一定是二分图。

1. 从前往后遍历每一个点

2. 如果某个点还没分好组，我们就分到左边去

1. 分完该点之后，遍历这个点和所有的邻点（和这个点连通的所有点），若该点属于左边，则它的所有相邻点属于右边，若该点属于右边，则它的所有相邻点属于左边，下图中第一个点属于左边，1代表左边，2代表右边，按照这种方式，将所有点可表示出来，其实就是在染色。

2. 染色的时候要注意：一条边的俩端点的颜色，一定不同，即这俩个点不能属于同一个集合。通过这种方式将图中的所有点进行染色。由于图中不存在奇数环，所以染色过程中不会出现矛盾。

只要一个图可以用染色法染完，而且不出现矛盾，它就是一个二分图。

这里采用深度优先遍历进行染色。

#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 1010, M = 2010;//由于是无向图，每条无向边要存俩条有向边，即俩个点之间有俩条线,所以M是N的2倍
int n, m;
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;//记录一下当前点颜色是C
    for (int i = h[u]; i != -1; i = ne[i])
    {
        int j = e[i];//存储当前点的编号
        if (!color[j])//如果当前点没有颜色
        {
            if (!dfs(j,3-c)) //如果当前是颜色1，就染成2，如果是颜色2，就染成1，这里写成3-c
                return false;//如果染色不成功就返回false,成功返回true
        }
        else if (color[j]==c)//如果这俩个相等，代表一条边的俩端点颜色一样，也有矛盾
        {
            return false;
        }
    }
    return true;
}
int main()
{
    scanf("%d %d", &n, &m);
    memset(h, -1, sizeof h);
    while (m--)
    {
        int a, b;
        scanf("%d %d", &a, &b);
        add(a, b);
        add(b, a);//由于是无向边，所以这里加俩条边，一条a到b，一条b到a

    }
    //开始染色
    bool flag = true;//标志位，表示染色成功或失败
    for(int i=1;i<=n;++i)
        if (!color[i])
        {
            if (!dfs(i, 1))//如果有矛盾dfs返回fasle,这里传1是想把有矛盾的地方染成第一种颜色，
            {
                flag = false;
                break;
            }
        }
    if (flag)
        puts("Yes");
    else
        puts("No");//有矛盾发生
    return 0;

}

匈牙利算法

假设有这样的一个二分图

匈牙利算法可以在一个比较块的时间内，算出左边和右边最大的匹配成功（不存在俩条边共用一个点）数是多少。

1.如果a和b点已经匹配好了，c点要和b匹配，此时c点找一下有没有其它能和a匹配的点，如果有，就让a和其它点去匹配，然后c自己和b匹配。

#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 510, M = 100010;
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]))//0表示还没匹配,或者说已经匹配到了左边的某点，但左边的这个点可以去匹配别的点
            {
                match[j] = x;
                return true;
            }
        }
    }
    return false;
}
int main()
{
    scanf("%d %d %d", &n1, &n2, &m);
    memset(h, -1, sizeof h);
    while (m--)
    {
        int a, b;
        scanf("%d %d", &a, &b);
        add(a, b);
    }
    int res = 0;//当前匹配的数量
    for (int i = 1; i <= n1; ++i)
    {
        memset(st, false, sizeof st);
        if (find(i)) res++;//如果左边右边匹配成功，直接++
    }
    printf("%d\n", res);
    return 0;
}