FZU 2173 Nostop（矩阵快速幂+floyd）

题目链接：【FZU 2173】

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

简化一下就是求k个floyd后城市1能否到达城市n

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <string>
using namespace std;
#define ll __int64

const ll inf=1e18;
int t, n, h, k;
struct mat
{
	ll a[55][55];
};
mat floyd(mat a, mat b)
{
	mat s;
	for(int i=1; i<=n; i++)
	{
		for(int j=1; j<=n; j++)
		{
			s.a[i][j] = inf;
		}
	}
	for(int l=1; l<=n; l++)
	{
		for(int i=1; i<=n; i++)
		{
			for(int j=1; j<=n; j++)
			{
				s.a[i][j] = min(s.a[i][j], a.a[i][l]+b.a[l][j]);
			}
		}
	}
	return s;
}
mat mpo(mat g)
{
	mat I=g;
	k--;
	while(k)
	{
		if(k&1) I = floyd(I, g);
		g = floyd(g, g);
		k>>=1;
	}
	return I;
}
int main()
{
	scanf("%d", &t);
	while(t--)
	{
		scanf("%d%d%d", &n, &h, &k);
		int u, v;
		ll c;
		mat dis;
		for(int i=1; i<=n; i++)
		{
			for(int j=1; j<=n; j++)
			{
				dis.a[i][j]=inf;
			}
		}
		for(int i=0; i<h; i++)
		{
			cin>>u>>v>>c;
			dis.a[u][v] = min(dis.a[u][v], c);
		}
		dis = mpo(dis);
		if(dis.a[1][n]==inf) printf("-1\n");
		else printf("%I64d\n", dis.a[1][n]);
	}
	return 0;
}