状态:=选择(proc(x)(x[1]=x[2]或x[5]=1)和(x[2]=x[3]或x[6]=1)和(x[3]=x[4]或x[7]=1)结束proc,[seq(seq(seq(seq(seq(seq(seq([a,b,c,d,e,f,g],g=0..1),f=0..1),e=0..1),d=0..1),c=0..1),b=0..1),a=0..1)]):
T: =矩阵(54,54,proc(i,j)局部k;
如果add(states[j,k]-states[i,k],k=1..4)>0,则返回0fi;
如果状态[j,5]>状态[i,5]或状态[j、6]>状态[i,6]或状态[j,7]>状态[1,7],则返回0 fi;
如果状态[i,1]>=状态[i,2]和状态[j,5]<>状态[i,5],则返回0fi;
如果状态[i,2]>=状态[i、3]和状态[j,6]<>状态[i和6],则返回0 fi;
如果状态[i,3]>=状态[i、4]和状态[j,7]<>状态[i和7],则返回0 fi;
1
结束进程):
U: =矢量(54,1):
E[0]:=向量(54):E[0][1]:=1:
对于从1到25的k,做E[k]:=T。E[k-1]od(E[k-1]od):
seq(U^%T.E[j],j=1..25); #罗伯特·伊斯雷尔2019年9月8日