蛙蛙推荐：算法练习：最大间隙问题

问题描述：
最大间隙问题：给定的n个实数x1,x2...,xn，求这N个数在实轴上相邻两个数之间最大差值。假设对任何实数的下去整耗时是O(1),设计最大间隙问题的线性时间算法。
数据输入：
输入数据由文件名为input.txt的文本文件提供，文件的第一行有一个正整数N，接下来的一行有N个实数x1,x2...,xn
数据输出：
程序运行结束时，将找到的最大间隙输出到output.txt里

这是《算法设计与实验题解》里的一道题，原文中的算法分析是用c写的，我转换成了c#语言，如下

using  System;
using  System.Diagnostics;
using  System.Collections;
using  System.IO;
using  System.Runtime.InteropServices;

namespace  ZuiDaJianXi {
     class  Program {
         static   int  MyFloor( double  d) {
            cast_struct c  =   new  cast_struct { d  =  d  -   0.499999999999   +   6755399441055744.0  };
             return  c.i;
        }

        [StructLayout(LayoutKind.Explicit)]
         struct  cast_struct {
            [FieldOffset( 0 )]
             public   double  d;
            [FieldOffset( 0 )]
             public   int  i;
        }
         static   unsafe   int  MyFloor2( double  dval) {
            dval  =  dval  -   0.499999999999 ;
             const   double  magic  =   6755399441055744.0 ;
            dval  +=  magic;
             return   * ( int * ) & dval;
        }
         static   void  Main( string [] args) {
            Random rnd  =   new  Random();
             double [] x  =   new   double [ 10000000 ];
             for  ( int  i  =   0 ; i  <  x.Length; i ++ )
                x[i]  =  rnd.NextDouble()  *   100000000 ;

            Stopwatch stop  =  Stopwatch.StartNew();
             // Array.Sort(x);
             // x = new double[] {1,2,3,4,5};
             // Array.ForEach(x, d => Console.WriteLine("{0}",d));


            Console.WriteLine( " 最大间隙为:{0},time={1} " ,
                maxgap(x.Length, x), stop.ElapsedTicks);
            stop.Reset();
            stop.Start();
            Console.WriteLine( " 最大间隙为:{0},time={1} " ,
                maxgap2(x), stop.ElapsedTicks);
            Console.ReadKey();
        }
         static   double  maxgap2( double [] x) {
            Array.Sort(x);
             double  tmp  =   0 , left  =  x[ 0 ];
             for  ( int  i  =   1 ; i  <  x.Length; i ++ ) {
                 double  thisgap  =  x[i]  -  left;
                left  =  x[i];
                 if  (thisgap  >  tmp) tmp  =  thisgap;
            }
             return  tmp;
        }
         static   double  maxgap( int  n,  double [] x) {
             double  minx  =  x[mini(n, x)], maxx  =  x[maxi(n, x)];

             // 用n-2个等间距的点分割区间[minx,maxx]
             // 产生n-1个桶，每个桶i中用high[i]和low[i]
             // 分别存储桶i的数中的最大数和最小数
             int [] count  =   new   int [n  -   1 ];
             double [] low  =   new   double [n  -   1 ];
             double [] high  =   new   double [n  -   1 ];

             // 桶初始化
             for  ( int  i  =   0 ; i  <  n  -   1 ; i ++ ) {
                count[i]  =   0 ;
                low[i]  =  maxx;
                high[i]  =  minx;
            }

            Stopwatch stop  =   new  Stopwatch();
            stop.Start();
             // 将n个数放进n-1个桶中
             for  ( int  i  =   0 ; i  <  n; i ++ ) {
                 // 底下这句不知道是基于什么数学原理能把一个数
                 // 确定的放入它自己所属的桶里？？
                 // int bucket = MyFloor(((n - 2) * (x[i] - minx) / (maxx - minx)));
                 int  bucket  =  MyFloor2(((n  -   2 )  *  (x[i]  -  minx)  /  (maxx  -  minx)));
                 // int bucket = (int)((n - 2) * (x[i] - minx) / (maxx - minx));
                count[bucket] ++ ;
                 if  (x[i]  <  low[bucket]) low[bucket]  =  x[i];
                 if  (x[i]  >  high[bucket]) high[bucket]  =  x[i];
            }
            Console.WriteLine( " 入桶时间：{0} " ,stop.ElapsedTicks);
             // 此时除了maxx和minx外的n-2个数放在n-1个桶中
             // 由鸽舍原理可知，至少有一个桶是空的
             // 这意味着最大间隙不会出现在同一个桶的两个数之间
             // 对每个桶做一次线性扫描即可找出最大间隙
             double  tmp  =   0 , left  =  high[ 0 ];
             for  ( int  i  =   0 ; i  <  n  -   1 ; i ++ ) {
                 if  (count[i]  !=   0 ) {
                     double  thisgap  =  low[i]  -  left;
                     if  (thisgap  >  tmp) {
                        tmp  =  thisgap;
                        Console.WriteLine( " left = {0},thisgap = {1} " , left, tmp);
                    }
                    left  =  high[i];
                }
            }

             return  tmp;
        }

         static   int  mini( int  n,  double [] x) {
             double  temp  =  x[ 0 ];
             int  k  =   0 ;
             for  ( int  i  =   0 ; i  <  n; i ++ ) {
                 if  (x[i]  <  temp) {
                    temp  =  x[i];
                    k  =  i;
                }
            }
             return  k;
        }
         static   int  maxi( int  n,  double [] x) {
             double  temp  =  x[ 0 ];
             int  k  =   0 ;
             for  ( int  i  =   0 ; i  <  n; i ++ ) {
                 if  (x[i]  >  temp) {
                    temp  =  x[i];
                    k  =  i;
                }
            }
             return  k;
        }
    }
}

经测试，在release模式下，输入数据在100w时，maxgap只比maxgap2快1秒，虽然只是快1秒，甚至debug模式下还要比maxgap2还慢，但算法的复杂度比maxgap2要简单很多。后来查了下原因，主要是因为c#的下取整比较耗时，因此造成入桶的时间占总体耗时的百分之八九十，后来在网上找了一些浮点截断成整形的快速算法，发现和c#自带的int强转double的性能差不了多少，汗了，谁要有更好的Floor算法给提供下。

参考链接：
float与double的范围和精度
http://www.cnblogs.com/tekson/archive/2009/07/16/1524604.html
浮点数到整数的快速转换
http://www.gamergroup.cn/html/38/n-4338.html
代码优化－之－优化浮点数取整
http://blog.csdn.net/housisong/archive/2007/05/19/1616026.aspx
其中MyFloor是脑袋提供的，原题出自网上下载的一个《算法设计与实验题解.pdf》文件