Генерация перестановок

using System;
using System.Collections.Generic;
 
namespace Problem
{
    class Program
    {
        static void Main()
        {
            int[] perm = {1, 2, 3};
            do { Console.WriteLine(string.Join(" ", perm)); } while (NextPermutation(perm)); //лексикографический перебор
 
            Console.WriteLine("************");
            
            int[] perm2 = {1, 2, 3};
            Console.WriteLine(string.Join(" ", perm2));
            foreach (var p in Permutation(3)) //код Грэя
            {
                for (int i = 0; i < perm2.Length; ++i)
                {
                    Console.Write("{0} ", perm2[p[i]]);
                }
                Console.WriteLine();
            }
            Console.ReadKey();
        }
 
        static public bool NextPermutation(int[] perm)
        {
            int n = perm.Length;
            int k = -1;
            for (int i = 1; i < n; i++)
                if (perm[i - 1] < perm[i])
                    k = i - 1;
            if (k == -1)
            {
                for (int i = 0; i < n; i++)
                    perm[i] = i;
                return false;
            }
 
            int l = k + 1;
            for (int i = l; i < n; i++)
                if (perm[k] < perm[i])
                    l = i;
 
            int t = perm[k];
            perm[k] = perm[l];
            perm[l] = t;
 
            Array.Reverse(perm, k + 1, perm.Length - (k + 1));
 
            return true;
        }
 
        static IEnumerable<int[]> Permutation(int k)
        {
            int[] a = new int[k + 1];
            int[] c = new int[k + 1];
            int[] o = new int[k + 1];
            for (int i = 0; i < o.Length; ++i)
            {
                a[i] = i - 1;
                c[i] = 0;
                o[i] = 1;
            }
            for (int j = k, s = 0; j != 1; )
            {
                int q = c[j] + o[j];
                if (q >= 0 && q != j)
                {
                    int i1 = j - c[j] + s;
                    int i2 = j - q + s;
                    Swap(ref a[i1], ref a[i2]);
                    c[j] = q;
                    int[] t = new int[a.Length - 1];
                    Array.Copy(a, 1, t, 0, t.Length);
                    yield return t;
                    j = k;
                    s = 0;
                    continue;
                }
                if(q >=0)
                    ++s;
                o[j] = -o[j];
                --j;
            }
        }
 
        private static void Swap(ref int a, ref int b)
        {
            int t = a;
            a = b;
            b = t;
        }
    }
}