Codeforces 51E Pentagon

最新推荐文章于 2021-02-05 09:44:22 发布

setoy

最新推荐文章于 2021-02-05 09:44:22 发布

阅读量466

点赞数

CC 4.0 BY-SA版权

分类专栏： Algorithms C-C++

本文链接：https://blog.csdn.net/setoy/article/details/104251536

Algorithms 同时被 2 个专栏收录

22 篇文章

订阅专栏

C-C++

13 篇文章

订阅专栏

本文介绍了一种在无向图中统计五元环数量的算法，通过枚举和矩阵运算实现。首先计算长度为2和3的路径数量，然后利用这些路径计数五元环，同时去除重复和无效计数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意

https://codeforces.com/problemset/problem/51/E

在无向图中统计有多少个五元环，顶点数 $\le 700$

算法：枚举

枚举五元环的两个顶点，找到两条各自包含2条和3条的不重复边的环。

$\to v$ 的长度是2的路径条数如何计数呢？可以类似Folyd的传递闭包方法，枚举一个中间点 $x$ 即可，根据乘法原理， $\sum_{u}^{x} = \sum_{u}^{x} \times \sum_{x}^{v}$ ，因此可以用二维矩阵保存 $u, v$ 的路径条数。若 $A$ 是原先的邻接矩阵，那么 $B=A\times A$ 就是长度是2的路径条数矩阵，而 $C=B\times A$ 就是长度是3的路径条数计数矩阵。

再根据乘法原理，包含路径长度分别是2和3的不重复边的环个数是 $\times C$ ，然而，其中有很多是重复的和无效的：

重复：对于一个包含5条边的环而言，上面的乘法会重复计算10次：一个点有两条边，顺时针逆时针构造是一样的，而一个环上共五个顶点
无效：由于长度是2的路径包含了诸如 $\to v \to u$ 这样的重复路径，因此一个三元环加上这样的重复路径是无效的，需要减掉这样的路径，路径条数是三个顶点度数之和-3（自身边重复减了3次）

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;

const int maxn = 700+2;
int n, m;
int deg[maxn], A[maxn][maxn], B[maxn][maxn], C[maxn][maxn];

void multiply(int A[][maxn], int B[][maxn], int C[][maxn])
{
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j <= i; j++)
        {
            for (int k = 0; k < n; k++)
            {
                C[i][j] += A[i][k] * B[k][j];
            }
            C[j][i] = C[i][j];
        }
    }
}


int main()
{
    int u, v;
    scanf("%d%d", &n, &m);
    for (int i = 0; i < m; i++)
    {
        scanf("%d%d", &u, &v);
        deg[u-1]++;
        deg[v-1]++;
        A[u-1][v-1] = A[v-1][u-1] = 1;
    }
    multiply(A, A, B);
    multiply(A, B, C);
    long long ans = 0;
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            ans += (long long)B[i][j] * C[i][j];
        }
    }
    ans /= 10;
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < i; j++)
        {
            if (!A[i][j]) continue;
            for (int k = 0; k < j; k++)
            {
                if (!A[i][k]) continue;
                if (!A[j][k]) continue;
                ans -= (deg[i] + deg[j] + deg[k] - 3);
            }
        }
    }
    printf("%lld\n", ans);
    
    return 0;
}