Mã C / C ++ mã lõi chỉ có một dòng!
static unsigned int gi = 0;
int rand7()
{
return (((rand() % 5 + 1) + (gi++ % 7)) % 7) + 1;
}
//call this seed before rand7
//maybe it's not best seed, if yo have any good idea tell me please
//and thanks JiminP again, he remind me to do this
void srand7()
{
int i, n = time(0);
for (i = 0; i < n % 7; i++)
rand7();
}
Srand7 () là hạt giống của rand7, phải gọi hàm này trước rand7, giống như gọi srand trước rand trong C.
Đây là một thứ rất tốt, vì nó chỉ gọi rand () một lần duy nhất và không có điều lặp, không có thêm ký ức.
Hãy để tôi giải thích nó: xem xét một mảng số nguyên với kích thước 5:
1st get one number from 1 2 3 4 5 by rand5
2nd get one number from 2 3 4 5 6
3rd get one number from 3 4 5 6 7
4th get one number from 4 5 6 7 1
5th get one number from 5 6 7 1 2
5th get one number from 6 7 1 2 3
7th get one number from 7 1 2 3 4
Vì vậy, chúng tôi đã có BẢNG, mỗi một trong số 1-7 xuất hiện 5 lần trong đó và có tất cả 35 số, vì vậy xác suất của mỗi số là 5/35 = 1/7. Và lần sau,
8th get one number from 1 2 3 4 5
9th get one number from 2 3 4 5 6
......
Sau khi đủ thời gian, chúng ta có thể nhận được phân phối đồng đều từ 1-7.
Vì vậy, chúng ta có thể phân bổ một mảng để khôi phục năm phần tử của 1-7 bằng vòng lặp trái-shift và nhận một số từ mảng mỗi lần bằng rand5. Thay vào đó, chúng ta có thể tạo ra tất cả bảy mảng trước đó và sử dụng chúng theo vòng tròn. Mã cũng đơn giản, có nhiều mã ngắn có thể làm điều này.
Nhưng, chúng ta có thể sử dụng các thuộc tính của% hoạt động, do đó, các hàng 1-7 của bảng tương đương với (rand5 + i)% 7, nghĩa là: a = rand ()% 5 + 1 là rand5 trong ngôn ngữ C, b = gi ++ % 7 tạo ra tất cả các hoán vị trong bảng trên và 0 - 6 thay thế 1 - 7 c = (a + b)% 7 + 1, tạo 1 - 7 đồng đều. Cuối cùng, chúng tôi đã nhận được mã này:
(((rand() % 5 + 1) + (gi++ % 7)) % 7) + 1
Nhưng, chúng tôi không thể nhận được 6 và 7 ở lần gọi đầu tiên, vì vậy chúng tôi cần một hạt giống, một số giống như srand cho rand trong C / C ++, để làm xáo trộn hoán vị cho cuộc gọi chính thức đầu tiên.
Đây là mã đầy đủ để thử nghiệm:
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
static unsigned int gi = 0;
//a = rand() % 5 + 1 is rand5 in C language,
//b = gi++ % 7 generates all permutations,
//c = (a + b) % 7 + 1, generates 1 - 7 uniformly.
//Dont forget call srand7 before rand7
int rand7()
{
return (((rand() % 5 + 1) + (gi++ % 7)) % 7) + 1;
}
//call this seed before rand7
//maybe it's not best seed, if yo have any good idea tell me please
//and thanks JiminP again, he remind me to do this
void srand7()
{
int i, n = time(0);
for (i = 0; i < n % 7; i++)
rand7();
}
void main(void)
{
unsigned int result[10] = {0};
int k;
srand((unsigned int)time(0)); //initialize the seed for rand
srand7() //initialize the rand7
for (k = 0; k < 100000; k++)
result[rand7() - 1]++;
for (k = 0; k < 7; k++)
printf("%d : %.05f\n", k + 1, (float)result[k]/100000);
}
