Warshall‘s algorithm 算法的实现及优化

本文介绍Warshall算法用于计算传递闭包的过程及优化方法,详细解释了该算法的时间复杂度,并通过经典代码片段展示了如何实现这一算法。此外还讨论了如何通过矩阵特性进一步优化算法效率。

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

用 Warshall’s 算法计算传递闭包

(1)时间复杂度为:O(nnn)
(2)用R的无穷闭包时间复杂度为O(nnn*(n - 1))

经典实现:

下面展示一些 经典代码片

void computeAPSP(const int n) {
    /* calculate shortest paths from every vertex to every vertex */
    for (int k = 0; k < n; k++) 
    {
        for (int i = 0; i < n; i++) 
        {
            for (int j = 0; j < n; j++) 
            {
                a[i][j] = min( a[i][j], a[i][k] + a[k][j] );
            }
        }
    }
}

利用矩阵的对称性优化:

void computeAPSP(const int n) 
{
    for (int k = 0; k < n; k++) 
    {
        for (int i = 0; i < n; i++) 
        {
            if (k != i) 
            {

                const int a_ki = (k < i) ? a[i][k] : a[k][i];

                for (int j = 0; j < min(k, i); j++)
                    a[i][j] = min( a[i][j], a_ki + a[k][j] );

                for (int j = k + 1; j < i; j++)
                    a[i][j] = min( a[i][j], a_ki + a[j][k] );

            }
        }
    }
}

只使用矩阵的下三角部分进行优化:

void computeAPSP(const int n) 
{
    for (int k = 0; k < n; k++) 
    {
        for (int i = 0; i < n; i++) 
        {
            if (k != i) 
            {

                const int a_ki = (k < i) ? a[i][k] : a[k][i];

                for (int j = 0; j < min(k, i); j++)
                    a[i][j] = min( a[i][j], a_ki + a[k][j] );

                for (int j = k + 1; j < i; j++)
                    a[i][j] = min( a[i][j], a_ki + a[j][k] );

            }
        }
    }
}

避免大量调用数学函数进行优化:

void computeAPSP(const int n) 
{
    for (int k = 0; k < n; k++) 
    {
        for (int i = 0; i < n; i++) 
        {
            if (k != i) 
            {

                const int a_ki = (k < i) ? a[i][k] : a[k][i];

                // skip if no path
                if (a_ki == POSITIVE_INFINITY)    continue;

                for (int j = 0; j < min(k, i); j++) 
                {
                    const int s_kj = a_ki + a[k][j];
                    if( s_kj < a[i][j] )    a[i][j] = s_kj;
                }

                for (int j = k + 1; j < i; j++) 
                {
                    const int s_jk = a_ki + a[j][k];
                    if( s_jk < a[i][j] )    a[i][j] = s_jk;
                }

            }
        }
    }
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值