poj3613Cow Relays——k边最短路

Posted Zinn

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了poj3613Cow Relays——k边最短路相关的知识,希望对你有一定的参考价值。

题目:http://poj.org/problem?id=3613

题意就是求从起点到终点的一条恰好经过k条边的最短路;

floyd+矩阵快速幂,矩阵中的第i行第j列表示从i到j的最短路,矩阵本身代表一个边数状态;

所以矩阵相乘就是floyd算法,两个矩阵相乘就得到它们所代表的边数相加边数的状态矩阵;

原始矩阵自乘k-1次,过程中取min,就得到答案;

因为只是自乘,所以可以使用快速幂。

代码如下:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
int n,cnt,s,t,k,nd[1005];
int rd()
{
    char dc;
    int x=0,f=1;dc=getchar();
    while(dc<0||dc>9)
    {
        if(dc==-)f=-1;
        dc=getchar();
    }
    while(dc>=0&&dc<=9)
    {
        x=x*10+(dc-0);
        dc=getchar();
    }
    return x*f;
}
struct Matrix{
    int a[205][205];
    Matrix operator * (const Matrix &y)  const
    {
        Matrix x;
        memset(x.a,0x3f,sizeof x.a);
        for(int i=1;i<=cnt;i++)
            for(int j=1;j<=cnt;j++)
                for(int k=1;k<=cnt;k++)
                    x.a[i][j]=min(x.a[i][j],a[i][k]+y.a[k][j]);
        return x;
    }
}sid,ans;
int main()
{
    k=rd();n=rd();s=rd();t=rd();
    memset(sid.a,0x3f,sizeof sid.a);//
    for(int i=1;i<=n;i++)
    {
        int x,y,z;
        z=rd();x=rd();y=rd();
        if(!nd[x])nd[x]=++cnt;
        if(!nd[y])nd[y]=++cnt;
        sid.a[nd[x]][nd[y]]=sid.a[nd[y]][nd[x]]=z;
    }
    k--;//已经有连了一条边的矩阵 
    ans=sid;//同上意义 
    while(k)
    {
//        if(k&1)ans+=sid.a[nd[s]][nd[t]];//不是普通加法
        if(k&1)ans=ans*sid;
        sid=sid*sid;
        k>>=1;
    }
    printf("%d",ans.a[nd[s]][nd[t]]);
    return 0;
}

 

以上是关于poj3613Cow Relays——k边最短路的主要内容,如果未能解决你的问题,请参考以下文章

poj 3613 Cow Relays(矩阵的图论意义)

POJ3613Cow Relays

「POJ3613」Cow Relays

poj 3613 Cow Relays矩阵快速幂+Floyd

[POJ3613] Cow Relays(Floyd+矩阵快速幂)

POJ3613 k边最短路