poj 3613(还是不懂，以后再看看)

最新推荐文章于 2021-09-27 10:55:36 发布

转载最新推荐文章于 2021-09-27 10:55:36 发布 · 247 阅读

收录于

图论

本文探讨了Floyd算法的原理及应用，特别是在求解带有限制条件的最短路径问题时的方法。通过矩阵快速幂的方式进行优化，解决了经过特定数量边的最短路径问题。

题意：给定一个T(2 <= T <= 100)条边的无向图，求S到E恰好经过N(2 <= N <= 1000000)条边的最短路。

以前一直没仔细想过floyd算法，觉得很简单，今天做这题的时候，看网上的报告都有一句：floyd是每次使用一个中间点k去更新i,j之间的距离，那么更新成功表示i,j之间恰有一个点k时的最短路，如果做N - 1次floyd那么不就是i,j之间借助N - 1 个点时的最短路了。看了很久不明白为什么。也对floyd的最外围的那个k<n产生了疑惑。后来突然想到了，floyd每次更新的都是本身，假如k=3时，dis[1][5]借助dis[1][3]和dis[3][5]成为最短路。不需要去知道1,3之间有多少个点，因为前面已经求出，，（dp?）只知道这一次求出的dis[1][5]多加了一个3这个点。

回到这题，floyd算法是对自身矩阵更新，而这道题却是更新到另一个矩阵上，所以不会出现刚更新过的值又来更新。。例如下面代码的b.mat[1][3] c.mat[3][5]就分别代表上面的dis[1][3],dis[3][5].我们不需要知道tmp.mat[1][5]已经有的点个数。（即已经更新的次数。）只知道，这次更新会加入一个3到他们中间。所以更新k-1次就行。

    #include <cstdio>  
    #include <cstring>  
    #include <map>  
    using namespace std;  
      
    int K, n, m, s, t, x, y, z;  
    map<int, int> mp;  
      
    struct Matrix {  
        int a[205][205];  
        Matrix operator * (Matrix &r) {  
            Matrix c;  
            memset(c.a, 0x3f, sizeof c.a);  
            for(int i = 1; i <= n; i++)  
            for(int j = 1; j <= n; j++)  
            for(int k = 1; k <= n; k++)  
                c.a[i][j] = min(c.a[i][j], a[i][k] + r.a[k][j]);  
            return c;  
        }  
    }st, ans;  
      
    void power() {  
        ans = st, K--;  
        while(K) {  
            if(K & 1) ans = ans * st;  
            st = st * st;  
            K >>= 1;  
        }  
    }  
      
    int main() {  
        scanf("%d%d%d%d", &K, &m, &s, &t);  
        memset(st.a, 0x3f, sizeof st.a);  
        while(m--) {  
            scanf("%d%d%d", &z, &x, &y);  
            if(mp[x]) x = mp[x]; else x = mp[x] = ++n;  
            if(mp[y]) y = mp[y]; else y = mp[y] = ++n;  
            st.a[x][y] = st.a[y][x] = z;  
        }  
        power();  
        printf("%d", ans.a[mp[s]][mp[t]]);  
        return 0;  
    }