// псевдо код с++
unsinged poo[lN]; // кольцевой буфер
usigned head, tail; // при старте == 0

void push_salt(unsigned salt)
{
    unsigned value = pool[head] + salt;
    pool[head] = value + (value < salt); // кольцевое сложение с переносом
    head = (head + 1) % N;
}

unsigned pop_salt()
{
    // y ~~ F(X**N, X**23 + X**17);
    unsigned result = rotate(pool[tail], 3);
    result += pool[tail + 17] ^ pool[tail + 23];
    pool[tail] = result + 1;
    tail = (tail + 1) % N;
    return result;
}

> // псевдо код с++
> unsinged poo[lN]; // кольцевой буфер
> usigned head, tail; // при старте == 0
> 
> void push_salt(unsigned salt)
> {
>     unsigned value = pool[head] + salt;
>     pool[head] = value + (value < salt); // кольцевое
> сложение с переносом
>     head = (head + 1) % N;
> }
> 
> unsigned pop_salt()
> {
>     // y ~~ F(X**N, X**23 + X**17);
>     unsigned result = rotate(pool[tail], 3);

>     result += pool[tail + 17] ^ pool[tail + 23];

Возможен выход индекса за границы массива.

>     pool[tail] = result + 1;
>     tail = (tail + 1) % N;
>     return result;
> }
> 

static struct
{
    unsigned Head, Tail;
    unsigned __int32 Pool[64];
} Random;

void LYSAP_RandomPushSalt(unsigned Salt)
{
    Random.Pool[Random.Head] += Salt;
    Random.Head = (Random.Head + 1) & 63;
}

void LYSAP_RandomInit(unsigned Seed)
{
    Random.Tail = 0;
    Random.Head = 0;

    for(unsigned i = 0; i < 64; i++)
        // можно заполнить нулями, но тогда первые
        // 4096 сгенерированных чисел лучше пропустить
        Random.Pool[i] = (Seed = 1664525ul * Seed + 1013904223ul);
}

unsigned LYSAP_GetRandom()
{
    unsigned __int32 Value =
        Random.Pool[Random.Tail]                // X ** 1
        ^ Random.Pool[(Random.Tail + 1) & 63]   // X ** 64
        ^ Random.Pool[(Random.Tail + 62) & 63]  // X ** 4
        ^ Random.Pool[(Random.Tail + 61) & 63]; // X ** 3

    unsigned __int32 Salt =
        1 // may be any non-zero
        + ROT_L(Random.Pool[(Random.Tail + 11) & 63], 5) 
        + ROT_L(Random.Pool[(Random.Tail + 23) & 63], 7)
        + ROT_L(Random.Pool[(Random.Tail + 37) & 63], 1);

    Value += Salt; Value += Value < Salt; // cyclic addition with carry
    Random.Pool[Random.Tail = (Random.Tail + 1) & 63] = Value; 
    return Value;
}

> unsigned LYSAP_GetRandom()
> {
>     unsigned __int32 Value =
> 	Random.Pool[Random.Tail]		// X ** 1
> 	^ Random.Pool[(Random.Tail + 1) & 63]	// X ** 64
> 	^ Random.Pool[(Random.Tail + 62) & 63]	// X ** 4
> 	^ Random.Pool[(Random.Tail + 61) & 63]; // X ** 3
> 
>     unsigned __int32 Salt =
> 	1 // may be any non-zero
> 	+ ROT_L(Random.Pool[(Random.Tail + 11) & 63], 5) 
> 	+ ROT_L(Random.Pool[(Random.Tail + 23) & 63], 7)
> 	+ ROT_L(Random.Pool[(Random.Tail + 37) & 63], 1);
> 
>     Value += Salt; Value += Value < Salt; // cyclic
> addition with carry
>     Random.Pool[Random.Tail = (Random.Tail + 1) & 63] =
> Value; 
>     return Value;
> }