图论——最小环

图论——最小环

对于一个有向图或一个无向图，如何求它的最小环的大小呢，默认边的权都是1。

Floyd算法求最小环（性能差）

设节点$i$和$j$为遍历节点$k$为中点，那么当$i=j$的时候，如果存在$i\to k$和$k\to j$，那么就会存在环$i\to k\to j(i)$此时最短路就是环的大小。

#include <bits/stdc++.h>

using namespace std;

#define FR freopen("in.txt","r",stdin)

#define UNC 1000

typedef long long ll;
int d[100][100];
int n,m;


int flody()
{
    int ans = UNC;
    for(int k = 1; k<=n; k++)
        for(int i = 1; i<=n; i++)
            for(int j = 1; j<=n; j++)
            {
                d[i][j] = min(d[i][j],d[i][k] + d[k][j]);
                if(i == j)
                    ans = min(ans,d[i][k] + d[k][j]);

            }
    return ans;
}

int main()
{
    cin >> n >> m;
    for(int i = 1; i<=n; i++)
        for(int j = 1; j<=n; j++)
            d[i][j] = UNC;

    for(int i = 0; i<m; i++)
    {
        int u,v;
        cin >> u >> v;
        d[u][v] = 1;
    }
    cout << flody();
    return 0;
}

网上还流传一种做法是在$k$循环内做两个内层循环，具体原理应该是把一个环中的节点按照从小到大排序，最小的那个就是$i$，第二大的那就是$j$，最大的就是$k$，因此循环的时候$i<j<k$，避免了重复计算，但是我不知道为什么要套两个循环，像我上面这样不行吗？

注意有向图请注意自环。

拓扑排序（不通过）

注意：从入度为0的节点开始消除，剩下的都是环，这句话没错，但是留下的图一定有环，只能说明存在环，不能计算最小环的大小。

Tarjan算法（通过）

在DFS过程中，如果遇到了后向边，那么就说明遇到了环，利用时间戳计算环的大小，最小大小就是最小环的大小。

P2661

#include <bits/stdc++.h>

using namespace std;

#define FR freopen("in.txt","r",stdin)
#define FW freopen("out1.txt","w",stdout)

typedef long long ll;

struct Edge
{
    int to;
    int nxt;
} e[200005];
int head[200005];
int d[200005];
bool vis[200005];
int ti = 0;
int tot = 0;
int n;

inline void add(int u,int v)
{
    tot++;
    e[tot].to = v;
    e[tot].nxt = head[u];
    head[u] = tot;
}

int ans = INT_MAX;

void dfs(int curr)
{
    ti++;
    d[curr] = ti;

    for(int ne = head[curr]; ne != 0; ne = e[ne].nxt)
    {
        int v = e[ne].to;

        if(d[v] != 0)
        {
            if(!vis[v])
                ans = min(ans,d[curr] - d[v] + 1);
        }
        else
        {
            dfs(v);
        }
    }
    vis[curr] = true;
}

int main()
{
    cin >> n;
    for(int i = 1; i<=n; i++)
    {
        int v;
        cin >> v;
        add(i,v);
    }
    for(int i = 1; i<=n; i++)
    {
        if(!vis[i]) dfs(i);
    }
    cout << ans;
    return 0;
}

并查集（通过）

代码取自anyway。只针对于本题，本题没有环套环，因此b一定是根节点。

#include<cstdio>
#include<iostream>
using namespace std;
int f[200002],d[200002],n,minn,last;   //f保存祖先节点，d保存到其祖先节点的路径长。 
int fa(int x)
{
    if (f[x]!=x)                       //查找时沿途更新祖先节点和路径长。 
    {
        int last=f[x];                 //记录父节点（会在递归中被更新）。 
        f[x]=fa(f[x]);                 //更新祖先节点。 
        d[x]+=d[last];                 //更新路径长（原来连在父节点上）。 
    }
    return f[x];
}
void check(int a,int b)
{
    int x=fa(a),y=fa(b);               //查找祖先节点。 
    if (x!=y) {f[x]=y; d[a]=d[b]+1;}   //若不相连，则连接两点，更新父节点和路径长。 
    else minn=min(minn,d[a]+d[b]+1);   //若已连接，则更新最小环长度。 
    return;
}
int main()
{
    int i,t;
    scanf("%d",&n);
    for (i=1;i<=n;i++) f[i]=i;         //祖先节点初始化为自己，路径长为0。 
    minn=0x7777777;
    for (i=1;i<=n;i++)
    {
        scanf("%d",&t);
        check(i,t);                    //检查当前两点是否已有边相连接。 
    }
    printf("%d",minn);
    return 0;
}