37 lines
620 B
C
37 lines
620 B
C
#include <stdlib.h>
|
|
#include <stdint.h>
|
|
#include <sys/types.h>
|
|
|
|
static unsigned int seed=1;
|
|
|
|
/* Knuth's TAOCP section 3.6 */
|
|
#define M ((1U<<31) -1)
|
|
#define A 48271
|
|
#define Q 44488 // M/A
|
|
#define R 3399 // M%A; R < Q !!!
|
|
|
|
// FIXME: ISO C/SuS want a longer period
|
|
int rand_r(unsigned int* seed)
|
|
{ int32_t X;
|
|
|
|
X = *seed;
|
|
X = A*(X%Q) - R * (int32_t) (X/Q);
|
|
if (X < 0)
|
|
X += M;
|
|
|
|
*seed = X;
|
|
return X;
|
|
}
|
|
|
|
int rand(void) {
|
|
return rand_r(&seed);
|
|
}
|
|
|
|
void srand(unsigned int i)
|
|
{
|
|
seed=i;
|
|
}
|
|
|
|
int random(void) __attribute__((alias("rand")));
|
|
void srandom(unsigned int i) __attribute__((alias("srand")));
|