fzuoj 2173(矩阵快速幂)

思路：用邻接矩阵存储图，然后矩阵的k次方即为答案。只需要修改矩阵乘法c[i][j] = min(c[i][j], a[i][k] + b[k][j])即可。并不难写关键是思路。

代码如下：

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;

const int LEN = 51;
typedef long long ll;
const ll LINF = 0x3f3f3f3f3f3f3f3fLL;
int n, h, k;


void debug(ll Mix[][LEN]){
    for(int i=0; i<n; i++){
        for(int j=0; j<n; j++){
            if(Mix[i][j] != LINF)cout << Mix[i][j] << ' ' ;
            else cout << "- "; 
        }
        cout << endl;
    }cout << endl;
}

void Mul(ll a[][LEN], ll b[][LEN]){
    ll c[LEN][LEN];
    memset(c, 0x3f, sizeof c);
    for(int i=0; i<n; i++){
        for(int j=0; j<n; j++){
            for(int k=0; k<n; k++){
                if(a[i][k] == LINF || b[k][j] == LINF) continue;
                c[i][j] = min(c[i][j], a[i][k] + b[k][j]);
            }
        }
    }
//    debug(a); debug(b);debug(c);
    for(int i=0; i<n; i++){
        for(int j=0; j<n; j++){
            a[i][j] = c[i][j];
        }
    }
}

void Mksm(ll t[][LEN], int k){
    ll Mix[LEN][LEN];
    for(int i=0; i<n; i++)
        for(int j=0; j<n; j++)
            Mix[i][j] = (i==j ? 0 : LINF);
    while(k){
        if(k & 1) Mul(Mix, t);
        Mul(t, t);
        k >>= 1;
    }
    for(int i=0; i<n; i++)
        for(int j=0; j<n; j++)
            t[i][j] = Mix[i][j];
}


int main()
{
   // freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    
    int T, a, b, val;
    ll Mix[LEN][LEN];
    cin >> T;
    while(T--){
        memset(Mix, 0x3f, sizeof Mix);
        cin >> n >> h >> k;
        for(int i=0; i<h; i++){
            cin >> a >> b >> val;
            a--, b--;
            Mix[a][b] = val;        
        }
        Mksm(Mix, k);
        if(Mix[0][n-1] != LINF) cout << Mix[0][n-1] << endl;
        else cout << -1 << endl;
    }
    return 0;
}