foj 2173 （floyd+矩阵快速幂 ）

Nostop

Time Limit: 3000MS

Memory Limit: 32768KB

64bit IO Format: %I64d & %I64u

Submit Status

Description

M国有N个城市，H条单向的道路，AekdyCoin从编号为1的城市出发，每经过一条道路要花一个单位的时间。假设他出发的时刻为0，他需要在K时刻到达编号为N的城市。并且，AekdyCoin不会在一个城市停留，每到一个城市他要立刻往下一个城市出发，最后在K时刻时他必须在城市N。虽然AekdyCoin经过任意一条道路的花费的时间都是1，但是每条道路的过路费不一定相同。现给出每条道路的过路费，问AekdyCoin从编号为1的城市出发，在K时刻到达编号为N的城市最小需要花费多少钱？注意AekdyCoin可以经过同一个城市任意多次，包括城市N。

Input

第一行输入一个整数T表示数据组数，接下来输入T组数据。对于每组数据，第一行输入三个整数N，H,K(1<=N<=50,1<=H<=3000,1<=K<=1000000000)，接下来输入H行，每行三个整数u、v、cost（1<=u,v<=n,1<=cost<=1000000），表示从u到v过路费为cost的一条单行道。

Output

对于每组数据输出一行一个整数表示最小花费，若无法在K时刻到达城市N，则输出-1。

Sample Input

15 5 31 2 12 5 11 3 103 4 104 5 10

Sample Output

30

分析：

转载

FOJ 2173 Nostop 从1点到n点恰好走了k次的最短路

题目链接：http://acm.fzu.edu.cn/problem.PHP?pid=2173

思路：

类似于传递闭包的性质

用矩阵mp[i][j] 表示i点到j点 走1次的最短路

--------------

若我们用 mp[i][j] 表示从i点到j点 走了k次的最短路距离

那么我们要通过 矩阵mp 得到 矩阵 ret[u][v] 表示 u->v 走了2*k次的最短路

就是：

mp[u][i] + mp[i][v]; i为任意点（即1-n）

显然我们转换一下上式就是：

ret[u][v] = inf;  
for(int i = 1; i <= n; i++)  
ret[u][v] = min(ret[u][v], mp[u][i]+mp[i][v]);  

然后求出整个的ret矩阵就是：

for(int u = 1; u<=n; u++)  
    for(int v = 1; v<=n; v++){  
        ret[u][v] = inf;  
        for(int i = 1; i <= n; i++)  
            ret[u][v] = min(ret[u][v], mp[u][i]+mp[i][v]);  
}  

显然就是 ret = mp*mp;
然后套个矩阵快速幂：

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <queue>
#define mem(p,k) memset(p,k,sizeof(p));
#define rep(a,b,c) for(int a=b;a<c;a++)
#define pb push_back
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define inf 0x6fffffff
#define ll long long
using namespace std;
int m,n,minn,maxx,x,y,cur;
ll h,le,k;
struct mat{
    ll ret[52][52];
};
void init(mat &A){
       for(int i=0;i<n;i++){
          for(int j=0;j<n;j++)
             A.ret[i][j]=1e17;
       }
       return;
}

mat mul(mat A,mat B){
    mat C;
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            C.ret[i][j]=1e17;
            for(int k=0;k<n;k++){
                C.ret[i][j]=min(C.ret[i][j],A.ret[i][k]+B.ret[k][j]);
            }
        }
    }
    return C;
}
void solve(mat A,int k){
    mat B;
    B=A;
    k--;
    while(k){

        if(k&1){
            B=mul(B,A);
        }
        A=mul(A,A);
        k>>=1;
    }
    if(B.ret[0][n-1]==1e17)cout<<-1<<endl;
    else cout<<B.ret[0][n-1]<<endl;
}
int main(){
    int t;
    cin>>t;
    while(t--){
        //ll le,h,k;
        cin>>n>>h>>k;
        mat A;
        init(A);
        for(int i=0;i<h;i++){
            cin>>x>>y>>le;
            x--,y--;
            A.ret[x][y]=min(A.ret[x][y],le);
        }
        solve(A,k);
    }

    return 0;
}