奶牛接力 矩阵乘法

题意/Description：

       为增强体质，约翰决定举办一场奶牛接力跑比赛。比赛现场有一些接力位置，这些位置间有T条路连接，第i条路的长度为Li。 有N头奶牛需要参加比赛，领头的奶牛从位置S出发，她会按照你的指示沿着一条路跑到下个位置，把接力棒交给等在那里的下一头奶牛，就休息去了。每头奶牛重复这个过程这条接力路线的终点必须在位置E上。



       奶牛数量较多，允许一些奶牛等候在同一位置，一条路也可以供多头奶牛奔跑。 奶牛们对接力跑兴趣不大，拜托你敷衍地设计一条总长度最短的路线。所以请设计一条由N段路组成的，起点在S，终点在E上的最短路线吧。

 
读入/Input：



       第一行：四个用空格分开的整数：表示N，T，S和E，2 ≤ N ≤ 106，2 ≤ T ≤ 100，1
≤ S, E ≤ 1

       第二行到T + 1行： 第i + 1行首先有一个正整数Li，表示第i条路的长度， 1 ≤ Li ≤ 1，其次是Ui和Vi，表示第i条路连接的两个位置，1 ≤ Ui, Vi ≤ 1

输出/Output：

       第一行：单个整数，表示起点为S，终点为E，且恰好经过N段路的最短路线长度



题解/solution：


      题目似乎十分简单，用一个邻接矩阵A就OK，T也不大，但N却十分吓人

，所以先离散化(1—1000)，再建造。


       这题求的是恰好经过N段路的最短路线长度，即矩阵A[i,j]表示从点i到点j的最小值,显然不能和往常那样写c[i,j]=c[i,j]+a[i,k]*b[k,j]了。所以在矩阵乘法那改一下，c[i,j]:=min(c[i,j],a[i,k]+b[k,j])且(0<a[i,j],b[k,j])即可。

       最后在程序中注意一下矩阵的值的大小，否则很容易出错。


       

代码/Code：

type
arr=array [0..101,0..101] of longint;
var
n,t,q,z,m:longint;
d,c:arr;
o:array [0..1001] of longint;
procedure mi(a,b:arr);
var
i,j,k:longint;
begin
for i:=1 to m do
for j:=1 to m do
c[i,j]:=maxlongint div 3;
for i:=1 to m do
for j:=1 to m do
for k:=1 to m do
if (a[i,k]<maxlongint) and (b[k,j]<maxlongint) then
if a[i,k]+b[k,j]<c[i,j] then
c[i,j]:=a[i,k]+b[k,j];
end;

procedure main(n:longint);
begin
if n<=1 then exit;
main(n div 2);
mi(c,c);
if odd(n) then mi(c,d);
end;

procedure init;
var
x,y,t,l,i,j:longint;
begin
readln(n,t,q,z);
m:=0;
for i:=1 to t do
begin
readln(l,x,y);
if o[x]=0 then
begin
inc(m);
o[x]:=m;
end;
if o[y]=0 then
begin
inc(m);
o[y]:=m;
end;
d[o[x],o[y]]:=l;
d[o[y],o[x]]:=l;
end;
for i:=1 to m do
for j:=1 to m do
if d[i,j]=0 then d[i,j]:=maxlongint;
c:=d;
end;

begin
init;
main(n);
write(c[o[q],o[z]]);
end.