Some new GP features permcycles permcycles returns the cycle decomposition of a permutation. Fixed points are returned as length-1 cycles. ? permcycles(Vecsmall([2,7,1,8,4,5,9,10,3,6])) %12 = [Vecsmall([1,2,7,9,3]) % ,Vecsmall([4,8,10,6,5])] ? permcycles(Vecsmall([3,1,4,5,9,2,6,8,7])) %13 = [Vecsmall([1,3,4,5,9,7,6,2]),Vecsmall([8])]