vijos1049送给圣诞夜的礼品

这题犯了两个sb错误，写下来，为以后做个警告

一、mul过程中将k作为了循环变量

二、看错了题……

题目中说是数到k行，而我却以为数k遍……

做矩阵乘法，只要记住一句话：置换一定可以写成矩阵的形式！

并且：矩阵满足结合律，不满足交换率！

代码：

 type matrix=array[..,..] of longint;
 var a:array[..] of matrix;
     b,c:matrix;
     i,j,n,m,k,x:longint;
     ans,d:array[..] of longint;
 procedure tiao(x:matrix);
  var i,j:longint;
  begin
  for i:= to n do
   begin
    for j:= to n do write(x[i,j]);
    writeln;
   end;
  end;
 procedure init;
  begin
    readln(n,m,k);
    fillchar(a,sizeof(a),);
    fillchar(b,sizeof(b),);
    fillchar(c,sizeof(c),);
    for i:= to m do
     begin
       for j:= to n do
         begin
           read(x);
           a[i,x,j]:=;
         end;
       readln;
     end;
    for i:= to n do begin b[i,i]:=;c[i,i]:=;end;
  end;
 procedure mul(var x,y,z:matrix);
  var t:matrix;
      i,j,l:longint;
  begin
    fillchar(t,sizeof(t),);
    for i:= to n do
      for j:= to n do
        for l:= to n do
          inc(t[i,j],x[i,l]*y[l,j]);
    z:=t;
  end;
 procedure ksm(cs:longint);
  begin
    while cs<> do
     begin
       if cs and = then mul(b,c,c);
       cs:=cs>>;
       mul(b,b,b);
     end;
  end;
 procedure main;
  begin
    for i:= to m do mul(b,a[i],b);
    ksm(k div m);
    for i:= to k mod m do mul(c,a[i],c);
    for i:= to n do d[i]:=i;
    fillchar(ans,sizeof(ans),);
    for i:= to n do
      for j:= to n do
       inc(ans[i],d[j]*c[j,i]);
    for i:= to n do write(ans[i],' ');
  end;
 begin
   init;
   main;
 end.