我创建了一个使用 rand() 在 C 中生成重复数字的程序。
然而,重复的数字不遵循 Central Limit Theorem
任何人都可以解决这个 rand() 错误问题,或者除了使用 rand() C 库来生成更好的随机数之外,还有其他选择吗?
这是屏幕截图:
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <time.h>
#include <windows.h>
#define TOTAL_THROW 10000000
typedef enum _COINTOSS {
UNDEFINED = 0,
HEAD = 1,
TAIL = 2
} COINTOSS;
COINTOSS toss_coin () {
int x = rand() % 2;
if (x == 0) return HEAD;
else if (x == 1) return TAIL;
}
void main () {
int x, i, j, v1 = 0, v2 = 200, total = 0;
int head_range[25] = {0};
int tail_range[25] = {0};
int no_range = 0;
int count = 0;
int repeated = 0;
COINTOSS previos_toss = UNDEFINED;
COINTOSS current_toss;
srand(time(NULL));
for (i=0; i<TOTAL_THROW; i++) {
current_toss = toss_coin(); // current toss
if (previos_toss == current_toss) {
count++;
} else {
if (current_toss == HEAD) head_range[count] += 1;
else if (current_toss == TAIL) tail_range[count] += 1;
previos_toss = current_toss;
count = 0;
}
}
for (i=24; i>=0; i--) {
printf("+%d = %d\n", i+1, head_range[i]);
}
puts("________________\n");
for (i=0; i<25; i++) {
printf("-%d = %d\n", i+1, tail_range[i]);
}
printf("\nTOTAL_THROW: %d\n", TOTAL_THROW);
printf("\nPress [ENTER] to exit. . .");
getchar();
}
最佳答案
您的问题是使用模数使您的随机数进入所需的范围,它使用较低的位(这是一个经典的问题):
int x = rand() % 2;
rand()
的低位(a linear congruential generator (LCG))不像高位那样随机。这适用于所有 LCG,无论库或语言如何。对于 [0..N) 的范围,您应该执行以下操作(使用高位):
int r = rand() / ( RAND_MAX / N + 1 );
关于c - rand() 不遵循高斯分布和中心极限定理,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/26008137/