解决方法
Here是GNU C库的random()函数中使用的算法的详细描述.
基本上:
- #include <stdio.h>
- #define MAX 1000
- #define seed 1
- main() {
- int r[MAX];
- int i;
- r[0] = seed;
- for (i=1; i<31; i++) {
- r[i] = (16807LL * r[i-1]) % 2147483647;
- if (r[i] < 0) {
- r[i] += 2147483647;
- }
- }
- for (i=31; i<34; i++) {
- r[i] = r[i-31];
- }
- for (i=34; i<344; i++) {
- r[i] = r[i-31] + r[i-3];
- }
- for (i=344; i<MAX; i++) {
- r[i] = r[i-31] + r[i-3];
- printf("%d\n",((unsigned int)r[i]) >> 1);
- }
- }