传说中的Stooge排序

using System; class Stooge {     public void sort(int[] data, int i ,int j)     {         if (data[i]>data[j])         {             int t = data[i];              data[i] = data[j];              data[j] = t;          }         if (i+1 >= j) return;         int k = (j-i+1)/3;          sort(data, i, j-k);          sort(data, i+k, j);          sort(data, i, j-k);      } } class Program {     public static void Main()     {          Random rnd = new Random();         int[] data = new int[10];         for (int i=0; i<10; i++)         {              data[i] = rnd.Next(100);              Console.Write("{0}\t",data[i]);          }          Console.WriteLine();          Stooge s = new Stooge();          s.sort(data, 0, 9);         for (int i=0; i<10; i++)              Console.Write("{0}\t",data[i]);          Console.Read();      } } 算法导论里面曾介绍过这个由Howard、Fine等教授提出的声称很厉害的排序算算法，其实说它厉害并不是因为它很快，实际上它比插入排序还要慢，之所以厉害是因为一般的脑子想想是无法想象出它的正确性的，但用循环不变式可以帮助理解算法的正确性：证明如下：  循环不变式：在每次 STOOGE-SORT(A, i, j) 返回时，数组 A[i..j] 是有序的。  初始化：当 i+1 ≥ j 时，不再进行递归，函数会立即返回。由于执行了第1、2行，可保证 A[i] ≤ A[j]。因为 i+1=j，所以数组A[i..j]有序。循环不变式成（）立。  保持：假设 STOOGE-SORT(A, i, j) 内部对STOOGE-SORT的所有递归调用都满足循环不变式，即STOOGE-SORT(A, i, j-k) 可使 A[i..j-k] 有序，STOOGE-SORT(A, i+k, j) 可使 A[i+k..j] 有序。  ∵ k =(j-i+1) div 3  ∴ j-i+1 ≥ 3k  ∴ (j-i+1)-2k ≥ k  ∵ A[i..j-k] 与 A[i+k..j] 重叠的部分为 A[i+k..j-k]，共(j-k)-(i+k)+1 = (j-k+1) - 2k ≥ k 个元素  而 A[i..j-k] 与 A[i+k..j] 不重叠的部分为 A[j-k+1..j] 共 j-(j-k+1)+1 = k 个元素。即，A[i..j-k] 与 A[i+k..j] 重叠部分的元素个数大于等于不重叠部分的元素个数。  ∴ 在执行了第 6、7 行的STOOGE-SORT(A, i, j-k) 和 STOOGE-SORT(A, i+k, j)，分别使A[i..j-k] 和 A[i+k..j]有序后，可保证 A[j-k+1..j] 中的元素是 A[i..j] 中最大的 k 个元素，且是有序的。  ∵ 在执行了第8行的STOOGE-SORT(A, i, j-k)后，可保证 A[i..j-k] 有序。  ∴ 综上，可保证在 STOOGE-SORT(A, i, j) 返回后，A[i..j] 是有序的。循环不变式成立。  终止：在最外层的 STOOGE-SORT(A, 1, n) 返回后，可使数组 A 有序，算法是正确的。  上面说了通废话，接下来切入正题，分析一下它的实践性能很容易写出递归式：  T(n)=Θ(1)当n=1  T(n)=3T(2n/3)+Θ(1) 当n>1  显然根据主定理可以得到它的复杂性为Θ(N^(log(3/2,3)))，大约就是Θ(n^2.7)，将近Θ(n^3)的复杂度，显然比插入排序还要慢！