RosettaCodeData/Task/Matrix-arithmetic/PARI-GP/matrix-arithmetic-2.pari

2 lines
73 B
Plaintext

matperm(M)=my(n=#M,t);sum(i=1,n!,t=numtoperm(n,i);prod(j=1,n,M[j,t[j]]))