vijos1049送给圣诞夜的礼品

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

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

二、看错了题……

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

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

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

代码：

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