邵珠庆の日记 生命只有一次，你可以用它来做很多伟大的事情–Make the world a little better and easier

北京市小汽车摇号程序的反编译、算法及存在的问题浅析-不重复随机数列生成算法

给定一个正整数n，需要输出一个长度为n的数组，数组元素是随机数，范围为0 – n-1，且元素不能重复。比如 n = 3 时，需要获取一个长度为3的数组，元素范围为0-2，

比如 0,2,1。

这个问题的通常解决方案就是设计一个 hashtable ，然后循环获取随机数，再到 hashtable 中找，如果hashtable 中没有这个数，则输出。下面给出这种算法的代码

        public static int[] GetRandomSequence0(int total)

        {

            int[] hashtable = new int[total];

            int[] output = new int[total];

            Random random = new Random();

            for (int i = 0; i < total; i++)

            {

                int num = random.Next(0, total);

                while (hashtable[num] > 0)

                {

                    num = random.Next(0, total);

                }

                output[i] = num;

                hashtable[num] = 1;

            }

            return output;

        }

代码很简单，从 0 到 total - 1 循环获取随机数，再去hashtable 中尝试匹配，如果这个数在hashtable中不存在，则输出，并把这个数在hashtable 中置1，否则循环尝试获取随机数，直到找到一个不在hashtable 中的数为止。这个算法的问题在于需要不断尝试获取随机数，在hashtable 接近满时，这个尝试失败的概率会越来越高。

那么有没有什么算法，不需要这样反复尝试吗？答案是肯定的。

如上图所示，我们设计一个顺序的数组，假设n = 4

第一轮，我们取 0 – 3 之间的随机数，假设为2，这时，我们把数组位置为2的数取出来输出，并把这个数从数组中删除，这时这个数组变成了

第二轮，我们再取 0-2 之间的随机数，假设为1，并把这个位置的数输出，同时把这个数从数组删除，以此类推，直到这个数组的长度为0。这时我们就可以得到一个随机的不重复的序列。

这个算法的好处是不需要用一个hashtable 来存储已获取的数字，不需要反复尝试。算法代码如下：

        public static int[] GetRandomSequence1(int total)

        {

            List<int> input = new List<int>();

            for (int i = 0; i < total; i++)

            {

                input.Add(i);

            }

            List<int> output = new List<int>();

            Random random = new Random();

            int end = total;

            for (int i = 0; i < total; i++)

            {

                int num = random.Next(0, end);

                output.Add(input[num]);

                input.RemoveAt(num);

                end--;

            }

            return output.ToArray();

        }

这个算法把两个循环改成了一个循环，算法复杂度大大降低了，按说速度应该比第一个算法要快才对，然而现实往往超出我们的想象，当total = 100000 时，测试下来，第一个算法用时 44ms, 第二个用时 1038 ms ，慢了很多！这是为什么呢？问题的关键就在这个 input.RemoveAt 上了，我们知道如果要删除一个数组元素，我们需要把这个数组元素后面的所有元素都向前移动1，这个移动操作是非常耗时的，这个算法慢就慢在这里。到这里，可能有人要说了，那我们不用数组，用链表，那删除不就很快了吗？没错，链表是能解决删除元素的效率问题，但查找的速度又大大降低了，无法像数组那样根据数组元素下标直接定位到元素。所以用链表也是不行的。到这里似乎我们已经走到了死胡同，难道我们只能用hashtable  反复尝试来做吗？在看下面内容之前，请各位读者先思考5分钟。

…… 思考5分钟

算法就像一层窗户纸，隔着窗户纸，你永远无法知道里面是什么，一旦捅穿，又觉得非常简单。这个算法对于我，只用了2分钟时间想出来，因为我经常实现算法，脑子里有一些模式，如果你的大脑还没有完成这种经验的积累，也许你要花比我长很多的时间来考虑这个问题，也许永远也找不到捅穿它的方法。不过不要紧，我把这个方法公布出来，有了这个方法，你只需轻轻一动，一个完全不同的世界便出现在你的眼前。原来就这么简单……。

还是上面那个例子，假设 n = 4

第一轮，我们随机获得2时，我们不将 2 从数组中移除，而是将数组的最后一个元素移动到2的位置

这时数组变成了

第二轮我们对 0-2 取随机数，这时数组可用的最后一个元素位置已经变成了2，而不是3。假设这时取到随机数为1

我们再把下标为2 的元素移动到下标1，这时数组变成了

以此类推，直到取出n个元素为止。

这个算法的优点是不需要用一个hashtable 来存储已获取的数字，不需要反复尝试，也不用像上一个算法那样删除数组元素，要做的只是每次把数组有效位置的最后一个元素移动到当前位置就可以了，这样算法的复杂度就降低为 O(n) ，速度大大提高。

经测试，在 n= 100000 时，这个算法的用时仅为7ms。

下面给出这个算法的实现代码

        /// <summary>

        /// Designed by eaglet

        /// </summary>

        /// <param name="total"></param>

        /// <returns></returns>

        public static int[] GetRandomSequence2(int total)

        {

            int[] sequence = new int[total];

            int[] output = new int[total];

            for (int i = 0; i < total; i++)

            {

                sequence[i] = i;

            }

            Random random = new Random();

            int end = total - 1;

            for (int i = 0; i < total; i++)

            {

                int num = random.Next(0, end + 1);

                output[i] = sequence[num];

                sequence[num] = sequence[end];

                end--;

            }

            return output;

        }

下面是n 等于1万，10万和100万时的测试数据，时间单位为毫秒。从测试数据看GetRandomSequence2的用时和n基本成正比，线性增长的，这个和理论上的算法复杂度O(n)也是一致的，另外两个算法则随着n的增大，用时超过了线性增长。在1百万时，我的算法比用hashtable的算法要快10倍以上。

现在摇号的程序及数据都可以在官网中查看，目的就是通过信息的透明度来堵住那些说摇号系统有猫腻之类的传言，我做为一个程序员，下意识的就想看看到底是否有猫腻。

通过下载程序，导入数据，将种子数写入后，的确没有什么猫腻，但还是不死心，想研究一下算法，结果发现了惊人的事情。
此算法使用的是伪随机，简单点将就是当种子数一定，摇号数据一定，每次随机结果都是一样的。举个例子，借用某申请编码2245102443992，摇号基数序号500015，取种子数范围在100000至120000之间能中签的种子数，结果显示共132个种子数，其概率为0.66%。随机找了一个种子数100575进行计算，得出结果在左侧表中，查询2245102443992是否中签，结果显示中签，在第15916行摇中；所以大家都应该知道怎么回事了吧，只要摇号池数据确定，查询出摇号基数序号，就能算出中签的种子数都有哪些，所以公布的数据和程序又有何用？关键在于种子数的确定。