{
vectoperm(v)=matrix(length(v),length(v),j,k,if(k==v[j],1,0))
}
{
ma(a,b,v,w,h,q)=if(q,,q=length(v));
a11=matrix(q,q,j,k,if(j==k,1+b,0));
a12=matrix(q,q,j,k,if(j==k,-1-a,0));
a21=vectoperm(v)*matrix(q,q,j,k,if(j==k,-a*(1+b),0));
a22=vectoperm(w)*matrix(q,q,j,k,if(j==k,b*(1+a),0));
m=matrix(2*q,2*q,j,k,if(j<=q,if(k<=q,a11[j,k],a12[j,k-q]),
if(k<=q,a21[j-q,k],a22[j-q,k-q])))+(a-b)*vectoperm(h);m
}
{
ls(a,c,x,v,p,q)=if(q,,q=length(v));
m=matrix(p+q,p+q,j,k,if(k<p+1,if(j<p+1,fu(j,k,a,x,p),(x[j-p]+(1+a)*(k-1))%p),
if(j<p+1,(-a*x[v[k-p]]+(1+a)*(j-1))%p,((j-1-p)+c*(k-1-p))%q+p)));m
}
{
fu(j,k,a,x,p)=e=((j-1)+a*(k-1))%p;for(n=1,q,if((j-k)%p == x[n],e=p-1+n));e
}
{
test(z,zv,r,k)=zv=[];k=matsize(z)[1];for(n=1,k,zv=concat(zv,z[n,]));
if(vecsort(vector(k^2,j,zv[j][1]+k*zv[j][2]))==vector(k^2,j,j-1),r=1,r=0);r
}