poj 3613Cow Relays (矩阵乘法 快速幂)

本题的大意就是问从S 到 T 经过边得个数恰为k的最短路是多少。

解题思路：

    首先看到要过k条边，就想到了矩阵的应用中的一道题：01邻接矩阵A的K次方C=A^K,C[i][j]表示i点到j点正好经过K条边的路径数

这道题与那道题其实差不多。只是C[I][J]变成了i到j的最短路。 所以要进行K次，就用矩阵乘法，但要稍微不同。

注意：矩阵快速幂并不只是用于矩阵的乘法，只要涉及多个矩阵进行运算，就可以使用，运算由我们自己选定。

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define INF ((1<<30)-1)
int n,t,s,e;
struct ss
 {
  int messi[211][211];
  ss(){  
        for(int i=0;i<211;i++)  
            for(int j=0;j<211;j++)messi[i][j]=INF;  
    }  
 };
 ss a,now;
int zan;
int sum[10000];


ss work(ss a1,ss b1)
 {
  ss c;
  for (int i=1;i<=zan;++i)
  for (int j=1;j<=zan;++j)
   for (int z=1;z<=zan;++z)
      {
     c.messi[i][j]=min(c.messi[i][j],a1.messi[i][z]+b1.messi[z][j]);
}
return c;
 }


int main()
{
scanf("%d %d %d %d",&n,&t,&s,&e);
memset(sum,0,sizeof(sum));
zan=0;
for (int i=1;i<=t;++i)
{
int length,x,y;
scanf("%d %d %d",&length,&x,&y);
if (sum[x]==0) sum[x]=++zan; 
if (sum[y]==0) sum[y]=++zan;
a.messi[sum[x]][sum[y]]=length; a.messi[sum[y]][sum[x]]=length;
}
for (int i=1;i<=zan;++i) now.messi[i][i]=0;
int u=n;
while (u>0)
{
if (u%2==1)
{
now=work(now,a);
 };
u=u/2;
a=work(a,a);
}
cout<<now.messi[sum[s]][sum[e]];
}