restart:with(LinearAlgebra): L:=407: > L2:=convert(L,binary): > N:=length(L2): > T:=2^N-L: > TB:=convert(T,binary): > k:=length(TB): > T2:=convert(T,base,2): > T2V:=Vector[row]([T2,0$(N-k)]): > T2VC:=Vector[row]([1$N]): > T2V1:=VectorScalarMultiply(T2V,-1): > T2VA1:=VectorAdd(T2VC,T2V1): > ADDT:=Vector[row]([1,T2VA1]): > B_N:=Array(Hessenberg[lower],1..N+1,1..N+1,(a,b)->1): > C:=convert(B_N, Matrix): > Permanent(C): > delB:=DeleteRow(C,N+1..N+1): > B(L):=Matrix(N+1,N+1,[[delB],[ADDT]]): > Permanent(B(L));

Please Wait...