随机系列生成算法（随机数生成）

1、问题描述

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

2、准备

          首先，让我们先构造如下函数：

[cpp]  view plain copy

1. //随机数生成函数，生成[a,b]之间的一个随机数字  
2. //其中0<=a<b<=RAND_MAX(即2^15-1=32767)  
3. int random(int a,int b)  
4. {  
5.     return a+rand()%(b-a+1);  
6. }  

3、解法一

          先给序列A[0...n-1]赋值，A[i]=i;(0=<i<n)。然后为每个元素A[i]赋一个优先级P[i]，然后依据优先级对序列A进行排序。这里要保证序列P中没有相同的元素，所以P[i]的取值范围可以选[0,(n-1)^3]。这样就尽可能保证P中元素无重复。伪代码如下：

[cpp]  view plain copy

1. n=length[A]  
2. for(i=0;i<n;i++)  
3.    P[i]=random(0,n^3);  
4. sort A,using p as sort keys  

          分析：显然解法一无法准确的确保P中元素无重复，且当n较大时，n的三次方会很大，如当n=10000时，n的三次方已经超出了无符号长整形的表示范围。该算法时间复杂度依赖排序时间，且这个排序不好实现。总之，这不是一个很好的随机序列生成算法。

4、解法二

        先给序列A[0...n-1]赋值，A[i]=i;(0=<i<n)。然后原地排列给定序列。代码如下：

[cpp]  view plain copy

1. //产生随机序列  
2. void randomSequence1(int *a,int n)  
3. {  
4.     int i,r,temp;  
5.     for(i=0;i<n;i++)  
6.         a[i]=i;  
7.     srand(time(NULL));  
8.     for(i=0;i<n;i++){  
9.         r=random(i,n-1);  
10.         //交换  
11.         temp=a[i];  
12.         a[i]=a[r];  
13.         a[r]=temp;  
14.     }    
15. }  

         分析：这个算法的时间复杂度为o（n）。

5、解法三

        这是在网上看的一种高效随机生成序列算法，详细请看原文：http://www.cnblogs.com/eaglet/archive/2011/01/17/1937083.html

        这里给出自己写的代码：

[cpp]  view plain copy

1. //产生随机序列  
2. void randomSequence2(int *a,int n)  
3. {  
4.     int i,end,r;  
5.     int *q=(int *)malloc(sizeof(int)*n);  
6.     for(i=0;i<n;i++)  
7.         q[i]=i;  
8.     srand(time(NULL));  
9.     end=n-1;  
10.     for(i=0;i<n;i++){  
11.         r=random(0,end);  
12.         a[i]=r;  
13.         q[r]=q[end];  
14.         end--;  
15.     }    
16. }  

6、附录

      参考书籍：《算法导论》     解法三原文：http://www.cnblogs.com/eaglet/archive/2011/01/17/1937083.html