RosettaCodeData/Task/Perfect-shuffle/PARI-GP/perfect-shuffle.parigp
2023-07-01 13:44:08 -04:00

4 lines
180 B
Text

magic(v)=vector(#v,i,v[if(i%2,1,#v/2)+i\2]);
shuffles_slow(n)=my(v=[1..n],o=v,s=1);while((v=magic(v))!=o,s++);s;
shuffles(n)=znorder(Mod(2,n-1));
vector(5000,n,shuffles_slow(2*n))