的确是伪随机~

最新推荐文章于 2023-06-25 09:59:30 发布

午后红茶馆

最新推荐文章于 2023-06-25 09:59:30 发布

阅读量560

点赞数 1

分类专栏： C++基础知识

本文链接：https://blog.csdn.net/u013351484/article/details/50603539

版权

C++基础知识专栏收录该内容

3 篇文章 0 订阅

订阅专栏

C语言的 <stdlib.h> 里 rand() 函数可以产生 0 ~ RAND_MAX （包括）之间的随机数，通常是经过 srand() 函数进行初始化后再使用。

RAND_MAX 的值默认是 32767 即 2^15-1 ，在 <stdlib.h> 里的定义是 #define RAND_MAX 0x7fff;

所以如果需要求 [a, b] （显然要有 b <= RAND_MAX）之间的随机数就可以这样：

value = rand() % (b - a + 1) + a;

然而这样做会出现一个小问题，虽然影响不大，但还是值得谈谈。

举个例子，假如需要产生 [1, 6] 之间的随机数，用上面的方法计算就是

ans = rand() % 6 + 1;

即把一个随机数对 6 取模，产生一个 0~5 的值，将这个值加上 1 并返回。但是，如果随机数生成函数所返回的最大值是 32767，那么这些值就不是“概率相等”。从 0~32765 返回的值所产生的 0~5 之间各个值的概率相等，都是 5461/32768。但是，最后两个值，32766 和 32767 的返回值将分别是 0 和 1，这使它们的出现概率有所增加（是 5462/32768）。由于我们需要的答案的范围很窄，所以这个差别是非常小的。

如果这个函数试图产生一个范围在 [1, 30000] 之间的随机数时，那么前 2768 个值的出现概率将是后面那些值的两倍！做个实验，我们要产生 10000000 个 [1, 30000] 之间的随机数；按照之前的说法，1~2768 之间的数出现的概率是 2/32768，而 2769~30000 之间的数出现的概率是 1/32768；因此 1~2768 之间的数出现的次数约为 10000000 * (2 / 32768) = 610 次，而 2768~30000 之间的数出现的次数约为 10000000 * (1 / 32768) = 305 次。

源代码：

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <string.h>
#define MAXN   30000   /* 随机数范围 */
#define N      10000000  /* 试验次数 */
int count[MAXN+2];

/* 设置随机种子 */
void Seed()
{
	static int is_seed = 0;
	if(!is_seed){
		srand((unsigned int)time(NULL));
		is_seed = 1;
	}
}

/* 计算 beg~end 的随机数*/
int Cal_rand1(int beg, int end)
{
	return rand() % (end - beg + 1) + beg;
}

int main()
{
	int i;
	freopen("out.txt", "w", stdout);
	Seed();
	for(i = 1; i <= N; i++){
		int num = Cal_rand1(1, MAXN);
		count[num]++;
	}
	for(i = 1; i <= MAXN; i++)
		printf("%d: %d\n", i, count[i]);
	return 0;
}

结果为：（显示一部分）

的确如此！那如果需要每个数出现的概率必须相同怎么办？有一个简单的方法：

比如要产生的最大随机数是 b，则令 t = ((RAND_MAX + 1) % b) * b - 1，那只要 rand() 函数一旦出现大于 t 的数就重新产生随机数（其实就是把多出来的那部分去掉~），直到产生的数不大于 t。

修改后的程序：

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <string.h>
#define MAXN   30000   /* 随机数范围 */
#define N      10000000  /* 试验次数 */
#define MAX_OK_RAND \
	(int)((((long)RAND_MAX +1) / MAXN) * MAXN - 1) /* 允许产生的随机数最大值 */
int count[MAXN+2];

/* 设置随机种子 */
void Seed()
{
	static int is_seed = 0;
	if(!is_seed){
		srand((unsigned int)time(NULL));
		is_seed = 1;
	}
}

/* 计算 beg~end 的随机数*/
int Cal_rand2(int beg, int end)
{
	int value;
	do{
		value = rand();
	}while(value > MAX_OK_RAND);
	return value % (end - beg + 1) + beg;
}

int main()
{
	int i;
	freopen("out.txt", "w", stdout);
	Seed();
	for(i = 1; i <= N; i++){
		int num = Cal_rand2(1, MAXN);
		count[num]++;
	}
	for(i = 1; i <= MAXN; i++)
		printf("%d: %d\n", i, count[i]);
	return 0;
}

理论上产生 10000000 个 [1, 30000] 之间的随机数，每个数出现的次数约为 10000000 * (1 / 30000) = 333 次，结果为：（显示一部分）