签名算法Rainbow的magma程序

 

//加扰动
//生成密钥的过程不变，签名的时候做相应改变
q:=2^8;
G<a>:=GF(q);
layer:=4;
m:=0;
o:=[6,8,9,1];
v:=[];
for i in [1..layer] do
 if i eq 1 then 
  v[i]:=8;
 else
  v[i]:=v[i-1]+o[i-1];
 end if;
 m:=m+o[i];
end for;
v[layer+1]:=v[layer]+o[layer];
n:=v[layer+1];
//Generate a random invertible matrix B
genRandInvMat:=function(N)
 repeat B:=Matrix(N,N,[Random(G):i in [1..N*N]]);
 until(Rank(B) eq N);
 return(B);
end function;//
procedure genMatrix1(r,k,~M)
 for j:=1 to k do
  L:=[Random(G):i in[1..Floor(r*(r+1)/2)]];
  Z:=ZeroMatrix(G,n,n);

  Mj:=InsertBlock(Z,UpperTriangularMatrix(L),1,1);
  Z:=ZeroMatrix(G,k,k);
  Mj:=InsertBlock(Mj,Z,r-k+1,r-k+1);
                N:=[Random(G):i in[1..Floor(o[layer]*(o[layer]+1)/2)]];
                Mj:=InsertBlock(Mj,UpperTriangularMatrix(N),v[layer+1]-o[layer]+1,v[layer+1]-o[layer]+1);
  M:=M cat [Mj];
 end for;
end procedure;
procedure genMatrix(r,k,~M)
 for j:=1 to k do
  L:=[Random(G):i in[1..Floor(r*(r+1)/2)]];
  Z:=ZeroMatrix(G,n,n);
  Mj:=InsertBlock(Z,UpperTriangularMatrix(L),1,1);
  Z:=ZeroMatrix(G,k,k);
  Mj:=InsertBlock(Mj,Z,r-k+1,r-k+1);
  M:=M cat [Mj];
 end for;
end procedure;
//生成各层私钥矩阵
Sys:=[];
for i in [1..layer] do
 if (i eq layer) or (i eq layer-1) then
   genMatrix1(v[i+1],o[i],~Sys);
 else
  genMatrix(v[i+1],o[i],~Sys);
 end if;
end for;