fzu 2173 floyd+矩阵快速幂

#include<stdio.h>

#define inf 1000000000000000

#define N 100

long long tmp[N][N],ma[N][N];

int n;

long long min(long long a,long long b)

{

    return a<b?a:b;

}

void floy(long long  a[N][N])

{

    int i,j,k;

    for(i=1;i<=n;i++)

        for(j=1;j<=n;j++)

            tmp[i][j]=inf;

      for(k=1;k<=n;k++)

        for(i=1;i<=n;i++)

          for(j=1;j<=n;j++)

            tmp[i][j]=min(tmp[i][j],a[i][k]+ma[k][j]);

      for(i=1;i<=n;i++)

        for(j=1;j<=n;j++)

          a[i][j]=tmp[i][j];

}

void f(int k,long long dis[N][N])

{

    k--;

    while(k)

    {

        if(k&1)

            floy(dis);

        floy(ma);

        k=k/2;

    }

}

int main()

{

    int t,m,k,i,j,u,v;

    long long cost;

    long long dis[N][N];

    scanf("%d",&t);

    while(t--)

    {

        scanf("%d%d%d",&n,&m,&k);

        for(i=1;i<=n;i++)

            for(j=1;j<=n;j++)

           {

                ma[i][j]=inf; dis[i][j]=inf;

            }

        for(i=1;i<=m;i++)

        {

            scanf("%d%d%I64d",&u,&v,&cost);

            if(ma[u][v]>cost){  ma[u][v]=cost;

            dis[u][v]=ma[u][v];

            }

        }

        f(k,dis);

        if(dis[1][n]>=inf)  printf("-1\n");

        else   printf("%I64d\n",dis[1][n]);

    }

    return 0;

}